DiffBind custom normalization factor
2
0
Entering edit mode
C T ▴ 140
@c-t-5858
Last seen 18 months ago
United States

I've searched through the manual and bioconductor support, there doesn't seem to be a way to specify a custom normalization factor in DiffBind when analyzing using DESeq2. I'm wondering if there is a relatively easy "hacky" way to do this? Thanks in advance for your help.

diffbind • 865 views
ADD COMMENT
0
Entering edit mode
Rory Stark ★ 5.2k
@rory-stark-5741
Last seen 28 days ago
Cambridge, UK

I'm expecting to add a straightforward way to supply normalization factors explicitly before the next release.

In the meantime, there is a hacky way to do this. When you are using the defaults method=DBA_DESEQ2 and bFullLibrarySize=TRUE, the normalization factors are set as follows:

libsize <- as.numeric(DBA$class['Reads',])
sizeFactors <- libsize/min(libsize)

You can change the values of DBA$class['Reads',] to whatever you like so that the final sizeFactors end up how you want them.

ADD COMMENT
0
Entering edit mode
Rory Stark ★ 5.2k
@rory-stark-5741
Last seen 28 days ago
Cambridge, UK

The latest version of DiffBind has a new interface function, dba.normalize(), that allows vectors of normalization factors and library sizes (or a matrix of normalization offsets) to be specified instead of using the internal normalization calculations.

ADD COMMENT

Login before adding your answer.

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