get_subset {remiod}

Extract specific parameters from MCMC samples

Description

Extract specific parameters from MCMC samples

Usage

get_subset(object, subset, warn = TRUE, mess = TRUE)

Arguments

object

an object of class MCMC.

subset

subset of parameters (columns of the mcmc object) to be used. See https://nerler.github.io/JointAI/articles/SelectingParameters.html for key-words of subseting parameters. Besides, selected_parms and selected_vars are new key-words for arbitrarily selecting parameters.

warn

logical, should warning messages be displayed?

mess

logical, should messages be displayed?

Examples

data(schizow)

test = remiod(formula = y6 ~ tx + y0 + y1 + y3, data = schizow,
              trtvar = 'tx', algorithm = 'jags', method="MAR",
              ord_cov_dummy = FALSE, n.adapt = 10, n.chains = 1,
              n.iter = 0, thin = 1, warn = FALSE, seed = 1234)

pms = c("beta[2]","alpha[2]","alpha[6]","alpha[9]")
mcsub = get_subset(object = test$mc.mar, subset=c(selected_parms = list(pms)))