ex=read.table("http://www.stat.unm.edu/~ghuerta/stat574/collett1.2.txt",col.names=c("time","event"))
library(survival)
fit=survfit( Surv(time,event)~1, data=ex)
summary(fit)
par(mfrow=c(2,2))
plot(fit)
plot(fit,fun="cumhaz")
plot(fit,fun="cloglog")
fit=survfit( Surv(time,event)~1, data=ex,type="fleming-harrington")
summary(fit)
par(mfrow=c(1,1))
plot(fit)
# no censoring
fit=survfit( Surv(ex$time, rep(1,13))~1)
par(mfrow=c(1,1))
plot(fit)
summary(fit)
lam=1/mean(ex$time)
d=exp(-ex$time*lam)
lines(ex$time,d)
tt=seq(0,250,by=0.1)
s=exp(-tt*lam)
lines(tt,s)