pbc {risksetROC} | R Documentation |
This is Mayo PBC data as obtained from the website: http://lib.stat.cmu.edu/datasets/pbc
A data frame with 418 observations and 20 variables: \itemidpatient id \itemfudaysfollow-up days, number of days between registration and the earlier of death, transplantation, or study analysis time in July, 1986 \itemstatussurvival status \itemdrug1 = D-penicillamine, 2 = placebo \itemageage in days \itemsex0 = male, 1 = female \itemascitespresence of asictes: 0=no 1=yes \itemhepatompresence of hepatomegaly: 0=no 1=yes \itemspiderspresence of spiders: 0=no 1=yes \itemedemapresence of edema: 0=no edema and no diuretic therapy for edema; .5 = edema present without diuretics, or edema resolved by diuretics; 1 = edema despite diuretic therapy \itembiliserum bilirubin in mg/dl \itemcholserum cholesterol in mg/dl \itemalbuminalbumin in gm/dl \itemcopperurine copper in ug/day \itemalkphosalkaline phosphatase in U/liter \itemsgotSGOT in U/ml \itemtrigtriglicerides in mg/dl \itemplateletplatelets per cubic ml / 1000 \itemprotimeprothrombin time in seconds \itemstagehistologic stage of disease
Patrick J. Heagerty
Heagerty, P.J., Zheng Y. (2005) Survival Model Predictive Accuracy and ROC curves Biometrics, 61, 92 – 105
library(MASS) data(VA) ## need to order the data in ascending order of survival time new.VA=VA[order(VA$stime),] risket.VA=riskset(new.VA)