PS.MA function

Calculate model probabilities for the propensity score model using a pseudo-MC3 algorithm