Gibbs sampler for linear regression

gibbs_lm(y, x, steps = 1000L)

Arguments

y

a vector of response values

x

a matrix of explanatory variables

steps

number of samples to obtain