Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
| Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior | ||
|
playground:playground [2011/05/13 19:41] walmes [section 1] |
playground:playground [2018/05/03 03:39] (atual) |
||
|---|---|---|---|
| Linha 1: | Linha 1: | ||
| ====== PlayGround ====== | ====== PlayGround ====== | ||
| - | |||
| - | <code R> | ||
| - | #------------------------------------------------------------------------------------------ | ||
| - | |||
| - | 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")) | ||
| - | |||
| - | #------------------------------------------------------------------------------------------ | ||
| - | </code> | ||
