Update Conjugate_priors_Bayes.R

Conjugate_priors_Bayes.R

@@ -27,13 +27,14 @@ beta.prior.mean <- c(0,1)
 beta.prior.precision <- 1*diag(p)
 
 ## Emulate the original Cauchy prior with some specific parameters of InvGamma
-prior.shape <- 0.5
-prior.rate <- 0.5
+prior.shape <- 2
+prior.rate <- 2
 
 bmod <- bayesLMConjugate(y ~ x, data=model_data, n.samples,beta.prior.mean, beta.prior.precision, prior.shape, prior.rate)
 beta <- as.vector(bmod$p.beta.tauSq.samples[,1])
 alpha <- as.vector(bmod$p.beta.tauSq.samples[,2])
-sigma <- as.vector(bmod$p.beta.tauSq.samples[,3])
+tauSq <- as.vector(bmod$p.beta.tauSq.samples[,3])
+sigma <- sqrt(tauSq)
 
 hist(beta)
 hist(alpha)