MuMIn软件包R中基于tweedie分布的模型选择

Question

我试图在R中用tweedie (复合泊松)分布数据进行AICc模型选择和模型平均。

我正在使用AICcmodavg R包，但没有成功，然后当我发现这里的建议(https://stats.stackexchange.com/questions/141806/glm-model-selection-using-aicc-with-tweedie-distribution)时，我决定尝试MuMIn包。

“您可以在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家族与这个包不兼容。

谢谢您抽时间见我。

Answer 1

职级= tweedie::AICtweedie

只是为了任何遇到这种情况的人