Added piece "%>% sapply(min, 0.999)" in M3 and M4 to avoid survival parameter > 1.
Also added this in 'init_calib' function.

All models tested and work.