PHeval-package | R Documentation |
This package provides tools for the evaluation of the goodness of fit of the proportional hazards model.
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.
Cecile Chauvel <cecile.chauvel@univ-lyon1.fr>
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. |
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.