Changeset 1abcc87 in tol

Timestamp:

Oct 26, 2008, 10:45:11 AM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

0691acd

Parents:

8349bbb

Message:

Preparing for common input missing

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

File:

• tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/dynhlm/_build.node.obs.tol (modified) (12 diffs)

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/dynhlm/_build.node.obs.tol

@@ -20 +20 @@
   //The last date to be analyzed
   Date last, 
+  //Missing label is used as prefix of missing parameters
+  Text missing.label,
   //If there are missing values it will be added prior information about it
   //multiplying the standard deviation of the serie by this factor. 
@@ -38 +40 @@
   {
     Real aux = AvrS(_.serie);
-    If(IsUnknown(aux), WriteLn("Serie "+_.name+" is fully unknown","W"));
+    If(IsUnknown(aux), WriteLn("Serie "+_.id_node+"::"+_.name+
+      " is fully unknown","W"));
     aux
   };
   Real _.stds = StDsS(_.serie);
+  Text _.missing.label       = missing.label;
+  Text _.missing.name        = _.missing.label+"_"+_.name;
   Real _.missing.sigmaFactor = missing.sigmaFactor; 
   Real _.missing.lowerBound  = missing.lowerBound; 
@@ -80 +85 @@
   {
     Text txt = FormatDate(dte,"%cy%Ym%md%dh%hi%is%s");
-    "MissingBlk::"+_.id_node+"::"+_.name+"::"+txt
+    "MissingBlk::"+_.missing.label+"::"+_.name+"::"+txt
   };
   //Returns the expression of the t-th regression coefficient that could be a
@@ -126 +131 @@
   //Time serie data
   Matrix serie, 
+  //Missing label is used as prefix of missing parameters
+  Text missing.label,
   //If there are missing values it will be added prior information about it
   //multiplying the standard deviation of the serie by this factor. 
@@ -148 +155 @@
   };
   Real _.stds = MatStDs(_.serie);
+  Text _.missing.label = missing.label;
+  Text _.missing.name  = _.missing.label+"_"+_.name;
   Real _.missing.sigmaFactor = missing.sigmaFactor; 
   Real _.missing.lowerBound  = missing.lowerBound; 
@@ -158 +167 @@
   {
     Text txt = IntText(num);
-    "MissingBlk::"+_.id_node+"::"+_.name+"::"+txt
+    "MissingBlk::"+_.missing.label+"::"+_.name+"::"+txt
   };
   //Returns the expression of the t-th regression coefficient that could be a
@@ -203 +212 @@
   //The last date to be analyzed or UnknownDate if it's a vector
   Date last, 
+  //Missing label is used as prefix of missing parameters
+  Text missing.label,
   //If there are missing values it will be added prior information about it
   //multiplying the standard deviation of the serie by this factor. 
@@ -221 +232 @@
   //WriteLn("TRACE [Obs.Data.Info] 1");
     Obs.Serie.Info(
-      id_node, name, serie, numIniVal, first, last, 
+      id_node, name, serie, numIniVal, first, last, missing.label,
       missing.sigmaFactor, missing.lowerBound, missing.upperBound)
   },
@@ -228 +239 @@
   //WriteLn("TRACE [Obs.Data.Info] 2");
     Obs.Vector.Info(
-      id_node, name, serie, 
+      id_node, name, serie, missing.label,
       missing.sigmaFactor, missing.lowerBound, missing.upperBound)
   },
@@ -335 +346 @@
     (
       id_node+"::Noise",
-      _.transferFunction,
+      Copy(_.transferFunction),
       EvalSet(_.omegaParam, Text(Bsr.Param.Info pi) { pi->Name}),
-      _.omega, 
-      _.delta,
+      Copy(_.omega), 
+      Copy(_.delta),
       If(!_.p,Rand(0,1,0,0), Tra(SerSetMat([[ser]], fst0, lst0))),
       Tra(SerSetMat([[ser]], fst,lst))
@@ -370 +381 @@
   //A set where each item is the result of calling Obs.TransFun
   Set transFun,
+  //Acumulate input.missing.labels to avoid repetitions
+  Set input.missing.labels,
   //The result of calling BysMcmc::Bsr::Import::Constraints
   NameBlock constraints,
@@ -416 +429 @@
     BinGroup("<<", EvalSet(_.input, Set(NameBlock inp)
     {
-      inp::_.serie.info::_.missing
+      Set label = [[inp::_.serie.info::_.missing.name]];
+      Text PutName(label[1],label[1]);
+      Real idx = FindIndexByName(input.missing.labels,label[1]);
+    //WriteLn("TRACE Obs Missing Label Name : "+label[1]+" -> "<<idx);
+      If(idx, Copy(Empty),
+      {
+        Set Append(input.missing.labels, label);
+        Real SetIndexByName(input.missing.labels);
+        inp::_.serie.info::_.missing
+      })
     }));