Compare Ratio Models

Compare multiple ratio models using LOO-CV or WAIC.

Usage

ratiod_compare(..., criterion = c("loo", "waic"))

Arguments

...

Multiple ratiod_fit objects to compare

criterion

Comparison criterion: "loo" (default) or "waic"

Value

A loo_compare object

Examples

# \donttest{
sim <- sim_ratiod(n = 100, family = ratiod_poisson_gamma(), seed = 1)
fit <- tratio(y_num | y_denom ~ x1 + (1 | group), data = sim$data,
              family = ratiod_poisson_gamma(),
              control = list(iter = 1000, chains = 2, seed = 1))
fit0 <- tratio(y_num | y_denom ~ 1 + (1 | group), data = sim$data,
               family = ratiod_poisson_gamma(),
               control = list(iter = 1000, chains = 2, seed = 1))
ratiod_compare(fit0, fit)
#> Warning: Some Pareto k diagnostic values are too high. See help('pareto-k-diagnostic') for details.
#>  model elpd_diff se_diff p_worse diag_diff       diag_elpd
#>    fit       0.0     0.0      NA                          
#>   fit0     -39.3    12.1    1.00           1 k_psis > 0.67
#> 
#> Diagnostic flags present.
#> See ?`loo-glossary` (sections `diag_diff` and `diag_elpd`)
#> or https://mc-stan.org/loo/reference/loo-glossary.html.
# }