coverage_plot {eesim}

Plot coverage of empirical confidence intervals

Description

Plots the relative risk point estimates and their confidence intervals for model fit results for each simulation, compared to the true relative risk. This gives a visualization of the coverage of the specified method for the relative risk. The confidence intervals which do not contain the true relative risk appear in red. The input to this function should be either the output of fit_mods or the second element of the output of eesim.

Usage

coverage_plot(summarystats, true_param)

Arguments

summarystats

A list or data frame of summary statistics from many repetitions of a simulation. Must include columns titled Estimate, lower_ci, and upper_ci. This could be the second object from the output of eesim, specified by using the format eesim_output[[2]].

true_param

The true value of the relative risk used to simulate the data.

Value

A plot displaying the coverage for the true value of the parameter by the confidence intervals resulting from each repetition of the simulation.

Examples

ex_sim <- eesim(n_reps = 100, n = 1000, central = 100, sd = 10,
          exposure_type = "continuous", average_outcome = 20, rr = 1.02,
          custom_model = spline_mod, custom_model_args = list(df_year = 1))
coverage_plot(ex_sim[[2]], true_param = 1.02)