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:
-
Nov 4, 2010, 10:56:57 AM (14 years ago)
- Author:
-
Víctor de Buen Remiro
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v9
|
v10
|
|
99 | 99 | [[LatexEquation( \alpha = \min \left\{ 1, \frac{ \pi\left( y \right) }{ \pi\left( x \right)} \right\} $$)]] |
100 | 100 | |
| 101 | Es decir, si el candidato es más probable que el punto anterior siempre es aceptado. |
| 102 | |
101 | 103 | En este caso se tiene una fórmula que da el óptimo ratio de aceptación propuesta en el artículo [http://projecteuclid.org/DPubS/Repository/1.0/Disseminate?view=body&id=pdf_1&handle=euclid.ss/1015346320 Optimal scaling for various Metropolis-Hastings algorithms (''Gareth O. Roberts and Jeffrey S. Rosenthal'')] |
102 | 104 | |
103 | | [[LatexEquation( \alpha^{MH}_n = 0.234 + 0.266*e^{1-n} $$)]] |
| 105 | [[LatexEquation( \alpha^{MH}_n = 0.234 + 0.266 e^{1-n} $$)]] |
104 | 106 | |
105 | 107 | |
… |
… |
|
123 | 125 | |
124 | 126 | ==== {{{Class @MulTryMet.RandWalk}}} ==== |
125 | | La clase {{{@MulTryMet.RandWalk : @MulTryMet, @RandWalk}}} implementa la modalidad más simple de [http://en.wikipedia.org/wiki/Multiple-try_Metropolis Multiple Try Metropolis]: el método Random Walk Multiple Try Metropolis (RWMTM). Por tratarse de una función de transición simétrica se puede tomar simplemente [[LatexEquation( \lambda \left( x,y \right) = 1/Q\left( x,y \right) $$)]] resultando que la función de pesos es la propia densidad objetivo: |
| 127 | La clase {{{@MulTryMet.RandWalk : @MulTryMet, @RandWalk}}} implementa la modalidad más simple de [http://en.wikipedia.org/wiki/Multiple-try_Metropolis Multiple Try Metropolis]: el método Random Walk Multiple Try Metropolis (RWMTM). Por tratarse de una función de transición simétrica se puede tomar simplemente |
126 | 128 | |
127 | | [[LatexEquation( w\left( x,y \right) = \pi\left( y \right) $$)]] |
| 129 | [[LatexEquation( \lambda \left( x,y \right) = \frac{1}{Q\left( x,y \right)} $$)]] |
| 130 | |
| 131 | resultando que la función de pesos es la propia densidad objetivo: |
| 132 | |
| 133 | [[LatexEquation( w\left( x,y \right) = \pi\left( x \right) $$)]] |
128 | 134 | |
129 | 135 | Hasta el momento no se ha podido encontrar una fórmula del ratio de aceptación óptimo análoga a la de Roberts-Rosenthal para el método Metrópolis-Hastings, o al menos el que escribe no tiene noticias de su existencia, por lo que se ha recurrido a una fórmula que heurísticamente da buenos resultados |