在 R 中使用 Bioconductor 進行 ChIP-seq
Peter Humburg
Statistician, Macquarie University
載入 peak 呼叫結果
peaks <- import.bed(peak_file, genome="hg19")
建立 BamViews 物件
bam_views <- BamViews(bam_file, bamRanges=peaks)
載入讀段
reads <- readGAlignments(bam_views)



在 R 中使用 Bioconductor 進行 ChIP-seq