Тест отношения правдоподобия в R, возвращающий странные значения при оценке линейной модели смешанного эффекта - PullRequest
Новая
       33

Тест отношения правдоподобия в R, возвращающий странные значения при оценке линейной модели смешанного эффекта

0 голосов
/ 15 апреля 2019

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

model <- lmer(
  DV ~ I(predictor^2) + predictor + (predictor | rand_effect),
  data = data, 
  REML=FALSE
)

, а сравнение: 

modelcomp <- lmer(
  DV ~ (predictor | rand_effect),
  data = data, 
  REML=FALSE
)

Когда я использую тест отношения правдоподобия с использованием функции anova:

lrt <- anova(model, modelcomp)

Я получаю очень высокий хи-квадрат (57779,45) и более относительное значение p ровно 0.У меня сложилось впечатление, что значения p обычно не могут быть ровно 0. Почему это может быть?Как мне исправить мой код, чтобы он вычислял разумное значение p?

...