Walmes,<div>Problema resolvido, valeu..</div><div>Ate mais<br><br><div class="gmail_quote">Em 4 de junho de 2011 19:46, Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><font face="trebuchet ms,sans-serif">Éder,<br><br>Se você quer obter essa matriz saturada, então é só remover os interceptos de cada porção e parear as matrizes (como o Benilton sugeriu), ó<br>
<br><span style="font-family:courier new,monospace">terms <- c(~-1+BLOCO, ~-1+CLONE, ~-1+LOCAL:CLONE)</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">str(terms)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">X <- cbind("(Intercept)"=1, do.call(cbind, lapply(terms, model.matrix, data=dados)))<br>
<br><span style="font-family:trebuchet ms,sans-serif">Para extrair a matriz de ajustes da lmer tem use abaixo</span><br><br>fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)<br>summary(fm1)<br>model.matrix(fm1) # efeitos fixos, X<br>
str(fm1)<br>str(fm1@Zt)<br>t(as.matrix(fm1@Zt)) # efeitos aleatórios, Z<br></span><br>À disposição.<br>Walmes.<br><br clear="all"></font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
<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>
<br></blockquote></div><br></div>