在R中使用Ordinal ModV运行调解分析

如何解决在R中使用Ordinal ModV运行调解分析

我从一项调查中收集了一些数据，希望以此来建立一个仲裁中介模型。我的四个变量是：

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

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

Bullying_c <- c(scale(Bullying,center = TRUE,scale = FALSE))
SelfEsteem_c <- c(scale(SelfEsteem,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个类。

解决方法

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