id	summary	reporter	owner	description	type	status	priority	milestone	component	resolution	keywords	cc
932	Resultado no esperado en estimación de parámetros	Ines Miranda	Pedro Gea	"Hola MMS,

para evitar el error de Choleski cuando tenemos poco histórico y no estamos midiendo bien los parámetros, hemos introducido en la definición de los parámetros, información  a priori ""fijando"" la media y dejando la sigma libre ""pero controlada"", de la siguiente manera:

{{{
// adiciona os ExpTer
  Text vble = ""Esp_SexVesDiaDosPais"";
  NameBlock expterm = submodel::CreateExpTerm_TransferFunction([[
    Text _.name = vble;
    NameBlock _.input = [[
      Text _.name = vble;
      Text _.variableIdentifier = vble
    ]];    
    Polyn _.transferFunction = 1+B+B^2+B^3
  ]]);

  // restrição de sinal. Efeito positivo
  Set param = expterm::GetParametersLinear(?);
  Set EvalSet(param, Real(NameBlock par) {
    Real par::SetInitialValue(0.01);
    Real par::SetConstraint([[
      Real _.inferiorValue = 0;
      Real _.superiorValue = ?
    ]])
  });  
  //Info a priori com sigma variavel e meia fixa
  Set EvalSet(param, NameBlock(NameBlock par) {
    MMS::@Hierarchy priorConSigmaVariable = model::CreateHierarchy([[
      Text _.name = Replace(par::GetIdentifier(?), ""__"", ""_"")<<""_Prior"";
      Set _.mElements = [[ par ]];
      NameBlock _.noise = [[
        Text _.type = ""Normal"";
        Real _.sigma = 10;
        Real _.sigmaFixed = 0.1
      ]]
    ]]);
    MMS::@HierarchyTerm  priorConSigmaVariable::CreateHierarchyTerm([[
      Text _.name = ""Mean"";
      Real _.initialParameter = 0.01; // Media
      Real _.isFixed = True
    ]])
  }); 
  1
}}}

Para este ejemplo en particular, el resultado del parámetro en algún caso (para algún output) está saliendo en torno a 0.073, pero en otros casos, nos están saliendo en torno a 0.5.

Puesto que estamos ""fijando"" la media en 0.01, a pesar de que la sigma esté libre y bastante amplia, no esperaría que me salieran parámetros tan alejados de la media definida.

Estamos definiendo bien nuestro objetivo o hay algun error en el proceso de estimación?

Si necesitáis hacer alguna prueba, puedo abrir un ticket en el trac privado y colocar un objeto forecast para que lo veais.

Muchas gracias!!
"	doubt	closed	blocker	Release 0.6	Parameters	fixed	sigma variable, prior, parametros	fortra.vivbra@…