polyaMLE {fabPrediction}

Obtain MLE of marginal Dirichlet-multinomial likelihood

Description

This function retuns the MLE of the prior concentration from a marginal Dirichlet-multinomial likelihood. Default method iterates a Newton-Raphson algorithm until convergence.

Usage

polyaMLE(
  D,
  init = NA,
  method = "Newton_Raphson",
  epsilon = 1e-04,
  print_progress = FALSE
)

Arguments

D

matrix (JxK) of counts; each row is a sample from a MN distribution with K categories

init

If NA, use method moment matching procedure to obtain good init values

method

"Newton_Raphson", "fixed_point", "separate", "precision_only"

epsilon

convergence diagnostic

print_progress

if TRUE, print progress to screen

Value

mle of prior concentration from marginal Dirichlet-multinomial likelihood