Why is there such a mismatch between the Model's predicted probability and theoretical probability in logistic regression?
It's probably due to the effect coding of categorical predictors. Eg, the regression coefficient for CPR = 7 is not zero but -(sum of regression coefficients for the other 6 levels). I guess EM should have an option to switch it to reference coding, then your way of computing the predicted probability should work.