Question

我正在尝试对使用lmer获得的混合效果模型进行简单的斜率分析。

该模型类似于以下内容：

data(Orthodont,package="nlme")
mod <- lme4::lmer(distance ~ age*Sex + (1|Subject), data=Orthodont)

当尝试使用软件包sim_slopes中的函数interactions进行简单的斜率分析时，我得到以下错误消息。

interactions::sim_slopes(model=mod, pred=age, modx=Sex)

错误：所请求的列之一不存在。
  回溯：
  1.互动：: sim_slopes（model = mod，pred = age，modx = Sex）
  5. jtools :::: summ.merMod（...）
  6. jtools ::: create_table（...）
  另外：警告消息：
  Johnson-Neyman间隔不适用于因子主持人。

有人可能对这个问题有任何想法吗？

请注意，它与软件包reghelper一起很好地工作，但是这些测试我需要p值，reghelper没有提供。

reghelper::simple_slopes(mod)

sessionInfo()的输出是：

R版本3.6.0（2019-04-26）   平台：x86_64-pc-linux-gnu（64位）在以下环境下运行：Ubuntu 18.04.2   LTS

矩阵产品：默认BLAS：
  /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.7.1 LAPACK：   /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.7.1

语言环境：[1] LC_CTYPE = en_CA.UTF-8 LC_NUMERIC = C
  LC_TIME = zh_CN.UTF-8 [4] LC_COLLATE = zh_CA.UTF-8
  LC_MONETARY = zh_CA.UTF-8 LC_MESSAGES = zh_CA.UTF-8 [7]   LC_PAPER = zh_CA.UTF-8 LC_NAME = C LC_ADDRESS = C
  [10] LC_TELEPHONE = C LC_MEASUREMENT = zh_CA.UTF-8   LC_IDENTIFICATION = C

附加的基本软件包：[1]统计图形grDevices utils
  数据集方法基础

通过名称空间（未附加）加载：
  [1] Rcpp_1.0.1 magrittr_1.5花键_3.6.0 MASS_7.3-51.1
  [5] munsell_0.5.0 colorspace_1.4-1格子_0.20-38 rlang_0.3.4
  [9] minqa_1.2.4 plyr_1.8.4 tools_3.6.0 grid_3.6.0
  [13] gtable_0.3.0 nlme_3.1-140 cli_1.1.0断言那个_0.2.1
  [17]摘要_0.6.19 lme4_1.1-21 lazyeval_0.2.2 tibble_2.1.2
  [21] crayon_1.3.4 Matrix_1.2-17 reghelper_0.3.4 nloptr_1.2.1

[25] ggplot2_3.1.1互动_1.1.0 jtools_2.0.1 pander_0.6.3
  [29]编译器_3.6.0支柱_1.4.1缩放_1.0.0引导_1.3-20
  [33] pkgconfig_2.0.2

Answer 1

（不是答案，但评论太久了。以后将酌情删除。）

在干净的R会话中，这对我有用：

data(Orthodont,package="nlme")
mod <- lme4::lmer(distance ~ age*Sex + (1|Subject), data=Orthodont)
interactions::sim_slopes(model=mod, pred=age, modx=Sex)
sessionInfo()

简单边坡分析
  性别=女性时的年龄斜率：
   美东时间。 S.E.吨价p
  ------ ------ -------- ------
   0.48 0.09 5.13 0.00
  性别=男性时的年龄斜率：
   美东时间。 S.E.吨价p
  ------ ------ -------- ------
    0.78 0.08 10.12 0.00
  警告信息：   Johnson-Neyman间隔不适用于因子主持人。

这是我的sessionInfo()：重要内容与您的内容相同（lme4 1.1-21，interactions 1.1.0，jtools 2.0.1），但是当然不完全相同... ??

R Under development (unstable) (2019-05-21 r76566)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.6 LTS

[matrix product and locale info deleted]   

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

loaded via a namespace (and not attached):
 [1] Rcpp_1.0.1         magrittr_1.5       splines_3.7.0      MASS_7.3-51.4     
 [5] tidyselect_0.2.5   munsell_0.5.0      colorspace_1.4-1   lattice_0.20-38   
 [9] R6_2.4.0           rlang_0.3.4        minqa_1.2.4        plyr_1.8.4        
[13] dplyr_0.8.1        tools_3.7.0        grid_3.7.0         gtable_0.3.0      
[17] nlme_3.1-140       cli_1.1.0          digest_0.6.19      assertthat_0.2.1  
[21] lme4_1.1-21        lazyeval_0.2.2     tibble_2.1.2       numDeriv_2016.8-1 
[25] crayon_1.3.4       Matrix_1.2-17      purrr_0.3.2        nloptr_1.2.1      
[29] ggplot2_3.1.1      jtools_2.0.1       interactions_1.1.0 glue_1.3.1        
[33] pander_0.6.3       compiler_3.7.0     pillar_1.4.1       scales_1.0.0      
[37] lmerTest_3.1-0     boot_1.3-22        pkgconfig_2.0.2

Answer 2

我通过拆下lmerTest解决了这个问题。

具有R包交互（jtools）的merMod对象的简单斜率测试

2 个答案: