RSTr (version 1.1.4)

load_samples: Load MCMC samples

Description

load_samples() gathers samples saved for model RSTr_obj. By default, loads the rate estimate samples lambda, but any model parameters can be loaded. Users can also specify a burn-in period.

Usage

load_samples(RSTr_obj, param = "lambda", burn = 2000)

Value

An array of samples from model RSTr_obj.

Arguments

RSTr_obj

RSTr model object to load in samples from.

param

Which parameter samples to load.

burn

Number of burn-in samples to discard.

Examples

data_min <- lapply(miheart, \(x) x[1:2, 1:3, 1:3])
adj_min <- list(2, 1)
on.exit(unlink(file.path(tempdir(), "test"), recursive = TRUE), add = TRUE)
mod_mst <- mstcar("test", data_min, adj_min, tempdir(), show_plots = FALSE, verbose = FALSE)
samples <- load_samples(mod_mst) * 1e5