PH Diagnosis Plot of Log Cumulative Hazard Rate versus time or log-time

This plot is also known as log negative log survival rate.

Usage

ph_diagplot_lch(
  km_fit,
  time_scale,
  log_time = TRUE,
  endpoint_name = "",
  subtitle = "",
  exclude_censor = TRUE
)

Arguments

km_fit

returned object from survival::survfit

time_scale

a character string, 'years', 'months', 'weeks' or 'days', time unit of median survival time

log_time

logical, TRUE (default) or FALSE

endpoint_name

a character string, name of the endpoint

subtitle

a character string, subtitle of the plot

exclude_censor

logical, should censored data point be plotted

Value

a plot of log cumulative hazard rate

Details

a diagnosis plot for proportional hazard assumption, versus log-time (default) or time

Examples

library(survival)
data(adtte_sat)
data(pseudo_ipd_sat)
combined_data <- rbind(adtte_sat[, c("TIME", "EVENT", "ARM")], pseudo_ipd_sat)
kmobj <- survfit(Surv(TIME, EVENT) ~ ARM, combined_data, conf.type = "log-log")
ph_diagplot_lch(kmobj,
  time_scale = "month", log_time = TRUE,
  endpoint_name = "OS", subtitle = "(Before Matching)"
)