ChIP-seq with Bioconductor in R
Peter Humburg
Statistician, Macquarie University

คำถามที่ต้องการหาคำตอบ:




qc_result <- ChIPQC("sample_info.csv", "hg19")
counts <- dba.count(qc_results, summits=250)
plotPrincomp(counts)
คำนวณระยะห่างระหว่างตัวอย่าง
distance <- dist(t(coverage))
สร้าง Dendrogram
dendro <- hclust(distance)
พล็อต Dendrogram
plot(dendro)
สร้าง Heatmap จากข้อมูล Coverage ของ Peak
dba.plotHeatmap(peaks, maxSites = peak_count, correlations = FALSE)

ChIP-seq with Bioconductor in R