Reads simulation summary tables from multiple output directories (one per CI method) and produces PDF figures comparing CI methods across sample sizes. For each (model type x variance estimator x table) combination two PDFs are saved: an estimator comparison (Bias, MSE) and a CI comparison (SE calibration, coverage, width).
Usage
simCompareMethodsFigures(
output.dirs = c(Bootstrap = "resiBootSim", Normal = "resiAsympNormalSim", QF =
"resiAsympQFSim", CF = "resiAsympCFSim"),
figures.dir = NULL,
alpha = 0.05,
fixed.knots = FALSE
)Arguments
- output.dirs
Named character vector mapping CI method labels to their simulation output directories. Default:
c(boot = "resiBootSim", normal = "resiAsympNormalSim", qf = "resiAsympQFSim", cf = "resiAsympCFSim"). Directories that do not exist are silently skipped.- figures.dir
Character, directory where comparison figures are saved. Default:
file.path(output.dirs[1], "figures", "method_comparison").- alpha
Numeric, nominal CI level used for coverage reference lines. Default 0.05.
- fixed.knots
Logical. Must match the value used in the original
insurancePlasmodeSimcall. DefaultFALSE.
Value
Invisibly returns the combined summary data.frame. Saves
estimator_<model>_<vcov>_<table>.pdf,
compare_<model>_<vcov>_<table>.pdf, and
covquant_<model>_<vcov>_<table>.pdf to figures.dir.
