[ Por favor, preencha este local (:wiki:site_notice) ]

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- playground:playground [2011/05/13 19:41] – [section 1] walmes
+++ playground:playground [2018/05/03 03:39] (atual) – edição externa 127.0.0.1
@@ Linha 1: / Linha 1: @@
 ====== 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>

Diferenças

Navegação

Busca

Ferramentas

QR Code