计算lm4 / lmerTest中空随机效应模型的p值

时间:2018-01-23 09:18:14

标签: r lme4 mixed-models lmertest

Edit1:正如@RolandASc报道的那样,CultureInfo.CurrentCulture.DateTimeFormat.ShortDatePattern中似乎有一个错误。我已经向报告该问题的软件包的维护者写了一封电子邮件。
Edit2:维护者响应:"我们正在进行更新,希望这些问题能够以更好的方式处理......"

我试图使用lmerTest / lme4获取空随机效果模型的p.values,但无法弄清楚为什么不为{{1}计算它们模型。

使用lmerTest数据我按如下方式定义模型:

null

我希望调用sleepstudy会为固定效果中的拦截打印一个p.value - 但library(lmerTest) lmer0 <- lmer(Reaction ~ 1 + (1|Subject), data = sleepstudy) 无法执行此操作并实际调用{{1}的摘要}:

summary(lmer0)

如果我在同一模型中使用lmerTest,那么一切看起来都是正确的:

lme4

P.S。如果这个问题更适合CrossValidated,请告诉我,我会把它移到那里。谢谢。

2 个答案:

答案 0 :(得分:2)

作为lmerTest开发人员之一,我可以确认这确实是一个错误。

它也可以在GitHub(https://github.com/runehaubo/lmerTest)上的开发版本中修复,你可以用它安装

library("devtools")
install_github("runehaubo/lmerTest")

请在GitHub上提出潜在的未来错误。

溴 符文

答案 1 :(得分:1)

这是lmerTest:::calcSummary中的错误。

预计resultmatrix,但在您的情况下会减少为vector

看起来像一个相当简单的修复。

相关问题