[R-br] Analise variância fatorial com tratamento adicional 3x2+1

[Fernando Souza]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140113/ea146324/attachment.html>