Changeset 10e4c70 in tol for doc/modeling

Timestamp:

Jul 30, 2009, 3:55:44 PM (17 years ago)

Author:

Jorge Pérez <josp.jorge@…>

Branches:

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

Children:

87e95d3

Parents:

176c3d90

Message:

updating to new interface based on Class

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

File:

• doc/modeling/BSR/bsr.import.exa.tol (modified) (4 diffs)

doc/modeling/BSR/bsr.import.exa.tol

@@ -150 +150 @@
 };
 
+// Este seria el objeto config para la 1.1.7
 /*
 NameBlock bsr.config = [[
@@ -174 +175 @@
   Real do.linear.effects  =  True
 ]];
-
+*/
+
+// El objeto config para la 2.0.1 es de tipo Config
+BysMcmc::Config bsr.config = BysMcmc::Config.Default( 0 );
 
 Serie WN = SubSer(Gaussian(0,1),y2000,y2008);
@@ -195 +199 @@
   {Serie BETA[_k]*X[_k]})) + SubSer(Z,y2007,y2008);
 
+// Construye el NameBlock que implementa el interfaz Import
 NameBlock bsr.import = Build.Example.BSR("/tmp/", Text "session01", Serie Y, Set X);
 
+// Si quiero escribir el archivo
 Text file.bsr = BysMcmc::Bsr::Import::Write(bsr.import);
 
+// si quiero parserarlo
 Set bsr.parsed = BSR.Parse(file.bsr);
 
+//  Esta seria la forma de estimar en la version 1.1.7
+/*
 Set bsr.estim = { BysMcmc::Bsr::Gibbs::Estim
                     (
@@ -207 +216 @@
                     )};
 */
+
+// Esta seria la forma de invocar Estim en la version 2.0.1
+
+// Creo un Cycler
+
+// El partial handler es el objeto que implementara el
+// condicionamiento, para este ejemplo usamos uno vacio.
+BysMcmc::McmcPartialHandler hdPartial = BysMcmc::McmcPartialHandler.Default( 0 );
+
+// El notifiar permite emitir notificaciones personalizadas de las
+// etapas del gibbs.
+BysMcmc::Notifier notifier = BysMcmc::NoNotifier( 0 );
+
+// El cycler se puede construir de 3 maneras
+
+BysMcmc::Cycler cycler = BysMcmc::Bsr::Gibbs::GetCycler( bsr.import, hdPartial, bsr.config, notifier );
+//Cycler cycler = GetCycler( file.bsr, hdPartial, bsr.config, notifier );
+//Cycler cycler = GetCycler( bsr.parsed, hdPartial, bsr.config, notifier );
+
+// ahora la estimacion
+BysMcmc::Estimation EE = BysMcmc::Estim( cycler );