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.