這篇文章主要為大家展示了“大數(shù)據(jù)中如何繪制曼哈頓圖”,內(nèi)容簡而易懂,條理清晰,希望能夠幫助大家解決疑惑,下面讓小編帶領大家一起研究并學習一下“大數(shù)據(jù)中如何繪制曼哈頓圖”這篇文章吧。
站在用戶的角度思考問題,與客戶深入溝通,找到婺源網(wǎng)站設計與婺源網(wǎng)站推廣的解決方案,憑借多年的經(jīng)驗,讓設計與互聯(lián)網(wǎng)技術結合,創(chuàng)造個性化、用戶體驗好的作品,建站類型包括:成都網(wǎng)站建設、做網(wǎng)站、企業(yè)官網(wǎng)、英文網(wǎng)站、手機端網(wǎng)站、網(wǎng)站推廣、主機域名、網(wǎng)頁空間、企業(yè)郵箱。業(yè)務覆蓋婺源地區(qū)。
曼哈頓圖的命名得益于其形狀,和紐約市曼哈頓區(qū)鱗次櫛比的大樓非常相近,曼哈頓區(qū)是摩天大樓最多的城市,標志性的景觀如下
曼哈頓圖示意如下
每條染色體可以看做是一座高樓,整體看上去形似曼哈頓區(qū)的摩天大樓,所以稱之為曼哈頓圖。
了解了曼哈頓圖的命名,再來看下它所展示的信息。從上圖可以看出,曼哈頓圖的x軸為snp位點在染色體上的位置,y軸為SNP位點對應的p值。從本質上來看,屬于散點圖。有x軸和y軸的坐標,人人都可以畫散點圖,那為何曼哈頓圖的繪制看上去并非如此簡單呢?
秘密就在于SNP位點的染色體位置,對于不同的染色體,在標記位置時都是從1開始計數(shù)的,比如1號和2號染色體的12bp對應的是兩個snp位點。通過染色體和位置兩列信息相結合, 可以唯一表征snp位點的位置。
然而在繪制散點圖時,x軸只能用1個數(shù)值來表示。在曼哈頓圖中,將染色體線性排列,1號染色體的位置不變,2號染色體的位置在原來的基礎上再加上的1號染色體的總長度,然后依次類推,通過這樣一種形式,將原始輸入文件中染色體和位置兩列信息,轉換為繪圖所用的x軸坐標信息。
在實際分析中,通過qqman
這個R包可以來實現(xiàn)曼哈頓圖的繪制,用法如下
輸出結果如下所示
可以看到,只需要準備好同樣格式的輸入文件,繪制曼哈頓圖就是幾秒鐘的事情。曼哈頓圖展示了每個SNP位點的p值信息,圖中的兩個水平橫線對應篩選顯著位點的閾值,一個對應的p值是1e-5,另一個對應的p值是5e-8。
當然該R包還提供了很多的功能,比如高亮標記指定的snp位點等,改變?nèi)旧w顯示的名稱等等,能夠滿足基本需求,但是還是由許多限制,比如不能展示性染色體。該包的源代碼中會自動對數(shù)據(jù)過濾,只保留常染色體上的SNP位點。
使用別人的R包就是這個樣子,別人給你什么,你用什么,想要個性化,要么通知R包的開發(fā)者新增功能,當然人家不一定會買賬,要么自己來。
將qqman中的核心源代碼部分截取出來,如下
這里省略了顏色的賦值,可以看到,大部分的工作是在計算snp位點對應的x軸坐標,畫圖的代碼只有寥寥幾句。
理解了曼哈頓圖的本質,就可以自已用R或者熟悉的軟件來定制曼哈頓圖。
以上是“大數(shù)據(jù)中如何繪制曼哈頓圖”這篇文章的所有內(nèi)容,感謝各位的閱讀!相信大家都有了一定的了解,希望分享的內(nèi)容對大家有所幫助,如果還想學習更多知識,歡迎關注創(chuàng)新互聯(lián)行業(yè)資訊頻道!
本文標題:大數(shù)據(jù)中如何繪制曼哈頓圖
URL標題:http://jinyejixie.com/article20/pgeejo.html
成都網(wǎng)站建設公司_創(chuàng)新互聯(lián),為您提供品牌網(wǎng)站設計、網(wǎng)站收錄、企業(yè)建站、網(wǎng)站內(nèi)鏈、移動網(wǎng)站建設、服務器托管
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉載內(nèi)容為主,如果涉及侵權請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉載,或轉載時需注明來源: 創(chuàng)新互聯(lián)