Método de simulación en estadística
En estadística computacional, el método Monte Carlo de cadena de Markov con salto reversible es una extensión de la metodología estándar de Monte Carlo de cadena de Markov (MCMC), introducida por Peter Green , que permite la simulación (la creación de muestras ) de la distribución posterior en espacios de dimensiones variables . [1]
Por lo tanto, la simulación es posible incluso si no se conoce el número de parámetros en el modelo . El "salto" se refiere al cambio de un espacio de parámetros a otro durante la ejecución de la cadena. El método RJMCMC es útil para comparar modelos de diferentes dimensiones para ver cuál se ajusta mejor a los datos. También es útil para las predicciones de nuevos puntos de datos, porque no necesitamos elegir y corregir un modelo, el método RJMCMC puede predecir directamente los nuevos valores para todos los modelos al mismo tiempo. Los modelos que se adaptan mejor a los datos se elegirán con mayor frecuencia que los peores.
Detalles sobre el proceso RJMCMC
Sea un indicador de modelo y el espacio de parámetros cuyo número de dimensiones depende del modelo . La indicación del modelo no necesita ser finita. La distribución estacionaria es la distribución posterior conjunta de que toma los valores .
La propuesta se puede construir con una aplicación de y , donde se extrae de un componente aleatorio con densidad en . El movimiento al estado se puede formular así como
La función
debe ser uno a uno y diferenciable, y tener un soporte distinto de cero:
de modo que existe una función inversa
que es diferenciable. Por lo tanto, y deben ser de igual dimensión, lo cual es el caso si el criterio de dimensión
se cumple donde es la dimensión de . Esto se conoce como coincidencia de dimensión .
Si entonces la condición de coincidencia dimensional se puede reducir a
con
La probabilidad de aceptación estará dada por
donde denota el valor absoluto y es la probabilidad posterior conjunta
donde es la constante normalizadora.
Paquetes de software
Hay una herramienta experimental RJ-MCMC disponible para el paquete BUGs de código abierto .
El sistema de programación probabilística Gen automatiza el cálculo de la probabilidad de aceptación de los núcleos MCMC de salto reversible definidos por el usuario como parte de su función Involution MCMC.
Referencias
- ^ Green, PJ (1995). "Computación de Monte Carlo con cadena de Markov de salto reversible y determinación de modelos bayesianos". Biometrika . 82 (4): 711–732. CiteSeerX 10.1.1.407.8942 . doi :10.1093/biomet/82.4.711. JSTOR 2337340. MR 1380810.