nlfh (version 0.1.0)

posterior_draws: Extract Posterior Draws from a Fay-Herriot Model

Description

Extract posterior draws stored in fitted Bayesian Fay-Herriot model objects. Draws are returned with one row per retained MCMC draw and one column per requested parameter or area-level quantity.

Usage

posterior_draws(object, variable = c("theta", "beta", "u", "A"), ...)
# S3 method for nlfh_fit
posterior_draws(object, variable = c("theta", "beta", "u", "A"), ...)

Value

A posterior::draws_df object when the posterior package is installed. Otherwise, a data frame with one row per posterior draw.

Arguments

object

A fitted model object.

variable

Character string naming the posterior quantity to extract. Supported values are "theta" for area-level estimates, "beta" for regression or hidden-layer coefficients when available, "u" for area-level random effects when available, and "A" for the random-effect variance parameter.

...

Additional arguments passed to methods.

Examples

data(acs_dat)
acs_small <- as.data.frame(acs_dat[1:500, ])
fit <- fit_fh(
  MedInc ~ SNAPRate + PovRate + White + Black + Hispanic + Asian,
  sampling_variance = MedIncSE^2,
  data = acs_small,
  method = "linear",
  control = list(n_iter = 50, burn_in = 25, progress = FALSE)
)

theta_draws <- posterior_draws(fit, variable = "theta")
beta_draws <- posterior_draws(fit, variable = "beta")