interval (version 0.7-5.5)

plot.icfit: 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.

See Also

icfit

Examples

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