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

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,