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)