set.seed(123)
n <- 200
X <- rnorm(n)
M1 <- 0.5 * X + rnorm(n)
M2 <- 0.3 * X + rnorm(n)
M3 <- 0.1 * X + rnorm(n)
Y <- 0.3 * X + 0.4 * M1 + 0.2 * M2 + rnorm(n)
data <- data.frame(X, Y)
mediators <- data.frame(M1, M2, M3)
if (requireNamespace("mediation")) {
results <- batch_mediate(data, mediators, nsims = 99)
print(results)
}
Run the code above in your browser using DataLab