0

我正在使用 R 包“vegan”进行冗余分析。我想获得用于分析的 p 值（或 t 值），以评估我的预测变量的重要性。我可以使用 anova 函数生成 p 值，分别为每个 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

但是，我真正需要的是每个 RDA 轴的每个预测变量的 p 值。所以像这样的表：

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

有谁知道如何以及是否有可能获得？任何帮助将非常感激！

4

0 回答 0