## Load BaSTA output:
data("bastaCMRout", package = "BaSTA")
## Plot traces for survival parameters:
plot(bastaCMRout)
## Plot posterior densities of survival parameters:
plot(bastaCMRout, densities = TRUE)
## Plot traces for proportional hazards parameter:
plot(bastaCMRout, trace.name = "gamma")
## Plot survival and mortality curves:
plot(bastaCMRout, plot.type = "demorates")
Run the code above in your browser using DataLab