基因组分析植物群体遗传学科研绘图

IBS距离矩阵

2022-02-24  本文已影响0人  千万别加香菜

IBS是指两个个体中共有的等位基因序列相同,IBS的遗传距离可以衡量样本间的相似性,评价其亲缘关系。

/home/software/plink --allow-extra-chr --chr-set 27 --file ld.QC.75_noinclude0-502502-geno02-maf03 --distance-matrix

#利用pheatmap可视化IBS距离矩阵

所需的数据格式为plink.mdist中第一行加入样本ID

library(pheatmap)
library(RColorBrewer)

###读取matrix数据
my_data <- read.delim(file.choose("plink.mdist"))      
rownames(my_data) = colnames(my_data)
pheatmap(my_data, scale = "none",       
         mycolor <- brewer.pal(9,"PiYG"),       ########修改颜色    
      show_rownames = T, show_colnames = T, main = "Heatmap", border_color = NA,)
上一篇下一篇

猜你喜欢

热点阅读