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

View source: R/icmle.R

plot.icsurvR Documentation

A plot method for the estimates produced by the estimation methods in Icens.

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)

Bioconductor/Icens documentation built on Oct. 29, 2023, 5:01 p.m.