Calculate enrichment values of Chip-seq peaks for hierarchical clustering
1
0
Entering edit mode
4.3 years ago

Hi everyone, I have Chip-seq samples + input controls (triplicate) for three different cell types. I want to study the changes of the Chip-seq samples between the three different cell types. I called already peaks using MACS2 and generated a global peak file (if a genomic regions pops up in at least one of the three cell types its in this file). No I wondering how I can calculate for each sample (should I merge the replicates, they show a very high pearson correlation) a enrichment value which considers also the input control. Of course the enrichment value should be also normalized by the total readcount to be comparable between the different samples. I would like to feed this enrichment value into pheatmap for hierarchical clustering.

Thanks a lot.

Flo

ChIP-Seq • 2.1k views
ADD COMMENT
0
Entering edit mode
4.3 years ago
igor 12k

You should look into DiffBind or csaw. There are a few tools for differential binding that will merge, quantify, and normalize peaks across multiple samples/groups.

ADD COMMENT

Login before adding your answer.

Traffic: 1793 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6