survobrien             package:survival             R Documentation

_O'_B_r_i_e_n'_s _T_e_s_t _f_o_r _A_s_s_o_c_i_a_t_i_o_n _o_f _a _S_i_n_g_l_e _V_a_r_i_a_b_l_e _w_i_t_h _S_u_r_v_i_v_a_l

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

     Peter O'Brien's test for association of a single variable with
     survival  This test is proposed in Biometrics, June 1978.

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

     survobrien(formula, data)

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

 formula: a valid formula for a cox model, without time dependent
          covariates.  

    data: a data frame.  

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

     a new data frame.  The original time and status variables are
     removed,  and have been replaced with 'start', 'stop', and
     'event'.  If a  predictor variable is a factor or is protected
     with 'I()', it is  retained as is.  Other predictor variables have
     been replaced with  time-dependent logit scores. 

     Because of the time dependent variables, the new data frame will
     have many  more rows that the original data, approximately #rows *
     #deaths /2.

_M_e_t_h_o_d:

     A time-dependent cox model can now be fit to the new data.  The
     univariate statistic, as originally proposed, is equivalent to 
     single variable score tests from the time-dependent model.  This
     equivalence is the rationale for using the time dependent model as
     a  multivariate extension of the original paper. 

     In O'Brien's method, the x variables are re-ranked at each death
     time.  A  simpler method, proposed by Prentice, ranks the data
     only once at the  start. The results are usually similar.

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

     O'Brien, Peter, "A Nonparametric Test for Association with
     Censored Data",  _Biometrics_ 34: 243-250, 1978.

_S_e_e _A_l_s_o:

     'survdiff'

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

     xx <- survobrien(Surv(futime, fustat) ~ age + factor(rx) + I(ecog.ps), 
                                    data=ovarian) 
     coxph(Surv(start, stop, event) ~ age, data=xx) 
     coxph(Surv(start, stop, event) ~ age + rx + ecog.ps, data=xx) 

