#------------------------------------------------------------ # Exemplo apresentado no slide da Patologia Experimental - # sobre a duração da visita domiciliar, idade do enfermeiro - # e da tipo de doença do paciente. Livro do Wayne W. Daniel - # Biostatistics: A foundation for analysis in the health - # sciences, pág. 349 - # - # Os dados originais no livro não atendem os pressupostos - # para validação da análise de variância. Assim, alterei - # alguns dados, forçando a validação dos pressupostos. - #------------------------------------------------------------ rm(list=ls()) Dados = read.csv2('C:/Aulas/Ano2013/Patologia/R/Fatorial/pag349.csv') attach(Dados) str(Dados) head(Dados) #------------- # Descritiva - #------------- (média.g = mean(tempo)) (variância.g = var(tempo)) (desvio.g = sd(tempo)) (cv.g = desvio.g / média.g * 100) (médias = tapply(tempo, idade, mean)) (variâncias = tapply(tempo, idade, var)) (desvios = tapply(tempo, idade, sd)) (cv = desvios / médias * 100) #-------------------------- # Análise de forma direta - #-------------------------- anava.fat = aov(tempo ~ trat) summary(anava.fat) bartlett.test(tempo ~ trat) shapiro.test(anava.fat$res) anava.fat = aov(tempo ~ tipo + idade + tipo * idade, data=Dados) summary(anava.fat) #-------------------------- # Usando o pacote do Eric - #-------------------------- require(ExpDes.pt) fat2.dic(tipo, idade, tempo, quali=c(TRUE, TRUE), fac.names=c("tipo", "Idade")) #------------------ # Página 355 - # Danos cerebrais - #------------------ rm(list=ls()) (fis = factor(rep(c('I','II','III','IV','V','VI'), each=12))) (psiq = factor(rep(LETTERS[1:4], 18))) (Trat = factor(interaction(fis,psiq))) (meses = c(11.0,9.4,12.5,13.2, 9.6,9.6,11.5,13.2, 10.8,9.6,10.5,13.5, 10.5,10.8,10.5,15.0, 11.5,10.5,11.8,14.6, 12.0,10.5,11.5,14.0, 12.0,11.5,11.8,12.8, 11.5,11.5,11.8,13.7, 11.8,12.3,12.3,13.1, 11.5,9.4,13.7,14.0, 11.8,9.1,13.5,15.0, 10.5,10.8,12.5,14.0, 11.0,11.2,14.4,13.0, 11.2,11.8,14.2,14.2, 10.0,10.2,13.5,13.7, 11.2,10.8,11.5,11.8, 10.8,11.5,10.2,12.8, 11.8,10.2,11.5,12.0)) (danos = data.frame(Trat, fis, psiq, meses)) attach(danos) names(danos) str(danos) #-------------------- # Para uso no Latex - #-------------------- # d = matrix(meses, 18, 4, byrow=T) # require(xtable) # xtable(d) anava.danos = aov(meses~Trat) summary(anava.danos) shapiro.test(anava.danos$res) bartlett.test(meses~Trat) anava.d.fat = aov(meses~fis+psiq+fis*psiq, data=danos) summary(anava.d.fat) require(ExpDes.pt) fat2.dic(fis, psiq, meses, quali=c(TRUE,TRUE), fac.names=c('fis','psiq'))