adjusted_KM: Calculate the adjusted survival probability

In Lesly1031/AdjKM.CIF: A package for generating covariate-adjusted Kaplan-Meier and cumulative incidence functions

Calculate the adjusted survival probability

Description

Use input data, time, status, grouping variables, adjusted covariates, whether to use stratified model, and reference group as inputs

Usage

adjusted_KM(
  data,
  time,
  status,
  group,
  covlist,
  stratified_cox = "Yes",
  reference_group = "G&B"
)

Arguments

data	the input dataset
time	column name of time variable
status	column name of event status
group	grouping variable
covlist	list of covariates that should be included in the model
stratified_cox	"Yes" refers to use stratified model, "No" refers to use Coxph
reference_group	NULL- unstratified cox when stratified = No; "G&B"- G&B when stratified = Yes; Otherwise, Storer's approach will be performed when using a self-defined reference

Value

Output is a dataframe with adjusted survival probabilities. If the PH assumption is invalid or if practitioners need a method by which the event time points of the adjusted function match those of the unadjusted function, the stratified model should be used (Gail and Byar and Storer et. al), otherwise, unstratified FG model can be used.

Examples

# Data preparation
library(KMsurv)
data(bmt)
bmt$arm <- bmt$group
bmt$arm = factor(as.character(bmt$arm), levels = c("2", "1", "3"))
bmt$z3 = as.character(bmt$z3)
bmt$t2 = bmt$t2 * 12/365.25

# Cox model
result1 = adjusted_KM(data = bmt, time = "t2", status = "d3", group = "arm", covlist = c("z1",
"z3"), stratified_cox = "No", reference_group = NULL)

# Stratified cox: Gail&Byar's approach
result2 = adjusted_KM(data = bmt, time = "t2", status = "d3", group = "arm", covlist = c("z1",
"z3"), stratified_cox = "Yes", reference_group = "G&B")

# Stratified cox: Storer's approach
result3 = adjusted_KM(data = bmt, time = "t2", status = "d3", group = "arm", covlist = c("z1",
"z3"), stratified_cox = "Yes", reference_group = "arm:2")

Lesly1031/AdjKM.CIF documentation built on March 10, 2023, 12:02 p.m.