riskset              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 function creates risk set at each unique failure time from a
     survival data set.

_U_s_a_g_e:

     riskset(dat, entry=FALSE)

_A_r_g_u_m_e_n_t_s:

     dat: survival dataset with at least three variables:
          survival.times, survival.status and marker, in that order.
          The survival data set may have additional variables. In case
          of interval censored data, the first four columns are: entry
          time, exit time, status at exit and marker 

   entry: default is FALSE indicating right censored data. TRUE if left
          truncated data

_D_e_t_a_i_l_s:

     This function creates risk set at each unique failure time from a
     survival data set and is needed for llCoxReg(). The function can
     handle both right censored and interval censored data.

_V_a_l_u_e:

     Returns a new data set with columns as follows: start, finish,
     newStatus and other variables from the original dataset except
     survival time and status. The first two columns correspond to the
     start and end of time intervals considered and the newStatus
     corresponds to the survival status of the patient corresponding to
     this interval, i.e. the status is 1 if the patient had event
     during this interval (start, finish] and 0 otherwise. Note that
     the survival time need to be in ascending order.

_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)

