library(vmsae)
example_model <- readRDS(system.file("example", "example_model.Rds", package = "vmsae"))
confint(example_model) # Get credible intervals for predicted values
confint(example_model, field = "beta_samples") # For fixed effects
Run the code above in your browser using DataLab