id	summary	reporter	owner	description	type	status	priority	milestone	component	version	severity	resolution	keywords	cc
951	No se satisfacen restricciones de dominio  con parámetros fijos en BSR.	mafernandez	Víctor de Buen Remiro	"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(?);
}}}

  "	defect	closed	highest		BSR	head	major	fixed