lines.survfit            package:survival            R Documentation

_A_d_d _L_i_n_e_s _o_r _P_o_i_n_t_s _t_o _a _S_u_r_v_i_v_a_l _P_l_o_t

_D_e_s_c_r_i_p_t_i_o_n:

     Often used to add the expected survival curve(s) to a Kaplan-Meier
     plot  generated with 'plot.survfit'.

_U_s_a_g_e:

     ## S3 method for class 'survfit':
     lines(x, type="s", mark=3, col=1, lty=1,
     lwd=1, mark.time=TRUE, 
             xscale=1,  firstx=0, firsty=1, xmax, fun, conf.int=FALSE,  ...)
     ## S3 method for class 'survfit':
     points(x, ...)

_A_r_g_u_m_e_n_t_s:

       x: a survival object, generated from the 'survfit' or 'survexp'
          functions.  

    type: the line type, as described in 'lines'.  The default is a
          step function  for 'survfit' objects, and a connected line
          for 'survexp' objects.  

mark, col, lty, lwd: vectors giving the mark symbol, color, line type
          and line width for the  added curves.  

     ...: other graphical parameters

mark.time: controls the labeling of the curves.    If 'FALSE', no
          labeling is done.    If 'TRUE', then curves are marked at
          each censoring time.    If 'mark.time' is a numeric vector,
          then curves are marked at   the specified time points.  

  xscale: a number used to divide the x values.  If time was originally
          in days, a  value of 365.25 would give a plotted scale in
          years.  

firstx, firsty: the starting point for the survival curves.  If either
          of these is set to  'NA' or < blank > the plot will start at
          the first time point of the curve.  

    xmax: the maximum horizontal plot coordinate.      This shortens
          the curve before plotting it, so unlike using the  'xlim'
          graphical parameter, warning messages about out of bounds
          points are  not generated.  

     fun: an arbitrary function defining a transformation of the
          survival curve.  For example 'fun=log' is an alternative way
          to draw a log-survival curve  (but with the axis labeled with
          log(S) values).  Four often used transformations can be
          specified with a character  argument instead: "log" is the
          same as using the 'log=T' option,  "event" plots cumulative
          events (f(y) = 1-y),  "cumhaz" plots the cumulative hazard
          function (f(y) = -log(y))  and "cloglog" creates a
          complimentary log-log survival plot   (f(y) = log(-log(y)
          along with log scale for the x-axis).  

conf.int: if 'TRUE', confidence bands for the curves are also plotted. 
          If set to '"only"', then only the CI bands are plotted, and
          the curve  itself is left off.    This can be useful for fine
          control over the colors or line types of a plot.  

_V_a_l_u_e:

     a list with components 'x' and 'y', containing the coordinates of
     the  last point on each of the curves (but not of the confidence
     limits).  This may be useful for labeling.

_S_i_d_e _E_f_f_e_c_t_s:

     one or more curves are added to the current plot.

_S_e_e _A_l_s_o:

     'lines',  'par',  'plot.survfit',  'survfit',  'survexp'.

_E_x_a_m_p_l_e_s:

     fit <- survfit(Surv(time, status) ~ sex, pbc,subset=1:312)
     plot(fit, mark.time=FALSE, xscale=365.25,
             xlab='Years', ylab='Survival')
     lines(fit[1], lwd=2, xscale=365.24)    #darken the first curve and add marks

     # Add expected survival curves for the two groups,
     #   based on the US census data
     # The data set does not have entry date, use the midpoint of the study
     efit <- survexp(~ ratetable(sex=sex,age=age*365.35,year=as.Date('1979/1/1')) +
                 sex, data=pbc,  times=(0:24)*182)
     temp <- lines(efit, lty=2, xscale=365.24, lwd=2:1)
     text(temp, c("Male", "Female"), adj= -.1) #labels just past the ends
     title(main="Primary Biliary Cirrhosis, Observed and Expected")