estou com experimento de três fatores e gostaria da opinião de vocês para modelar uma análise de variância.
De forma mais direta eu tenho dois ambientes e em cada um deles tenho um experimento em arranjo fatorial em DIC com um
tratamento adicional.
#CMR
set.seed(102)
resp <- rnorm(100,50,13)
resp
# O primeiro fator se refere a dois ambientes (temperaturas diferentes).
temp <- gl(2,50)
# O segundo fator se refere a duas substâncias aplicadas e um controle com a ausência das substâncias dentro de cada um dos
# ambientes.
subs <- factor(c(gl(2,20),rep("Adicional1",10),gl(2,20),rep("Adicional2",10)))
# O terceiro fator se refere a duas doses de cada uma das substâncias e um controle (dose 0 - já citado a cima)
# dentro de cada um dos ambientes.
dose <- factor(c(gl(2,10),gl(2,10),rep("Adicional1",10),gl(2,10),gl(2,10),rep("Adicional2",10)))
O arquivo de dados seria o seguinte:
dat <- data.frame(temp,subs,dose,resp)
Por fim montei um quadro para a provável ordem das fontes de variação que consegui enxergar para o modelo. Vocês concordam
com a ordem e com as fontes de variação? Alguém tem ideia de como modelar este problema via função aov() ou lm()?
| FV |
| temp |
| resíduo 1 |
| subs |
| dose |
| subs x dose |
| subs x temp |
| dose x temp |
| subs x dose x temp |
| temp1 / adicional1 x fatorial1 |
| temp2 / adicional2 x fatorial2 |
| resíduo 2 |
Agradeço desde já por qualquer ajuda.
Att.
Tiago.
#################################################################
Tiago de Souza Marçal - Engenheiro Agrônomo pelo CCA-UFES
Mestrando em Genética e Melhoramento de Plantas pelo CCA-UFES
#################################################################