R语言 LM型号的穿刺针交互-错误类别

fdbelqdn  于 2022-12-20  发布在  其他
关注(0)|答案(1)|浏览(174)

我是R的新手,一直在尝试获得一个包含一个因子(2个水平)和一个连续变量的双向交互作用的简单斜率。尝试使用R中交互作用包中的probe_interactions,会出现以下关于该类的错误。Sim_slopes也会出现相同的错误。我知道如何在Stata/SPSS上执行此操作,但我需要在R上运行它。
有人有什么信息可以解决这个问题吗?
提前感谢您的帮助!

> model1 <- lm(y ~  dummy * cont, data = data_new)
> probe_interaction(model = model1, pred = cont, modx = dummy, data = data_new)
Error:
! `class` must be a character vector, not a list.
Backtrace:
  1. interactions::probe_interaction(...)
  3. interactions::sim_slopes(...)
  5. jtools:::summ.lm(...)
  6. jtools:::do_robust(model, robust, cluster, data, vcov)
  7. jtools::get_robust_se(...)
  8. jtools::stop_wrap(...)
 10. rlang (local) `<fn>`(class = <named list>)

sim_斜率(

  • 模型=模型1,
  • 前药=对照品,
  • modx =虚拟,
  • 数据=新数据)错误:!class必须是字符向量,而不是列表。反向跟踪:

1.交互作用::sim_斜率(...)

  1. jtools:::总和(...)
  2. jtools:::do_robust(模型、健壮、群集、数据、vcov)
  3. jtools::获取健壮的代码(...)
  4. jtools::停止 Package (...)
  5. rlang(本地)<fn>(类=)
yvt65v4c

yvt65v4c1#

您可能想尝试the marginaleffects package(免责声明:我是作者)。它可以做到这一点,甚至更多。该网站包括一个简短的入门指南和数千字的教程:
https://vincentarelbundock.github.io/marginaleffects/
下面是一个非常简单的可复制的例子,其中涉及一个二元变量和一个连续变量之间的交互:

library(marginaleffects)

mod <- lm(mpg ~ am * hp, data = mtcars)

mfx <- marginaleffects(mod)
summary(mfx)
#>   Term   Effect Std. Error z value   Pr(>|z|)    2.5 %   97.5 %
#> 1   am  5.27675   1.098721   4.803 1.5660e-06  3.12330  7.43021
#> 2   hp -0.05897   0.008726  -6.759 1.3922e-11 -0.07608 -0.04187
#> 
#> Model type:  lm 
#> Prediction type:  response

plot_cap(mod, condition = c("hp", "am"))

相关问题