0

我正在尝试使用exactLRT 来查看我是否应该在我的模型中包含一个随机截距。下面是一个可重现的例子：

图书馆（lmer）

m.intercept <- lmer(Reaction ~ Days+(1|Subject), data=sleepstudy)

m0 <- lm（反应〜天，数据=睡眠研究）

精确轻轨(m=m.intercept, m0=m0)

我收到以下错误消息：“exactLRT 中的错误（m = m.intercept，m0 = m0）：模型中的多个随机效应 - exactLRT 需要 'm' 且只有一个随机效应。”

但是，据我了解，我在 m.intercept 中包含了一个随机效应（即随机截距）。

关于问题的任何想法？

如果有用：

R 版本 3.1.2 (2014-10-31) 平台：i686-pc-linux-gnu (32-bit)

语言环境：[1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C
LC_TIME=en_CA.UTF-8 LC_COLLATE=en_CA.UTF-8 [5] LC_MONETARY=en_CA.UTF-8 LC_MESSAGES=en_CA.UTF-8
LC_PAPER=en_CA。 UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C
LC_MEASUREMENT=en_CA.UTF-8 LC_IDENTIFICATION=C

附加的基础包：[1] 样条线网格统计图形 grDevices 实用程序数据集方法基础

其他附加软件包：[1] RLRsim_3.0 MuMIn_1.10.0
psych_1.4.5 WriteXLS_3.5.0 reshape2_1.2.2 languageR_1.4.1 [7] e1071_1.6-3 HLMdiag_0.2.5 gdata_2.13.3 psy_1.1
car_2.0-19 ggplot2_0.9.3 .1 [13] date_1.2-34 plyr_1.8.1
Hmisc_3.14-4 Formula_1.1-1survival_2.37-7 lattice_0.20-29 [19] foreign_0.8-62 lmerTest_2.0-6 lme4_1.1- 6 Rcpp_0.11.2
矩阵_1.1-5

通过命名空间加载（未附加）：[1] bitops_1.0-6
caTools_1.17 class_7.3-11 cluster_1.15.3
colorspace_1.2-4 [6] dichromat_2.0-0 digest_0.6.4
gplots_2.13.0 gtable_0。 1.2 gtools_3.4.0 [11] KernSmooth_2.23-13 labeling_0.2 latticeExtra_0.6-24 MASS_7.3-37 minqa_1.2.3 [16] munsell_0.4.2
nlme_3.1-119 nnet_7.3-8 numDeriv_2012.9-1
pbkrtest_0 .3-8 [21] proto_0.3-10 RColorBrewer_1.0-5 RcppEigen_0.3.2.1.2 scales_0.2.3 stats4_3.1.2 [26] stringr_0.6.2 tools_3.1.2

4

1 回答 1

0

感谢@alexforrence 为我指明了正确的方向。lmerTest如果加载（我猜exactLRT不接受 class 的对象merModLmerTest，但只接受 class 的对象merMod）或用于运行混合效果模型，则上面的代码不起作用。

所以以下工作并允许lmerTest加载：

图书馆（lmer）

库（RLRsim）

数据（睡眠研究）

m.intercept <- lme4::lmer（反应〜天+（1|主题），数据=睡眠研究）

m0 <- lm（反应〜天，数据=睡眠研究）

精确轻轨(m=m.intercept, m0=m0)

抱歉，如果回答我自己的问题不在这个论坛的规则之内，我只是不想让它挂在我想通之后。

于 2015-01-27T21:08:01.510 回答