“用R进行有调节的中介路径分析”的版本间的差异
来自OBHRM百科
Lichaoping(讨论 | 贡献) (→脚本与注释) |
Lichaoping(讨论 | 贡献) (→脚本与注释) |
||
第7行: | 第7行: | ||
names(pathdata) <-c("y","m","x","z") # 给变量命名 | names(pathdata) <-c("y","m","x","z") # 给变量命名 | ||
attach(pathdata) # 激活数据 | attach(pathdata) # 激活数据 | ||
− | pathdata$xz <- (x-mean(x))*(z-mean(z)) # | + | pathdata$xz <- (x-mean(x))*(z-mean(z)) # 变量去中心化,计算得到调节变量xz |
pathmodel <- ' m ~ x +z + xz # 设置模型,pathmodel还是自己根据情况命名。因为是路径分析,用的~ | pathmodel <- ' m ~ x +z + xz # 设置模型,pathmodel还是自己根据情况命名。因为是路径分析,用的~ | ||
y ~ x +z + m+ xz' | y ~ x +z + m+ xz' |
2017年3月20日 (一) 16:23的版本
示意图
脚本与注释
library(lavaan) # 调用lavaan包,如果没有安装,需要先安装 curwd=setwd("c:\\program files\\mplus\\mplus examples\\user's guide examples") # 设定数据文件所在的目录为工作目录,如果脚本与数据文件在同一目录,则不需要设定 pathdata <- read.table("ex3.18.dat") # 读取数据文件到pathdata,pathdata是自己命名,可以随便定 names(pathdata) <-c("y","m","x","z") # 给变量命名 attach(pathdata) # 激活数据 pathdata$xz <- (x-mean(x))*(z-mean(z)) # 变量去中心化,计算得到调节变量xz pathmodel <- ' m ~ x +z + xz # 设置模型,pathmodel还是自己根据情况命名。因为是路径分析,用的~ y ~ x +z + m+ xz' pathfit <- sem(pathmodel, pathdata) # 进行SEM,pathfit为自己的名称,方便后面调用 summary(pathfit,fit.measures="TRUE") # 显示SEM的总体结果 fitMeasures(pathfit,fit.measures="all", baseline.model=NULL) # 显示所有拟合指数 standardizedSolution(pathfit) # 显示标准化的结果