“用R进行结构方程模型分析”的版本间的差异
来自OBHRM百科
Lichaoping（讨论 | 贡献） （→脚本与注释） |
Lichaoping（讨论 | 贡献） （→脚本与注释） |
||
| 第2行： | 第2行： | ||
<pre>library(lavaan) # 调用lavaan包，如果没有安装，需要先安装 | <pre>library(lavaan) # 调用lavaan包，如果没有安装，需要先安装 | ||
curwd=setwd("c:\\program files\\mplus\\mplus examples\\user's guide examples") # 设定数据文件所在的目录为工作目录，如果脚本与数据文件在同一目录，则不需要设定 | curwd=setwd("c:\\program files\\mplus\\mplus examples\\user's guide examples") # 设定数据文件所在的目录为工作目录，如果脚本与数据文件在同一目录，则不需要设定 | ||
| − | semdata <- read.table("ex5.11.dat") | + | semdata <- read.table("ex5.11.dat") # 读取数据文件到cfadata，cfadata是自己命名，可以随便定 |
| − | names(semdata) <-c(paste("y", 1:12, sep="")) | + | names(semdata) <-c(paste("y", 1:12, sep="")) # 给变量命名 |
| − | semmodel <- ' | + | semmodel <- ' # 设置模型，semmodel还是自己命名 |
| − | f1=~ y1+y2+y3 | + | f1=~ y1+y2+y3 # 测量模型，用=~ |
f2=~ y4+y5+y6 | f2=~ y4+y5+y6 | ||
f3=~ y7+y8+y9 | f3=~ y7+y8+y9 | ||
f4=~ y10+y11+y12 | f4=~ y10+y11+y12 | ||
| − | f3 ~ f1+f2 | + | f3 ~ f1+f2 # 路径 |
f4 ~ f3 | f4 ~ f3 | ||
| − | f1 ~~ f2 | + | f1 ~~ f2 # f1与f2相关 |
' | ' | ||
| − | semfit <- sem(semmodel, semdata) | + | semfit <- sem(semmodel, semdata) # 计算，semfit为自己的名称，方便后面调用 |
| − | summary(semfit,fit.measures="TRUE") | + | summary(semfit,fit.measures="TRUE") # 显示SEM的总体结果 |
| − | fitMeasures(semfit,fit.measures="all", baseline.model=NULL) | + | fitMeasures(semfit,fit.measures="all", baseline.model=NULL) # 显示所有拟合指数 |
| − | standardizedSolution(semfit) | + | standardizedSolution(semfit) # 显示标准化的结果</pre> |
==结果== | ==结果== | ||
2017年3月6日 (一) 20:57的版本
脚本与注释
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) # 计算，semfit为自己的名称，方便后面调用
summary(semfit,fit.measures="TRUE") # 显示SEM的总体结果
fitMeasures(semfit,fit.measures="all", baseline.model=NULL) # 显示所有拟合指数
standardizedSolution(semfit) # 显示标准化的结果
