Add support functions for Mixture IS LOO

We have the Mix-IS-LOO vignette https://mc-stan.org/loo/articles/loo2-mixis.html, but it would be useful to add support for creating a loo object from mixture posterior log_lik. Assuming the Stan model code implements the mixture distribution and generates log_lik variable, Mix-IS-LOO elpd_loo estimate is given by

log_lik_mix <- extract(fit_mix)$log_lik
l_common_mix <- rowLogSumExps(-log_lik_mix)
log_weights <- -log_lik_mix - l_common_mix
elpd_mixis <- logSumExp(-l_common_mix) - rowLogSumExps(t(log_weights))

We can also compute Pareto-k's for each pointwise elpd, but there is a difference that the mixture posterior elpd itself is already different from the regular posterior, so the error of Mix-IS-LOO tend to be smaller than regular PSIS-LOO.

Read the Mix-IS-LOO vignette at https://mc-stan.org/loo/articles/loo2-mixis.html and use the issue's log_lik_mix, l_common_mix, log_weights, and elpd_mixis calculations as the behavioral reference. Done should provide support for creating a loo object from mixture posterior log_lik and compute Pareto-k values where applicable.