[R-br] Automatizar a inclusão de váriaveis

[Nicolay Cunha]

Prezados,



Gostaria de ajuda para automatizar a inclusão de variáveis dependentes em
um modelo de RDA.



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.



Basicamente o que preciso fazer é calcular RDA
​s​
com o número crescente de Ys.



Abaixo um CMR.





# library(vegan)

y <- matrix(sample(rnorm(3540)+10), 60)

x <- rep(c("A", "B", "C"), each = 20)





# Aqui gero o modelo com somente as colunas 1 e 2 da matriz y

m1 <- capscale((y[,1:2]) ~ x, dist="bray")

anova(m1)


# agora incluindo a terceira coluna

m2 <- capscale((y[,1:3]) ~ x, dist="bray")

anova(m2)


# adicionando a quarta coluna

m3 <- capscale((y[,1:4]) ~ x, dist="bray")

anova(m3)



# e assim até incluir as 59 colunas na RDA

m59 <- capscale((y[,1:59]) ~ x, dist="bray")

anova(m59)



Feito isso, preciso acessar os valores de F gerados para cada modelo para
fazer a comparação.



Alguém saberia me ajudar com isso?



Muito Obrigado,


Nicolay​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150921/f12bba70/attachment.html>