Плохая совместимость MuMin с мультиномными (nnet) объектами, имеющими 3-х уровневые переменные ответа?

Question

Я сталкиваюсь с проблемой при попытке использовать функцию model.avg из MuMin пакет на объектах, полученных в результате функции multinom (пакет nnet). Хотя ручной список как multinom совместим, model.avg Функция не возвращает средние коэффициенты модели для объясняющих переменных, как я ожидаю.

Я знаю, что проблема не существует, когда multinom Переменная reponse только в 2 уровня, так что, я думаю, это зависит от структуры результирующего объекта с множеством объектов. Однако я не знаю, как решить эту проблему.

Вот пример кода, где функция model.avg работает как положено:

#### MuMIn example  
library(MuMIn)  
data(Cement)  
fm1 <- lm(y ~ ., data = Cement)  
dd <- dredge(fm1)  
summary(model.avg(dd, subset = delta < 4))$coefmat

Результаты:

               Estimate Std. Error Adjusted SE   z value  Pr(>|z|)
(Intercept) 64.69312754 22.2354794  22.4624137 2.8800613 0.0039760
X1           1.45579791  0.2036676   0.2193043 6.6382554 0.0000000
X2           0.62502598  0.1202570   0.1291705 4.8387682 0.0000013
X4          -0.47600710  0.2215167   0.2309395 2.0611767 0.0392862
X3          -0.02153196  0.3767120   0.3924356 0.0548675 0.9562440

В этом случае функция model.avg возвращает Model-averaged coefficients (колонка Estimate) для объясняющих переменных (то есть X1-X4), как и ожидалось.

Теперь, если я хочу применить это к multinom объект (переменная ответа минимум с 3 уровнями):

#### application to multinomial regression
library(foreign) ; library(nnet) ; library(MuMIn)
ml <- read.dta("http://www.ats.ucla.edu/stat/data/hsbdemo.dta")
ml$prog2 <- relevel(ml$prog, ref = "academic")
test <- multinom(prog2 ~ ses + write + schtyp, data = ml)
# [...] verbosity not displayed here
dd <- dredge(test)
# [...] verbosity not displayed here
summary(model.avg(dd, subset = delta < 4))$coefmat

Результаты

         Estimate Std. Error  z value    Pr(>|z|)
general  2.850316   1.169132 2.437977 0.014769732
vocation 5.134921   1.173878 4.374322 0.000012181

Здесь функция возвращает "средние коэффициенты модели" для разных уровней переменной отклика (т.е. general а также vocation), а не для модельных средних коэффициентов для объясняющих переменных.

Если это возможно, подскажите пожалуйста, как получить Model average coefficients для пояснительных переменных из multinom объекты?

Спасибо заранее.

Answer 1

Это действительно была ошибка с извлечением коэффициентов из мультиномных объектов. Это исправлено в MuMIn 1.9.13.