Kaplan-Meier多组比较

假设我有以下数据帧

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))

库（生存）
图书馆（multcomp）
数据（cml）
cml$groupmultcomp的generalsiminf.pdf中有一个例子。这里简化

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")))

库（multcomp）
图书馆（生存）
如果（！file.exists（“AML_Bullinger.rda”））
加载（url）http://www.stat.uni-muenchen.de/~hothorn/data/AML_Bullinger.rda“，open=“r”））
风险multcomp的generalsiminf.pdf中有一个例子。这里简化

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")))

库（multcomp）
图书馆（生存）
如果（！file.exists（“AML_Bullinger.rda”））
加载（url）http://www.stat.uni-muenchen.de/~hothorn/data/AML_Bullinger.rda“，open=“r”））
非常感谢。我有一个问题。survreg符合参数生存回归模型。我拥有的是简单的知识管理分析估算。我可以使用survreg和multicomp对KM估计值进行多重比较吗？很好。我不知道KM多个测试包，但如果这是严重的，我可能会去Bayes（如果您的老板/编辑允许）试试。这个问题当然是正确的。让我们看看会发生什么。非常感谢你。我有一个问题。survreg符合参数生存回归模型。我拥有的是简单的知识管理分析估算。我可以使用survreg和multicomp对KM估计值进行多重比较吗？很好。我不知道KM多个测试包，但如果这是严重的，我可能会去Bayes（如果您的老板/编辑允许）试试。这个问题当然是正确的。让我们看看会发生什么。