diffbind : Is it possible to analyze one sample only ?
1
0
Entering edit mode
ZheFrench ▴ 40
@zhefrench-11689
Last seen 15 months ago
France

I have a pair of chipseq. One is really bad. Mapping failed.

I was wondering if it's possible to analyse only one without the pair ? 

Other question :

I tried to do the analyze with both of the samples even is one is not very good.

I did two different things : the first time I set the rule that peaks must be in both replicates...that gaves me 1116 differential peaks.

dba_chip_consensus <- dba.peakset(dba_chip, consensus=c(type),minOverlap=0.99)

Then I apply the following rule : the peaks must be at least in one or the other replicate. It gave me 1000 differential peaks.

dba_chip_consensus <- dba.peakset(dba_chip, consensus=c(type),minOverlap=0.01)

Strange, no ? I was hoping that it gaves me more peaks using one or the other... Is that possible ?

 

diffbind • 689 views
ADD COMMENT
0
Entering edit mode
Rory Stark ★ 5.1k
@rory-stark-5741
Last seen 8 weeks ago
Cambridge, UK

When you include all the peaks (from either replicate), it is likely that most of the additional peaks are unique to one replicate, so they won't be identified as differential due to high variance (and inadequate power). At the same time, adding these peaks changes the distribution of p-values and hence the calculation of FDR, making it less likely that any specific peak will be below your threshold. So what you are seeing looks right to me!

-Rory

ADD COMMENT

Login before adding your answer.

Traffic: 784 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