带有exactLRT的错误消息,用于测试lmer模型中随机截距的适用性

时间:2015-01-27 15:40:25

标签: r random lmer

我试图使用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

1 个答案:

答案 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)

很抱歉,如果回答我自己的问题不在本论坛的规则中,我就不想让它暂停,因为我已经明白了。