EM using Matlab algorithm for matrix exponential in combination with Armadillo
EMstep_PADE(h, alpha, S, obs, weight, rcens, rcweight)
nuisance parameter
initial probalities
sub-intensity
the observations
the weights for the observations
censored observations
the weights for the censored observations