tests/regtest-survival.R

In multcomp: Simultaneous Inference in General Parametric Models

library("multcomp")
set.seed(29)

if (requireNamespace("survival")) {

bladder1 <- bladder[bladder$enum < 5, ] 
bladder1$number <- as.factor(bladder1$number)

cm <- coxph(Surv(stop, event) ~ rx + size + number + strata(enum) + 
            cluster(id), data = bladder1)
print(summary(glht(cm, mcp(number = "Dunnett"))))

wm <- survreg(Surv(stop, event) ~ rx + size + number + strata(enum) + 
              cluster(id), data = bladder1)
print(summary(glht(wm, mcp(number = "Dunnett"))))

}