# NOT RUN {
# Hierarchical model for responses in a disease with multiple subtypes
dat <- thallhierarchicalbinary_parameters_demo()
samp <- rstan::sampling(stanmodels$ThallHierarchicalBinary, data = dat)
rstan::plot(samp, pars = 'p') # Plot the modelled response rates in the subtypes
# }
Run the code above in your browser using DataLab