时间戳

首页 > 解决方案 > 在 R 中使用 Ordinal ModV 运行调节中介分析

r - 在 R 中使用 Ordinal ModV 运行调节中介分析

问题描述

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

我已经计算过了，发现Self-Esteem是一个简单的 Mediator 的效果Bullying -> Dissatisfaction。我现在想看看是否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 个类。

标签： rstatistics

解决方案

通过不断的研究，我偶然发现了rosetta包含 function 的包，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

推荐阅读