mml function

Approximate marginal model likelihood