IBS距離矩陣

IBS是指兩個(gè)個(gè)體中共有的等位基因序列相同,IBS的遺傳距離可以衡量樣本間的相似性,評(píng)價(jià)其親緣關(guān)系。

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

#利用pheatmap可視化IBS距離矩陣

所需的數(shù)據(jù)格式為plink.mdist中第一行加入樣本ID

library(pheatmap)
library(RColorBrewer)

###讀取matrix數(shù)據(jù)
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,)
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容