plot.rmtfit: Plot the estimated treatment effect curve

In rmt: Restricted Mean Time in Favor of Treatment

Description

Plot the estimated overall or stage-wise restricted mean times in favor of treatment as a function of follow-up time.

Usage

## S3 method for class 'rmtfit'
plot(
  x,
  k = NULL,
  conf = FALSE,
  main = NULL,
  xlim = NULL,
  ylim = NULL,
  xlab = "Follow-up time",
  ylab = "Restricted mean time in favor",
  conf.col = "black",
  conf.lty = 3,
  ...
)

Arguments

x	An object returned by rmtfit.
k	If specified, μ_k(τ) is plotted; otherwise, μ(τ) is plotted.
conf	If TRUE, 95% confidence limits for the target curve are overlaid.
main	A main title for the plot
xlim	The x limits of the plot.
ylim	The y limits of the plot.
xlab	A label for the x axis, defaults to a description of x.
ylab	A label for the y axis, defaults to a description of y.
conf.col	Color for the confidence limits if conf=TRUE.
conf.lty	Line type for the confidence limits if conf=TRUE.
...	Other arguments that can be passed to the underlying plot method.

Value

No return value, called for side effects.

See Also

rmtfit, summary.rmtfit, bouquet.

Examples

# load the colon cancer trial data
library(rmt)
head(colon_lev)
# fit the data
obj=rmtfit(ms(id,time,status)~rx,data=colon_lev)
# plot overal effect mu(tau)
plot(obj)
# set-up plot parameters
oldpar <- par(mfrow = par("mfrow"))
par(mfrow=c(1,2))
# Plot of component-wise RMT in favor of treatment over time
plot(obj,k=2,conf=TRUE,col='red',conf.col='blue', xlab="Follow-up time (years)",
    ylab="RMT in favor of treatment (years)",main="Survival")
plot(obj,k=1,conf=TRUE,col='red',conf.col='blue', xlab="Follow-up time (years)",
    ylab="RMT in favor of treatment (years)",main="Pre-relapse")
par(oldpar)

rmt documentation built on May 25, 2021, 9:06 a.m.