plot.icfit {interval} | R Documentation |
Creates survival curves for an icfit
object. If there is more than one strata,
all strata will be fit. Note that for interval censored data usually the changes in the
NPMLE of the survival function do not occur at unique points but occur within some interval
where any of an infinite number of curves will maximize the likelihood. We show those intervals
were the NPMLE is indeterminate as a gray rectangle.
## S3 method for class 'icfit': plot(x,XLAB="time",YLAB="Survival",COL=gray((8:1)*.1),LTY=1:9,LEGEND=NULL, XLEG=NULL,YLEG=.1,...)
x |
an icfit object, see icfit |
XLAB |
x label |
YLAB |
y label |
COL |
a vector representing color of rectangles of indeterminate NPMLE, COL[i] used for ith strata |
LTY |
a vector for lty values for lines, LTY[i] used for ith strata |
LEGEND |
logical value, include legend or not, if NULL set to TRUE only if number of strata>1 |
XLEG |
x location for legend, if NULL then gives maximum of 0 and minimum time from intmap |
YLEG |
y location for legend |
... |
other arguments passed to the plot function |
Returns a list of arguments for the legend. Values are x,y, legend, fill, lty. See legend
help.
An object of class 'icsurv' from the Icens
package can use this plot function by redefining its class to 'icfit' and 'plot.icfit' will work on it.
data(bcos) fit1<-icfit(Surv(left,right,type="interval2")~treatment,data=bcos) summary(fit1) plot(fit1)