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 [2013/08/31 17:51] fmoro [section 1] |
||
---|---|---|---|
Linha 1: | Linha 1: | ||
====== PlayGround ====== | ====== PlayGround ====== | ||
- | + | [[pessoais:fernando|Fernando G. Moro]] | |
- | <code R> | + | [[pessoais:alguem|Alguem Silva]] |
- | #------------------------------------------------------------------------------------------ | + | |
- | + | ||
- | 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> | + | |