The BGLIMM Procedure

Conjugate Sampling

Conjugate prior distributions are a type of prior distribution in which the prior and posterior distributions are in the same family of distributions. For example, if you model an independently and identically distributed random variable y Subscript i by using a normal likelihood with known variance sigma squared,

y Subscript i Baseline tilde normal left-parenthesis bold x prime Subscript i Baseline bold-italic beta comma sigma squared right-parenthesis

then a normal prior on bold-italic beta,

bold-italic beta tilde normal left-parenthesis bold-italic beta 0 comma bold upper Sigma 0 right-parenthesis

is a conjugate prior, because the posterior distribution of bold-italic beta is also a normal distribution, where the covariance is bold upper Sigma 0 Superscript negative 1 Baseline plus StartFraction 1 Over sigma squared EndFraction sigma-summation Underscript i equals 1 Overscript upper N Endscripts bold x Subscript i Baseline bold x prime Subscript i and the mean is left-parenthesis bold upper Sigma 0 Superscript negative 1 Baseline plus StartFraction 1 Over sigma squared EndFraction sigma-summation Underscript i equals 1 Overscript upper N Endscripts bold x Subscript i Baseline bold x prime Subscript i right-parenthesis Superscript negative 1 Baseline left-parenthesis bold upper Sigma 0 Superscript negative 1 Baseline bold-italic beta 0 plus StartFraction 1 Over sigma squared EndFraction sigma-summation Underscript i equals 1 Overscript upper N Endscripts bold x Subscript i Baseline y Subscript i Baseline right-parenthesis

PROC BGLIMM uses conjugate samplers in the normal and multivariate normal cases, as shown in Table 14.

Table 14: Conjugate Sampling in PROC BGLIMM

Family Parameter Prior
Normal bold-italic beta Multivariate normal
Normal Variance sigma squared Inverse gamma
Multivariate normal bold-italic beta Multivariate normal
Multivariate normal Covariance bold upper Sigma Inverse Wishart


The fixed-effects parameters bold-italic beta and the covariances bold upper G and bold upper R are sampled when applicable.

Last updated: December 09, 2022