欢迎您访问 最编程 本站为您分享编程语言代码,编程技术文章!
您现在的位置是: 首页

在 R 代码数据中实现用于结构方程建模和路径图可视化的 sem

最编程 2024-06-07 11:05:41
...

原文链接:tecdat.cn/?p=23312

最近我们被客户要求撰写关于结构方程建模的研究报告,包括一些图形和统计输出。

结构方程模型是一个线性模型框架,它对潜变量同时进行回归方程建模

引言

 诸如线性回归、多元回归、路径分析、确认性因子分析和结构回归等模型都可以被认为是SEM的特例。在SEM中可能存在以下关系。

  • 观察到的变量与观察到的变量之间的关系(γ,如回归)。
  • 潜变量与观察变量(λ,如确认性因子分析)。
  • 潜变量与潜变量(γ,β,如结构回归)。

SEM独特地包含了测量和结构模型。测量模型将观测变量与潜变量联系起来,结构模型将潜变量与潜变量联系起来。目前有多种软件处理SEM模型,包括Mplus、EQS、SAS PROC CALIS、Stata的sem和最近的R的lavaan。R的好处是它是开源的,可以免费使用,而且相对容易使用。

本文将介绍属于SEM框架的最常见的模型,包括

  • 简单回归
  • 多元回归
  • 多变量回归
  • 路径分析
  • 确认性因素分析
  • 结构回归

目的是在每个模型中介绍其

  • 矩阵表述
  • 路径图
  • lavaan语法
  • 参数和输出

在这次训练结束时,你应该能够理解这些概念,足以正确识别模型,认识矩阵表述中的每个参数,并解释每个模型的输出。

语法简介

语法一:f3~f1+f2(路径模型)

结构方程模型的路径部分可以看作是一个回归方程。而在R中,回归方程可以表示为yax1+bx2+c,“”的左边的因变量,右边是自变量,“+”把多个自变量组合在一起。那么把y看作是内生潜变量,把x看作是外生潜变量,略去截距,就构成了语法一。

语法二:f1 =~ item1 + item2 + item3(测量模型)

"=~"的左边是潜变量,右边是观测变量,整句理解为潜变量f1由观测变量item1、item2和item3表现。

语法三:item1 item1 , item1 item2

"~~"的两边相同,表示该变量的方差,不同的话表示两者的协方差

语法四:f1 ~ 1

表示截距

基础知识

加载数据

在这种情况下,我们将模拟数据。

y ~ .5*f  #有外部标准的回归强度


f =~ .8*x1 + .8*x2 + .8*x3 + .8*x4 + .8*x5  #定义因子f,在5个项目上的载荷。


x1 ~~ (1-.8^2)*x1 #残差。请注意,通过使用1平方的载荷,我们在每个指标中实现了1.0的总变异性(标准化的)。
......

#产生数据;注意,标准化的lv是默认的
simData 

#看一下数据
describe(simData)[,1:4]

图片

指定模型

y ~ f # "~回归"
f =~ x1+ x2 + x3 + x4 + x5 # "=~被测量的是"
x1 ~~ x1 # 方差
x2 ~~ x2 #方差
x3~~x3 #变量
x4~~x4 #变量
x5~~x5 #变量
#x4~~x5将是协方差的一个例子

拟合模型

summary(model_m)

图片

图片

图片

inspect(model_m)

图片

Paths

图片

路径分析

与上述步骤相同,但主要侧重于回归路径。值得注意的是这种方法对调节分析的效用。


点击标题查阅往期内容

图片

R语言:结构方程模型、潜变量分析

图片

左右滑动查看更多

图片

01

图片

02

图片

03

图片

04

图片

##加载数据
set.seed(1234)

Data <- data.frame(X = X, Y = Y, M = M)

指定模型

 # 直接效应
             Y ~ c*#使用字符来命名回归路径
           # 调节变量
             M ~ a*X
             Y ~ b*M
           # 间接效应(a*b)
             ab := a*#定义新参数
           # 总效应
             total := c + (a*b) #使用":="定义新参数

