Changeset 12df85b in tol

Timestamp:

Nov 28, 2008, 2:02:53 PM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

d29988b

Parents:

fdc8a04

Message:

Handling with normal prior information about parameters of non linear filters

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@236 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol (modified) (3 diffs)

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/gibbs/_nonLinBlock.tol

@@ -81 +81 @@
     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 +164 @@
       Set bounds = _.filterEval::get.bounds(i,u);
       Code fltEv = _.filterEval::eval;
+      Code get.prior = _.get.prior;
       Matrix Z = _.Z;
       Real x0 = bounds[1];
@@ -177 +181 @@
         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)