Нулевая инфляция NB - объект 'model_count' не найден - PullRequest
Новая
       47

Нулевая инфляция NB - объект 'model_count' не найден

1 голос
/ 22 апреля 2020

Я моделирую опосредованную нулевую модель с отрицательным биномиальным индексом (ZINB). Я следую инструкциям O'Rourke & Vazquez (2019) -> https://www.sciencedirect.com/science/article/abs/pii/S0306460319301078

Несколько дней go я управлял другой моделью ZINB, которая работала совершенно нормально , Однако сегодня я написал другую модель, те же переменные, но почему-то она больше не работает и выдает странную ошибку. Когда я пробую свою предыдущую модель, я внезапно получаю ту же ошибку, а именно: 

Error in zeroinfl(Y1 ~ X1 + M1 | X1 +  : object 'model_count' not found

Остальная часть моего кода: 

#loading required packages
library(psych)
library(foreign)
library(ggplot2)
library(MASS)
library(pscl)
library(nonnest2)
library(lmtest)
library(boot)

#Import data
mydata

#Fit a ZINB model
ex1zinb <- zeroinfl(X1 ~ Y1 + M1 | Y1 + M1, data = mydata, dist="negbin", EM= TRUE)

Я уже проверил допущение для статистического анализа et c., И особенно странно, что аналогичная модель работала несколько дней go и больше. Я попытался установить несколько пакетов сегодня, но столкнулся с ненулевым состоянием выхода. После некоторых комментариев к stackoverflow я установил пакет, добавив зависимости = TRUE, но он запустился. После этого начались проблемы. Может быть, что-то не так с моими пакетами?

При загрузке библиотек я получаю следующие сообщения: 

> library(psych)
> library(foreign)
> library(ggplot2)
Need help? Try Stackoverflow: https://stackoverflow.com/tags/ggplot2

Attaching package: ‘ggplot2’

The following objects are masked from ‘package:psych’:

    %+%, alpha

> library(MASS)
> library(pscl)
Classes and Methods for R developed in the
Political Science Computational Laboratory
Department of Political Science
Stanford University
Simon Jackman
hurdle and zeroinfl functions by Achim Zeileis
> library(nonnest2)
This is nonnest2 0.5-3.
nonnest2 has not been tested with all combinations of model classes.
> library(lmtest)
Loading required package: zoo

Attaching package: ‘zoo’

The following objects are masked from ‘package:base’:

    as.Date, as.Date.numeric

> library(boot)

Attaching package: ‘boot’

The following object is masked from ‘package:psych’:

    logit

У кого-нибудь есть предложения или идеи? Заранее спасибо!

Редактировать: я получаю нормальный вывод при работе ненулевой раздутой отрицательной биномиальной модели. Следующий код работает гладко: 

summary(ex1nb <- glm.nb(Y1~ X1 + M1, data = mydata))

Так что я думаю, что мои данные в порядке?

...