plot.icfit {interval}R Documentation

Plot icfit object

Description

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.

Usage

## 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,...)

Arguments

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

Value

Returns a list of arguments for the legend. Values are x,y, legend, fill, lty. See legend help.

Note

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.

See Also

icfit

Examples

data(bcos)
fit1<-icfit(Surv(left,right,type="interval2")~treatment,data=bcos)
summary(fit1)
plot(fit1)
 

[Package interval version 0.7-5.5 Index]