ピークの統合
ピークコールへのアノテーション方法を探る前に、すべてのサンプルのコールを取り込んだマージ済みピークセットを取得する方法を確認しておきましょう。パッケージ ChIPpeakAnno には、そのための関数 findOverlapsOfPeaks() が用意されています。これは、最大5サンプルまでの GenomicRanges で表現されたピークセットをマージできます。
この関数を使う前に、先に作成した ChIPQCexperiment オブジェクトからピークコール情報を抽出する必要があります。ピークの位置は peaks() 関数で取得できます。これは、QC を通過しなかったものも含め、すべてのサンプルのピークコールを返す点に注意してください。QC の段階で今後の解析対象として選択したサンプルのインデックスを含むベクターは qc_pass として利用できます。
この演習はコースの一部です
Rで学ぶBioconductorによるChIP-seq
演習の手順
- ChIPQCexperiment オブジェクトからピークを抽出します。
- QC に不合格のサンプルをすべて除外します。
findOverlapsOfPeaks()関数でピークセット間のオーバーラップを探索します。mergedPeaksエントリとして利用できるマージ済みピークセットを確認します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Extract peaks from ChIPQCexperiment object
peak_calls <- ___(ar_calls)
# Only keep samples that passed QC
peak_passed <- ___[qc_pass]
# Find overlaps between peak sets
peaks_combined <- ___(peak_passed[[1]], peak_passed[[2]], peak_passed[[3]], peak_passed[[4]], maxgap=50)
# Examine merged peak set
print(___)