基于R语言实现中介效应检验以及sobel检验代码

数据格式

随机数据，不一定好
y是因变量，x是自变量，m是中介变量

基本原理

M ~ X
Y ~ X
Y ~ X + M

直接上代码

library(mediation)
library(bda)

# 加载readxl包
library(readxl)  # 读取Excel表格


# 读取数据  我是从剪切板读取的
data = readClipboard(sep = '\t')

# 读取Excel文件的数据
excel_data <- readxl::read_excel("C:\\Users\\12810\\Desktop\\中介效应代码\\test.xlsx")

# 查看读取的数据
print(excel_data)

## sobel 检验
mediation.test(excel_data$m, excel_data$x, excel_data$y)



### 三步检验法

model1 <- lm(m~x, data = excel_data)
model2 <- lm(y~x, data = excel_data)
model3 <- lm(y~m+x, data = excel_data)
med.out <- mediation::mediate(model2, model3, treat ="x", mediator="m",
                              boot=TRUE, sims=500)
summary(med.out)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30