偏最小二乘法判别分析(PLS_DA)是一种用于判别分析的多变量统计分析方法,一种根据观察或测量到的若干变量值,来判断研究对象如何分类的常用统计分析方法。对不同处理样本(如观测样本、对照样本)的特性分别进行训练,产生训练集,并检验训练集的可信度。以下是PLS_DA分析绘图的一个小示例。(本文仅作绘图示例,不涉及分析说明)

【绘图代码】

install.packages("BiocManager")
## install mixOmics
BiocManager::install('mixOmics')
library(mixOmics)data1=read.table("test.txt", header=T, row.names=1,sep="\t")
group=read.table("group.txt",sep="\t",header=T,row.names=1)
X=t(data1[, rownames(group)])
Y=group$group
##PLS-DA分析,选取2个主成分进行分析
plsda.datatm <-plsda(X, Y, ncomp = 3)
plotIndiv(plsda.datatm, ind.names = FALSE, legend=TRUE,ellipse = TRUE, title="sPLS-DA - final result")
background <- background.predict(plsda.datatm, comp.predicted=2,dist = "max.dist")

#plotVar(plsda.datatm)
plotIndiv(plsda.datatm, comp = 1:2, ind.names = FALSE, title = "Maximum distance",legend = TRUE,  background = background,ellipse = TRUE)
#auc.plsda <- auroc(plsda.datatm)

软文同步于“作图帮”公众号,感兴趣的小伙伴可以关注一下哦~

R语言|PLS_DA分析绘图示例相关推荐

  1. 科研绘图之R语言生存分析KM曲线累计风险表放在图片内部

    科研绘图之R语言生存分析KM曲线和累计风险表 KM估计 R语言展示KM估计的生存函数曲线 1.最简单的方法 2.利用survminer包绘制 3.进一步美化,添加累计风险表格.图例.文本注释 KM估计 ...

  2. R语言生存分析之COX比例风险模型构建及列线图(nomogram)、校准曲线(calibration curve)绘制示例

    R语言生存分析之COX比例风险模型构建及列线图(nomogram).校准曲线(calibration curve)绘制示例 列线图(Alignment Diagram),又称诺莫图(Nomogram图 ...

  3. 技巧 | 如何使用R语言的基础绘图系统的拼图功能

    我们知道ggplot2工具包有很多方便的拼图拓展包,如cowplot.patchwork等,而本篇就来介绍在使用R语言的基础绘图系统如何进行拼图.需要明确的是,基础绘图系统的拼图功能不需要借助任何其他 ...

  4. R语言生存分析可视化分析

    生存分析指的是一系列用来探究所感兴趣的事件的发生的时间的统计方法. 生存分析被用于各种领域,例如: 癌症研究为患者生存时间分析, "事件历史分析"的社会学 在工程的"故障 ...

  5. 视频|分类模型评估:精确率、召回率、ROC曲线、AUC与R语言生存分析时间依赖性ROC实现

    最近我们被客户要求撰写关于分类模型的研究报告,包括一些图形和统计输出. 本文将帮助您回答以下问题: ROC曲线是什么? 曲线下的面积是多少? 二元分类的决策阈值是多少? 分类模型可接受的 AUC值是多 ...

  6. [R语言] R语言PCA分析教程 Principal Component Methods in R

    R语言PCA分析教程 Principal Component Methods in R(代码下载) 主成分分析Principal Component Methods(PCA)允许我们总结和可视化包含由 ...

  7. R语言生存分析数据分析可视化案例

    目标 本文的目的是对如何在R中进行生存分析进行简短而全面的评估.关于该主题的文献很广泛,仅涉及有限数量的(常见)问题.最近我们被客户要求撰写关于生存分析的研究报告,包括一些图形和统计输出. 可用的R包 ...

  8. R语言可以直接从数据库读取数据么?R语言读取数据库数据示例

    R语言可以直接从数据库读取数据么?R语言读取数据库数据示例 目录 R语言可以直接从数据库读取数据么?R语言读取数据库数据示例 R语言是解决什么问题的? R语言可以直接从数据库读取数据么?R语言读取数据 ...

  9. R语言效用分析 ( 效能分析、Power analysis)、除了pwr包之外还有其它包、例如、基因研究中的效能分析、MBESS包可用于各种形式的效能分析和最少样本量确定、其他效用分析包的简要介绍

    R语言效用分析 ( 效能分析.Power analysis).除了pwr包之外还有其它包.例如.基因研究中的效能分析(power analysis).MBESS包可用于各种形式的效能分析(power ...

最新文章

  1. api怎么写_API数据加密框架monkeyapiencrypt
  2. JavaScrip入门-浏览器里的js
  3. 怎么设置ppt页面的长度和宽度_将PPT里的字弄很小,PPT就有逼格吗?
  4. jQuery实战之 attr() 和 prop() 的区别
  5. C语言 迷宫问题求解(顺序栈应用示例)
  6. gtest的单独编译和使用
  7. integer是值传递还是引用传递_值传递与引用传递
  8. java的md5盐值加密_MD5盐值加密
  9. JSON和JavaScript对象互转
  10. c++已知2点求中垂线_呆哥数学函数合集——函数的概念【2】
  11. 树莓派蜜罐节点部署实战
  12. 无法在PC上找到,android app 创建的文件 问题.
  13. 博文视点“阿里云全系列技术图书”隆重亮相2021云栖大会
  14. ufw防火墙配置命令
  15. SageMath密码学密码体制及加解密
  16. 俞敏洪沉默,新东方落泪
  17. 前花椒直播CEO胡震生:我如何不花钱引爆的花椒直播
  18. python爬取58同城的兼职信息
  19. 【汇正财经】扬帆起航,医美市场行业规模扩张
  20. 【Add_two_numbers】

热门文章

  1. 用太极拳讲分布式理论,真舒服!
  2. 肝了3版才满意:分布式系统之CAP理论,我们对它的理解和误解
  3. 你知道Redis可以实现延迟队列吗?
  4. Linux服务端最大并发数是多少?
  5. HashMap 为什么会导致 CPU 100%?文章看不懂?
  6. Redis AOF 持久化详解
  7. 《Java EE企业级应用开发》,《分布式爬虫》等书包邮送50本!企业开发利器!...
  8. 关于远程办公,微软MVP 15年研发团队的经验分享
  9. 如何确定敏捷是否适合你的团队?
  10. 从《翔谈》说起,看美团