评估线性混合效应模型时，R中的似然比测试返回奇异值

Question

我已经创建了一些数据的模型，我相信这些数据会沿着二次曲线下降，我将完整模型写为：

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

，比较方式为：

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

当我使用方差分析功能进行似然比检验时：

lrt <- anova(model, modelcomp)

我得到一个非常高的卡方（57779.45），而更令人关注的p值恰好为0。我的印象是p值通常不能恰好为0。为什么会这样呢？我该如何纠正我的代码，以便计算出合理的p值？