原始数据为遗传距离矩阵,使用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变换后的数据读入数据框</