A plot method for the estimates produced by the estimation methods in Icens.
Produces nice plots of the estimated NPMLE.
## 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, ...)
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. |
No value is returned. A plot of the NPMLE is made on the active graphics device.
Alain Vandal and Robert Gentleman.
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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.