gibbs_lm_R: Gibbs sampler for linear regression

In iNZightVIT/iNZightBayes: Bayesian Inference with 'iNZight'

Description

Gibbs sampler for linear regression

Usage

gibbs_lm_R(y, X, steps)

Arguments

y	a vector of the response variable
X	a matrix whose columns are the explanatory variables
steps	number of iterations to run Gibbs sampler for

Details

Assumes Normal likelihood and joint prior proportional to 1 / sigma2

Value

a list of posterior samples for regression coefficients

Examples

## Not run: 
set.seed(1234)
n = 100
sigma = 40
beta = c(150, 5, 10)
x = seq(-4, 10, length = n)
y = beta[1] + beta[2] * x + beta[3] * x ^ 2 + rnorm(n, 0, sigma)
X = cbind(x, x ^ 2)
mcmc = gibbs_lm_R(y, X, 10000)

## End(Not run)

iNZightVIT/iNZightBayes documentation built on Dec. 20, 2021, 5:56 p.m.