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 |