讀取 BED 檔
在這個練習中,你會從 BED 檔載入峰值呼叫,並利用峰值位置的資訊,從 BAM 檔中擷取與這些峰值重疊的讀取。
本練習屬於課程
在 R 中使用 Bioconductor 進行 ChIP-seq
練習說明
- 使用
import.bed函式,從「chr20_peaks」載入峰值呼叫。 - 使用這些峰值呼叫為「chr20_bam」建立
BamViews物件。 - 載入與峰值呼叫重疊的讀取。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Load peak calls from chr20_peaks
peaks <- ___(chr20_peaks)
# Create a BamViews object
bam_views <- ___(___, bamRanges=peaks)
# Load the reads
reads <- ___(___)