Plots VPC/ICC across the models compared by compare_maihda.
Dispatched via plot() on the classed result.
Usage
# S3 method for class 'maihda_comparison'
plot(x, ...)Arguments
- x
A
maihda_comparisonobject fromcompare_maihda.- ...
Additional arguments (not used).
Examples
# \donttest{
strata <- make_strata(maihda_sim_data, vars = c("gender", "race"))
null_model <- fit_maihda(health_outcome ~ 1 + (1 | stratum), data = strata$data)
adj_model <- fit_maihda(health_outcome ~ age + (1 | stratum), data = strata$data)
comparison <- compare_maihda(null_model, adj_model, bootstrap = TRUE)
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
#> boundary (singular) fit: see help('isSingular')
plot(comparison)
# }
