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 < [email protected]> 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 <[email protected]> > *Para:* [email protected] > *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: [email protected] > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > R-br mailing list > [email protected] > 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 > [email protected] > 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 [email protected] 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.
