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 9 and Version 10 of OfficialTolArchiveNetworkBysSampler


Ignore:
Timestamp:
Nov 4, 2010, 10:56:57 AM (14 years ago)
Author:
Víctor de Buen Remiro
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • OfficialTolArchiveNetworkBysSampler

    v9 v10  
    9999[[LatexEquation( \alpha = \min \left\{ 1, \frac{ \pi\left( y \right)  }{ \pi\left( x \right)}  \right\} $$)]]
    100100
     101Es decir, si el candidato es más probable que el punto anterior siempre es aceptado.
     102
    101103En 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'')]
    102104
    103 [[LatexEquation( \alpha^{MH}_n = 0.234 + 0.266*e^{1-n} $$)]]
     105[[LatexEquation( \alpha^{MH}_n = 0.234 + 0.266 e^{1-n} $$)]]
    104106
    105107
     
    123125
    124126==== {{{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:
     127La 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
    126128
    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
     131resultando que la función de pesos es la propia densidad objetivo:
     132
     133[[LatexEquation( w\left( x,y \right) = \pi\left( x \right) $$)]]
    128134
    129135Hasta 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