Question: Calculate enrichment values of Chip-seq peaks for hierarchical clustering
gravatar for florian.noack
3.7 years ago by
florian.noack20 wrote:

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.


chip-seq • 1.8k views
ADD COMMENTlink modified 3.6 years ago by igor11k • written 3.7 years ago by florian.noack20
gravatar for igor
3.6 years ago by
United States
igor11k wrote:

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 COMMENTlink written 3.6 years ago by igor11k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1704 users visited in the last hour