Выбор модели с помощью твидовых распределений в пакете MuMIn R

Question

Я пытаюсь сделать выбор модели AICc и усреднение модели с помощью распределенных данных Твиди (составной Пуассон) в R.

Я безуспешно работал с пакетом 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 несовместимо с этим пакетом.

Спасибо за ваше время.

Answer 1

rank = tweedie::AICtweedie

Просто для тех, кто сталкивается с этим