疏导功能错误-R package MuMln

Question

我必须对数据集进行统计分析。我想创建所有可能的模型，并使用疏导功能测试它们，但它不起作用。实际上，当我输入时：

glm1<-glm(presabs~dca1+dca2+se1+se2, family=binomial(logit))
dredge(glm1)

我得到了这个错误：

Erreur in dredge(glm1) : 
'global.model''s 'na.action' argument is not set and options('na.action') is "na.omit"

有人能帮我吗？

Answer 1

使用options(na.action = "na.fail")的问题是它改变了R的全局设置。如果你有一个大的脚本，改变全局设置可能会影响你的代码的其他部分，你隐含地依赖于R的默认设置。有两种方法可以避免这种情况：

通过options(na.action = "na.omit").将设置更改回使用dredge后的

或者更好的方式...

1. 利用回归函数的功能来设置na.action参数。在您的案例中：

glm1 <- glm(presabs ~ dca1+dca2+se1+se2,
            family=binomial(logit),
            na.action = "na.fail")

Answer 2

请参阅?dredge

# Example from Burnham and Anderson (2002), page 100:
data(Cement)
options(na.action = "na.fail")   #  prevent fitting models to different datasets

fm1 <- lm(y ~ ., data = Cement)
dd <- dredge(fm1)

如果您跳过第二行，则会弹出所描述的错误，因为这些模型适合不同的数据集(由于删除了NAs)。