我试图使用exactLRT来查看我是否应该在模型中包含随机拦截。 以下是一个可重复的例子:
库(11聚物)
m.intercept< - lmer(反应〜天+(1 |主题),数据= sleepstudy)
m0< - lm(反应〜天,数据=睡眠研究)
exactLRT(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位)
区域设置:[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 utils数据集方法基础
其他附件包:[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-1 survival_2.37-7 lattice_0.20-29 [19] foreign_0.8-62 lmerTest_2.0-6 lme4_1.1-6 Rcpp_0.11.2
Matrix_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
答案 0 :(得分:0)
归功于@alexforrence指出我正确的方向。如果加载了lmerTest
,上面的代码就不起作用(我猜exactLRT
不接受类merModLmerTest
的对象,但只接受类merMod
的对象或者习惯了运行混合效应模型。
因此以下工作并允许加载lmerTest
:
库(11聚物)
库(RLRsim)
数据(sleepstudy)
m.intercept< - lme4 :: lmer(反应〜天+(1 |主题),数据= sleepstudy)
m0< - lm(反应〜天,数据=睡眠研究)
exactLRT(m = m.intercept,m0 = m0)
很抱歉,如果回答我自己的问题不在本论坛的规则中,我就不想让它暂停,因为我已经明白了。