plot.icsurv: A plot method for the estimates produced by the estimation...

In Icens: NPMLE for Censored and Truncated Data

Description

Produces nice plots of the estimated NPMLE.

Usage

## S3 method for class 'icsurv'
plot(x, type="eq", surv=FALSE, bounds=FALSE, shade=3, density=30,
angle=45, lty=1, new=TRUE, xlab="Time", ylab="Probability", main="GMLE",
ltybnds=2, ...)

Arguments

x	The estimate of the NPMLE.
type	Three options, "eq" for equivalence call, "gw" for the Groeneboom-Wellner estimate, and "lc" for the left-continuous estimate.
surv	Logical indicating whether or not to plot the survival curve.
bounds	Logical indicating whether or not to include bounds around the estimate.
shade	An integer in 1, 2, or 3 denoting what component of the plot to shade.
density	The density of shading lines, in lines per inch.
angle	The slope of shading lines, given as an angle in degrees (counter-clockwise).
lty	The line type for the estimates.
new	Logical indicating whether or not to create a new plot.
xlab	The x-axis label.
ylab	The y-axis label.
main	The main title for the plot.
ltybnds	The line type for the bounds on the estimates.
...	Additional arguments passed to the plot function.

Value

No value is returned. A plot of the NPMLE is made on the active graphics device.

Author(s)

Alain Vandal and Robert Gentleman.

See Also

VEM, ISDM, EMICM, PGM

Examples

    data(cosmesis)
    csub1 <- subset(cosmesis, subset=Trt==0, select=c(L,R))
    e1 <- VEM(csub1)
    par(mfrow=c(2,2))
    plot(e1)
    data(pruitt)
    e2 <- EM(csub1)
    plot(e2)
    e3 <- PGM(csub1)
    plot(e3)
    e4 <- EMICM(csub1)
    plot(e4)

Example output

Loading required package: survival

Icens documentation built on Nov. 8, 2020, 8:12 p.m.