本文主要介绍基于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)相关推荐

  1. 基于R语言的Meta分析【全流程、不确定性分析】方法与Meta机器学习技术应用

    Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...

  2. 基于R语言的Meta分析【全流程、不确定性分析】方法与Meta机器学习应用

    Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...

  3. 基于R语言的关联规则分析项目

    转自http://blog.163.com/dm_team/blog/static/2379750132014891084989/ 摘要 抱着成为一名优秀的数据分析师/数据挖掘师的理想,我们组成了一个 ...

  4. R数据分析:潜在转化分析LTA的做法和解释(一)

    之前给大家写了很多潜在类别分析的教程Mplus教程:如何做潜在类别分析LCA R数据分析:用R语言做潜类别分析LCA Mplus数据分析:潜在类别分析(LCA)流程(详细版) R数据分析:再写潜在类别 ...

  5. 基于R语言的Copula变量相关性分析及应用

    在工程.水文和金融等各学科的研究中,总是会遇到很多变量,研究这些相互纠缠的变量间的相关关系是各学科的研究的重点.虽然皮尔逊相关.秩相关等相关系数提供了变量间相关关系的粗略结果,但这些系数都存在着无法克 ...

  6. Mplus教程:如何做潜在类别分析LCA

    之前潜在类别分析已经给大家写的很详细了,不过都是基于R软件的,我发现一个很奇怪的地方哦,找我做潜类别的基本都不用R,找我解决R语言问题的基本都不知道潜类别是干啥的,看样子做潜变量分析还是mplus受众 ...

  7. 基于R语言极值统计学及其在相关领域中的应用

    受到气候变化.温室效应以及人类活动等因素的影响,自然界中极端高温.极端环境污染.大洪水和大暴雨等现象的发生日益频繁:在人类社会中,股市崩溃.金融危机等极端情况也时有发生:今年的新冠疫情就是非常典型的极 ...

  8. 【问题导向】利用R语言进行情感分析

    目录 前言 问题分析 R语言的安装 情感分析 分类的原理 未来学习计划 引文 前言 R语言是统计领域的一个重要工具,其中有很多独特的包能够实现特定功能,在进行GIS分析时有非常大的用处. 笔者最近在帮 ...

  9. 基于R语言的代理模型(高斯过程、贝叶斯优化、敏感性分析、异方差性等)高级技术应用

    基于R语言的代理模型(高斯过程.贝叶斯优化.敏感性分析.异方差性等)高级技术应用 直播时间:10月30日-10月31日.11月6日-7日(4天+1周辅导练习) (上午9:30-12:00  下午14: ...

最新文章

  1. python3里的pillow怎么安装_“python安装pillow教程“python3.4怎么安装pil
  2. (25)2-9-9-12分页(下)
  3. linux无法打开共享对象文件或目录,linux - libmodbus.so.5:无法打开共享对象文件:没有这样的文件或目录 - 堆栈内存溢出...
  4. 中石油在意大利加油站0.5分/升,会不会赔死啊?
  5. oracle暂停索引,Oracle索引被抑制情况
  6. -32767转化为二进制_谁说前端不需要懂二进制
  7. 微信公众号发送小程序卡片_微信公众号里怎么添加小程序-如何在微信[[公众号]]添加小程序卡片-微信关联小程序...
  8. ikbc poker2
  9. 横版过关游戏开发-游戏流程
  10. threeJS 创建逼真海洋落日
  11. 英飞凌TC264之双核互斥
  12. docker一键部署springboot项目(三)
  13. OpenCL基础(一)
  14. (转载自)章文嵩博士和他背后的负载均衡帝国
  15. 墨客科技执行董事袁英:MOAC区块链赋能实体产业的方案与实践
  16. 【拆解】Apple Watch Series 6 ,电池更大、陶瓷和蓝宝石外壳更薄,更强硬,更耐磨!...
  17. 实习技术员的基本功(二)
  18. Linux QT 4G发送HTTP POST请求发送JSON格式的数据
  19. excel账龄计算:两个经典公式快速制作账龄统计表
  20. git--“HEAD detached from”的解决方法

热门文章

  1. 不借助 matlab 内置函数,生撸均值方差模型
  2. Java培训班靠谱吗?
  3. 复旦计算机教师,复旦教师:我有无人能及的相声天赋
  4. matlab双纵轴刻度覆盖问题,求助: matlab双纵轴换图问题
  5. 中国房价不可能下降的19个理由
  6. swoole-面向连接(tcp)的心跳检测和合包分包
  7. HTML5期末大作业:动漫人物介绍网站设计——柯南(5页) HTML+CSS+JavaScript 学生DW网页设计作业成品 学生动漫网页设计模板下载
  8. 基于ruoyi+vue+elementUI实现列表,新增,附件上传,tab+springBoot+mybatis+oracle序列+批量新增
  9. keil c语言 1602时钟代码,基于单片机的LCD1602电子时钟设计.doc
  10. 无线wifi-PJ-之在开启WPS下使用reaver