Question: sicer chip seq analysis
gravatar for robertacar
5.2 years ago by
robertacar0 wrote:

Hi everyone,

I'm new in ChiP-seq analysis. I'm using SICER to analyze histone modification data with a control library. I have some doubts regarding the output file. In the output file (summary of filtered islands with an FDR < 0,01) there are some peaks which appear to be more enriched in the control library respect to the chip one. Is it normal? I thought I would have peaks only enriched in chip library. How is the fold change value calculated? Does it exist a best used threashold to filter this value?


Thank you

chip-seq • 3.4k views
ADD COMMENTlink modified 3.3 years ago by ariel.balter140 • written 5.2 years ago by robertacar0
gravatar for ayyappakumar.s
5.2 years ago by
ayyappakumar.s20 wrote:

General behaviour is we see control samples are not much enriched when compared to the ChiP peak regions. Ensure that you have equal number of replicates (biological replicate if you are merging them before peak detection). If you have unequal biological replicates :

for ex: Control sample has two biological replicates and Test sample has single replicate then it may cause such differences. So, in this case try tp run peak detection by considering individual biological replicate.

In such case, then check the FDR values for the enriched peaks in the control samples.

General way of Fold change calculation is by :  Control Sample divided by ChIP Sample (or) vice versa

ADD COMMENTlink written 5.2 years ago by ayyappakumar.s20
gravatar for ariel.balter
3.3 years ago by
ariel.balter140 wrote:

You might want to look into SICER's new incarnation in EPIC. You will find it much more accessible. And the creator is very responsive to issues on GitHub and answering questions if they are posted here or on Stackoverflow.

ADD COMMENTlink written 3.3 years ago by ariel.balter140
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: 879 users visited in the last hour