[R-br] Parcelas sub-divididas com fatorial nas parcelas

[Tiago Souza Marçal]

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