原始数据为遗传距离矩阵,使用ape包的pcoa函数进行主坐标分析,然后使用ggplot2进行绘图。

欢迎大家批评指正

转载请标明出处:https://www.ivistang.com/articles/278

简述

主坐标分析(Principal Coordinates Analysis,PCoA),可呈现研究数据相似性或差异性的可视化坐标,是一种非约束性的数据降维分析方法,可用来研究样本群落组成的相似性或相异性。PCoA通过一系列的特征值和特征向量进行排序后,选择主要排在前几位的特征值,找到距离矩阵中最主要的坐标,结果是数据矩阵的一个旋转,它没有改变样本点之间的相互位置关系,只是改变了坐标系统。PCoA的目的是在尽力保证原本的距离关系不发生改变的前提下对高维数据进行低维投影,也就是使得原始数据间点的距离与投影中即结果中各点之间的距离尽可能相关。

代码

library(ape)

library(ggplot2)

ind = read.table("905_ind.txt",skip = 1,row.names = 1, stringsAsFactors = FALSE) ####读取遗传距离矩阵

library(readr)

group = read_csv("905_group.csv", col_names = FALSE) ####读取分组数据

colnames(ind)=rownames(ind)

pcoa_ind = pcoa(ind) ####pcoa计算

pc = as.data.frame(pcoa_ind$vectors) #### 将pcoa变换后的数据读入数据框

colname

r语言绘制精美pcoa图_[R语言 学习笔记]用R做主坐标分析(PCoA)相关推荐

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

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

  2. R语言多层桑基图_绘制炫酷桑基图,R语言(networkD3包)能搞定?

    你知道如何使用R语言绘制桑基图吗? 小锐当然知道啦,并且教你 绘制酷炫的桑基图. 开课啦 桑基图(Sankey)是一种特定类型的流程图,图中延伸的分支的宽度对应数据流量的大小,其最明显的特征就是,始末 ...

  3. R语言绘制核密度图实战(Kernel Density Plot)

    R语言绘制核密度图实战(Kernel Density Plot) 目录 R语言绘制核密度图实战(Kernel Density Plot) #仿真数据

  4. R语言绘制气泡矩阵图

    R语言绘制气泡矩阵图 示例图 之前在一些文章中看到过气泡矩阵的表达方法,该图形表达的意思就是不同样本中不同物种的丰度分布情况,气泡越大则是代表物种的相对丰度(或者说16S得到的绝对丰度)越大,在这个例 ...

  5. python中数据用折线图表示_使用PyQtGraph进行Python数据可视化:绘制精美折线图(以 上证指数走势为例)...

    在前两篇文章中,我们介绍了: 在了解了基本的PyQtGraph模块绘制图形功能之后,我们通过几个常用常见的数据可视化图形来演示使用PyQtGraph进行Python数据可视化. 本篇,我们介绍使用Py ...

  6. R语言小白学习笔记3—R语言读取数据

    R语言小白学习笔记3-R语言读取数据 笔记链接 想说的话 学习笔记3-R语言读取数据 3.1 读取CSV文件 3.1.1 read_delim函数 3.1.2 fread函数 3.2 读取Excel数 ...

  7. Python语言入门这一篇就够了-学习笔记(十二万字)

    Python语言入门这一篇就够了-学习笔记(十二万字) 友情提示:先关注收藏,再查看,12万字保姆级 Python语言从入门到精通教程. 文章目录 Python语言入门这一篇就够了-学习笔记(十二万字 ...

  8. 0.《新概念51单片机C语言教程》(郭天祥)学习笔记

    看到之前的阅读量,有标题党的嫌疑,再接再厉,加油坚持自己的文风,相信会写的越来越好:) 继续更(一个人,一个记事本,一个PDF阅读器,装备全靠打) <新概念51单片机C语言教程>(郭天祥) ...

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

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

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

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

最新文章

  1. 景德镇陶瓷学院计算机协会,景德镇陶瓷学院计算机协会换届大会总结报告
  2. 独家 | 不同机器学习模型的决策边界(附代码)
  3. cordova 不打开浏览器_[Cordova inAppBrowser 在App内打开浏览器]
  4. 使用ANT打包Android应用
  5. AJAX之表单即时验证
  6. linux还原系统_怎么成为一名合格的Linux运维工程师
  7. 诗与远方:无题(六十二)
  8. cbc系统是指_制动EBD,CBC是什么意思?
  9. 记录数据库内一条记录的更新时间
  10. MOE 分子动力学模拟操作
  11. 开心电视助手v3.8最新绿色版,安卓设备远程管理工具神器
  12. android手机蓝牙连接扫码枪,android 扫码枪解惑
  13. 12 款适用于开发人员的最佳 Web 开发软件
  14. 英语发音规则之26个字母发音规则(A字母)
  15. ZZULIOJ:1132: 数字字符统计
  16. 简历中使用STAR法则
  17. 国家二级计算机office知识点,2017年全国计算机二级office常见知识点
  18. LABS1000-01空盒气压计检定系统
  19. 读《IBM蓝色基因百年智慧》
  20. Java库:Jansi - 彩色日志输出体验

热门文章

  1. CASAIM自动化精密尺寸测量设备全尺寸检测铸件自动化检测铸件
  2. 计算机毕业设计Java基于的药品销售系统(源码+系统+mysql数据库+lW文档)
  3. 阿里云服务器购买及使用流程(新手入门图文教程)
  4. Transact-SQL语言基础语句
  5. Revit出图问题:打印机中新建纸张尺寸?批量导出图纸?
  6. 3D效果海报怎么制作设计?如何用PS制作出来!
  7. 酷睿i7 1165g7相当于什么水平 i71165g7属于哪个档次
  8. 好有作为的Java 程序员,年薪40W其实并不高。
  9. LVM Linear vs Striped Logical Volumes
  10. Afterlogic Aurora Corporate Crack