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
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
