Fitting Discrete Distribution Models to Count Data

Provides functions for fitting discrete distribution models to count data. Included are the Poisson, the negative binomial and, most importantly, a new implementation of the Poisson-beta distribution (density, distribution and quantile functions, and random number generator) together with a needed new implementation of Kummer's function (also: confluent hypergeometric function of the first kind). Three different implementations of the Gillespie algorithm allow data simulation based on the basic, switching or bursting mRNA generating processes. Moreover, likelihood functions for four variants of each of the three aforementioned distributions are also available. The variants include one population and two population mixtures, both with and without zero-inflation. The package depends on the 'MPFR' libraries (<>) which need to be installed separately (see description at <>). This package is supplement to the paper "A mechanistic model for the negative binomial distribution of single-cell mRNA counts" by Lisa Amrhein, Kumar Harsha and Christiane Fuchs (2019) available on bioRxiv.


Reference manual

It appears you don't have a PDF plugin for this browser. You can click here to download the reference manual.


1.0.1 by Lisa Amrhein, a year ago

Browse source code at

Authors: Lisa Amrhein [aut, cre] , Kumar Harsha [aut] , Christiane Fuchs [aut] , Pavel Holoborodko [ctb] (Author and copyright holder of 'mpreal.h')

Documentation:   PDF Manual  

Task views: Probability Distributions

GPL-3 license

Imports Rcpp

Suggests knitr, rmarkdown, testthat

Linking to Rcpp

System requirements: gmp (>= 4.2.3), mpfr (>= 3.0.0)

See at CRAN