| plot.icsurv {Icens} | R Documentation |
Procduces nice plots of the estimated NPMLE.
plot.icsurv(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 |
~~Describe surv here~~ |
bounds |
~~Describe bounds here~~ |
shade |
~~Describe shade here~~ |
density |
~~Describe density here~~ |
angle |
~~Describe angle here~~ |
lty |
~~Describe lty here~~ |
new |
~~Describe new here~~ |
xlab |
~~Describe xlab here~~ |
ylab |
~~Describe ylab here~~ |
main |
~~Describe main here~~ |
ltybnds |
~~Describe ltybnds here~~ |
... |
~~Describe ... here~~ |
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)