Changeset 3e558bb in tol for tol_tests

Timestamp:

Apr 12, 2010, 2:09:09 PM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

a30612b

Parents:

e06a79e

Message:

accept_reject debe llamar siempre a doAfterDraw, se acepte o no el candidato, y debe hacerlo después de actualizar los indicadores de aceptacón.

No es necesario recalcular el Set _.arima ni el NameBlock _.almagro si no se ha aceptado el nuevo candidato, pero sí se deben generar nuevos valores iniciales, puesto que en realidad se trata de un bloque de Gibss aparte

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

File:

• tol_tests/tol/ARIMA/mcmc/mtrymh.tol (modified) (6 diffs)

tol_tests/tol/ARIMA/mcmc/mtrymh.tol

@@ -134 +134 @@
     {
       VMatrix _.x := Copy(y);
-      Real doAfterDraw(0);
       Real _.num_accept := _.num_accept +1
     });
@@ -147 +146 @@
       Real _.s := _.s*_.sFactor
     });
+  //Real If(_.accept, { Real doAfterDraw(0) });
+    Real doAfterDraw(0);
     _.accept
   };
@@ -507 +508 @@
   ////////////////////////////////////////////////////////////////////////////
   {
-    Polyn  _.ar := ARIMAGetAR(_.arima);
-    Polyn  _.ma := ARIMAGetMA(_.arima);
-    NameBlock _.almagro := ARMAProcess::Eval.Almagro(_.ar,_.ma,_.z,_.sigma);
+    Real If(_.accept,
+    {
+      Polyn  _.ar := ARIMAGetAR(_.arima);
+      Polyn  _.ma := ARIMAGetMA(_.arima);
+      NameBlock _.almagro := ARMAProcess::Eval.Almagro(_.ar,_.ma,_.z,_.sigma);
+      True
+    });
     VMatrix _.u0 := _.almagro::Draw.U_cond_Z(0);
     VMatrix _.z0 := _.almagro::Get.Z0(_.u0);
@@ -525 +530 @@
       SetOfPolyn(1-f->AR, 1-f->MA)
     });
-        VMatrix _.x := arima2col(_.arima);
+    VMatrix _.x := arima2col(_.arima);
+    Real _.accept := True;
     Real ok = setArima(0);
     Real _.n := VRows(_.x);
@@ -557 +563 @@
   ////////////////////////////////////////////////////////////////////////////
   {
-    Set _.arima := col2arima(_.x);
-    Real setArima(void);
-    Real doAfterDraw_cum(void)
+    Set If(_.accept, _.arima := col2arima(_.x));
+    Real setArima(void)
+//  Real doAfterDraw_cum(void)
   };
 
@@ -698 +704 @@
 Real TolOprProfiler.Enabled := True;
 While(numSim<2500, append_draw(mtm_arima::draw_multiple_try(0)));
-Real TolOprProfiler.Dump("draw_multiple_try.03.csv");
+Real TolOprProfiler.Dump("draw_multiple_try.csv");
 Real mtm.time = (Copy(Time) - mtm.t0)/numSim;
 WriteLn("Time by MTM simulation "<<mtm.time);
 /* */
 
-Matrix real_param = Group("ConcatColumns", For(1,Columns(PARAM),Matrix(Real j) 
+Matrix cmp = mcmc | Group("ConcatColumns", For(1,Rows(PARAM),Matrix(Real j) 
 {
-  Constant(Rows(PARAM),1,MatDat(PARAM,j,1))
+  Constant(Rows(mcmc),1,MatDat(PARAM,j,1))
 }));
 
 /* */
 
+
+