[R-br] Ajuda com modelos mistos (medidas repetidas)

Olympio Neto olympiotneto em gmail.com
Terça Março 13 09:36:35 BRT 2012


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



Mais detalhes sobre a lista de discussão R-br