1

我正在尝试计算涉及三向交互的简单斜率（apoe4:time:sleep_eff_10）。见下文

library(nlme)
m1  <-  lme(EF ~ sex_recoded + race_recoded + educ_year_centered + BMI_masterfile_centered + apoe4 + sleep_meds + CVD + CESD_withoutslp + baseline_age_masterfile_centered + baseline_age_masterfile_centered:time + time + sleep_eff_10c + time:apoe4 +  sleep_eff_10c:time + apoe4:sleep_eff_10c + apoe4:time:sleep_eff_10c, random = ~1|idno, 
            data = df2, method = "REML", na.action = na.omit)
summary(m1)

我正在使用 reghelper 包来执行此操作（请参见下面的代码）。

simple_slopes(m1, levels=list(apoe4=c(0, 1, 'sstest'),
                              time=c(0, 1, 2, 'sstest'), 
                              sleep_eff_10c=c(0.02, -4.67, 1.81, 'sstest')))

当我这样做时，我收到以下错误和警告消息：

Error in summary(new_model)$tTable[dummy_names, ] : 
  subscript out of bounds
In addition: Warning message:
Using formula(x) is deprecated when x is a character vector of length > 1.
  Consider formula(paste(x, collapse = " ")) instead.

我相信我正在遵循包的要求，所以我不明白为什么我会收到这些消息。有没有人有任何建议让这个工作？

4

0 回答 0