r - 如何在 R 中使用 OpenMx 或 lavaan 包构建这样的 SEM 框架

Question

嗨，我是临床科学家，对 R 和编码相对较新。在研究医学文献时，我经常遇到 SEM 框架图，它在一个有趣的图中解释了多元回归或相关性。我也在努力制定一个但无法解决。我如何从以下的总结（拟合）中制定这个图表。您也可以解释摘要（适合）吗？

library("lavaan.survey")
install.packages("lavaan.survey")

data("pisa.be.2003")

model.pisa <- "
    math =~ PV1MATH1 + PV1MATH2 + PV1MATH3 + PV1MATH4
    neg.efficacy =~ ST31Q01 + ST31Q02 + ST31Q03 + ST31Q04 + 
                    ST31Q05 + ST31Q06 + ST31Q07 + ST31Q08
    neg.selfconcept =~ ST32Q02 + ST32Q04 + ST32Q06 + ST32Q07 + ST32Q09

    neg.selfconcept ~ neg.efficacy + ESCS + male
    neg.efficacy ~ neg.selfconcept + school.type + ESCS + male
    math ~ neg.selfconcept + neg.efficacy + school.type + ESCS + male
"

fit <- lavaan(model.pisa, data = pisa.be.2003, auto.var = TRUE, std.lv = TRUE,
  meanstructure = TRUE, int.ov.free = TRUE, estimator = "MLM")
summary(fit)

![:ČukićI,BatesTC(2015)The Association between NeuroticismandHeartRateVariabilityIsNotfullyExplainedbyCardiovascularDisease and Depression.PLoSONE10(5):e0125882. doi:10.1371/journal.pone.0125882] 1

Answer 1

这看起来像您所描述的：

library(lavaanPlot)
lavaanPlot(model = fit, coefs = TRUE)

Answer 2

另一种选择是来自包 semPlot 的 semPaths：

library(semPlot)
semPaths(fit, "std")

链接中有几个情节选择