When bs(varname) has been used in the formula of a coxph model, this function can be used to reconstruct the predicted relative risk of that parameter over time.
# S3 method for coxph
reconstruct.bs(
mod,
varname,
fun = NULL,
probs = c(0.025, 0.975),
bw = FALSE,
xlab = NULL,
ylab = NULL,
plot = TRUE,
...
)median, upper and lower confidence bands for the effect of varname over time; the funciton also produces a plot.
model output, created by function survspat
name of the variable modelled by a B-spline
optional function to feed in. Default is to plot relative risk against the covariate of interest. Useful choices include "identity" (but with no quotes), which plots the non-linear effect on the scale of the linear predictor.
upper and lower quantiles for confidence regions to plot> The default is c(0.025,0.975).
Logical. Plot in black/white/greyscale? Default is to produce a colour plot. Useful for producing plots for journals that do not accept colour plots.
label for x axis, there is a sensible default
label for y axis, there is a sensible default
logical, whether to plot the effect of varname over time
other arguments to be passed to the plotting function.