Multiple Comparisons of Survival Curves
Calculate pairwise comparisons between group levels with corrections for multiple testing.
pairwise_survdiff(formula, data, p.adjust.method = "BH", na.action, rho = 0)
formula |
a formula expression as for other survival models, of the form Surv(time, status) ~ predictors. |
data |
a data frame in which to interpret the variables occurring in the formula. |
p.adjust.method |
method for adjusting p values (see
|
na.action |
a missing-data filter function. Default is options()$na.action. |
rho |
a scalar parameter that controls the type of test. Allowed values include 0 (for Log-Rank test) and 1 (for peto & peto test). |
Returns an object of class "pairwise.htest", which is a list containing the p values.
Alboukadel Kassambara, alboukadel.kassambara@gmail.com
survival::survdiff
library(survival)
library(survminer)
data(myeloma)
# Pairwise survdiff
res <- pairwise_survdiff(Surv(time, event) ~ molecular_group,
data = myeloma)
res
# Symbolic number coding
symnum(res$p.value, cutpoints = c(0, 0.0001, 0.001, 0.01, 0.05, 0.1, 1),
symbols = c("****", "***", "**", "*", "+", " "),
abbr.colnames = FALSE, na = "")Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.