Changeset 86179b6 in tol

Timestamp:

Nov 4, 2008, 3:42:56 PM (18 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

aee0562

Parents:

fb0fadd

Message:

Allowing fixed sigma in observational nodes

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

Location:

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/dynhlm

Files:

• _db_api.load.tol (modified) (1 diff)
• _db_api.tools.modses.obs.tol (modified) (4 diffs)

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/dynhlm/_db_api.load.tol

@@ -223 +223 @@
       db_output->missingPriorSigmaFactor, 
       db_output->missingPriorLowerBound, 
-      db_output->missingPriorUpperBound), 
+      db_output->missingPriorUpperBound),
+    db_sigma,
     arima)};
 //WriteLn("TRACE LoadObsNode 6");

tol/stdlib/math/stat/models/bayesian/bysMcmc/bsr/dynhlm/_db_api.tools.modses.obs.tol

@@ -17 +17 @@
                               Real mis_pri_max,
                               Date start,
-                              Date end)
+                              Date end,
+                              Real sigma)
 //////////////////////////////////////////////////////////////////////////////
 {
@@ -28 +29 @@
   SqlFormatReal(mis_pri_max)+","+
   SqlFormatDate(start)+","+
-  SqlFormatDate(end));
+  SqlFormatDate(end)+","+
+  SqlFormatReal(sigma)+")");
   Eval("ModSes.Obs.Output.Create.ok=ok")
 };
@@ -82 +84 @@
   SqlFormatReal(mis_pri_max));
   Eval("ModSes.Obs.Input.Create.ok=ok")
+};
+
+//////////////////////////////////////////////////////////////////////////////
+//Inserts a new latent sigma block register
+Real ModSes.Obs.SigBlk.Create(Text dbName, 
+                              Text model, 
+                              Text session,
+                              Text node, 
+                              Real sigma)
+//////////////////////////////////////////////////////////////////////////////
+{
+  NameBlock _table = BysMcmc::Bsr::DynHlm::DBApi::_.table;
+  Real ok = BysMcmc::Bsr::DynHlm::DBApi::SqlExecTrace(
+  "INSERT INTO "+dbName+"."+_table::ObsSigBlk+"\n"+
+  "VALUES('"+model+"','"+session+"','"+node+"',"+
+  If(IsUnknown(sigma),"NULL",BysMcmc::Bsr::Import::RealTextS(sigma))+")");
+  Eval("ModSes.Lat.SigBlk.Create.ok=ok")
 };
 
@@ -176 +195 @@
       out::_.mis_pri_max,
       out::_.firstDate,
-      out::_.lastDate
+      out::_.lastDate,
+      out::_.sigma
     )
   },