Boa noite pessoal, 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 #################################################################