我正在尝试使用 R 中的 tweedie（复合泊松）分布式数据进行 AICc 模型选择和模型平均。
我正在使用 AICcmodavg R 包但没有成功，然后当我在这里遇到建议时决定尝试 MuMIn 包（https://stats.stackexchange.com/questions/141806/glm-model-selection-using- aicc-with-tweedie-distribution）
“您可以直接在 MuMIn 的函数中使用 AICtweedie，只需将其指定为排名参数即可。”
我按如下方式设置我的模型我的响应变量 (NVIR) 是东部蝾螈成虫的每单位努力捕获量，我的解释变量是我的采样点的各种栖息地特征。
m1<- glm(NVIR~Water_T+cond+DO+ORP+pH+max_depth+type,
family = tweedie(link.power=0, var.power=1.3), data = cpue)
m2<- glm(NVIR~Water_T+cond+DO+ORP+pH+littoral_slope+type,
family = tweedie(link.power=0, var.power=1.3), data = cpue)
m3<- glm(NVIR~pH+DO+cond+max_depth+type,
family = tweedie(link.power=0, var.power=1.3), data = cpue)
m4<- glm(NVIR~pH+DO+cond+littoral_slope+type,
family = tweedie(link.power=0, var.power=1.3), data = cpue)
m5<- glm(NVIR~cond+type+pH+max_depth,
family = tweedie(link.power=0, var.power=1.3), data = cpue)
然后尝试了这条线
model.sel(m1, m2, m3, m4, m5, rank = AICc, rank.args = AICtweedie)
并收到错误
Error in UseMethod("logLik") :
no applicable method for 'logLik' applied to an object of class "function"
In addition: Warning message:
In model.sel.default(m1, m2, m3, m4, m5, rank = AICc, rank.args = AICtweedie) :
models are not all fitted to the same data
或者，我也试过这条线
model.sel(m1,m2,m3,m4,m5, rank.args=AICtweedie)
并得到这个错误：
Error in get(x) : object 'Tweedie' not found
In addition: Warning message:
In model.sel.default(m1, m2, m3, m4, m5, rank.args = AICtweedie) :
models are not all fitted to the same data
我想知道问题是否出在我的代码上，或者 tweedie 系列是否与此软件包不兼容。
感谢您的时间。