在 R 中使用 Ordinal ModV 运行调节中介分析答案

【问题标题】：Running a Moderated Mediation analysis with Ordinal ModV in R在 R 中使用 Ordinal ModV 运行调节中介分析
【发布时间】：2020-09-13 09:14:19
【问题描述】：

我从一项调查中收集了数据，我想用这些数据建立一个有调节的中介模型。我的四个变量是：

Bullying = 个人对多个李克特量表问题的持续平均回答，范围 1-5 (IV)
Self-Esteem = 连续，如上计算（中介）
Dissatisfaction = 连续，如上计算 (DV)
Year Level = 离散序数，范围为 1-10（主持人）

我已经计算过，发现Self-Esteem是效果Bullying -> Dissatisfaction的简单Mediator。我现在想看看Year level 是否是此中介模型的主持人，但是我能找到的唯一指南依赖于将主持人分成两个二分组，我不想这样做。我想我通过代码发现Year Level并不是Bullying -> Dissatisfaction效果的简单版主：

Bullying_c <- c(scale(Bullying, center = TRUE, scale = FALSE))
SelfEsteem_c <- c(scale(SelfEsteem, center = TRUE, scale = FALSE))

fitMod <- lm(Dissatisfaction ~ Bullying_c + SelfEsteem_c + Bullying_c*SelfEsteem_c)
fitModB <- Boot(fitMod, R = 1000)
summary (fitModB)

哪些输出：

Residuals:
    Min      1Q  Median      3Q     Max 
-1.5143 -0.6560 -0.2014  0.5426  3.0809 

Coefficients: 
                         Estimate Std. Error t value Pr(>|t|)    
(Intercept)              1.853848   0.032137  57.686  < 2e-16 ***
Bullying_c               0.203157   0.027922   7.276  8.5e-13 ***
SelfEsteem_c             -0.036391   0.015824  -2.300   0.0217 *  
Bullying_c:SelfEsteem_c   0.007795   0.012430   0.627   0.5308

这表明Self-esteem 不是一个整体的调解人，因此我正在查看调解调解，而不是调解调解。

我已尝试查找包：psych、mediation 和 lavaan，但我一直无法找到使用序号主持人运行分析的方法。大多数指南希望我选择主持人的两个值，但我想包括所有 10 个类。

【问题讨论】：

标签： r statistics

【解决方案1】：

通过不断的研究，我偶然发现了rosetta 包，其中包含函数gemm()，它可以使用以下代码创建此模型：

result <- gemm(data = mydata,
xvar = "Bullying_c",
mvar = "SelfEsteem_c",
yvar = "Dissatisfaction",
xmmod = "YearLevel",
nboot = 5000)

print(result)

这个函数也允许协变量。 gemm() 函数在下面几页中有清楚的解释：

函数详情：
https://rdrr.io/github/psytext/rosetta/man/gemm.html

运行这些分析和解释结果的教程的 PDF 下载链接： https://psyarxiv.com/mj2ug/download

【讨论】：