I updated my DESeq recently, but got problem when assigning values (No. of total mapped read) to the size factor.
here is my code:
The R has warnings like this:
1: In log(ifelse(y == 0, 1, y/mu)) : NaNs produced
2: step size truncated due to divergence
Does anyone encounter the same problem or know the reason.
Thanks a lot!
here is my code:
Code:
libsizes<-c( wt_rep1=18303917, wt_rep2=18947532, mu_rep1=20295971, mu_rep2=18256798) sizeFactors(cds)<-libsizes cds <- estimateDispersions( cds )
1: In log(ifelse(y == 0, 1, y/mu)) : NaNs produced
2: step size truncated due to divergence
Does anyone encounter the same problem or know the reason.
Thanks a lot!
Comment