Function reference

All functions

ddirmnom()

calculates the PDF of a dirichlet-multinomial compound distribution for the data x

dirmnom_fit()

calculates parameter-optimized fit (including LL) using a dirichlet multinomial model given a matrix of data (k), a model matrix with the features (mm), and a concentration parameter.

estimate_overdispersion()

Estimates overdispersion in multinomial models

rough_dirmnom_ll()

calculates a parameter-non-optimized log-likelihood using a dirichlet multinomial model given a matrix of data (k), a model matrix with the features (mm), and a concentration parameter. This is useful for data sets that fail to converge the optimization process. e.g. those with not so many samples.