No se satisfacen restricciones de dominio con parámetros fijos en BSR.

[mafernandez]

Expongo un modelo sencillo que consiste en un output y dos inputs.

El modelo posee una restricción de dominio que obliga al parámetro lineal del segundo input a ser < 0. Hasta aquí todo bien, ese parámetro estimado resulta ser negativo.

Sin embargo cuando se pide que el parámetro del primer input quede fijo a un valor dado entonces el parámetro del segundo input deja de satisfacer la restricción que se le pide.

Dentro del ejemplo que adjunto, los parámetros son:

OutputD__ExpTermD1__Deg0.Coef;
OutputD__ExpTermD2__Deg0.Coef;

La restricción es la siguiente:

-10<=+1*OutputD__ExpTermD2__Deg0.Coef;
+0>=+1*OutputD__ExpTermD2__Deg0.Coef;

El parámetro que se fija es OutputD__ExpTermD1__Deg0.Coef a través de BysMcmc::@McmcPartialHandler

Este es el código que uso para realizar la prueba:

Text ascii = "C:/y2010m07d27h12i35s14.00/master.bsr";
Text moduleType = "master";

Set parsed = BSR.Parse(ascii,moduleType);
Set fixedNames  = [[ "OutputD__ExpTermD1__Deg0.Coef" ]];

VMatrix  fixedValues = Mat2VMat(Col(0.12));

NameBlock _.hdPartial  = BysMcmc::@McmcPartialHandler::New(fixedNames, fixedValues, Set Empty, BysMcmc::@McmcReloader::Null(?)); 
NameBlock _.bsr.config   = BysMcmc::@Config::Default(?);
NameBlock _.bsr.notifier = BysMcmc::@Notifier::Null(?);

BysMcmc::@Cycler cycler = BysMcmc::Bsr::Gibbs::GetCyclerFromParsed(
parsed, 
_.hdPartial, 
_.bsr.config, 
_.bsr.notifier
);

BysMcmc::@Estimation estim := BysMcmc::@Estimation::New(cycler);
Real estim::Run(?);

[Víctor de Buen Remiro]

Falta definir el fechado DatingD

[mafernandez]

Replying to vdebuen:

Falta definir el fechado DatingD

En efecto, TimeSet DatingD = D(1);

[Víctor de Buen Remiro]

(In [2414]) Fixes #951
In method @BlockStdLin::matchConstrains, line

Set _.r2 = Select(nz_idx,Real(Real j) { j>0 });

should use := instead of =

[Víctor de Buen Remiro]

(In [2417]) New package version: 1.2
Solved problem with constrained linear block in partial simulation (Refs #951)

[Víctor de Buen Remiro]

(In [2541]) refs #951

[Víctor de Buen Remiro]

(In [3091]) Refs #951
BSR can show warnings

[Víctor de Buen Remiro]

(In [3145]) Refs #951

[Víctor de Buen Remiro]

(In [4092]) refs #951
upgrading test