R: permlmer，lmer，lme4，谓词均值中混合模型似然比检验的置换检验误差

Question

我想用置换检验来检验分类变量对似然比检验的主要影响。我有一个连续的结果和一个二分法的分组预测器和一个分类时间预测器(日，5级)。

数据通过这个驱动器链接以rda格式暂时可用。

library(lme4)
lmer1 <- lmer(outcome ~ Group*Day + (1 | ID), data = data, REML = F, na.action=na.exclude)
lmer2 <- lmer(outcome ~ Group + (1 | ID), data = data, REML = F, na.action=na.exclude)

library(predictmeans)
permlmer(lmer2,lmer1)

但是，这段代码给出了以下错误：

Density.default中的错误(c(lrtest1，lrtest)，内核= "epanechnikov")：至少需要两个点才能自动选择带宽

下面的代码确实有效，但实际上并没有给出我相信的经过变异的LR测试的结果：

library(nlme)
lme1 <- lme(outcome ~ Genotype*Day,
                  random = ~1 | ID,
                  data = data,
                  na.action = na.exclude)

library(pgirmess)
PermTest(lme1)

谁能指出为什么我在使用permlmer 函数？时会出现"epanechnikov“错误？

谢谢!

Answer 1

问题是NANs，从数据集中删除所有nans并重新运行模型。我也有同样的问题，这解决了它。