Deviance Information Criterion (DIC)
Calculate the DIC for a model fitted using the nma() function.
dic(x, ...)
x |
A fitted model object, inheriting class stan_nma |
... |
Other arguments (not used) |
A nma_dic object.
print.nma_dic() for printing details, plot.nma_dic() for
producing plots of residual deviance contributions.
## Smoking cessation
# Run smoking FE NMA example if not already available
if (!exists("smk_fit_FE")) example("example_smk_fe", run.donttest = TRUE)
# Run smoking RE NMA example if not already available
if (!exists("smk_fit_RE")) example("example_smk_re", run.donttest = TRUE)
# Compare DIC of FE and RE models
(smk_dic_FE <- dic(smk_fit_FE))
(smk_dic_RE <- dic(smk_fit_RE)) # substantially better fit
# Plot residual deviance contributions under RE model
plot(smk_dic_RE)
# Check for inconsistency using UME model
# Run smoking UME NMA example if not already available
if (!exists("smk_fit_RE_UME")) example("example_smk_ume", run.donttest = TRUE)
# Compare DIC
smk_dic_RE
(smk_dic_RE_UME <- dic(smk_fit_RE_UME)) # no difference in fit
# Compare residual deviance contributions
plot(smk_dic_RE, smk_dic_RE_UME, show_uncertainty = FALSE)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.