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.