拟合模型

summary(model_m)

图片

图片

图片

Paths(model)

图片

间接效应的Bootstrapping置信区间

除了指定对5000个样本的标准误差进行bootstrapping外,下面的语法还指出标准误差应进行偏差校正(但不是accelearted)。这种方法将产生与SPSS中的PROCESS宏程序类似的结果,即对标准误差进行偏差修正。

sem(medmodel,se = "bootstrap")

图片

图片

图片

确认性因素分析

加载数据

我们将使用例子中的相同数据

指定模型

'
f =~ x1 + x2 + x3 +x4 + x5
x1~~x1
x2~~x2
x3~~x3
x4~~x4
x5~~x5
'

拟合模型

sem(fit, simData)

图片

图片

Paths(fit)

图片

anova

图片

正如各模型的LRT所示,sem()和cfa()是具有相同默认值的软件包。CFA可以很容易地使用cfa()或sem()完成 结构方程模型

加载数据

在这种情况下,我将模拟数据。

#结构成分
y ~ .5*f1 + .7*f2 #用外部标准回归的强度


#测量部分
f1 =~ .8*x1 + .6*x2 + .7*x3 + .8*x4 + .75*x5 #定义因子f,在5个项目上的载荷。

x1 ~~ (1-.8^2)*x1 #残差。注意,通过使用1平方的载荷,我们实现了每个指标的总变异性为1.0(标准化)。
...

#生成数据;注意,标准化的lv是默认的
sim <- sim(tosim)

#看一下数据
describe(sim )

 图片 

指定模型

测试正确的模型

#结构性
y ~ f1+ f2
#测量
f1 =~ x1 + x2 + x3 + x4 + x5 
f2 =~ x6 + x7

测试不正确的模型。假设我们错误地认为X4和X5负载于因子2。

incorrect
#结构性
y ~ f1+ f2
#测量
f1 =~ x1 + x2 + x3 
f2 =~ x6 + x7 + x4 + x5

拟合模型

正确的模型

图片

summary(model_m)

图片

图片

图片

不正确的模型

图片

summary(incorrectmodel_m, fit.measures = TRUE)

图片

图片

图片

比较模型

正确模型

图片

 不正确模型

Paths(incorrec)

图片

anova

图片

除了不正确模型的整体拟合指数较差--如CFI<0.95,RMSEA>0.06,SRMR>0.08和Chi-square test<0.05所示,正确模型也优于不正确模型,如正确模型的AIC和BIC低得多所示。

图片

点击文末 “阅读原文”

获取全文完整资料。

本文选自《R语言中实现sem进行结构方程建模和路径图可视化》。

点击标题查阅往期内容

结构方程模型SEM分析心理学和营销研究数据路径图可视化结构方程模型SEM、路径分析房价和犯罪率数据、预测智力影响因素可视化2案例
结构方程模型 SEM 多元回归和模型诊断分析学生测试成绩数据与可视化
在R语言中实现sem进行结构方程建模和路径图可视化
R语言结构方程SEM中的power analysis 效能检验分析stata如何处理结构方程模型(SEM)中具有缺失值的协变量
R语言基于协方差的结构方程拟合的卡方检验
R语言基于协方差的SEM结构方程模型中的拟合指数
R语言:结构方程模型、潜变量分析
R语言对用电负荷时间序列数据进行K-medoids聚类建模和GAM回归
R语言进行支持向量机回归SVR和网格搜索超参数优化
R语言贝叶斯MCMC:用rstan建立线性回归模型分析汽车数据和可视化诊断
R语言深度学习:用keras神经网络回归模型预测时间序列数据
R语言贝叶斯MCMC:GLM逻辑回归、Rstan线性回归、Metropolis Hastings与Gibbs采样算法实例
R语言使用虚拟变量(Dummy Variables) 回归分析工资影响因素

推荐阅读