CauchyCP {CauchyCP}R Documentation

A robust test under non-proportional hazards using Cauchy combination of change-point Cox regressions.

Description

A robust test under non-proportional hazards using Cauchy combination of change-point Cox regressions.

Usage

CauchyCP(time, status, x, covar = rep(1, length(time)),
  cutpoints = c(0, quantile(time[status == 1])[2:4]))

Arguments

time

- Follow up time for right censored data.

status

- The event status indicator, 0=censored, 1=event.

x

- The variable of interest, e.g. a treatment indicator.

covar

- The matrix of covariates. If no covariates, a vector of ones should be used (default).

cutpoints

- The pre-specified change-points. The default choice is a vector of 0th, 25th, 50th and 75th percentiles of the event time.

Value

1. A matrix of estimated hazard ratios before and after the change-points. 2. the vector of p-values corresponding to the change-points. 3. a final p-value.

References

Hong Zhang, Qing Li, Devan Mehrotra and Judong Shen. "CauchyCP: a powerful test under non-proportional hazards using Cauchy combination of change-point Cox regressions", arXiv:2101.00059.

Examples

data(gast)
CauchyCP(time=gast$time, status=gast$status, x=gast$trt)

[Package CauchyCP version 0.1.1 Index]