This function computes the probability for each observed choice at the
(normalized, burned and thinned) samples from the posterior. These
probabilities are required to compute the `WAIC`

and the
marginal model likelihood `mml`

.

## Usage

`compute_p_si(x, ncores = parallel::detectCores() - 1, recompute = FALSE)`