pg_gibbs: Polya-Gamma Gibbs sampler

View source: R/logistic_regression.R

pg_gibbsR Documentation

Polya-Gamma Gibbs sampler

Description

This implements the sampler proposed in Nicholas G Polson, James G Scott, and Jesse Windle. Bayesian inference for logistic models using Polya–Gamma latent variables. Journal of the American statistical Association, 108(504):1339–1349, 2013. The arguments are:

  • niterations: the number of desired MCMC iterations,

  • logistic_setting: a list of precomputed quantities obtained via 'logistic_precomputation'.

Usage

pg_gibbs(niterations, logistic_setting)

Value

a matrix where each row corresponds to an iteration of the sampler, and contains the regression coefficient at that iteration.


pierrejacob/debiasedmcmc documentation built on Aug. 22, 2022, 12:41 a.m.