R：如何获得熔岩中介模型的诊断图？答案

【问题标题】：R: How to obtain diagnostic plots for a lavaan mediation model?R：如何获得熔岩中介模型的诊断图？
【发布时间】：2022-06-19 20:30:22
【问题描述】：

我不确定这是否更适合在这里询问或 CrossValidated，因为我特别询问使用 R / lavaan...

我不确定我是否完全误解了如何检查违反假设的情况。我知道我们可以通过以下方式获得线性模型的诊断图：

model <- lm(data$outcome ~ data$predictor)

plot(model, which = c(1:6))

但是我无法弄清楚我应该如何为这样的中介模型执行此操作：

model <- 'outcome ~ c*predictor + b*mediator
              mediator ~ a*predictor
              indirect_effect := a*b
              total_effect := c + (a*b)
              '

model.fit <- lavaan::sem(
  model = model,
  data = data,
  missing = "FIML",
  estimator = "ML")

然后，如果我尝试以相同的方式 (plot(model.fit, which = c(1:6))) 获取绘图，我会得到 Error in as.double(y) : cannot coerce type 'S4' to vector of type 'double'。

另外，为了检查是否违反了 Pearson 相关性的假设，我们是通过单独查看每个变量的结构，还是通过制作线性模型 (lm(data$outcome ~ data$predictor)) 或使用相关性本身 (cor.test(data$var1, data$var2) ) 以某种方式？

【问题讨论】：

标签： r r-lavaan

【解决方案1】：

试试：lavaanPlot::lavaanPlot(model = model.fit, coefs=T)

【讨论】：