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.

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.

Usage

dirmnom_fit(k, mm, conc)

Arguments

k

A matrix with the data to fit a dirichlet multinomial model.

mm

A model matrix with the features.

conc

The concentration parameter that accounts for overdispersion. To estimate this value use the estimate_overdispersion() function.

Value

A list with beta = a matrix with the best set of parameters found; and ll = the log-likelihood.