Получение p-значений для анализа RDA с использованием пакета R "Веган"

Question

Я провожу анализ избыточности с использованием пакета R "Веган". Я хотел бы получить p-значения (или альтернативно t-значения) для анализа, чтобы оценить значимость моих предикторных переменных. Я могу генерировать p-значения, используя функцию anova, отдельно для каждой оси RDA (оси) и каждого предиктора (терминов):

model<-rda(precip ~ SOI + IOD + AAO, predictors)

anova_axis<-anova(model,by="axis")
anova_terms<-anova<-anova(model,by="term")


#### Outputs ####

Model: rda(formula = precip ~ SOI + IOD + AAO, data = predictors)
          Df Var    F         N.Perm Pr(>F)   
RDA1      1  640315 6.1841    199    0.005 **
RDA2      1  255421 2.4669    199    0.010 **
RDA3      1   75211 0.7264     99    0.710   
Residual  25 2588537  


Model: rda(formula = precip ~ SOI + IOD + AAO, data = predictors)
          Df    Var      F     N.Perm Pr(>F)   
SOI       1  435080 4.2020     99     0.01 **
IOD       1  419880 4.0552     99     0.01 **
SAM_AAO   1  115987 1.1202     99     0.23   
Residual 25 2588537                        

Однако, что мне действительно нужно, это p-значения для каждого предиктора, для каждой оси RDA. Итак, таблица примерно такая:

Predictor    RDA1    RDA2    RDA3
SOI          xxx     xxx     xxx
IOD          xxx     xxx     xxx
AAO          xxx     xxx     xxx

Кто-нибудь знает как и возможно ли это получить? Любая помощь приветствуется!