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:

     data(ovarian)
     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)

