Nothing
#############################
# rmst1 (one-arm) -- hidden
#############################
rmst1=function(time, status, tau, alpha=0.05){
#-- time
#-- statuts
#-- tau -- truncation time
#-- alpha -- gives (1-alpha) confidence interval
ft= survfit(Surv(time, status)~1)
idx=ft$time<=tau
wk.time=sort(c(ft$time[idx],tau))
wk.surv=ft$surv[idx]
wk.n.risk =ft$n.risk[idx]
wk.n.event=ft$n.event[idx]
time.diff <- diff(c(0, wk.time))
areas <- time.diff * c(1, wk.surv)
rmst = sum(areas)
rmst
wk.var <- ifelse((wk.n.risk-wk.n.event)==0, 0,
wk.n.event /(wk.n.risk *(wk.n.risk - wk.n.event)))
wk.var =c(wk.var,0)
rmst.var = sum( cumsum(rev(areas[-1]))^2 * rev(wk.var)[-1])
rmst.se = sqrt(rmst.var)
#--- check ---
# print(ft, rmean=tau)
#--- output ---
out=matrix(0,2,4)
out[1,]=c(rmst, rmst.se, rmst-qnorm(1-alpha/2)*rmst.se, rmst+qnorm(1-alpha/2)*rmst.se)
out[2,]=c(tau-out[1,1], rmst.se, tau-out[1,4], tau-out[1,3])
rownames(out)=c("RMST","RMTL")
colnames(out)=c("Est.", "se", paste("lower .",round((1-alpha)*100, digits=0), sep=""), paste("upper .",round((1-alpha)*100, digits=0), sep=""))
Z=list()
Z$result=out
Z$rmst = out[1,]
Z$rmtl = out[2,]
Z$tau=tau
Z$rmst.var = rmst.var
Z$fit=ft
class(Z)="rmst1"
return(Z)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.