Fábio, Não tem bloco dentro do local, foi em DIC, apenas uma repetição dentro de cada local. Att Em 22 de outubro de 2012 07:05, Fabio Mathias Corrêa < fabio.ufla@yahoo.com.br> escreveu:
Não tem bloco dentro de local?
Foi inteiramente casualizado?
Valeu!
Fábio Mathias Corrêa
Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET
Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia
Tel.: 73-3680-5076 ------------------------------ *De:* Walmes Zeviani <walmeszeviani@gmail.com> *Para:* r-br@listas.c3sl.ufpr.br *Enviadas:* Domingo, 21 de Outubro de 2012 16:30 *Assunto:* Re: [R-br] Blocos com algumas repetições do material
Éder,
Agora entendi suas colocações. Bem, dessa maneira acredito que o modelo deva ser declarado dessa maneira.
#------------------------------------------------------------------------------------------ # modelo com efeito aleatório e local, cultivar, local:cultivar e resídual # só é possível estimar resídual se haver repetições dentro de locais:cultivar
# número de níveis nl <- 10; nc <- 8; nr <- 3 # efeitos aleatórios l <- rnorm(nl,0,2); # de local c <- rnorm(nc,0,1.5); # de cultivar lc <- rnorm(nl*nc,0,1); # de local:cultivar (interação genótipo:ambiente) e <- rnorm(nl*(nc+nr),0,0.5) # residual, var entre rep no mesmo local:cult
# dados artificiais da1 <- expand.grid(local=gl(nl,1), cultivar=gl(nc,1), r=1) da2 <- expand.grid(local=gl(nl,1), cultivar=gl(nr,1), r=2) da <- rbind(da1, da2)
# matriz de delineamento Z <- with(da, cbind(model.matrix(~-1+local), model.matrix(~-1+cultivar), model.matrix(~-1+local:cultivar))) dim(Z) ranef <- c(l, c, lc) length(ranef) y <- Z%*%ranef+e
require(lme4)
m0 <- lmer(y~(1|local)+(1|cultivar)+(1|local:cultivar), data=da) summary(m0)
#------------------------------------------------------------------------------------------
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.