plot_bas_hazard: Plot the Baseline Hazard Step-Function

In TimeDepFrail: Time Dependent Shared Frailty Cox Model

Plot the Baseline Hazard Step-Function

Description

This function plots the baseline hazard step-function based on the estimated parameters from the Adapted Paik et al.'s model.

Usage

plot_bas_hazard(
  result,
  xlim = c(min(result$TimeDomain), max(result$TimeDomain)),
  ylim = c(0, max(result$BaselineHazard)),
  xlab = "Time",
  ylab = "Values",
  main = "Baseline hazard step-function",
  color = "black",
  pch = 21,
  bg = "black",
  cex_points = 0.7
)

Arguments

result	S3 object of class 'AdPaik', returned by the method call 'AdPaikModel(...)'.
xlim	A numeric vector specifying the x-axis limits. Default is set to the interval min-max of the time-domain.
ylim	A numeric vector specifying the y-axis limits. Default is 0 to the maximum value of the baseline hazard.
xlab, ylab	String giving the x and y axis name. Default values are 'x' and 'y'.
main	Title of the plot. Default title is 'Baseline hazard step-function'.
color	Color used for plotting the horizontal segments of the step-function. Default one is 'black'.
pch	Symbol for marking the boundaries of each segment. Default is a dot (value 21).
bg	Color for the boundary symbols. Default matches the plot color ('black').
cex_points	Size of the boundary symbols. Default is 0.7.

Details

The function plots a horizontal segment for each interval of the time domain, representing the baseline hazard. The boundaries of each segment are marked with colored dots, and subsequent segments are intentionally left unconnected to reflect the discrete nature of the intervals.

Value

Plot of the baseline hazard step-function and value of the function in each interval.

Examples

# Import data
data(data_dropout)

# Define the variables needed for the model execution
eps_paik <- 1e-10
categories_range_min <- c(-8, -2, eps_paik, eps_paik, eps_paik)
categories_range_max <- c(-eps_paik, 0.4, 1 - eps_paik, 1, 10)
time_axis <- c(1.0, 1.4, 1.8, 2.3, 3.1, 3.8, 4.3, 5.0, 5.5, 5.8, 6.0)
formula <- time_to_event ~ Gender + CFUP + cluster(group)


# Call the main model function
result <- AdPaikModel(formula, data_dropout, time_axis, categories_range_min, categories_range_max)

plot_bas_hazard(result)

TimeDepFrail documentation built on April 11, 2025, 5:41 p.m.