IBS距离矩阵

作者: 千万别加香菜 | 来源:发表于2022-02-24 19:40 被阅读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,)
    

    相关文章

      网友评论

        本文标题:IBS距离矩阵

        本文链接:https://www.haomeiwen.com/subject/eipstrtx.html