“用R进行结构方程模型分析”的版本间的差异
来自OBHRM百科
Lichaoping(讨论 | 贡献) (创建页面,内容为“==脚本与注释== ==结果==”) |
Lichaoping(讨论 | 贡献) (→脚本与注释) |
||
第1行: | 第1行: | ||
==脚本与注释== | ==脚本与注释== | ||
+ | <pre>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) # 显示标准化的结果</pre> | ||
==结果== | ==结果== |
2017年3月5日 (日) 10:59的版本
脚本与注释
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) # 显示标准化的结果