I have count-data from RNA-seq of miRNAs. I am trying to use Deseq() in order to find the miRNAs that are DE between control group and patients, but no matter what I do, some rows do not converge. I even replaced the call to Deseq() with the following lines:
dds <- estimateSizeFactors(dds)
nc <- counts(dds, normalized=TRUE)
filter <- rowSums(nc >= 10) >= 2
dds <- dds[filter,]
dds <- estimateDispersions(dds)
dds <- nbinomWaldTest(dds, maxit=500000)
sessionInfo( )
still, I get this error message: "6 rows did not converge in beta, labelled in mcols(object)$betaConv. Use larger maxit argument with nbinomWaldTest". The run is really slow as it is, should I really choose a larger maxit?
