This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DiffBind normalization returned NULL

Hello,

I am trying to normalize some ChIPseq data in DiffBind, and having some difficulty with the command dba.normalize. When I try to examine the normalization methods and library size, I simply get NULL returned.

> B.DBcount
6 Samples, 1173 sites in matrix:
    ID Condition Treatment Replicate    Reads FRiP
1 wt1b        WT      none         1 16902874 0.20
2 mt1b        MT      none         1 14594474 0.13
3 wt2b        WT      none         2 16878697 0.22
4 mt2b        MT      none         2 14289354 0.15
5 wt3b        WT      none         3 18699283 0.24
6 mt3b        MT      none         3 15424246 0.20
> norm <- dba.normalize(B.DBcount, bRetrieve=TRUE)
> norm
NULL
> norm$norm.method
NULL

Any help greatly appreciated! Thanks

diffbind dba.normalize normalization chipseq

1 answer

You can't normalize and retrieve the result in the same call to dba.normalize(). In order to retrieve information about how normalization has been computed, you need to normalize the data first.

You can get the default normalization as follows:

B.DBcount <- dba.normalize(B.DBcount)
norm <- dba.normalize(B.DBcount, bRetrieve=TRUE)
norm

That makes sense - this worked.

Thank you!

Log in to answer this question.