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.