普氏分析在生信中的应用

普氏分析(Procrustes Analysis)

在微生物群落研究的过程中，我们经常需要评估微生物群落结构与环境因子整体之间是否具有显著的相关性，此时，通常使用的方式是Mantel test和普氏分析。

普鲁克分析(Procrustes Analysis) 又名普氏分析，是一种用来分析形状分布的统计方法。应用在数据分析中，可以理解为比较两组数据一致性的方法，主要用于表示样品不同方面的数据关联度。在宏基因组测序中， Procrustes分析常用于解释细菌组成与耐药基因的相关性，细菌与功能基因的相关性；在转录组测序中， Procrustes分析可用于解释基因与表型的相关性等；在物种分类中，普氏分析可以解释不同基因对物种鉴定的一致性。

在这里，我用Procrustes Analysis分析了地表水和地下水微生物群落的相似性和一致性。

library(vegan)
groundwater<-read.table("clipboard",header = T)
surface<-read.table("clipboard",header=T)
groundwater.dist<-vegdist(groundwater,method = "bray")
surface.dist<-vegdist(surface,method = "bray")

mds.ground<-monoMDS(groundwater.dist)
mds.surface<-monoMDS(surface.dist)

pro.g.s<-procrustes(mds.ground,mds.surface)
pro.g.s
protest(mds.ground,mds.surface)

Y <- cbind(data.frame(pro.g.s$Yrot), data.frame(pro.g.s$X))
X <- data.frame(pro.g.s$rotation)

Y$ID <- rownames(Y)

p <- ggplot(Y) +
  geom_segment