EMstep_mdph function

EM for multivariate discrete phase-type