0

我正在使用 glmulti 来评估几个预测变量对计数的相对重要性。当我使用 glm.nb 时，我正在使用包装器。我想添加一个抵消项，但我不确定它应该去哪里，或者是否可能。我尝试将其直接添加到公式中，但当时不使用偏移量。我还尝试添加它类似于我如何使用随机效果来添加它，但随后偏移量只是作为另一个预测器出现。然后我尝试将其作为重量放入，尽管偏移量不是重量，但它给了我极高的 AIC 值（2000 与我运行的相同“常规”glm.nb 模型给出的 740 比较） .

请看下面的代码。我非常感谢您对此提出一些建议！谢谢！

  1. 版本：输入公式时忽略偏移量
glmnb.glmulti<-function(formula,data,...){
glm.nb(paste(deparse(formula)),data=data)
}
test<-glmulti(counts~a+b+c+offset(log(d)),data=mydata,level=1,fitfunction=glmnb.glmulti,crit="aic")
  1. 版本：当添加类似于随机效应时，偏移量成为另一个预测器
glmnb.glmulti<-function(formula,data,offset,...){
glm.nb(paste(deparse(formula),offset),data=data)
}
selection<-glmulti(counts~a+b+c,data=mydata,offset="+(log(d))",level=1,fitfunction=glmnb.glmulti,crit="aic")
4

0 回答 0