sample_posterior: Sample the posterior distribution

In JasperHG90/blm: Bayesian linear model

Description

Calling the sample_posterior() function executes a blm sampling plan. The user should set priors, initial values etc. before executing this function, because by executing it the user will 'lock' the sampling plan. That is, no changes may be made except for changing the 'burn' parameter.

Usage

sample_posterior(x)

Arguments

x	blm object

Details

For details on the implementation of the Gibbs sampler, see https://tinyurl.com/y5vk9x35. For details on the implementation of the Metropolis-Hastings sampler, see https://tinyurl.com/y3zp3l3w.

Value

blm object containing sampled posterior

Examples

data("directors")
fit <- blm("Compensation ~ Age", data=directors) %>%
   sample_posterior()

JasperHG90/blm documentation built on Sept. 4, 2019, 11:16 a.m.