Функция absoluteRLRT () `в` RLRsim` не может распознать объект `lmerMod` - PullRequest
Новая
       6

Функция absoluteRLRT () `в` RLRsim` не может распознать объект `lmerMod`

1 голос
/ 02 февраля 2020

У меня проблемы с получением функции exactRLRT() в пакете RLRsim. 

library(faraway)
library(RLRsim)
data(eggs)
cmods <- lmer(Fat ~ 1 + (1|Lab:Technician:Sample), data=eggs)
exactRLRT(cmods)

Это приводит к ошибке 

Error in exactRLRT(cmods) : Invalid <m> specified.

В справке для RLRsim::exactRLRT он говорит 

The fitted model under the alternative or, for testing in models with multiple variance components, the reduced model containing only the random effect to be tested (see Details), an lme, lmerMod or spm object

Но cmods был lmerMod объектом, так почему же это не работает?

1 Ответ

1 голос
/ 02 февраля 2020

у меня работает, я выкладываю здесь, чтобы вы могли проверить и посмотреть, что вызывает вашу ошибку: 

library(faraway)
library(RLRsim)
library(lme4)
data(eggs)
cmods <- lmer(Fat ~ 1 + (1|Lab:Technician:Sample), data=eggs)
exactRLRT(cmods)


    simulated finite sample distribution of RLRT.

    (p-value based on 10000 simulated values)

data:  
RLRT = 14.382, p-value < 2.2e-16

Мой sessionInfo (): 

R version 3.6.1 (2019-07-05)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: OS X El Capitan 10.11.6

Matrix products: default
BLAS:   /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRlapack.dylib

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] lme4_1.1-21   Matrix_1.2-17 RLRsim_3.1-3  faraway_1.0.7

loaded via a namespace (and not attached):
 [1] minqa_1.2.4     MASS_7.3-51.4   compiler_3.6.1  tools_3.6.1    
 [5] mgcv_1.8-28     Rcpp_1.0.2      splines_3.6.1   nlme_3.1-140   
 [9] grid_3.6.1      nloptr_1.2.1    boot_1.3-23     lattice_0.20-38
...