wagner_plot function

Produces a Wagner plot for the specified choice of covariates and/or time slices.