Predict survival from a stan_surv object
1 2 |
1 2 3 4 5 6 7 8 9 10 | pbc2 <- survival::pbc
pbc2 <- pbc2[!is.na(pbc2$trt),]
pbc2$status <- as.integer(pbc2$status > 0)
m1 <- stan_surv(survival::Surv(time, status) ~ trt, data = pbc2)
df <- flexsurv::bc
m2 <- stan_surv(survival::Surv(rectime, censrec) ~ group,
data = df, cores = 1, chains = 1, iter = 2000,
basehaz = "fpm", iknots = c(6.594869, 7.285963 ),
degree = 2, prior_aux = normal(0, 2, autoscale = F))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.