diff --git a/test/numerical/bayesian/LinearRegression.jl b/test/numerical/bayesian/LinearRegression.jl
index f49defe..8758c25 100644
--- a/test/numerical/bayesian/LinearRegression.jl
+++ b/test/numerical/bayesian/LinearRegression.jl
@@ -4,7 +4,7 @@ tests = [
     (Prior_Ridge(), 20.080877893580514),
     (Prior_Laplace(), 20.070783434589128),
     (Prior_Cauchy(), 20.019759144845644),
-    (Prior_TDist(), 20.042614331921428),
+    (Prior_TDist(), 20.08147561106022),
     (Prior_HorseShoe(), 20.042984550677183),
 ]
 
diff --git a/test/numerical/bayesian/LogisticRegression.jl b/test/numerical/bayesian/LogisticRegression.jl
index ac9ad61..2708cfa 100644
--- a/test/numerical/bayesian/LogisticRegression.jl
+++ b/test/numerical/bayesian/LogisticRegression.jl
@@ -32,9 +32,9 @@ tests = [
         Prior_TDist(),
         (
             (Logit(), 0.588835403024102),
-            (Probit(), 0.7635642627091132),
-            (Cloglog(), 0.7609943137312546),
-            (Cauchit(), 0.772095066757767)
+            (Probit(), 0.7642595382152266),
+            (Cloglog(), 0.7571826865012884),
+            (Cauchit(), 0.7713699945746971)
         )
     ),
     (
diff --git a/test/numerical/bayesian/NegBinomialRegression.jl b/test/numerical/bayesian/NegBinomialRegression.jl
index 58ed951..0d75a20 100644
--- a/test/numerical/bayesian/NegBinomialRegression.jl
+++ b/test/numerical/bayesian/NegBinomialRegression.jl
@@ -1,11 +1,11 @@
 sanction = dataset("Zelig", "sanction")
 
 tests = [
-    (Prior_Ridge(), 6.8753100988051274),
-    (Prior_Laplace(), 6.908332048475347),
-    (Prior_Cauchy(), 6.9829255933233645),
-    (Prior_TDist(), 6.915515248823249),
-    (Prior_HorseShoe(), 6.703023191644206),
+    (Prior_Ridge(), 6.999865486088317),
+    (Prior_Laplace(), 6.886529206600885),
+    (Prior_Cauchy(), 6.900001819752649),
+    (Prior_TDist(), 6.876415480722939),
+    (Prior_HorseShoe(), 6.902138507950901),
 ]
 
 for (prior, test_mean) in tests
diff --git a/test/numerical/bayesian/PoissonRegression.jl b/test/numerical/bayesian/PoissonRegression.jl
index aa7b6f0..d80b490 100644
--- a/test/numerical/bayesian/PoissonRegression.jl
+++ b/test/numerical/bayesian/PoissonRegression.jl
@@ -1,11 +1,11 @@
 sanction = dataset("Zelig", "sanction")
 
 tests = [
-    (Prior_Ridge(), 7.177578002644547),
-    (Prior_Laplace(), 7.1454141602741785),
-    (Prior_Cauchy(), 7.148699646242317),
-    (Prior_TDist(), 7.165968828611132),
-    (Prior_HorseShoe(), 7.144190707091213),
+    (Prior_Ridge(), 7.163048138457556),
+    (Prior_Laplace(), 7.164837449702468),
+    (Prior_Cauchy(), 7.166326185314563),
+    (Prior_TDist(), 7.167147727917408),
+    (Prior_HorseShoe(), 7.158818008027834),
 ]
 
 for (prior, test_mean) in tests