get_subset: Extract specific parameters from MCMC samples

In xsswang/remiod: Reference-Based Multiple Imputation for Ordinal/Binary Response

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)))

xsswang/remiod documentation built on Feb. 20, 2023, 9:29 a.m.