NBumiFitDispVsMean {M3Drop} | R Documentation |
Fits a power-law relationship between mean expression and dispersion.
NBumiFitDispVsMean(fit, suppress.plot=TRUE)
fit |
output from NBumiFitModel or NBumiFitBasicModel. |
suppress.plot |
Whether to plot the calculated fit. |
Fits a power-law relationship between mean expression and fitted gene-specific dispersions, using linear regression on ln transformed values. Lowly expressed genes, mean expression < 2^4, are excluded as long as at least 2000 genes remain for fitting.
Coefficients of linear regression.
library(M3DExampleData) counts <- as.matrix(Mmus_example_list$data); counts <- counts[rowSums(counts) > 0,]; fit <- NBumiFitModel(counts); coeffs <- NBumiFitDispVsMean(fit, suppress.plot=TRUE);