# NOT RUN {
if(require("survival") & require("TH.data")) {
## Load data
data(GBSG2, package = "TH.data")
## Weibull model
bmod <- survreg(Surv(time, cens) ~ horTh, data = GBSG2, model = TRUE)
## Coefficient table
grid.newpage()
coeftable.survreg(bmod)
## partitioned model
tr <- pmtree(bmod)
## plot
plot(tr, terminal_panel = node_pmterminal(tr, plotfun = survreg_plot,
confint = TRUE, coeftable = coeftable.survreg))
}
# }
Run the code above in your browser using DataLab