PHeval (version 1.1)

PHeval-package: Evaluation of the proportional hazards assumption with a standardized score process

Description

This package provides tools for the evaluation of the goodness of fit of the proportional hazards model.

Author

Cecile Chauvel <[e-mail]>

Details

Package:	PHeval
Type:	Package
Version:	1.1
Date:	2024-09-02
License:	GPL (>=2.0)

This package provides functions to compute and plot the standardized score process of OQuigley (2008). The mathematical properties of this process are studied in Chauvel (2014). The restrained adaptive test introduced in Chauvel and OQuigley (2014) is implemented. The plot of this process over the ranks of the failure times gives an indication of the validity of the proportional hazards assumption. A function to evaluate the \(R^2\) coefficient of OQuigley and Flandre (1994) is provided. This coefficient is a measure of the predictive ability of the proportional hazards model.

References

Chauvel, C (2014). PhD thesis (in French): Processus empiriques pour l'inférence dans le

modèle de survie à risques non proportionnels.

Université Pierre et Marie Curie - Paris VI.

Chauvel, C, OQuigley, J (2014). Tests for comparing estimated survival functions.

Biometrika 101, 3, 535 – 552.

OQuigley, J (2008). Proportional hazards regression. Springer New-York.

OQuigley, J, Flandre, P (1994). Predictive capability of proportional hazards regression. PNAS 91, 2310-2314.

Examples

library(survival)
S=standscore(futime+fustat~resid.ds+rx,data=ovarian)
plotscore(S,printCB = TRUE)
testscore(futime+fustat~resid.ds+rx,data=ovarian)
R2(futime+fustat~resid.ds+rx,data=ovarian)