R语言复现 | 潜在剖面模型(LPA)解读+R语言代码

培训通知

Nhanes数据库数据挖掘，快速发表发文的利器，你来试试吧！欢迎报名郑老师团队统计课程，4.20直播。

潜在剖面模型

本文是潜变量文章第5篇。

大家是不是手头有的数据都是一堆定量数据（年龄、血压、血脂、血糖、BMI……），如果按照正常写文章的套路，只要这些不是自变量因变量，就只能出现在表1里统计描述里，没有了其他价值。

所以，能不能换一个思路，重新利用起来这些数据呢？

假如我根据这些数据将人群分类了呢？

那我就可以得到年轻亚健康组（年纪小但有三高）；身体健康组（年纪大但没有三高）；老年亚健康组（年纪大有三高）。根据这个分类，是不是又可以产出一篇文章啦！

按照常规套路，应该先介绍一下潜变量、显变量之类的，但是弯弯绕绕一大堆也不太理解，

本文就去掉定义，你就只需要知道：

潜在剖面模型：针对横断面数据，把多个定量数据变为几个分类，也就是把人群重新分组，然后接下去该用logsitic/线性等统计方法就该用哪种统计方法来分析。

至于具体是几分类，数据放进去，代码跑出来的结果自然会告诉你！

一、文献案例介绍

今天要介绍简单的中文文献。

本公众号回复“ 沙龙 ”即可获得R语言代码，PPT，数据等资料

1.1 文章介绍

我们略过其他的，直接看关于潜在剖面模型的重点。文章是根据短视频成瘾量表的四个维度得分，也就是4个定量数据在来把人群分类。

1.2 统计学方法

这里讲了每个指标的情况，也就是用这些指标来判断具体分几个剖面（分几类）

①赤池信息准则（AIC）、贝叶斯信息准则（BIC）和校正后aBIC的值越小，则模型的拟合越好;

②LMR和BLRT两个指标的P<0.05时，表明k个类别的模型显著优于k-1个类别的模型;

③Entropy的值达到0.80以上时表明分类的精准性较高，越接近于1表明分类的可信度越高。

1.3 文章结果展示

研究结果表明，以1分类模型为基础，逐渐增加剖面数量，探索1～6分类的潜在剖面，各模型的适配度拟合指标见表1。

• AIC、BIC、aBIC值随分类数增加而逐渐减小

• BLRT所有模型均提示有统计学意义（P值均<0.01）

• LMR在2,3,5分类模型中有统计学意义（P值均<0.01），4和6分类模型LMR无统计学意义（P值均>0.05），故排除;

• 2和5分类模型的Entropy值<0.8，也予以排除。

综上，选择3分类模型（C1、C2、C3）为最佳剖面模型。

这部分介绍就写的很详细啦，大家按照这样的流程，再结合自己专业分组情况综合判断哦！

二、R语言演示

本次演示包括的统计学方法有：

• 查看数据有无缺失

• 查看数据应该分几个剖面（分几类）

• 拟合模型

• 画图

• 整合数据

因为潜在剖面模型的公共数据库的文章比较少，所以这次就没办法复现，只能用普通的数据来做啦，这次的数据跑出来的结果不是特别有代表性，但是说不定换上你们自己的数据各个结果都很好看哦！

（注：本公众号回复“沙龙”即可获取R语言代码）

2.1 首先，加载R包和导入数据。

首先，安装需要的包，接下去导入从qgraph包自带五大人格量表的部分数据，数据一共500行。

2.2 查看数据有无缺失

潜在剖面模型要求数据无确失，所以先查看一下数据是不是有缺失的，可以看到整个数据都是没有缺失的。

2.3 探索一下对于我们的数据我们应该分为几个剖面（几类）

①根据BIC

结果就是下图这样了，是不是有点开始看不懂了，其实初学者不用管什么EEV还是VEI等，你可以理解为这些都是不同模型罢了。

只需要关注逗号后面是9，说明它建议你分为9分类。

②根据ICL

运行结果如下：

结果ICL的结果跟前面BIC做出来的一致，都是建议我们把数据分为9类

③根据BLRT

前面都是EEV的方式做出来的9分类，所以这个红框框部分就写“EEV”

运行结果如下：

建立这里就需要解释了，1分类和2分类比有差异，2分类和3分类比没有差异，所以其实分为2分类是BLRT建议我们的最优选择。

但是前面也说了这个数据不典型，如果数据好的话，这里会跟前面BIC和ICL结果一致的。

具体分几类是综合考虑的结果，这些模型拟合告诉你的数据互相比较，结合自己的专业领域的分类，选出最佳的就可以。那我这里就先选9分类啦，我们接着往下。

（注：本公众号回复“沙龙”即可获取R语言代码）

2.4 拟合模型

这里的9和EEV是前面已经解释过的啦，LPA_md可以自己命名

运行结果如下：

然后，这里就是可以看到9分类，每一类里面的具体样本量。

2.5 画图

提取均数是因为要根据均数来画图，其他所有都不用动，把你的数据集名称填在这就可以啦

运行结果如下：

2.6 数据集整合

最后，把分类放进原来数据集里，你就能看到每行数据在哪一类了。

本公众号回复“ 沙龙 ”即可获得R语言代码，PPT，数据等资料

本公众提供各种科研服务了！

一、课程培训 2022年以来，我们召集了一批富有经验的高校专业队伍，着手举行短期统计课程培训班，包括R语言、meta分析、临床预测模型、真实世界临床研究、问卷与量表分析、医学统计与SPSS、临床试验数据分析、重复测量资料分析、nhanes、孟德尔随机化等10余门课。如果您有需求，不妨点击查看： 发文后退款：2024-2025年科研统计课程介绍 二、数据分析服务 浙江中医药大学郑老师团队接单各项医学研究数据分析的服务，提供高质量统计分析报告。有兴趣了解一下详情： 课题、论文、毕业数据分析   临床试验设计与分析 、公共数据库挖掘与统计