Index: tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_basicMaster.tol
===================================================================
--- tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_basicMaster.tol	(revision 0d654d62735a661c578b2f49da57bdd035be7184)
+++ tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_basicMaster.tol	(revision ff5d8e229af5a6432915082010bf1211b3b0e8c5)
@@ -75,5 +75,7 @@
       Real filter.blk(BysMcmc::@BlockSampler extFlt)
       {
-        extFlt::do.filter(_.Y)
+        Real ok = extFlt::do.filter(_.Y);
+      //WriteLn("TRACE"+_MID+"DoNonLinOutputFilter 4 _.Y=\n"<<VMat2Mat(_.Y,1));
+        ok
       };
       SetMin(EvalSet(_.nonLinOutputFilter,filter.blk))
@@ -149,4 +151,5 @@
   ////////////////////////////////////////////////////////////////////////////
   {
+  //WriteLn("TRACE"+_MID+" define.all 0");
     Set _.nonLinBlk := 
     {
@@ -202,6 +205,11 @@
   {
   //WriteLn("TRACE"+_MID+" initialize 0");
+    Matrix givenBeta=
+    SetCol(EvalSet(_.modelDef->LinearBlock,Real(@BSR.LinearBlock b)
+    { b->InitialValue }));
+    VMatrix beta0 = findMatchingBeta(Mat2VMat(givenBeta));
+    VMatrix _.filter := _.modelDef->X*beta0;
     Real initialize.blocks(unused);
-
+  //WriteLn("TRACE"+_MID+"initialize _.filter=\n"<<Matrix VMat2Mat(_.filter,1));
     Real sel.extFlt (BysMcmc::@BlockSampler nlb) { !(nlb::_.isInputFilter) };
     Set _.nonLinOutputFilter := Select(_.nonLinBlk, sel.extFlt);
@@ -215,8 +223,4 @@
     VMatrix _.XplusU := _.mis.blk::filter.input(_.modelDef->X);
   //WriteLn("TRACE"+_MID+" initialize 1.2 _.XplusU("<<VRows(_.XplusU)+"x"<<VColumns(_.XplusU)+")");
-    Matrix givenBeta=
-    SetCol(EvalSet(_.modelDef->LinearBlock,Real(@BSR.LinearBlock b)
-    { b->InitialValue }));
-    VMatrix beta0 = findMatchingBeta(Mat2VMat(givenBeta));
   //WriteLn("TRACE"+_MID+" initialize 3 beta0("<<VRows(beta0)+"x"<<VColumns(beta0)+")");
     Real ok = _.mln.blk::matchConstrains(beta0, _.modelDef->A, _.modelDef->a);
Index: tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol
===================================================================
--- tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol	(revision 0d654d62735a661c578b2f49da57bdd035be7184)
+++ tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol	(revision ff5d8e229af5a6432915082010bf1211b3b0e8c5)
@@ -243,8 +243,9 @@
   //WriteLn("TRACE "+_MID+"::setStore 1 values("<<Rows(values)+"x"<<Columns(values)+")");
     Real putStore(values);
+  //WriteLn("TRACE "+_MID+"::setStore 3");
+    Matrix _.Fi := _.filter::eval(VMat2Mat(_.store));
+  //WriteLn("TRACE "+_MID+"::setStore 3.1");
     Real If(_.N==0, True,
     {
-    //WriteLn("TRACE "+_MID+"::setStore 3");
-      Matrix _.Fi := _.filter::eval(VMat2Mat(_.store));
     //WriteLn("TRACE "+_MID+"::setStore 4");
       Real ok = If(!_.isInputFilter, 0,