Errores en el ejemplo: Basics A

[josp]

En samples/Basics/A2_B_03_estimation_bsr.tol veo mensajes como

<W>
Warning: [3] [@BlockSampler @BSR.MainLinBlock.A2_BSR__1.0] [_.freeBlock] Cholesky solving had rounding error in interval [6.000000000000001e-08,16 for 66 times of 300
</W>
...
el error va creciendo
....
<W>
Warning: [20] [@BlockSampler @BSR.MainLinBlock.A2_BSR__1.0] [_.freeBlock] Cholesky solving had rounding error in interval [6.000000000000001e-08,201326592 for 742 times of 2000
</W>

Ese termina

Pero este samples/Basics/A3_B_03_estimation_bsr.tol deja esta traza

Incluyendo /home/jsperez/TOL/MMS/MMS/samples/Basics/A3_B_03_estimation_bsr.tol
Model - 0     (D: 00.00s     T: 00.00s)
Model - 1     (D: 00.01s     T: 00.01s)
DataSet - 0     (D: 00.00s     T: 00.00s)
DataSet - 1     (D: 00.00s     T: 00.00s)
DataSet - 2     (D: 00.00s     T: 00.01s)
DataSet - 3     (D: 00.02s     T: 00.03s)
DataSet - X     (D: 00.00s     T: 00.03s)
Model - 2     (D: 00.04s     T: 00.05s)
Model - 3     (D: 00.01s     T: 00.06s)
Model - 4     (D: 00.02s     T: 00.08s)
Model - 5     (D: 00.00s     T: 00.09s)
Model - 6     (D: 00.00s     T: 00.09s)
Model - 7     (D: 00.00s     T: 00.09s)
Model - X     (D: 00.00s     T: 00.10s)
[Bsr::Import] Write ...
[Bsr::Import] Write.Module master ...
[Bsr::Import]  Writing master module
  "/home/jsperez/.tol/MMS/BSR/A3_BSR__1.0/2016.02.18_14.05.20/master.bsr"
[Bsr::Import] Write.Module primary ...
[Bsr::Import]  Writing primary module
  "/home/jsperez/.tol/MMS/BSR/A3_BSR__1.0/2016.02.18_14.05.20/Submodels/OutputA/primary.bsr"
<E>
ERROR: [1] Intento de acceso al elemento 1-esimo de un conjunto de 0 elementos</E>

[Call stack]
  [19] Set estimationA3::estimation::model::submodel::submodel::submodel::SubData_Degrees. (Set data., Real minDegree, Real maxDegree)
  [18] Set estimationA3::estimation::model::submodel::expTerm::expTerm::expTermRatio::SubData. (Set data.)
  [17] Set estimationA3::estimation::model::submodel::expTerm::expTerm::expTermRatio::GetInitialFilter. (Real void)
  [16] Anything estimationA3::estimation::model::submodel::expTerm::expTerm::expTermRatio::GetInitialFilter (Real void)
  [15] Matrix cycler::cycler::segment::submodelImport::_Get.Matrix.X (Real void)
  [14] Anything cycler::cycler::segment::submodelImport::_UseCache (Code method)
  [13] Matrix cycler::cycler::segment::submodelImport::Get.Matrix.X (Real void)
  [12] VMatrix cycler::cycler::segment::submodelImport::Get.InputVMatrix (Real void)
  [11] Real BysMcmc::Bsr::Import::Write.Data.RegressionMatrices (NameBlock segment, Real fileHandler, Text path.prefix)
  [10] Text BysMcmc::Bsr::Import::Write.Primary (NameBlock model)
  [9] Text BysMcmc::Bsr::Import::Write.Module (NameBlock model, Text moduleType)
  [8] Text BysMcmc::Bsr::Import::Write.Master (NameBlock model)
  [7] Text BysMcmc::Bsr::Import::Write.Module (NameBlock model, Text moduleType)
  [6] Text BysMcmc::Bsr::Import::Write (NameBlock model)
  [5] NameBlock BysMcmc::Bsr::Gibbs::GetCycler (Anything modelDef, NameBlock mcmcPartialHandler, NameBlock config, NameBlock notifier)
  [4] NameBlock estimationA3::estimation::strategy::strategyBSR::_BuildCycler (Anything void)
  [3] Real estimationA3::estimation::strategy::strategyBSR::Prepare (Real void)
  [2] Real estimationA3::estimation::strategy::strategyBSR::Execute (Real void)
  [1] Real space::_E::Execute (Real void)
<W>
Warning: [22] La función Dating ha fallado en : 
Dating(data.[[1])</W>
<E>
ERROR: [2] Argumentos erróneos para Date Succ (Date fecha1, TimeSet fechado, [Real n=1])</E>

luego termina con SEGFAULT

[Pedro Gea]

El primer caso (A2) es un problema de mal diseño del modelo, un parámetro de omitido no converge y genera esos errores.
Si se ponen prior a los omitidos, por ejemplo:

Set EvalSet(modelA2::GetMVariables(?), Real (NameBlock mv) {
  Real mvAvr = AvrS(mv::GetData(100));
  Real mvStd = StDsS(mv::GetData(100));
  mv::SetBaseParameterMissing([[
    NameBlock _.prior = [[
      Real _.mean = mvAvr;
      Real _.sigma = mvStd
    ]]
  ]]);
1});

se soluciona.

[Pedro Gea]

(In [4792]) Refs #1226
Se añaden priors para los omitidos que se construyan.

[Pedro Gea]

El segundo error corresponde a un desarrollo sin completar. Véase [4723].
En general el tratamiento de los términos de tipo Ratio tiene algunos problemas de ambigüedad en el tratamiento de los parámetros de los polinomios así como de los parámetros para la integración del input.

Los términos explicativos no lineales, más generales, disponen de dos variantes de ya prediseñadas:

• MMS::@NonLinearFilter.Delta: que permite estimar ratios con ganancia unitaria del tipo: (1-delta)/(1-delta*B) y con un parámetro de integración conocido como "delta-missing".
• MMS::@NonLinearFilter.SDelta0: que permite estimar ratios con ganancia unitaria y estacionalidad del tipo: (1-delta)/(1-delta*B^period) y sin parámetros de integración (se suponen ceros).

[Pedro Gea]

(In [4793]) Refs #1226
Se reescribe el modelo para poder ser estimado con los desarrollos actuales.

[Pedro Gea]

Sobre la evolución del desarrollo de los términos explicativo de tipo Ratio véase #899.