# NOT RUN {
library("spsurv")
data("veteran") ## imports from survival package
fit_mle <- spbp(Surv(time, status) ~ karno + factor(celltype),
data = veteran, model = "po")
summary(fit_mle)
fit_bayes <- spbp(Surv(time, status) ~ karno + factor(celltype),
data = veteran, model = "po", approach = "bayes",
cores = 1, iter = 300, chains = 1,
priors = list(beta = c("normal(0,4)"),
gamma = "lognormal(0,4)"))
summary(fit_bayes)
# }
Run the code above in your browser using DataLab