Index: tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol
===================================================================
--- tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol	(revision fdc8a04705533397e975cd36f7665e7862cbb5a0)
+++ tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol	(revision 12df85b306e13d5485fdf41d91946f5aed84b6b0)
@@ -81,4 +81,7 @@
     0
   });
+  Code _.get.prior = If(FunctionExist("Matrix","_.filterEval::get.prior"),
+    _.filterEval::get.prior,
+    Matrix (Real paramIdx) {Rand(0,1,0,0) } );
   //Names of affected parameters in linear block 
   //(just for input filters)
@@ -161,4 +164,5 @@
       Set bounds = _.filterEval::get.bounds(i,u);
       Code fltEv = _.filterEval::eval;
+      Code get.prior = _.get.prior;
       Matrix Z = _.Z;
       Real x0 = bounds[1];
@@ -177,5 +181,6 @@
         Matrix Fi = fltEv(u);
         Real Do PutMatDat(u,i,1,old);
-        VMat2Mat(si*Mat2VMat(Z-Fi*omega))
+        Matrix fn = VMat2Mat(si*Mat2VMat(Z-Fi*omega));
+        fn << get.prior(i)
       };
       Real LogDens.WhiteNoise(Real x)