pbc                package:risksetROC                R Documentation

_I_n_c_i_d_e_n_t/_D_y_n_a_m_i_c (_I/_D) _R_O_C _c_u_r_v_e, _A_U_C _a_n_d _i_n_t_e_g_r_a_t_e_d _A_U_C (_i_A_U_C)
_e_s_t_i_m_a_t_i_o_n _o_f _c_e_n_s_o_r_e_d _s_u_r_v_i_v_a_l _d_a_t_a

_D_e_s_c_r_i_p_t_i_o_n:

     This is Mayo PBC data as obtained from the website:
     http://lib.stat.cmu.edu/datasets/pbc

_F_o_r_m_a_t:

     A data frame with 418 observations and 20 variables:

_i_d patient id

_f_u_d_a_y_s follow-up days, number of days between registration and the
     earlier of death, transplantation, or study analysis time in July,
     1986 

_s_t_a_t_u_s survival status

_d_r_u_g 1 = D-penicillamine, 2 = placebo

_a_g_e age in days

_s_e_x 0 = male, 1 = female

_a_s_c_i_t_e_s presence of asictes: 0=no 1=yes

_h_e_p_a_t_o_m presence of hepatomegaly: 0=no 1=yes

_s_p_i_d_e_r_s presence of spiders: 0=no 1=yes

_e_d_e_m_a presence 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

_b_i_l_i serum bilirubin in mg/dl

_c_h_o_l serum cholesterol in mg/dl

_a_l_b_u_m_i_n albumin in gm/dl

_c_o_p_p_e_r urine copper in ug/day

_a_l_k_p_h_o_s alkaline phosphatase in U/liter

_s_g_o_t SGOT in U/ml

_t_r_i_g triglicerides in mg/dl

_p_l_a_t_e_l_e_t platelets per cubic ml / 1000

_p_r_o_t_i_m_e prothrombin time in seconds

_s_t_a_g_e histologic stage of disease

_A_u_t_h_o_r(_s):

     Patrick J. Heagerty

_R_e_f_e_r_e_n_c_e_s:

     Heagerty, P.J., Zheng Y. (2005) Survival Model Predictive Accuracy
     and ROC curves _Biometrics_, *61*, 92 - 105

_E_x_a_m_p_l_e_s:

     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)

