PCoA分析,即主坐标分析(principal co-ordinates analysis),也是一种非约束性的数据降维分析方法,可用来研究样本的相似性或差异性,与PCA分析类似;主要区别在于,PCA基于欧氏距离,PCoA基于除欧氏距离以外的其它距离。PCoA分析,首先对一系列的特征值和特征向量进行排序,然后选择排在前几位的最主要特征值,并将其表现在坐标系里,结果相当于是距离矩阵的一个旋转,它没有改变样本点之间的相互位置关系,只是改变了坐标系统。

输入:

1、样本距离矩阵文件,可由分析模块“Generate distance matrix from fpkm matrix”生成。

示例:

T4    T5    T6    T7    T8    T9

T4    0       556.755292869315 548.144810678248 1840.00551404162 4595.29924214683 1500.48982692357

T5    556.755292869315 0       608.699577853257 1817.04805740656 4504.26458143164 1526.05178655135

T6    548.144810678248 608.699577853257 0       1768.22280076409 4633.56038886132 1354.17671707347

T7    1840.00551404162 1817.04805740656 1768.22280076409 0       3697.34366245185 1251.5445135536

T8    4595.29924214683 4504.26458143164 4633.56038886132 3697.34366245185 0       4612.87039733818

T9    1500.48982692357 1526.05178655135 1354.17671707347 1251.5445135536    4612.87039733818 0

2、样品分组信息表(可选)。

示例:

T4     group1

T5     group1

T6     group1

T7     group2

T8     group2

T9     group2

输出:

1、pcoa_sites.txt:记录了样本在各个维度上的位置,用于作图的数据。对PC1和PC2作散点图,则PC1为x轴,PC2为y轴。

2、pcoa_rotation.txt:记录了对各主成分的贡献度。

3、pcoa_importance.txt:记录了各维度解释结果的百分比。如果PC1值为50%,则表示PC1的差异可以解释全面分析结果的50%。

4、pcoa_plot.pdf: PCoA分析图

示例:

PCoA分析图(不提供分组文件)。

PCoA分析图(提供分组文件)。不同分组的样品以不同的颜色标识。

注:不同颜色或形状的点代表不同分组情况下的样品(由分组信息表提供)。横、纵坐标轴的刻度是相对距离,无实际意义。PC1、PC2分别代表对于两组样本基因表达模式的可能影响因素,需要结合样本特征信息归纳总结,例如A组(红色)和B组(蓝色)样本在PC2轴的方向上分离开来,则可分析为PC2是导致A组和B组分开的因素(可以是两种不同的处理或两个不同的品种),同时验证了这个因素有较高的可能性影响了基因的表达模式。

分析模块引用R语言(v3.2.1)中的PCoA统计和分析函数进行数据运算和做图。

