Caro Walmes,<div>Muito obrigada pela ajuda!! </div><div>Você tem algum livro do R sobre modelos mistos em pdf? Ou qual você indicaria?</div><div><br></div><div>Obrigada,</div><div>Catarina<br><br><div class="gmail_quote">2012/4/14 Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>></span><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">Bem, vamos primeiro recapitular para dar check list.<br>Você tem 11 unidades, 6 no nível A e 5 no nível B de capoeira (fator da parcela). Cada uma dessas 11 foi divida para receber os 2 níveis do fator abertura (0% e 35%), então são as 22 subparcelas. Em cada uma das 22 subparcelas foram colocados 45 exemplares dos dois níveis de espécie (I e II). A unidade experimental é o exemplar ou muda. Até aqui são 11*2*2*45=1980 registros. A covariável idade (age) é valor observado por parcela, portanto 11. A covariável luz é valor observado por muda, portanto, 1980. Assim nos temos que os termos de efeito aleatório são parcelas (11 níveis), subparcelas (ou capoeira:abertura, 22 níveis) e muda (1980 níveis). Este último é o desvio resídual. Os termos de efeito fixo são tipo de capoeira (2 níveis), abertura de capoeira (2 níveis), espécie (2 níveis) e as covariáveis idade e luz. Então na lme() poderia se declarar assim (veja que o CMR não precisa conter os dados verdadeiros, apenas representar o design corretamente),<br>
<br><span style="font-family:courier new,monospace">da <- expand.grid(uni=1:11, abr=c("0%","35%"), esp=c("I","II"), mud=1:45)</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">str(da)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">da$cap <- factor(ifelse(da$uni<=5, "A", "B"))</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">da$age <- runif(11)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">da$luz <- runif(nrow(da))</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">da$uni <- factor(da$uni, ordered=FALSE)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">da$sub <- with(da, interaction(uni, abr))</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">head(da, 20)</span><br style="font-family:courier new,monospace"><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">lapply(da, function(x) if(is.factor(x)) levels(x) else NULL)</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">da$y <- rnorm(nrow(da))</span><br style="font-family:courier new,monospace"><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">require(nlme)</span><br style="font-family:courier new,monospace">
<br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">m0 <- lme(y~age+luz+cap*abr*esp, random=~1|uni/abr, data=da)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">summary(m0) # estimativas sob uma particular restrição paramétrica</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">anova(m0) # teste de Wald para os termos de efeito fixo do modelo, similar a ANOVA.</span><br style="font-family:courier new,monospace"><br>O que você indica como erro que sair nomes como no summary() abaixo é padrão do R e acho que você estava pensando que isso era uma anova. São as estimativas dos parâmetros sob uma determinada restrição paramétrica, nesse caso a de zerar o efeito do primeiro nível (ordem alfanumérica) de cada fator. Isso está explicado nos livros de R que abordam uso da função lm(), glm() e demais.<br>
<br><span style="font-family:courier new,monospace"> Value Std.Error DF t-value p-value</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">(Intercept) -0.14754343 0.10744471 1953 -1.3732033 0.1698</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">age 0.22709384 0.11910516 8 1.9066667 0.0930</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">luz -0.00750802 0.08007717 1953 -0.0937598 0.9253</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">capB 0.09953498 0.09803749 8 1.0152746 0.3397</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">abr35% -0.06785446 0.09670281 9 -0.7016803 0.5006</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">espII 0.06385534 0.09657213 1953 0.6612192 0.5085</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">capB:abr35% 0.19696937 0.13092688 9 1.5044227 0.1667</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">capB:espII -0.12451915 0.13076225 1953 -0.9522561 0.3411</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">abr35%:espII 0.05027986 0.13669501 1953 0.3678252 0.7130</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">capB:abr35%:espII -0.20166323 0.18506593 1953 -1.0896832 0.2760</span><br style="font-family:courier new,monospace"><br>À disposição.<span class="HOEnZb"><font color="#888888"><br>
Walmes.<br><br><br clear="all"></font></span></font><div class="HOEnZb"><div class="h5"><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: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a></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>
</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><br clear="all"><div><br></div>-- <br>
<div>MSc. Ana Catarina C. Jakovac</div>PhD. Candidate<div>Forest Ecology and Forest Management Group</div><div><br></div><div>Centre for Ecosystem Studies</div><div>Wageningen University and Research Centre</div><div>P.O. Box 47 - Lumen Building room 1.212</div>
<div>6700 AA Wageningen, The Netherlands</div><div><br></div><br>
</div>