<div dir="ltr">Estou enviando um material aqui, leia o capitulo 3 e 4 que nele esta sua solução.</div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 19 de novembro de 2013 00:52, FHRB Toledo <span dir="ltr"><<a href="mailto:fernandohtoledo@gmail.com" target="_blank">fernandohtoledo@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra">Senhores,</div><div class="gmail_extra"><br></div><div class="gmail_extra">Seguindo no assunto proposto pelo Alisson, no thread sobre pvalor com a lme4, estou tentando parametrizar um modelo misto com o intuito de fazer a inferência pelo perfil da verossimilhança.</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Entretanto, estou tentando "traduzir" um modelo da lm() pela parametrização da lme4 e enfrento grandes dificuldades, assim gostaria da ajuda de vocês com esse intento...</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Seguem dois modelos da lm() que são equivalentes, visto que o fator "trat" tem 9 níveis e "eras" tem um certo confundimento com "trat", assim procedo:</div>
<div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_extra">mf1 <- lm(terms(resp ~ local/rep + trat + local:trat, keep = TRUE), data = dados)</div><div class="gmail_extra">mf2 <- lm(terms(resp ~ local/rep + eras/trat + local:eras/trat, keep = TRUE), data = dados)</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Esses dois modelo tem o mesmo erro e gostaria de reproduzi-los usando a lmer(). Quando faço o primeiro modelo, sem problemas:</div><div class="gmail_extra"><br>
</div>
<div class="gmail_extra">mm1 <- lmer(resp ~ ( 1 | local/rep ) + ( 1 | trat ) + ( 1 | local:trat ), data = dados)<br></div><div class="gmail_extra"><br></div><div class="gmail_extra">Recupero a mesma variância do resíduo da lm(), todavia não consgigo fazer o segundo modelo recuperar essa mesma variância... Alguém tem alguma sugestão?</div>
</div><div class="gmail_extra"><br></div><div class="gmail_extra">P.S.: Se, por acaso alguém tiver outra sugestão, que use a lme(), ou lme2() gostaria de saber... como disse, meu objetivo a obter os perfis da verossimilhança!</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Seja os dados de exemplo:</div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_extra">dados <- dput(dados)</div><div class="gmail_extra">
structure(list(local = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, </div><div class="gmail_extra">1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, </div><div class="gmail_extra">1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, </div>
<div class="gmail_extra">2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, </div><div class="gmail_extra">2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, </div><div class="gmail_extra">2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, </div>
<div class="gmail_extra">3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, </div><div class="gmail_extra">3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, </div><div class="gmail_extra">4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, </div>
<div class="gmail_extra">4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("1", "2", "3", </div><div class="gmail_extra">"4"), class = "factor"), rep = structure(c(1L, 1L, 1L, 1L, 1L, </div>
<div class="gmail_extra">1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, </div><div class="gmail_extra">3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 1L, </div><div class="gmail_extra">1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, </div>
<div class="gmail_extra">2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, </div><div class="gmail_extra">4L, 4L, 4L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, </div><div class="gmail_extra">2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, </div>
<div class="gmail_extra">4L, 4L, 4L, 4L, 4L, 4L, 4L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, </div><div class="gmail_extra">2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, </div><div class="gmail_extra">3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("1", </div>
<div class="gmail_extra">"2", "3", "4"), class = "factor"), trat = structure(c(1L, 6L, </div><div class="gmail_extra">2L, 3L, 5L, 8L, 7L, 4L, 9L, 1L, 3L, 7L, 6L, 5L, 4L, 2L, 8L, 9L, </div>
<div class="gmail_extra">1L, 8L, 4L, 6L, 3L, 9L, 2L, 5L, 7L, 1L, 5L, 9L, 6L, 8L, 7L, 2L, </div><div class="gmail_extra">3L, 4L, 9L, 8L, 3L, 2L, 7L, 4L, 1L, 6L, 5L, 9L, 2L, 1L, 8L, 7L, </div><div class="gmail_extra">6L, 3L, 4L, 5L, 9L, 4L, 6L, 8L, 2L, 5L, 3L, 7L, 1L, 9L, 7L, 5L, </div>
<div class="gmail_extra">8L, 4L, 1L, 3L, 2L, 6L, 4L, 9L, 2L, 3L, 8L, 6L, 5L, 1L, 7L, 4L, </div><div class="gmail_extra">3L, 5L, 9L, 8L, 1L, 2L, 6L, 7L, 4L, 6L, 1L, 9L, 3L, 7L, 2L, 8L, </div><div class="gmail_extra">5L, 4L, 8L, 7L, 9L, 6L, 5L, 2L, 3L, 1L, 9L, 2L, 7L, 8L, 4L, 5L, </div>
<div class="gmail_extra">6L, 3L, 1L, 9L, 8L, 6L, 2L, 4L, 3L, 7L, 5L, 1L, 9L, 5L, 3L, 2L, </div><div class="gmail_extra">8L, 1L, 7L, 4L, 6L, 9L, 4L, 1L, 2L, 5L, 6L, 7L, 8L, 3L), .Label = c("a", </div><div class="gmail_extra">
"b", "c", "d", "e", "f", "g", "h", "i"), class = "factor"), eras = structure(c(2L, </div><div class="gmail_extra">4L, 2L, 3L, 3L, 1L, 4L, 3L, 1L, 2L, 3L, 4L, 4L, 3L, 3L, 2L, 1L, </div>
<div class="gmail_extra">1L, 2L, 1L, 3L, 4L, 3L, 1L, 2L, 3L, 4L, 2L, 3L, 1L, 4L, 1L, 4L, </div><div class="gmail_extra">2L, 3L, 3L, 1L, 1L, 3L, 2L, 4L, 3L, 2L, 4L, 3L, 1L, 2L, 2L, 1L, </div><div class="gmail_extra">4L, 4L, 3L, 3L, 3L, 1L, 3L, 4L, 1L, 2L, 3L, 3L, 4L, 2L, 1L, 4L, </div>
<div class="gmail_extra">3L, 1L, 3L, 2L, 3L, 2L, 4L, 3L, 1L, 2L, 3L, 1L, 4L, 3L, 2L, 4L, </div><div class="gmail_extra">3L, 3L, 3L, 1L, 1L, 2L, 2L, 4L, 4L, 3L, 4L, 2L, 1L, 3L, 4L, 2L, </div><div class="gmail_extra">1L, 3L, 3L, 1L, 4L, 1L, 4L, 3L, 2L, 3L, 2L, 1L, 2L, 4L, 1L, 3L, </div>
<div class="gmail_extra">3L, 4L, 3L, 2L, 1L, 1L, 4L, 2L, 3L, 3L, 4L, 3L, 2L, 1L, 3L, 3L, </div><div class="gmail_extra">2L, 1L, 2L, 4L, 3L, 4L, 1L, 3L, 2L, 2L, 3L, 4L, 4L, 1L, 3L), .Label = c("0", </div><div class="gmail_extra">
"1", "2", "3"), class = "factor"), resp = c(1878L, 3417L, 3546L, </div><div class="gmail_extra">2284L, 3220L, 2692L, 2243L, 2088L, 1871L, 1793L, 1983L, 2087L, </div><div class="gmail_extra">
2628L, 2909L, 1928L, 3065L, 2650L, 2390L, 2428L, 3125L, 3010L, </div><div class="gmail_extra">2181L, 2478L, 2575L, 2348L, 2492L, 2295L, 1892L, 2257L, 2003L, </div><div class="gmail_extra">3123L, 2673L, 2449L, 3333L, 3040L, 3212L, 1767L, 2514L, 2503L, </div>
<div class="gmail_extra">2547L, 2055L, 1591L, 1898L, 2748L, 2376L, 1458L, 1709L, 1649L, </div><div class="gmail_extra">1784L, 1724L, 1956L, 1706L, 1043L, 1056L, 1346L, 1657L, 1970L, </div><div class="gmail_extra">1483L, 1103L, 985L, 779L, 915L, 800L, 1155L, 1074L, 1435L, 988L, </div>
<div class="gmail_extra">1052L, 932L, 1112L, 1728L, 2746L, 5150L, 3101L, 5020L, 5828L, </div><div class="gmail_extra">3927L, 5123L, 5765L, 3620L, 4901L, 5859L, 5050L, 5383L, 2920L, </div><div class="gmail_extra">4002L, 3850L, 4727L, 4998L, 5702L, 5154L, 5098L, 4039L, 3493L, </div>
<div class="gmail_extra">5417L, 5097L, 4954L, 3964L, 4371L, 5681L, 4090L, 4644L, 3370L, </div><div class="gmail_extra">5529L, 5491L, 5145L, 5158L, 3864L, 2777L, 2697L, 2739L, 1911L, </div><div class="gmail_extra">2103L, 2959L, 3241L, 2457L, 2321L, 2410L, 2085L, 3511L, 2657L, </div>
<div class="gmail_extra">3349L, 3157L, 2948L, 3491L, 2951L, 2626L, 2396L, 1997L, 2895L, </div><div class="gmail_extra">2244L, 2807L, 3408L, 3350L, 4104L, 2613L, 2589L, 2617L, 2791L, </div><div class="gmail_extra">3063L, 3665L, 2918L, 2092L, 2583L)), .Names = c("local", "rep", </div>
<div class="gmail_extra">"trat", "eras", "resp"), class = "data.frame", row.names = c(NA, </div><div class="gmail_extra">-144L))</div><div class="gmail_extra"><br></div><div class="gmail_extra">
att,</div><div class="gmail_extra">FH</div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>