[R-br] Blocos com algumas repetições do material

Eder David Borges da Silva eder em leg.ufpr.br
Segunda Outubro 22 09:14:51 BRST 2012


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 em 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 em gmail.com>
> *Para:* r-br em 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 em ufpr.br
> skype: walmeszeviani
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em 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 em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121022/6ea82637/attachment.html>


Mais detalhes sobre a lista de discussão R-br