library(survival)
lung$status <- ifelse(lung$status == 2, 1, 0)
lung$age <- arules::discretize(lung$age, breaks = 3, labels = c("Low","Medium","High"))
library(SurvivalTests)
# \donttest{
out = SANOVA(time~age+status, lung, alpha = 0.05)
Spaircomp(out, adjust.method = "bonferroni")
# to illustrate pairwise comparison, alpha is set to 0.5
out2 = SANOVA(time~age+status, lung, alpha = 0.5)
Spaircomp(out2, adjust.method = "bonferroni")
Spaircomp(out2, adjust.method = "hochberg")
Spaircomp(out2, adjust.method = "hommel")
Spaircomp(out2, adjust.method = "holm")
# }
Run the code above in your browser using DataLab