<div dir="ltr"><div class="gmail_default"><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 class="gmail_signature"><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 class="gmail_default" 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 class="gmail_default" 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>