Предупреждающее сообщение: модели не все приспособлены к одинаковому количеству наблюдений - PullRequest
Новая
       70

Предупреждающее сообщение: модели не все приспособлены к одинаковому количеству наблюдений

1 голос
/ 19 апреля 2019

Я подгоняю несколько моделей с другой структурой. Пример каждой структуры: 

> m1 <- glmer(X~Y1+Y2+Y3+(1|B/C),control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = dataset, family=binomial(link="logit"))
> 
> m2 <- glmer(X~Y1*A+Y2*A+Y3*A+(1|B/C),control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = dataset, family=binomial(link="logit"))
> 
> m3 <- glmer(X~(1+Y1|B)+Y2+Y3,control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = dataset, family=binomial(link="logit"))

Я использую один и тот же набор данных для всех них и хочу сравнить значения AIC, используя формулу AIC() в пакете stats. После запуска AIC(m1, m2, m3) я получаю следующее предупреждение: 

Warning message:
In AIC.default(m1, m2, m3) :
  models are not all fitted to the same number of observations

Я проверил количество наблюдений для всех них, и они действительно разные: 

> nobs(m1)
[1] 157278

> nobs(m2)
[1] 150164

> nobs(m3)
[1] 157278

Но я не понимаю, почему это значение отличается, поскольку я использую один и тот же набор данных (я запускал его дважды с самого начала, а также модели, чтобы убедиться, что он одинаков для всех моделей).

Может ли это быть связано с формулой модели?

...