Kaplan_Meier2

library(survival)
library(survminer)

## Loading required package: ggplot2

## Loading required package: ggpubr

## 
## Attaching package: 'survminer'

## The following object is masked from 'package:survival':
## 
##     myeloma

d1=data.frame(time=c(3,4,6,11,17,21,24,25,26,30),
              event=c(1,1,0,1,0,1,0,0,1,0))
kmc=with(d1,Surv(time,event));kmc

##  [1]  3   4   6+ 11  17+ 21  24+ 25+ 26  30+

plot(kmc)

plot(kmc,fun="cumhaz")

kmc2=surv_fit(Surv(time,event)~1,data = d1);kmc2

## Call: survfit(formula = Surv(time, event) ~ 1, data = structure(list(
##     time = c(3, 4, 6, 11, 17, 21, 24, 25, 26, 30), event = c(1, 
##     1, 0, 1, 0, 1, 0, 0, 1, 0)), class = "data.frame", row.names = c(NA, 
## -10L)))
## 
##       n events median 0.95LCL 0.95UCL
## [1,] 10      5     26      11      NA

summary(kmc2)

## Call: survfit(formula = Surv(time, event) ~ 1, data = structure(list(
##     time = c(3, 4, 6, 11, 17, 21, 24, 25, 26, 30), event = c(1, 
##     1, 0, 1, 0, 1, 0, 0, 1, 0)), class = "data.frame", row.names = c(NA, 
## -10L)))
## 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     3     10       1    0.900  0.0949       0.7320            1
##     4      9       1    0.800  0.1265       0.5868            1
##    11      7       1    0.686  0.1515       0.4447            1
##    21      5       1    0.549  0.1724       0.2963            1
##    26      2       1    0.274  0.2122       0.0602            1