Bayesian Regression: lm.bayes()

This is a set of functions that is similar to the lm() function in R for ordinary least squares (OLS). Except this is code for Bayesian regression with reference priors on the coefficients and the variance parameters.   Here is a .pdf file with the full explanation and derivation of the posterior full conditionals for the parameters.