R语言基于协方差的结构方程拟合的卡方检验
创始人
2024-04-23 02:34:09
0

在评估结构方程模型的拟合,很常见的应用是研究χ2进行测试,因为在给定足够大的样本量的情况下,它几乎总会检测出模型与数据之间的统计上的显着差异。因为,我们的模型几乎总是数据的近似值。如果我们的模型的协方差矩阵实际上匹配抽样变异中的样本协方差矩阵,该χ2 无论样本量多大,该检验在统计学上均无统计学意义。

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

 相关视频:结构方程模型SEM分析心理学营销数据路径图可视化|数据分享

什么是结构方程建模SEM和R语言心理学和营销研究数据路径图可视化

因为到大样本量,从业人员往往依赖于其他拟合指数,如RMSEACFITLI-所有这些都是基于χ 2。在lavaan中,您会自动使用置信区间和p值对RMSEA进行紧密拟合测试。这个测试实际上使用χ2分布。

RMSEA的公式为:

其中,χ2是χ2模型的检验统计量,dF是模型自由度,N是样本量。

如果你的模型拟合数据完美,分子为零;这是标准的假设χ 2χ2-test测试。如果我们在RMSEA进行测试中,使用χ 2参数对应于RMSEA为0.05的分布。Lavaan将测试结果报告为拟合统计之一。

那么这对我们有什么帮助呢?非中心参数(λ )在lavaan的RMSEA测试实际上是χ 2 - d ˚Fχ2-dF对应于RMSEA为0.05的值。

因此,对于测试,λ 是:

对于中等拟合的测试,λ 是:

请注意,lavaan的处理方式可能有所不同。

因此,给定模型的自由度和样本量,我们可以计算出非中心性参数(λ )。给定λ中,χ2 值和模型的自由度,我们可以计算p值进行测试。

R的语法是:

示范

# 具有HolzingerSwineford1939数据集的双因子模型的模型语法
writeLines(syntax <- paste(paste("g =~", paste0("x", 1:9, collapse = " + ")),# paste("visual =~", paste0("x", 1:3, collapse = " + ")),paste("textual =~", paste0("x", 4:6, collapse = " + ")),paste("speed =~", paste0("x", 7:9, collapse = " + ")),sep = "\n"
))g =~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9

运行模型并报告拟合度。仅报告统计信息:

 lavaan (0.5-23.1097) converged normally after  25 iterationsNumber of observations                           301Estimator                                         MLMinimum Function Test Statistic               42.291Degrees of freedom                                21P-value (Chi-square)                           0.004Root Mean Square Error of Approximation:RMSEA                                          0.05890 Percent Confidence Interval          0.032  0.083P-value RMSEA <= 0.05                          0.276

卡方统计意义显着,该完美拟合检验表明,由于样本的变异性。

默认的卡方检验:

pchisq [1] 0.003867178

使用上面的公式计算紧密度测试的非中心参数:.0025乘以模型自由度乘以样本大小-1

 ncp.close  [1] 15.75

计算紧密拟合的卡方检验:

pchisq [1] 0.2740353

紧密契合度测试的p值为.27,接近lavaan报告的值。

如果我们降低标准以进行中等拟合的卡方检验:.0064乘以模型自由度乘以样本大小-1

 ncp.med  [1] 40.32pchisq [1] 0.9199686

我们在模型中观察模型隐含的协方差矩阵的可能性为92%。非常好。


最后,SEM从业者通常报告χ 2-test,但通常希望该测试能够检测到模型规范错误,因此在实践中经常将其忽略。


PS:潜在变量建模的另一种方法是PLS路径建模。这是一种基于OLS回归的SEM方法。


  1. MacCallum, R. C., Browne, M. W., & Sugawara, H. M. (1996). Power analysis and determination of sample size for covariance structure modeling. Psychological Methods, 1(2), 130–149. https://doi.org/10.1037/1082-989X.1.2.130 ↩


 

相关内容

热门资讯

【MySQL】锁 锁 文章目录锁全局锁表级锁表锁元数据锁(MDL)意向锁AUTO-INC锁...
【内网安全】 隧道搭建穿透上线... 文章目录内网穿透-Ngrok-入门-上线1、服务端配置:2、客户端连接服务端ÿ...
GCN的几种模型复现笔记 引言 本篇笔记紧接上文,主要是上一篇看写了快2w字,再去接入代码感觉有点...
数据分页展示逻辑 import java.util.Arrays;import java.util.List;impo...
Redis为什么选择单线程?R... 目录专栏导读一、Redis版本迭代二、Redis4.0之前为什么一直采用单线程?三、R...
【已解决】ERROR: Cou... 正确指令: pip install pyyaml
关于测试,我发现了哪些新大陆 关于测试 平常也只是听说过一些关于测试的术语,但并没有使用过测试工具。偶然看到编程老师...
Lock 接口解读 前置知识点Synchronized synchronized 是 Java 中的关键字,...
Win7 专业版安装中文包、汉... 参考资料:http://www.metsky.com/archives/350.htm...
3 ROS1通讯编程提高(1) 3 ROS1通讯编程提高3.1 使用VS Code编译ROS13.1.1 VS Code的安装和配置...
大模型未来趋势 大模型是人工智能领域的重要发展趋势之一,未来有着广阔的应用前景和发展空间。以下是大模型未来的趋势和展...
python实战应用讲解-【n... 目录 如何在Python中计算残余的平方和 方法1:使用其Base公式 方法2:使用statsmod...
学习u-boot 需要了解的m... 一、常用函数 1. origin 函数 origin 函数的返回值就是变量来源。使用格式如下...
常用python爬虫库介绍与简... 通用 urllib -网络库(stdlib)。 requests -网络库。 grab – 网络库&...
药品批准文号查询|药融云-中国... 药品批文是国家食品药品监督管理局(NMPA)对药品的审评和批准的证明文件...
【2023-03-22】SRS... 【2023-03-22】SRS推流搭配FFmpeg实现目标检测 说明: 外侧测试使用SRS播放器测...
有限元三角形单元的等效节点力 文章目录前言一、重新复习一下有限元三角形单元的理论1、三角形单元的形函数(Nÿ...
初级算法-哈希表 主要记录算法和数据结构学习笔记,新的一年更上一层楼! 初级算法-哈希表...
进程间通信【Linux】 1. 进程间通信 1.1 什么是进程间通信 在 Linux 系统中,进程间通信...
【Docker】P3 Dock... Docker数据卷、宿主机与挂载数据卷的概念及作用挂载宿主机配置数据卷挂载操作示例一个容器挂载多个目...