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

Question

Hi well i am clinical scientist and relatively new to R and coding. 在研究医学文献的时候，我经常会遇到SEM框架工作图，用一张图就能解释多变量回归或相关性，这很有意思，可以理解。我也在努力制定一个但无法整理出来。我怎么能从下面的summary(fit)制定这个图。您也能解释摘要(fit)吗？

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)TheAssociation betweenNeuroticismandHeartRateVariabilityIsNotFullExplainedbyCardiovascularDiseaseandDepression.PLoSONE10(5):e0125882.doi:10.1371journal.pone.0125882] 。1

Answer 1

这看起来像你所描述的。

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

Answer 2

另一种选择是包semPlot中的semPaths。

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

有几个图的替代品在 联系