解释疏浚的输出

机器算法验证 r 广义线性模型 解释 生态
2022-03-17 17:14:47

我第一次在 R 中运行 GLM,但不知道如何解释结果。这是输入:

global.modelAcar <- lm(Acar ~ logNutrientsc*logNDSc*logNNNc, data = dat)
summary(global.modelAcar)
options(na.action=na.fail)
MAcar <- dredge(global.modelAcar)
MAcar

这是结果的前几行:

Global model call: lm(formula = Acar ~ logNutrientsc * logNDSc * logNNNc, data = dat)
---
Model selection table 
        (Int)     lND      lNN      lgN  ... lND:lNN:lgN df  logLik AICc delta weight
2   2.159e-17 -0.2590                    ...              3 -26.445 59.6  0.00  0.214
3   1.682e-17         -0.25420           ...              3 -26.497 59.7  0.10  0.203
1   7.778e-18                            ...              2 -27.799 59.9  0.36  0.179
5   1.220e-17                  -0.16580  ...              3 -27.256 61.2  1.62  0.095
6   2.249e-17 -0.2295          -0.09269  ...              4 -26.283 61.7  2.17  0.072

问题:我知道模型 2 是最好的模型,并且表明 lND 对多样性有负面影响。第二好的模型显示 lNN 具有负面影响。没有价值意味着没有效果。AIC 值表明这些模型的信息量不是很大。这种解释是正确的还是我遗漏了什么?

1个回答

该函数MuMIn::dredge只返回一个模型列表,其中包含预测变量的每种可能组合。至于你的结果,请允许我不同意你所说的:

我知道模型 2 是最好的模型,并且表明 lND 对多样性有负面影响。

部分正确,1ND 确实对多样性有负面影响,但从 delta (=delta AIC) 你无法将模型 2 与 3、1 和 5 区分开来,因为(使用通用拇指规则)它们的 dAIC < 2

没有价值意味着没有效果。

这是不正确的。dredge 返回一个包含所有可能的变量组合的列表,如果变量没有值,则表示它没有包含在模型中。比如模型3只有1NN,除了截距明显。

AIC 值表明这些模型的信息量不是很大。这种解释是正确的还是我遗漏了什么?

由于前四个模型具有相似的支持(还要注意它们的 Akaike 权重,从 0 到 1 不等,不是相对高)我强烈建议您使用模型平均,看看MuMIn::model.avg并阅读Burnham & Anderson的第 4 章(2002 年)我希望这已经足够清楚了,但请随时再次询问