<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Não sei se eu entendi bem mas a sua duvida é saber qual dos modelos abaixo usar?<div><br></div><div>Yijk = m + Gi + B/Ajk + Aj + GAij + Eijk</div><div><br></div><div>Yijk = m + Gi + Bk + Aj + GAij + Eijk</div><div><br></div><div>Att.</div><div><br></div><div>Tiago.<br><br><div><div id="SkyDrivePlaceholder"></div><hr id="stopSpelling">Date: Sat, 20 Oct 2012 15:36:31 -0300<br>From: eder@leg.ufpr.br<br>To: R-br@listas.c3sl.ufpr.br<br>Subject: Re: [R-br] Blocos com algumas repetições do material<br><br><div class="ecxgmail_quote"><br><blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex"><div>Boa tarde Pessoal,</div><div><br></div><div>Tenho experimentos em blocos(Locais) e dentro deles gostaria de por algumas de repetições de alguns tratamentos para estimar a variância dentro do bloco, que desconfio ser grande. Por exemplo, tenho 5 Cultivares gostaria de repetir 3 delas dentro do bloco, tendo ao todo 8 parcelas, o ideal seria ter uma repetição de tudo, porem na pratica isso não é possível. </div>
<div>Minhas duvidas são:</div><div>1) Estou definindo de maneira correta essas repetições dentro do modelo?</div><div>2) Como o método de calculo dos dois modelos abaixo é igual para os efeitos fixos, por que o intercepto esta com valores diferentes?</div>
<div>3) Neste dois casos(m0 e m1, abaixo) o BLUP seria a média mais o efeito da cultivar? </div><div>Obrigado pela ajuda</div><div><br></div><div>require(lme4)</div><div>### Simulando dados</div><div>dados <- expand.grid(Variedade=factor(c(letters[1:3],letters[1:5])),Local=factor(1:5))</div>
<div>dados$Rep <- factor(rep(c(2,2,2,1,1,1,1,1),t=5))</div><div>dados$Resp <- sort(rnorm(nrow(dados)))</div><div>dados</div><div><br></div><div>### Modelos</div><div># Sem efeito dentro do local</div><div>m0 <- lmer(Resp~1+(1|Local)+(1|Variedade),dados)</div>
<div>summary(m0)</div><div>ranef(m0)</div><div># Com efeito dentro do local</div><div>m1 <- lmer(Resp~1+(0+Rep|Local)+(1|Variedade),dados)</div><div>summary(m1)</div><div>ranef(m1)</div></blockquote></div>
<br>_______________________________________________
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.</div></div> </div></body>
</html>