Diffbind-weird Fold value
1
0
Entering edit mode
7 months ago
wei.zhang • 0

Hi all, I am new to Diffbind, when I used it to ran my ChIP-seq data, I got below result. From the dba.report if my understanding is correct Fold=Conc_GNP1-Conc_GNP2, but for my data is obviously not, anyone can help me with this? Thanks a lot!

result<-dba.report(Epigenome2,th=1)
result
GRanges object with 59140 ranges and 6 metadata columns:

seqnames            ranges strand |      Conc Conc_GNP1 Conc_GNP2         Fold     p-value       FDR
           <Rle>         <IRanges>  <Rle> | <numeric> <numeric> <numeric>    <numeric>   <numeric> <numeric>

chr20 16984425-16985425      * |   5.53053   5.85727   5.10743   0.00212363 0.000522680  0.999946
ChIP-Seq • 291 views
ADD COMMENT
2
Entering edit mode
7 months ago
Rory Stark ★ 1.2k

In the latest version of DiffBind, the default is to use fold changes calculated by the underlying analysis package. The default analysis package DESeq2 does a shrinkage estimation of fold change; the value in the report reflects this value, and not a simple calculation based on the sample group concentrations.

ADD COMMENT
0
Entering edit mode

Thanks a lot Rory!!!

ADD REPLY

Login before adding your answer.

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