PlayGround

#------------------------------------------------------------------------------------------
 
da <- expand.grid(trat=gl(5,1,la="tr"),
                  bloc=gl(4,1,la="bl"),
                  temp=1:7)
da$parc <- with(da, factor(paste(trat, bloc, sep="-"))) # fator parcela
trat.ef <- c(1,2,-1,0)                   # t-1 efeitos fixos para tratamento
bloc.ef <- c(1,0,-0.5)                   # b-1 efeitos fixos para bloco
parc.ef <- rnorm(nlevels(da$parc), 0, 2) # p efeitos aleatórios para parcela
temp.ef <- 0.5                           # efeito fixo do tempo
 
da$y <-
  model.matrix(~bloc+trat, da)%*%c(1, bloc.ef, trat.ef)+ # parte fixa bloc+trat
  model.matrix(~-1+parc, da)%*%parc.ef+ # parte aleatória de parcela
  model.matrix(~-1+temp, da)%*%temp.ef+ # parte fixa de tempo
  rnorm(nrow(da), 0, 0.5) # erro aleatório
 
require(lattice)
 
xyplot(y~temp|bloc, groups=trat, data=da)
 
require(nlme)
 
mm0 <- lme(y~bloc+trat*temp, # efeitos fixos
           random=~1|parc,   # efeito aleatório de parcela no intercepto
           data=da)
summary(mm0)                 # tabela de estimativas
ranef(mm0)                   # efeitos aleatórios preditos
 
p0 <- predict(mm0, level=0)  # predição só com efeitos fixos
p1 <- predict(mm0, level=1)  # predição com os efeitos aleatórios
 
xyplot(y+p0+p1~temp|bloc*trat, data=da,
       distribute.type=TRUE, type=c("p","l","l"))
 
#------------------------------------------------------------------------------------------