EMstep_bivdph function

EM for discrete bivariate phase-type