Plot method for nomatchfit objects — plot.nomatchfit • nomatch

Create a panel plot of cumulative incidence and effectiveness estimates across all evaluation time points.

Usage

# S3 method for class 'nomatchfit'
plot(
  x,
  effect = c("risk_ratio", "risk_difference", "relative_risk_reduction"),
  ci_type = NULL,
  color = "#0072B2",
  ...
)

Arguments

x: An object of class nomatchfit created by nomatch() or matching().
effect: The effect measure to plot next to the cumulative incidence plots. Either "risk_ratio"(default), "relative_risk_reduction" or "risk_difference".
ci_type: Character string specifying the type of confidence interval band to plot, one of "wald", "percentile", "simul", or "none". Must choose a ci_type whose lower and upper bounds are already computed in estimates component of x. Uses "wald" by default if available.
color: Aesthetic value to map data values to. Default: "#0072B2" (blue)
...: Additional arguments (currently ignored).

Value

a ggplot2 object with three faceted panels (for cumulative incidences and the chosen effect measure)

Details

The two cumulative incidence panels share y-axis limits. The effect measure panel has its own independent y-axis scaling.

Examples

fit <- nomatch(
 data = simdata,
 outcome_time = "Y",
 outcome_status = "event",
 exposure = "V",
 exposure_time = "D_obs",
 covariates = c("x1", "x2"),
 timepoints = seq(30, 180, by = 30),
 immune_lag = 14,
 boot_reps = 5
)
#> Bootstrapping 5 samples...
#> Bootstrap completed in 1.66 secs
plot(fit)