In the frequentist approach that rms used, it is usually straightforward to use the maximum likelihood (ML) estimator as a point estimator for nomogram - as there is no other choice/everything is Normal in frequentist world so there is no difference whatsoever.

However, in the Bayesian context, there are not just ML estimator: posterior mean, mode, and median, all are point estimators. I think the current approach in rmsb is the marginal posterior mode. However, as pointed out by Florian Hartig the joint mode might not correspond to the marginal one. That makes me wonder if it is really the optimal option.

That said, to minimize the error, shouldn’t we aim for the transformed mean of the mortality risk (i.e. the mean of inv_logit of the log-odd), rather than the MAP/joint Mode?