179 | | * ''División en bloques de Gibbs'': Debe haber un '''proceso maestro''' (master) que se encargue de disparar |
180 | | la generación de cada bloque según una estrategia dada (secuencial, aleatoria o arbitraria). Al |
181 | | contrario de lo que ocurre con el actual BSR, el proceso maestro no necesitará conocer la |
182 | | semántica de cada uno de los bloques, no hace distingos entre bloques lineales, ARIMA o de |
183 | | varianzas, sino que le basta con que cada uno tenga un método capaz de generar una muestra |
184 | | de forma condicionada al resto de bloques. Es necesario crear un |
185 | | [wiki:TolDocBaySimGibbsGenericConditioning protocolo básico de condicionamiento entre bloques] |
186 | | que sea totalmente general y no implique que ni el master |
187 | | ni los propios bloques tengan que conocer los intríngulis de los demás, sino tan solo una API |
188 | | mínima. Hay que tener en cuenta que la acción a llevar a cabo para condicionar |
189 | | puede suponer desde modificar simplemente un valor escalar (como la sigma de un bloque de |
190 | | regresión lineal normal), hasta modificar total o parcialmente un valor vectorial o matricial, |
191 | | (como ocurre en el caso de los omitidos que deben trasladarse a un subconjunto arbitrariamente |
192 | | salteado de entre las celdas de la matriz de input o de output). Además, los valores con los |
193 | | que hay que sustituir los valores a condicionar no tienen porqué ser directamente los valores |
194 | | condicionadores, sino que en general serán el resultado de una función dependiente de los mismos. |
195 | | Por estos motivos no resulta ni mucho menos trivial condicionar de una forma genérica y eficiente. |
196 | | También habría que copiar el sub-vector del bloque recién generado en la posición correspondiente |
197 | | dentro del vector conjunto ubicado en el master, que será el que se almacene en el disco como cadena |
198 | | MCMC. |
199 | | * Generación dentro de cada bloque tipo ''Metropolis-Hastings'' y derivados mediante {{{BysSampler}}} o |
200 | | cualquier otra metodología alternativa con tal cumpla una API mínima, pero siempre basado únicamente |
201 | | en el logaritmo, salvo una constante, de la función de densidad condicionada por el resto de bloques. |
202 | | La parte más complicada a la hora de aplicar el condicionamiento es la que se refiere a la información |
203 | | a priori y restricciones que afectan a variables de distintos bloques. |