Caros amigos Estou analisando uns dados de produção de gases onde foi avaliado o efeito da adição de 3 níveis de farelo de coco em uma ração base (feno+soja+Milho) sobre a produção de metano. Foi utilizado um DIC com arranjo fatorial 3 x 2 (nível de farelo adicionado x Industria de origem do farelo). Nesse experimento a produção de metano da dieta sem adição de farelo (nível controle tambem foi avaliada. Eu quero avaliar este experimento em relação a dieta controle? Como posso fazer isso, uma vez que não tenho um fatorial para o controle? Penso que este experimento deva ser avaliado como um fatorial com tratamento adicional 3 x 2 + 1, mas não sei como fazê-lo. Alguém pode me orientar? dados<-structure(list(Amostra = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), Industria = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), Metano = c(23.82046, 22.58261, 6.90982, 9.52078, 18.04296, 24.80488, 19.28557, 13.93736, 19.13685, 8.62221, 16.07414, 9.38481, 10.82267, 17.51822, 22.07641, 23.04412, 22.8202, 26.22151, 8.506, 19.75686, 2.87987, 14.44739, 13.23224, 2.10336, 9.61595, 1.94412, 8.66299, 4.49037, 14.61952, 15.54833, 17.3914, 1.89799, 15.31825, NA, 14.09188, 4.86393, NA, 10.15676, 4.80198, 14.74247, 18.07063, 8.45602, 15.9324, 8.86496, 13.61907, 4.17938, 16.34748, 14.99311, 12.23218)), .Names = c("Amostra", "Industria", "Metano"), class = "data.frame", row.names = c(NA, -49L)) Amostra= Nível de farelo adicionado (nível 0 é o controle) Industria= (Industria de origem do farelo), há duas industria avaliada dados<-transform(dados,Amostra=factor(Amostra),Industria=factor(Industria)) ANALISE<-aov(Metano~Amostra*Industria,data=dados[-c(32,3,21),]) drop1(ANALISE,~.,test="F") shapiro.test(resid(ANALISE) --- Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus está ativa. http://www.avast.com