library(survival)
library(survminer)
## Warning: package 'survminer' was built under R version 4.4.3
## Loading required package: ggplot2
## Loading required package: ggpubr
## Warning: package 'ggpubr' was built under R version 4.4.3
##
## Attaching package: 'survminer'
## The following object is masked from 'package:survival':
##
## myeloma
d1=data.frame(time=c(2,4,6,10),event=c(1,0,1,1))
kmc=with(d1,Surv(time,event));kmc
## [1] 2 4+ 6 10
kmc2=surv_fit(Surv(time,event)~1,data=d1);kmc2
## Call: survfit(formula = Surv(time, event) ~ 1, data = structure(list(
## time = c(2, 4, 6, 10), event = c(1, 0, 1, 1)), class = "data.frame", row.names = c(NA,
## -4L)))
##
## n events median 0.95LCL 0.95UCL
## [1,] 4 3 6 2 NA
summary(kmc2)
## Call: survfit(formula = Surv(time, event) ~ 1, data = structure(list(
## time = c(2, 4, 6, 10), event = c(1, 0, 1, 1)), class = "data.frame", row.names = c(NA,
## -4L)))
##
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 2 4 1 0.750 0.217 0.4259 1
## 6 2 1 0.375 0.286 0.0839 1
## 10 1 1 0.000 NaN NA NA
plot(kmc,col="red")

plot(kmc,fun="cumhaz")
