[统计模型] 基于R的潜在剖面分析(LPA)
本文主要介绍基于R进行潜在剖面分析(Latent Profile Aanalysis, LPA)。
关于潜在类别分析LCA/潜在剖面分析LPA的概念和区别,可以参考一下这篇文章:
潜在类别分析LCA \潜在剖面分析LPA介绍 https://www.jianshu.com/p/4c830094131e?utm_campaign=maleskine
但是此文对潜类别分析的阐释还比较生涩,而且没有具体的代码编写、模型优选等内容。本文旨在丰富LPA的内容,并希望和大家交流、共同进步。
1. 什么是潜类别分析、潜剖面分析,两者的区别在哪里?
通俗的说,就是我们想对我们数据中的个案进行分类,但是没有一个确切的指标可以告诉我们把他们分为几类。举个例子,比如我有1万个社会网络,社会网络有很多结构指标(比如网络的节点数、网络密度、网络的直径等),我想对这1万个网络按照网络指标进行分类,但是并不知道通过哪些变量进行分类,潜类别分析就派上了用场,潜类别分析就是用潜在的变量表征所有的结构指标变量,然后这1万个网络按照这些潜在变量就可以分成若干类别。
从上文中的简书内容我们可以看到,潜类别分析的变量是分类变量,而潜剖面分析纳入的变量是连续变量,这是两种分析方法最明显的区别。
2. 基于R的tidyLPA包进行LPA分析
假设数据已经准备好了,均是连续型变量,R的分析代码如下:
install.packages("tidyLPA")
install.packages("mclust")
install.packages("MplusAutomation")
library(MplusAutomation)
library(mclust)
library(tidyLPA)#读取数据
lpa_data<-read.csv('lpa_continuous_variable.csv',header = TRUE)#要分成几类,n_profiles就设置为几类
#关于models的取值,建议参考R包的文档,里面有解释,我后期会更新本文进行详解
m2<-estimate_profiles(df=lpa_data, n_profiles=2, models=1,select_vars=c('opioid', 'syn', 'idu','minority', 'female', 'unmarry','low_edu', 'west', 'mid_age', 'mmt','nodes', 'density', 'avg_deg', 'modula','diameter', 'avg_dis'))m2_fit<-get_fit(m2)
m2_class<-get_data(m2)
m2_estimates<-get_estimates(m2)
write.csv(m2_class,"m2_class.csv")
write.csv(m2_estimates,"m2_estimates.csv")
write.csv(m2_fit,"m2_fit.csv")
3.模型优选
我们将数据分为了2、3、4、5、6类,到底分为几类是最好的呢?这里建议参考两本书:
《潜变量建模与Mplus应用-基础篇》和 《潜变量建模与Mplus应用-进阶篇》,关注并点赞本博客,私信我邮箱给你发送哟~
其实大致意思就是把各类模型结果中的AIC、BIC、aBIC做成折线图,看拐点在几分类上,然后结合现实意义进行解释即可~
[统计模型] 基于R的潜在剖面分析(LPA)相关推荐
- 基于R语言的Meta分析【全流程、不确定性分析】方法与Meta机器学习技术应用
Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...
- 基于R语言的Meta分析【全流程、不确定性分析】方法与Meta机器学习应用
Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...
- 基于R语言的关联规则分析项目
转自http://blog.163.com/dm_team/blog/static/2379750132014891084989/ 摘要 抱着成为一名优秀的数据分析师/数据挖掘师的理想,我们组成了一个 ...
- R数据分析:潜在转化分析LTA的做法和解释(一)
之前给大家写了很多潜在类别分析的教程Mplus教程:如何做潜在类别分析LCA R数据分析:用R语言做潜类别分析LCA Mplus数据分析:潜在类别分析(LCA)流程(详细版) R数据分析:再写潜在类别 ...
- 基于R语言的Copula变量相关性分析及应用
在工程.水文和金融等各学科的研究中,总是会遇到很多变量,研究这些相互纠缠的变量间的相关关系是各学科的研究的重点.虽然皮尔逊相关.秩相关等相关系数提供了变量间相关关系的粗略结果,但这些系数都存在着无法克 ...
- Mplus教程:如何做潜在类别分析LCA
之前潜在类别分析已经给大家写的很详细了,不过都是基于R软件的,我发现一个很奇怪的地方哦,找我做潜类别的基本都不用R,找我解决R语言问题的基本都不知道潜类别是干啥的,看样子做潜变量分析还是mplus受众 ...
- 基于R语言极值统计学及其在相关领域中的应用
受到气候变化.温室效应以及人类活动等因素的影响,自然界中极端高温.极端环境污染.大洪水和大暴雨等现象的发生日益频繁:在人类社会中,股市崩溃.金融危机等极端情况也时有发生:今年的新冠疫情就是非常典型的极 ...
- 【问题导向】利用R语言进行情感分析
目录 前言 问题分析 R语言的安装 情感分析 分类的原理 未来学习计划 引文 前言 R语言是统计领域的一个重要工具,其中有很多独特的包能够实现特定功能,在进行GIS分析时有非常大的用处. 笔者最近在帮 ...
- 基于R语言的代理模型(高斯过程、贝叶斯优化、敏感性分析、异方差性等)高级技术应用
基于R语言的代理模型(高斯过程.贝叶斯优化.敏感性分析.异方差性等)高级技术应用 直播时间:10月30日-10月31日.11月6日-7日(4天+1周辅导练习) (上午9:30-12:00 下午14: ...
最新文章
- python3里的pillow怎么安装_“python安装pillow教程“python3.4怎么安装pil
- (25)2-9-9-12分页(下)
- linux无法打开共享对象文件或目录,linux - libmodbus.so.5:无法打开共享对象文件:没有这样的文件或目录 - 堆栈内存溢出...
- 中石油在意大利加油站0.5分/升,会不会赔死啊?
- oracle暂停索引,Oracle索引被抑制情况
- -32767转化为二进制_谁说前端不需要懂二进制
- 微信公众号发送小程序卡片_微信公众号里怎么添加小程序-如何在微信[[公众号]]添加小程序卡片-微信关联小程序...
- ikbc poker2
- 横版过关游戏开发-游戏流程
- threeJS 创建逼真海洋落日
- 英飞凌TC264之双核互斥
- docker一键部署springboot项目(三)
- OpenCL基础(一)
- (转载自)章文嵩博士和他背后的负载均衡帝国
- 墨客科技执行董事袁英:MOAC区块链赋能实体产业的方案与实践
- 【拆解】Apple Watch Series 6 ,电池更大、陶瓷和蓝宝石外壳更薄,更强硬,更耐磨!...
- 实习技术员的基本功(二)
- Linux QT 4G发送HTTP POST请求发送JSON格式的数据
- excel账龄计算:两个经典公式快速制作账龄统计表
- git--“HEAD detached from”的解决方法
热门文章
- 不借助 matlab 内置函数,生撸均值方差模型
- Java培训班靠谱吗?
- 复旦计算机教师,复旦教师:我有无人能及的相声天赋
- matlab双纵轴刻度覆盖问题,求助: matlab双纵轴换图问题
- 中国房价不可能下降的19个理由
- swoole-面向连接(tcp)的心跳检测和合包分包
- HTML5期末大作业:动漫人物介绍网站设计——柯南(5页) HTML+CSS+JavaScript 学生DW网页设计作业成品 学生动漫网页设计模板下载
- 基于ruoyi+vue+elementUI实现列表,新增,附件上传,tab+springBoot+mybatis+oracle序列+批量新增
- keil c语言 1602时钟代码,基于单片机的LCD1602电子时钟设计.doc
- 无线wifi-PJ-之在开启WPS下使用reaver