
13 Mar
2012
13 Mar
'12
12:36
Prof. Ivan, com certeza não... há grandes efeitos de interação na análise gráfica e meu sentimento é que não há esta independência. Pensando agora no problema, poderia tratar, também como ANCOVA... mas não sei se perco muita informação. ############## Como é meu dataframe ############# Id: Unidade Amostral #Estrutura dos dados: set.seed(13032012) d=expand.grid(Trat1=gl(3,18),Tempo=gl(2,1)) d$Trat2=gl(2,9,108) d$Id=as.factor(rep(1:54,2)) d$Y=abs(rnorm(108,c(10,30,8),10)) #Verificar estrutura para a unidade amostral i=3, por exemplo d[d$Id==3,]: Tratamento1 =1, Tratamento2=1 Tempo=1, 2. #Ajuste que fiz usando lme library(nlme) f=lme(Y~Trat1*Trat2*Tempo, random=~1|Id,d) Abs e muito obrigado. Olympio