peakcalling軟件MACS如何使用

peak calling軟件MACS如何使用，針對這個問題，這篇文章詳細介紹了相對應(yīng)的分析和解答，希望可以幫助更多想解決這個問題的小伙伴找到更簡單易行的方法。

創(chuàng)新互聯(lián)建站專注于化州網(wǎng)站建設(shè)服務(wù)及定制，我們擁有豐富的企業(yè)做網(wǎng)站經(jīng)驗。 熱誠為您提供化州營銷型網(wǎng)站建設(shè)，化州網(wǎng)站制作、化州網(wǎng)頁設(shè)計、化州網(wǎng)站官網(wǎng)定制、小程序定制開發(fā)服務(wù)，打造化州網(wǎng)絡(luò)公司原創(chuàng)品牌,更為您提供化州網(wǎng)站排名全網(wǎng)營銷落地服務(wù)。

MACS全稱是Model-based Analysis of ChIP-Seq，是使用的最廣泛的peak calling軟件之一，其基本原理簡介如下

1. 預(yù)測測序reads與peak的偏倚距離

在chip_seq等數(shù)據(jù)中科學(xué)家發(fā)現(xiàn)在真實的結(jié)合位點兩側(cè)，正負鏈的測序深度分布如下圖所示，對應(yīng)峰值的中心距離peak中心有一定的偏移

MACS首先通過一個模型來評估真實的peak中心和測序峰值的偏移距離，給定參數(shù)bandwidth和mfold, 采用一個大小為2倍bandwidth的滑動窗口，比較該窗口內(nèi)真實測序深度的分布與隨機測序的差異，如果二者的差異倍數(shù)超過了閾值mfold，則認為該窗口是一個peak區(qū)域。識別到初始的peak區(qū)域之后，隨機挑選1000個高可信度的peak區(qū)域，分別計算正鏈和負鏈的測序深度分布，示意如下

通過這種方式識別到正負鏈峰值之間的距離，定義為d。在后續(xù)peak calling時，會在初始計算結(jié)果的基礎(chǔ)上向3’端偏移d/2的距離。

2. peak 檢測

通常認為基因組上測序深度的分布是一個泊松分布，根據(jù)基因組測序深度的值，首先估算出整體泊松分布的參數(shù)

將初始peak區(qū)域偏移d/2的距離之后，再次已2d的滑動窗口進行滑動，尋找富集的區(qū)域，如果一個區(qū)域的測序分布的值大于整體泊松分布的值，則認為該區(qū)域是富集的。對于overlap的富集區(qū)域，進行合并，然后以peak中心區(qū)域開始，整體向3’端偏移d個堿基，其中測序深度最高的點作為peak的峰summit。

在比較input和IP樣本的差異時，MACS也是基于泊松分布來計算的 ，只不過考慮到測序錯誤，重復(fù)區(qū)域，CNV等因素的影響，MACS采用了一個動態(tài)的泊松分布模型，即不是認為基因組所有區(qū)域都符合同一個泊松分布，而是認為不同的區(qū)域泊松分布的參數(shù)可以不同。對于候選的peak區(qū)域，通過以下模型去估算泊松分布的參數(shù)

1k代表以peak中心為中心的1k窗口內(nèi)的數(shù)據(jù)估算出來的值，后面的5k, 10k類似。然后通過這個局部泊松分布的參數(shù)去計算每個peak區(qū)域的p值，如果小于指定的pvalue的閾值，則輸出該peak區(qū)域，該區(qū)域內(nèi)IP樣本的序列條數(shù)和局部泊松分布參數(shù)的比值作為該peak的fold enrichment。

關(guān)于peak calling軟件MACS如何使用問題的解答就分享到這里了，希望以上內(nèi)容可以對大家有一定的幫助，如果你還有很多疑惑沒有解開，可以關(guān)注創(chuàng)新互聯(lián)行業(yè)資訊頻道了解更多相關(guān)知識。

本文題目：peakcalling軟件MACS如何使用
URL地址：http://muchs.cn/article14/ghijde.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供微信小程序、App開發(fā)、建站公司、品牌網(wǎng)站設(shè)計、企業(yè)網(wǎng)站制作、營銷型網(wǎng)站建設(shè)

廣告

聲明：本網(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)