Distribución a priori semidefinida en los parámetros

[Pedro Gea]

El diseño de MMS_0.6 impide la construcción de priors semidefinidos (con media o varianza desconocida). Para ello se propone crear constructores de jerarquías específicos para esta situación.

[mafernandez]

(In [1641])
closes #313

Se ha creado el método de model CreatePriorFreeSigma2 para implementar los priors con sigma desconocida. El uso es el siguiente:

@MMS.Hierarchy hierarchyp = model::CreatePriorFreeSigma2(

Set _.element
,
Real mean);

Donde el Set _.element especifica el parámetro y Real mean es el valor promedio de la Normal del prior.

Como ejemplo para el modelo de matriculación de vehículos podríamos tener el siguiente prior semidefinido

@MMS.Hierarchy hierarchyp = model::CreatePriorFreeSigma2(

"Veh.Bus.Mat", "Pib.es.Men", [["Linear", 0? ]]
,
Real mean = 1.0 );

[Pedro Gea]

(In [1643]) Cambios menores en la implementación de priors semidefinidos.
Refs #313

[Pedro Gea]

Se solicita considerar la posibilidad de ponerle prior a la sigma.

[Pedro Gea]

Se incorpora un trozo de código que permite hacer esovía una jerarquía con hiperparámetro fijo:

MMS::@Model model = ...;
MMS::@Parameter parameter = ...;

MMS::@Hierarchy priorConSigmaVariable = model::CreateHierarchy([[
  Text _.name = Replace(parameter::GetIdentifier(?), "__", "_")<<"_Prior";
  Set _.mElements = [[ parameter ]];
  NameBlock _.noise = [[
    Text _.type = "Normal";
    Real _.sigma = 0.1;
    Real _.sigmaFixed = 0.5
  ]]
]]);
MMS::@HierarchyTerm mediFijaDelPrior = priorConSigmaVariable::CreateHierarchyTerm([[
  Text _.name = "Mean";
  Real _.initialParameter = 0.5; // Media
  Real _.isFixed = True
]]);