close
Warning:
Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.
- Timestamp:
-
Jul 5, 2011, 10:15:35 PM (14 years ago)
- Author:
-
Víctor de Buen Remiro
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
|
v3
|
v4
|
|
| 18 | 18 | |
| 19 | 19 | Si se conoce un método eficaz para simular cada bloque condicionado al resto entonces es aplicable |
| 20 | | el método de simulación de Gibbs. Si no es el caso pero sí es posible evaluar la función de |
| 21 | | log-verosimilitud salvo constante condicionada al resto de bloques |
| | 20 | el método de simulación de Gibbs por bloques. |
| | 21 | |
| | 22 | Normalmente resulta mucho más sencillo evaluar la función de log-verosimilitud salvo constante condicionada |
| | 23 | al resto de bloques, |
| 22 | 24 | |
| 23 | 25 | [[LatexEquation( \ln\pi_{i}\left(\left.\Omega_{i}\right|\Omega_{j\neq i}\right) )]] |
| 24 | 26 | |
| 25 | | entonces el método a aplicar se llama ''Metropolis within Gibbs'', que consiste en simular |
| 26 | | alternativamente los distintos bloques con el método de Metropolis-Hastigs (MH) o alguno similar |
| 27 | | como el Multiple Try Metroplis (MTM), dejando invariantes todos aquellos de los que depende. En cada |
| 28 | | simulación es necesario evaluar la log-verosimilitud un número de veces dependiente del método usado. |
| | 27 | lo cual permite aplicar el método Metropolis Within Gibbs, que consiste en simular |
| | 28 | alternativamente los distintos bloques con el método de Metropolis-Hastigs (MH) o alguno similar como el |
| | 29 | Multiple Try Metroplis (MTM), dejando invariantes todos aquellos de los que depende. En cada simulación |
| | 30 | es necesario evaluar la log-verosimilitud un número de veces dependiente del método usado. |
| | 31 | |
| | 32 | Este esquema tiene la ventaja de que sirve para simular cualquier clase de modelos implementando bloques de |
| | 33 | simulación para cada tipo de elemento que intervenga y definiendo aparte las relaciones de condicionamiento |
| | 34 | como meras reglas de asignación total o parcial de miembros de un bloque condicionado en función de los |
| | 35 | miembros del bloque condicionante. |
| | 36 | |
| | 37 | Si para un bloque concreto se dispone de un método de generación condicionada que resulte más sencillo o |
| | 38 | eficaz que el MH o sus derivados, es perfectamente lícito aplicarlo siempre y cuando se respeten las |
| | 39 | relaciones de condicionamiento. |
| 29 | 40 | |
| 30 | 41 | == Arquitectura del sistema == |
| … |
… |
|
| 67 | 78 | mucho menos trivial condicionar de una forma genérica y eficiente. |
| 68 | 79 | |
| 69 | | Existen pues tres entidades en el sistema, cada una de las cuales dará lugar a una jerarquía de clases |
| 70 | | que permita adaptarse a todas las situaciones con una API mínima. |
| | 80 | Existen varias entidades en el sistema, algunas de las cuales darán lugar a una jerarquía de clases |
| | 81 | que permitirán adaptarse a todas las situaciones con una API mínima, mientras que otros conceptos son |
| | 82 | relativos al uso y no a su estructura por lo que no precisan de una clase |
| 71 | 83 | |
| 72 | 84 | * {{{@RandVar}}}: Variable aleatoria generable mediante log-verosimilitud, sin perjuicio de que pueda haber |
| 73 | | otro método de generación alternativo más eficaz en alguna clase heredada. |
| | 85 | otro método de generación alternativo más eficaz en alguna clase heredada. |
| 74 | 86 | * {{{@Relation}}}: Relación entre un bloque condicionante y otro condicionado que especifica sin ambigüedad |
| 75 | 87 | posible todas las acciones a realizar. |
| 76 | | * {{{@Model}}}: Conjunto de variables aleatorias y de relaciones de condicioanmiento. Es un tipo especial |
| 77 | | de variale aleatoria, por lo que hereda de {{{@RandVar}}}. |
| | 88 | * {{{@Model}}}: Conjunto de bloques de variables aleatorias y de relaciones de condicioanmiento. Es un |
| | 89 | tipo especial de variale aleatoria, por lo que hereda de {{{@RandVar}}}. |
| | 90 | * Bloque: Variable aleatoria cuyo generador es llamado por un modelo según una estrategia de ejecución. |
| | 91 | de algún modo compatible con el algoritmo de Gibbs. Un modelo puede funcionar como un bloque de otro modelo. |
| | 92 | * Estrategia: Forma de alternar o simultanear la simulación de los distintos bloques de Gibbs, de forma |
| | 93 | compatible con la estructura de relaciones. |
| | 94 | * Master: Modelo principal formado por bloques que pueden ser a su vez modelos y que lanza las simulaciones |
| | 95 | de cada uno de ellos según la estrategia seleccionada y la va almacenando en disco. |
| 78 | 96 | |
| 79 | 97 | == Clase {{{@RandVar}}} == |
| 80 | | A la clase base de la que derivan todos los bloques la llamaremos {{{Class @RandVar}}}, pues lo que |
| 81 | | se está planteando no es otra cosa que un mecanismo de gestión de variables aleatorias. |
| | 98 | A la clase base de la que derivan todos las variables aleatorias (v.a.) la llamaremos {{{Class @RandVar}}}, pues |
| | 99 | lo que se está planteando no es otra cosa que un mecanismo de gestión de variables aleatorias por la |
| | 100 | vía de la generación de simulaciones por uno u otro método. |
| | 101 | |
| | 102 | No toda v.a. definida tiene porqué ser un bloque de Gibbs de un modelo. Muchas veces un bloque |
| | 103 | se puede definir mediante operaciones algebraicas de otras variables aleatorias las cuales no pertenecerán |
| | 104 | directamente al modelo sino indirectamente a través de uno de sus bloques. Otro tipo de v.a. de gran |
| | 105 | utilidad serán las que definan la información ''a priori'' relativa a un bloque o a un modelo, simplemente |
| | 106 | sumando su log-verosimilitud. |
| 82 | 107 | |
| 83 | 108 | En este nivel raíz estará disponible una batería de algoritmos de generación MCMC basados en la |
| … |
… |
|
| 104 | 129 | que exista otro método directo más eficaz de muestreo. |
| 105 | 130 | |
| 106 | | Cuando no existe dicho método directo la clase derivada final debe implementar obligatoriamente el |
| 107 | | método que devuelve el logaritmo de la verosimilitud para poder aplicar alguno de los métodos MCMC |
| | 131 | Cuando no existe dicho método directo la clase derivada final, ésta debe implementar obligatoriamente |
| | 132 | el método que devuelve el logaritmo de la verosimilitud para poder aplicar alguno de los métodos MCMC |
| 108 | 133 | {{{ |
| 109 | 134 | #!java |