r语言做绘制精美pcoa图_PCOA分析相关推荐

  1. r语言做绘制精美pcoa图_R语言进行PCoA分析

    相信大家在做微生物多样性研究时经常听到PCA分析.PCoA分析,NMDS分析,CCA分析,RDA分析. 它们对 物种(或基因.功能)的分析具有重要作用,因而频频出现在16S测序及宏基因组测序中. 那么 ...

  2. r语言做绘制精美pcoa图_科学网—R语言 PCA PCoA ggplot2 - 靳泽星的博文

    这一篇是衔接上一篇的,就是要用ggplot2程序包对PCA和PCoA进行可视化.代码我直接照搬过来了,只是绘图的时候用ggplot函数.ggplot2包实现了一个在R中基于全面一致的语法创建图形时的系 ...

  3. r语言做绘制精美pcoa图_三分钟绘制一张优美的PCoA图 | 云平台

    本文主要为您迅速绘制一张优美的PCoA图提供步骤指导(1-4要点),另附细节优化方法(5及之后的要点),可后期慢慢研究调整,同时提供相关分析方法说明和参考文献,为您深入了解提供支持. 1.了解输入数据 ...

  4. r语言做绘制精美pcoa图_R语言统计与绘图:绘制QQ图

    QQ图一般用来比较两组数字的分布,比如两个变量或两个向量. 在国外,顾客习惯性会给服用人员小费,给多少小费是一个经常被讨论的话题,今天使用的数据集为reshape2包的tips数据集,用来演示QQ图的 ...

  5. r语言做绘制精美pcoa图_如何绘制精美的PCoA图形

    匿名用户 1级 2017-06-20 回答 虽然一般的16S或者宏基因组等分析e69da5e6ba90e799bee5baa6e79fa5e9819331333363373762流程当中都会包含PCo ...

  6. r语言做绘制精美pcoa图_科学网-R语言 PCA PCoA ggplot2-靳泽星的博文

    这一篇是衔接上一篇的,就是要用ggplot2程序包对PCA和PCoA进行可视化.代码我直接照搬过来了,只是绘图的时候用ggplot函数.ggplot2包实现了一个在R中基于全面一致的语法创建图形时的系 ...

  7. r语言做绘制精美pcoa图_R语言:Bary-Curtis PCoA

    导读 Bary-Curtis PCoA降维分析菌群结构. 一.输入数据 1.1 菌属丰度矩阵 df = data.frame(abs(round(matrix(rnorm(729, 100, 50), ...

  8. r语言绘制精美pcoa图_如何绘制精美的PCoA图形?

    原标题:如何绘制精美的PCoA图形? 今天我们来分享干货--PCoA图形的代码.继PCA.火山图.热图等代码后,基迪奥的程序猿又整理出PCoA代码.具体往期我们分享过的代码贴,可以在文末查看哦. 什么 ...

  9. r语言绘制精美pcoa图_R语言绘制交互式热图

    热图 通过热图可以简单地聚合大量数据,并使用一种渐进的色带来优雅地表现,最终效果一般优于离散点的直接显示,可以很直观地展现空间数据的疏密程度或频率高低.但也由于很直观,热图在数据表现的准确性并不能保证 ...

最新文章

  1. 全网最细 | 21张图带你领略集合的线程不安全
  2. 如何模拟将CPU、IO打满?
  3. php mysql.x86 64_centos7.2yum安装php70w.x86_64
  4. RabbitMQ系列教程之二:工作队列(Work Queues)
  5. mlp神经网络_白天鹅黑天鹅灰天鹅?卷积神经网络帮你搞定识别
  6. UNIX(进程间通信):03---僵尸进程
  7. 使用MetalLB在您的Raspberry Pi家庭实验室上安装Kubernetes负载均衡器
  8. 【UML】用例图Use Case diagram(转)
  9. 大数据可视化的优点有哪些
  10. Java动态规划走金字塔_【动态规划基础】数字金字塔
  11. java学习软件_刚学习java,用哪些学习软件比较好?
  12. 苏大计算机考研专业课,苏州大学软件工程考研初试科目考什么?
  13. win10安装(U盘启动)
  14. matlab中Ctphi什么意思,Matlab在物理上的应用举例
  15. Vue不同用户权限显示不同内容
  16. 查看pcie总线上的设备ID以及带宽计算
  17. 联想的高清壁纸蛮不错的,速提!
  18. word图文混排复制到CuteEditor图片不显示
  19. 计算机网络胡工程施工税率,弱电项目增值税6%、9%、13%税率怎样区分?项目经理必知...
  20. GoogleMap的简单用法

热门文章

  1. 九章算法 | 苏州微软面试题:程序检查
  2. 第一章 第二节 数字化含义及作用——丁老师
  3. WPF全屏和退出全屏最简单的方法
  4. 多级评论回复功能(纯前端,未连数据库)
  5. 微信白名单服务器ip地址,微信服务器白名单IP列表
  6. Android面试常见问题汇总
  7. 一个高尚的人,一个纯粹的人,一个有道德的人,一个脱离了低级趣味的人,一个有益于人民的人。...
  8. CVPR 2022 | 百变发型!中科大/微软/港城大提出HairCLIP:基于文本和参考图像的头发编辑方法...
  9. 试验数据管理系统TDM6.0产品功能介绍之——签到管理
  10. 全连接层输入为什么是固定维度的(拉直/压扁Flatten成为列向量)