I have the data as follows:
library(survival)
library(survminer)
data(aml)
aml$x <- as.character(aml$x)
aml[10,3] <- 'SuperMaintained'
aml[11,3] <- 'SuperMaintained'
aml[22,3] <- 'SuperMaintained'
aml[23,3] <- 'SuperMaintained'
aml$x <- factor(aml$x, levels = c('Nonmaintained','Maintained','SuperMaintained'))
fit <- survfit(Surv(time, status) ~ x, data = aml)
ggsurvplot(fit, conf.int = FALSE, surv.median.line = c('hv'), data = aml, pval = TRUE, risk.table = FALSE)
[![Survivalplot][1]][1]
res=pairwise_survdiff(Surv(time, status) ~ x, data = aml)
res$p.value
Output look like this
Nonmaintained Maintained
Maintained 0.085155296 NA
SuperMaintained 0.007303126 0.0851553
How do I add this table to the main survival plot?

