“用R进行结构方程模型分析”的版本间的差异

来自OBHRM百科
跳转至: 导航搜索
(创建页面,内容为“==脚本与注释== ==结果==”)
 
脚本与注释
第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)                                                    # 显示标准化的结果

结果