ChIP-seq med Bioconductor i R
Peter Humburg
Statistician, Macquarie University

Frågor vi vill besvara:




qc_result <- ChIPQC("sample_info.csv", "hg19")
counts <- dba.count(qc_results, summits=250)
plotPrincomp(counts)
Beräkna avstånd mellan sampel
distance <- dist(t(coverage))
Skapa dendrogram
dendro <- hclust(distance)
Rita dendrogram
plot(dendro)
Skapa heatmap från täckningsdata för toppar
dba.plotHeatmap(peaks, maxSites = peak_count, correlations = FALSE)

ChIP-seq med Bioconductor i R