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.

Changes between Version 16 and Version 17 of OfficialTolArchiveNetworkBysSamplerRandWalkInPolytope


Ignore:
Timestamp:
Feb 1, 2011, 6:25:18 PM (14 years ago)
Author:
Víctor de Buen Remiro
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • OfficialTolArchiveNetworkBysSamplerRandWalkInPolytope

    v16 v17  
    33
    44== Descripción del problema ==
    5 Si tenemos que generar muestras que cumplan una serie de restricciones de desigualdad lineal
     5Tenemos que generar una cadena de Markov de muestras que se distribuyan con cierta función de la
     6que se conoce su log-densidad salvo una cosntante, a la par que cumplen un conjunto de restricciones
     7de desigualdad lineal
    68
    79[[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} )]]
    810
    9 caben dos posibilidades:
    10  1. utilizar un generador de candidatos que cumpla las restricciones por construcción
    11  1. usar uno libre y luego rechazar los candidatos no factibles.
     11Lo único que deben cumplir las restricciones es que el politopo generado por ellas, o sea,
     12el conjunto de puntos que las cumplen, tenga medida no nula en [[LatexEquation( \mathbb{R}^{n} )]],
     13para que tengan sentido los algoritmos que se van a plantear. La matriz [[LatexEquation(A)]] puede
     14ser singular o tener cualquier forma con tal cumpla esa condición. En particular puede haber
     15restricciones no activas, es decir, que se pueden quitar sin que cambie en absoluto la forma del
     16politopo generado por las restricciones.
     17
     18Nos caben dos posibilidades a la hora de diseñar el generador de candidatos:
     19 1. utilizar un generador de candidatos que cumpla las restricciones por construcción.
     20 1. usar un generador no restringido y luego rechazar los candidatos no factibles.
    1221 
    13 La ventaja de la segunda es que puede ser simétrico y evita el cálculo de la verosimilitud del
     22La ventaja de la segunda es que el generador puede ser simétrico y evitarse el cálculo de la verosimilitud del
    1423candidato, pero el problema es que puede ser que tarde mucho en encontrar uno factible si el punto
    1524actual está demasiado cerca de la frontera, lo cual será muy habitual si el punto de máxima
    1625verosimilitud se encuentra fuera del politopo definido por las anteriores inecuaciones.
    1726
    18 En la siguiente figura se observa un caso en el que el punto máximo verosímil sin restricciones
     27En la siguiente figura se observa un caso en el que el punto máximo-verosímil sin restricciones
    1928(en verde) se encuentra fuera del politopo (en gris), por lo que el punto máximo-verosimil
    2029restringido (en rojo) se encuentra en la frontera del politopo. Las elipses verdes son líneas
     
    150159aún a los vértices donde a veces podría concentrarse la verosimilitud.
    151160
    152 El paseo aleatorio radial asimétrico tiene la ventaja de que puede acercarse mucho más a la frontera
    153 pero resulta mucho más caro de evaluar, pues puede precisar resolver hasta dos problemas de optimización
    154 univariante para cada generación y para cada evaluación de la log-densidad.
     161El paseo aleatorio radial asimétrico tiene la ventaja de que puede acercarse mucho más a la frontera,
     162de hecho puede simular puntos fronterizos, pero resulta mucho más caro de evaluar, pues puede precisar
     163resolver hasta dos problemas de optimización univariante para cada generación y para cada evaluación de
     164la log-densidad.
    155165
    156166Si dos métodos de generar cadenas de Markov cumplen la propiedad del