2

Rでtweedie(複合ポアソン)分散データを使用して、AICcモデルの選択とモデルの平均化を試みています。

私は AICcmodavg R パッケージで作業していましたが成功しませんでしたが、ここで提案に出くわしたときに MuMIn パッケージを試すことにしました ( https://stats.stackexchange.com/questions/141806/glm-model-selection-using- aicc-with-tweedie-distribution ) その

「ランク引数として指定するだけで、AICtweedie を MuMIN の関数で直接使用できます。」

モデルを次のように設定しました。応答変数 (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 ファミリーがこのパッケージと互換性がないのか疑問に思っています。

お時間をいただきありがとうございます。

4

1 に答える 1