Kaplan-Meier meerdere groepsvergelijkingen

Laten we zeggen dat ik het volgende dataframe heb

library(survival)
library(multcomp)
data(cml)
cml$group<-sample(1:5, 507, replace=T)
plot(survfit(Surv(time=cml$time, cml$status)~cml$group))
(survdiff(Surv(time=cml$time, cml$status)~cml$group))

Hoe kan ik een vergelijkingstest uitvoeren om bijvoorbeeld groep0 te vergelijken met alle andere groepen? Of elke groep met elkaar?

Is er een leuke manier om deze meerdere vergelijkingen te plotten (zoals bijvoorbeeld in plot.TukeyHSD() in aov() ?

1

1 antwoord

Er is een voorbeeld in de generalalsiminf.pdf van multcomp. Hier vereenvoudigd

library(multcomp)
library(survival)
if (!file.exists("AML_Bullinger.rda"))
  load(url("http://www.stat.uni-muenchen.de/~hothorn/data/AML_Bullinger.rda", open = "r"))
risk <- rep(0, nrow(clinical))
rlev <- levels(clinical[, "Cytogenetic.group"])
risk[clinical[, "Cytogenetic.group"] %in% rlev[c(7,8,4)]] <- "low"
risk[clinical[, "Cytogenetic.group"] %in% rlev[c(5, 9)]] <- "intermediate"
risk[clinical[, "Cytogenetic.group"] %in% rlev[-c(4,5, 7,8,9)]] <- "high"
risk <- as.factor(risk)
names(clinical)[6] <- "FLT3"
save(clinical,file="AML_Bullinger.rda")
smod <- survreg(Surv(time, event) ~ Sex + Age + WBC+risk,
                 data = clinical)
summary(glht(smod, linfct = mcp(risk = "Tukey")))
1
toegevoegd
Dank u zeer. Ik heb een vraag. survreg past op een parametrisch overlevingsregressiemodel. Wat ik heb, is een eenvoudige analyse van de KM-analyse. Kan ik survreg en multicomp gebruiken om meerdere KM-schattingen te vergelijken?
toegevoegd de auteur ECII, de bron
Moet ik dit vragen in stats.stackexchange.com ?
toegevoegd de auteur ECII, de bron
Oke ik heb het gedaan. Laten we eens kijken wat er aan de hand is. stats.stackexchange.com/questions/31077/&hellip ;
toegevoegd de auteur ECII, de bron
Goed punt. Ik weet niets over KM-meervoudige testpakketten, maar als dit serieus is, zou ik waarschijnlijk naar Bayes gaan (als je baas/redacteur het toestaat) cran.r-project.org/web/views/Survival.html
toegevoegd de auteur Dieter Menne, de bron
Proberen. De vraag is zeker geldig.
toegevoegd de auteur Dieter Menne, de bron