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)