改变rda anova中的排列数

Question

我正在尝试使用anova（）函数测试rda模型的整体重要性，但我意识到，如果我多次运行anova，它不仅每次都会给出不同的p值，但它也给出了不同的排列数。

  

Rfpreds＆lt; -rda（Rorders~Cornyield + Respiration + Nmin + logNase + logBGase + logPase，data = Rfunctions，na.action =“na.omit”）

     

ANOVA（Rfpreds）

所以我想看看我是否会使用更多排列获得更一致的结果。我在使用R文档方面不是很有经验，但我的理解是我应该能够设置排列编号：

  

anova（Rfpreds，permutations = 999）

但这导致错误：

  

match.arg（model）中的错误：'arg'必须为NULL或字符向量

我无法弄清楚我在这里做错了什么，但如果每次运行代码时它都不同，我当然无法报告p值。

Answer 1

可能你有一个旧版本的素食主义者。可能比2.2-0还要大。

在素食主义者 2.2-0之前，我们使用了一种惰性策略，一旦我们确定排列 p -value低于临界阈值（通常是 p = 0.05），但在2.2-0中，我们总是使用相同数量的排列，如参数中所示。

您有两种选择：（1）升级纯素并使用参数permutations或（2）将参数step和perm.max设置为相同给出所需排列数的值。