用R进行结构方程模型分析

脚本与注释

library(lavaan)                                                                 # 调用lavaan包，如果没有安装，需要先安装
curwd=setwd("c:\\program files\\mplus\\mplus examples\\user's guide examples")  # 设定数据文件所在的目录为工作目录，如果脚本与数据文件在同一目录，则不需要设定
semdata <- read.table("ex5.11.dat")                                              # 读取数据文件到cfadata，cfadata是自己命名，可以随便定
names(semdata) <-c(paste("y", 1:12, sep=""))                                      # 给变量命名
semmodel <- '                                                                    # 设置模型，semmodel还是自己命名
              f1=~ y1+y2+y3                                                      # 测量模型，用=~
              f2=~ y4+y5+y6
              f3=~ y7+y8+y9
              f4=~ y10+y11+y12
              f3 ~ f1+f2                                                         # 路径
              f4 ~ f3
              f1 ~~ f2                                                           # f1与f2相关
 '
semfit <- sem(semmodel, semdata)                                                # 进行CFA，cfafit为自己的名称，方便后面调用
summary(semfit,fit.measures="TRUE")                                             # 显示SEM的总体结果
fitMeasures(semfit,fit.measures="all", baseline.model=NULL)                     # 显示所有拟合指数
standardizedSolution(semfit)                                                    # 显示标准化的结果