Marginal likelihood estimation of negative binomial parameters with applications to RNA-seq data.