Gibbs Sampling for Simple Linear Regression
devtools::install_github("https://github.com/jiang275/Gibbs", build_vignettes = TRUE)
library(Gibbs)
data(simulated_data)
draws <- Gibbs(x=simulated_data$x,y=simulated_data$y, sigma2_initial=1, beta0_initial=1,beta1_initial=1,method="deterministic",B=1000)
apply(draws, 2, mean)