Learn R
Icens (version 1.44.0)

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

Description

Produces nice plots of the estimated NPMLE.

Usage

"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.

See Also

VEM, ISDM, EMICM, PGM

Examples

Run this code
    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)

Run the code above in your browser using DataCamp Workspace