IBS距离矩阵

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,)
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容