plot_marker: Plots a Markers Mean Curve with Pointwise Quantiles

In VAJointSurv: Variational Approximation for Joint Survival and Marker Models

Plots a Markers Mean Curve with Pointwise Quantiles

Description

Plots a Markers Mean Curve with Pointwise Quantiles

Usage

plot_marker(
  time_fixef,
  time_rng,
  fixef_vary,
  x_range,
  vcov_vary,
  p = 0.95,
  xlab = "Time",
  ylab = "Marker",
  newdata = NULL,
  ...
)

Arguments

time_fixef	the time-varying fixed effects. See .e.g. poly_term.
time_rng	the time-varying random effects. See .e.g. poly_term.
fixef_vary	fixed effect coefficients for time_fixef.
x_range	2D numeric vector with start and end points.
vcov_vary	the covariance matrix for time_rng.
p	coverage of the two quantiles.
xlab, ylab, ...	arguments passed to plot.
newdata	data.frame with data for the weights if any.

Value

A list containing data for plotting.

Examples

# load in the data
library(survival)
data(pbc, package = "survival")

# re-scale by year
pbcseq <- transform(pbcseq, day_use = day / 365.25)
pbc <- transform(pbc, time_use = time / 365.25)

# create the marker terms
m1 <- marker_term(
  log(bili) ~ 1, id = id, data = pbcseq,
  time_fixef = bs_term(day_use, df = 5L),
  time_rng = poly_term(day_use, degree = 1L, raw = TRUE, intercept = TRUE))

fixef_vary <- c(-0.1048, 0.2583, 1.0578, 2.4006, 2.9734)
vcov_vary <- rbind(c(0.96580, 0.09543), c(0.09543,  0.03998))

# plot marker's trajectory
plot_marker(
  time_fixef = m1$time_fixef,
  time_rng = m1$time_rng,
  fixef_vary = fixef_vary,
  vcov_vary = vcov_vary, x_range = c(0,5))

VAJointSurv documentation built on Aug. 14, 2022, 9:05 a.m.