This is a function to compute posterior probability of the
crossover interaction. It takes the output of
do_bms as input. summary must be set to
FALSE when calling do_bms. By default, the
joint conditional probabilities are returend. If prob
is set to "conditional", the conditional probability given each
model is returned. If the option is set to "marginal", a scalar
value is returned.
get_co(fit, prob = "joint")A named vector or scalar of posterior probability.