haploview如何繪制曼哈頓圖

這篇文章給大家分享的是有關(guān)haploview如何繪制曼哈頓圖的內(nèi)容。小編覺得挺實用的，因此分享給大家做個參考，一起跟隨小編過來看看吧。

我們提供的服務(wù)有：成都網(wǎng)站制作、成都網(wǎng)站設(shè)計、外貿(mào)營銷網(wǎng)站建設(shè)、微信公眾號開發(fā)、網(wǎng)站優(yōu)化、網(wǎng)站認(rèn)證、隆安ssl等。為近千家企事業(yè)單位解決了網(wǎng)站和推廣的問題。提供周到的售前咨詢和貼心的售后服務(wù)，是有科學(xué)管理、有技術(shù)的隆安網(wǎng)站制作公司

對于GWAS分析的結(jié)果，最常見的可視化手段就是曼哈頓圖了。一個典型的曼哈頓圖示例如下

圖中的每個點代表一個SNP位點，橫坐標(biāo)是SNP位點在染色體上的位置，縱坐標(biāo)是關(guān)聯(lián)分析計算出的p值。為了美觀，不同染色體上的SNP位點用不同顏色標(biāo)識?？v軸上的虛線代表閾值，p值越小，對應(yīng)到縱軸的值越大。通常在閾值上方的SNP位點認(rèn)為是和疾病或者形狀相關(guān)聯(lián)的位點。

曼哈頓圖本質(zhì)上就是一個散點圖，繪制方法比較簡單，有很多的R包可以直接出圖，本篇主要介紹haploview軟件。 這個軟件可以直接讀取plink關(guān)聯(lián)分析的結(jié)果，然后繪制曼哈頓圖。

在之前介紹了plink關(guān)聯(lián)分析的方法，根據(jù)輸入的hapmap1.ped和hapmap1.map兩個文件，最終可以得到一個allel水平的關(guān)聯(lián)分析結(jié)果，文件名為as1.assoc。haploview需要讀取hapmap1.map和as1.assoc兩個文件。

1. 導(dǎo)入數(shù)據(jù)

啟動haploview之后，選擇PLINK Format選項，Results File 選擇關(guān)聯(lián)分析的結(jié)果文件，可以是allel關(guān)聯(lián)分析的后綴為.assoc的文件，也可以是包含了genotype關(guān)聯(lián)分析的后綴為.model的文件；Map file 選擇后綴為.map的文件，當(dāng)然也可以選擇后綴為.bim的二進(jìn)制的map文件，選擇好之后，點擊OK即可。

2. 篩選數(shù)據(jù)

導(dǎo)入成功之后，可以看到如下所示的表格，這個表格里就是導(dǎo)入的關(guān)聯(lián)分析的結(jié)果。紅色方框標(biāo)識的區(qū)域可以對關(guān)聯(lián)分析的結(jié)果進(jìn)行篩選，可以通過Chr, Start Kb, End kb對染色體區(qū)域進(jìn)行篩選，也可以通過Filter對符合要求的SNP位點進(jìn)行篩選。

3. 繪圖

篩選好之后，點擊Plot按鈕，設(shè)置繪圖的參數(shù)。X軸肯定是Chromosomes, Y軸通常選擇P值，如果你有進(jìn)行多重假設(shè)檢驗的校正，也可以選擇校正之后的P值。scale表示對y軸的值進(jìn)行處理，有一下3種處理

1. untransformed

2. -log10

3. ln

一般情況下選擇-log10。Export to SVG可以導(dǎo)出圖片，格式為svg。

設(shè)置好之后，點擊OK就可以出圖了。示例如下

由于測試數(shù)據(jù)的SNP位點比較少，最后的效果圖并不是特別美觀。

感謝各位的閱讀！關(guān)于“haploview如何繪制曼哈頓圖”這篇文章就分享到這里了，希望以上內(nèi)容可以對大家有一定的幫助，讓大家可以學(xué)到更多知識，如果覺得文章不錯，可以把它分享出去讓更多的人看到吧！

新聞標(biāo)題：haploview如何繪制曼哈頓圖
標(biāo)題路徑：http://www.muchs.cn/article6/gceiig.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供網(wǎng)站營銷、企業(yè)網(wǎng)站制作、虛擬主機(jī)、網(wǎng)站維護(hù)、網(wǎng)站改版、網(wǎng)站策劃

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時需注明來源： 創(chuàng)新互聯(lián)