marg_loglik function

Calculate marginal loglikelihood at high density points of a MAR model.