c - R 的 hier.part() 取决于预测器顺序

Question

我尝试使用 hier.part R 包(它是 hier.part() 函数)来估计变量的独立贡献。突然之间，我对预测顺序产生了强烈的依赖。

有一篇论文 ( http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0011698 ) 报告了超过 9 个预测变量的类似行为。包的文档还报告了 9-12 个预测变量的“较小舍入误差”。但即使在具有 4 个自变量的包示例中，我也经历过这种情况。

> library(hier.part)
> data(amphipod)
> env1 <- amphipod[,2:5]
> hier.part(amphipod$australis, env1, fam = "binomial",
+           gof = "logLik", barplot=F)$I.perc
              I
fimp   19.74685
fconn  18.24679
densep 25.24765
unseal 36.75871
> hier.part(amphipod$australis, env1[,c(4,3,2,1)], fam = "binomial",
+           gof = "logLik", barplot=F)$I.perc
              I
unseal 16.16072
densep 17.98755
fconn  32.00246
fimp   33.84927

此外，调用之间的确切值并不稳定。该函数使用C代码(在partition()调用中)，并且至少“调用之间的稳定性”问题可以归因于它，因为所有相关的R代码都是确定性的。

Answer 1

这次你的问题在我的R中没有显示出来，两个I$perc是一样的。-晨曦