Index: tol_pkg/BysVecLinReg/test/test_0001/test.tol
===================================================================
--- tol_pkg/BysVecLinReg/test/test_0001/test.tol	(revision 94219974c42c5f1bab43c13b538701a0afbedb9b)
+++ tol_pkg/BysVecLinReg/test/test_0001/test.tol	(revision db2b6d3d36f64d002b40ab13905683ff9256f717)
@@ -12,4 +12,5 @@
 Real numErr0 = Copy(NError);
 
+#Require BysVecLinReg;
 #Embed "../gen_multivariant_linear_regression.tol";
 
@@ -41,5 +42,4 @@
 
 VMatrix vec.A_ = Tra(Constant(1,sampleLength,1/sampleLength)*mcmc);
-VMatrix vec.A = BysVecLinReg::vmatrix2vec(A);
 VMatrix cmp.A = vec.A | vec.A_;
 
@@ -59,5 +59,5 @@
 Set partialResults_ = 
 [[
-  numErr0, numErr1, stdErr, stdErr_, quality.adjust
+  numErr0, numErr1, stdErr, stdErr_, quality.adjust,
   ok_.X.Li, ok_.IdxXL, ok_.IdxXLi, okLLt
 
Index: tol_pkg/BysVecLinReg/test/test_0002/test.tol
===================================================================
--- tol_pkg/BysVecLinReg/test/test_0002/test.tol	(revision 94219974c42c5f1bab43c13b538701a0afbedb9b)
+++ tol_pkg/BysVecLinReg/test/test_0002/test.tol	(revision db2b6d3d36f64d002b40ab13905683ff9256f717)
@@ -12,4 +12,5 @@
 Real numErr0 = Copy(NError);
 
+#Require BysVecLinReg;
 #Embed "../gen_multivariant_linear_regression.tol";
 
@@ -23,5 +24,5 @@
 
 //Real MV::set.prior.V.none(0);
-Real MV::set.prior.V.inverseWishart(0.90, V);
+Real MV::set.prior.V.inverseWishart(0.10, V);
 
 VMatrix mcmc = Group("ConcatRows", For(1,sampleLength,VMatrix(Real numSim)
Index: tol_pkg/BysVecLinReg/test/test_0003/test.tol
===================================================================
--- tol_pkg/BysVecLinReg/test/test_0003/test.tol	(revision 94219974c42c5f1bab43c13b538701a0afbedb9b)
+++ tol_pkg/BysVecLinReg/test/test_0003/test.tol	(revision db2b6d3d36f64d002b40ab13905683ff9256f717)
@@ -13,5 +13,8 @@
 Real numErr0 = Copy(NError);
 
+#Require BysVecLinReg;
+//Se generará un modelo con r.eq restricciones de igualdad lineales
 Real r.eq = 3;
+Real rndSeed = 792466945;
 #Embed "../gen_multivariant_linear_regression.tol";
 
@@ -42,4 +45,40 @@
 }));
 
+//////////////////////////////////////////////////////////////////////////////
+//Standard report
+//////////////////////////////////////////////////////////////////////////////
+
+Set report = { 
+  BysMcmc::@Config cfg = BysMcmc::@Config::Default(0);
+  BysMcmc::Inference::report(
+  MV::_.desc, 
+  Set colNames_ = SetConcat(For(1,MV::_.m,Set(Real j)
+  {
+    For(1,MV::_.d,Text(Real i)
+    {
+      "linblk_a_"<<i<<"_"<<j+"_"+MV::_.outputNames[i]+"_"+MV::_.inputNames[j]
+    })
+  })) << SetConcat(For(1,MV::_.d,Set(Real j)
+  {
+    For(1,MV::_.d,Text(Real i)
+    {
+      "sigblk_s_"<<i<<"_"<<j+"_"+MV::_.outputNames[i]+"_"+MV::_.outputNames[j]
+    })
+  })),
+  Set BysMcmc::@BysInf.Report.Config
+  (
+    cfg::report.raftery.diag.q,
+    cfg::report.raftery.diag.r,
+    cfg::report.raftery.diag.s,
+    cfg::report.raftery.diag.eps,
+    cfg::report.acf.lag,
+    cfg::report.histogram.parts,
+    cfg::report.kerDens.points
+  ),
+  Matrix VMat2Mat(mcmc))};
+
+//////////////////////////////////////////////////////////////////////////////
+//Reporting about average of simulations
+//////////////////////////////////////////////////////////////////////////////
 VMatrix avg = Tra(Constant(1,sampleLength,1/sampleLength)*mcmc);
 VMatrix vec.A_ = BysVecLinReg::vmatrix2vec(