<div dir="ltr"><div class="gmail_default" style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0)">Prezados, </div><div class="gmail_default" style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0)">O problema foi resolvido. Abaixo a solução:</div><div class="gmail_default" style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0)"><div class="gmail_default">fs=as.numeric()</div><div class="gmail_default">for(i in 2:ncol(y)){</div><div class="gmail_default">m1 <- capscale((y[,1:i]) ~ x, dist="bray")</div><div class="gmail_default">fs[i]=anova(m1)$F[1]</div><div class="gmail_default">}</div></div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_default" style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0)">​Nicolay​</div><br clear="all"><div><div class="gmail_signature"></div></div>
<br><div class="gmail_quote">2015-09-21 16:28 GMT-04:00 Nicolay Cunha <span dir="ltr"><<a href="mailto:nicolaycunha@gmail.com" target="_blank">nicolaycunha@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><font color="#000000" face="trebuchet ms, sans-serif">Prezados</font><span style="color:rgb(34,34,34);font-family:arial,sans-serif">, </span><span style="color:rgb(34,34,34);font-family:arial,sans-serif"> </span><br></div><div><div><div dir="ltr"><div><p class="MsoNormal"></p>

<p class="MsoNormal"> </p>

<p class="MsoNormal">Gostaria de ajuda para
automatizar a inclusão de variáveis dependentes em um modelo de RDA. </p>

<p class="MsoNormal"> </p>

<p class="MsoNormal">Usarei o valor de F
como critério para seleção do número de variáveis Y a serem usadas
posteriormente, entretanto, preciso de uma forma automatizada para incluir as
59 colunas que tenho em minha matriz.</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal">Basicamente o que
preciso fazer é calcular RDA</p><div style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0);display:inline">​s​</div> com o número crescente de Ys. <p></p>

<p class="MsoNormal"> </p>

<p class="MsoNormal">Abaixo um CMR.</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal"> </p>

<p class="MsoNormal"># library(vegan)</p>

<p class="MsoNormal">y <-
matrix(sample(rnorm(3540)+10), 60)</p>

<p class="MsoNormal">x <-
rep(c("A", "B", "C"), each = 20)</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal"> </p>

<p class="MsoNormal"># Aqui gero o modelo
com somente as colunas 1 e 2 da matriz y</p>

<p class="MsoNormal">m1 <- capscale((y[,1:2])
~ x, dist="bray")</p>

<p class="MsoNormal">anova(m1)</p><p class="MsoNormal"><br></p>

<p class="MsoNormal"># agora incluindo a
terceira coluna</p>

<p class="MsoNormal">m2 <-
capscale((y[,1:3]) ~ x, dist="bray")</p>

<p class="MsoNormal">anova(m2)</p><p class="MsoNormal"><br></p>

<p class="MsoNormal"># adicionando a quarta
coluna</p>

<p class="MsoNormal">m3 <-
capscale((y[,1:4]) ~ x, dist="bray")</p>

<p class="MsoNormal">anova(m3)</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal"># e assim até incluir
as 59 colunas na RDA</p>

<p class="MsoNormal">m59 <-
capscale((y[,1:59]) ~ x, dist="bray")</p>

<p class="MsoNormal">anova(m59)</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal">Feito isso, preciso acessar
os valores de F gerados para cada modelo para fazer a comparação.</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal">Alguém saberia me
ajudar com isso?</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal">Muito Obrigado,</p><p class="MsoNormal"><br></p>

<p class="MsoNormal"></p><div style="font-family:'trebuchet ms',sans-serif;color:rgb(0,0,0)"><span style="font-family:arial,sans-serif;color:rgb(34,34,34)">Nicolay</span>​</div><br></div><div><br></div><div><br></div><div><br></div></div></div></div>
</div>
</blockquote></div><br></div></div>