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:
-
Feb 1, 2011, 6:35:08 PM (14 years ago)
- Author:
-
Víctor de Buen Remiro
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v17
|
v18
|
|
4 | 4 | == Descripción del problema == |
5 | 5 | Tenemos que generar una cadena de Markov de muestras que se distribuyan con cierta función de la |
6 | | que se conoce su log-densidad salvo una cosntante, a la par que cumplen un conjunto de restricciones |
| 6 | que se conoce su log-densidad salvo una constante, a la par que cumplen un conjunto de restricciones |
7 | 7 | de desigualdad lineal |
8 | 8 | |
9 | | [[LatexEquation( A x \ge a \wedge x\in\mathbb{R}^{n} \wedge a\in\mathbb{R}^{r} \wedge A\in\mathbb{R}^{r\times n} )]] |
| 9 | [[LatexEquation( A x \ge a \wedge x\in\mathbb{R}^{n} )]] |
| 10 | |
| 11 | [[LatexEquation( x\in\mathbb{R}^{n} )]] |
| 12 | |
| 13 | [[LatexEquation( a\in\mathbb{R}^{r} )]] |
| 14 | |
| 15 | [[LatexEquation( A\in\mathbb{R}^{r\times n} )]] |
| 16 | |
10 | 17 | |
11 | 18 | Lo único que deben cumplir las restricciones es que el politopo generado por ellas, o sea, |
… |
… |
|
14 | 21 | ser singular o tener cualquier forma con tal cumpla esa condición. En particular puede haber |
15 | 22 | restricciones no activas, es decir, que se pueden quitar sin que cambie en absoluto la forma del |
16 | | politopo generado por las restricciones. |
| 23 | politopo generado. |
| 24 | |
| 25 | Para obtener un elemento de la cadena dado el anterior se utilizará algún método similar al |
| 26 | Metropolis-Hastings, como podría ser el Multiple Try Metropolis. Todos ellos necesitan un método |
| 27 | de generación aleatoria de candidatos o precandidatos que luego serán aceptados o no. |
17 | 28 | |
18 | 29 | Nos caben dos posibilidades a la hora de diseñar el generador de candidatos: |
… |
… |
|
46 | 57 | |
47 | 58 | [[Image(source:/tolp/OfficialTolArchiveNetwork/BysSampler/doc/image/RandWalk.InPolytope.chart.726488582.png)]] |
| 59 | |
| 60 | En la práctica se observa que cuando hay muchas restricciones la cadena tiende a colapsar en un |
| 61 | punto de la frontera, lo cual impide que haya una convergencia real en un intervalo de tiempo |
| 62 | razonable. |
48 | 63 | |
49 | 64 | == Diseños de generadores de candidatos factibles == |