36 votos

¿Cuáles son las propiedades de un medio de distribución de Cauchy?

Actualmente estoy trabajando en un problema, cuando tengo que desarrollar una cadena de Markov Monte Carlo (MCMC) algoritmo para un modelo de espacio de estado.

Para ser capaz de resolver el problema, se me ha dado la siguiente probabilidad de $\tau$: p($\tau$) = 2I($\tau$>0)/(1+$\tau^2$). $\tau$ siendo la desviación estándar de $x$.

Así que ahora sé que es un medio de Cauchy de distribución, porque he de reconocer que viendo ejemplos y, porque me lo dijeron así. Pero yo no entender completamente por qué es un "Medio de Cauchy" la distribución y las propiedades que vienen con él.

En términos de propiedades no estoy seguro de lo que quiero. Soy bastante nuevo en este tipo de econometría teoría. Así que es más para mí para entender la distribución y cómo la usamos en un espacio de estado del modelo de contexto. El modelo en sí se parece a esto: \begin{align} y_t &= x_t + e_t \\ x_{t+1} &= x_t + a_{t+1} \\[10pt] a_{t+1} &\sim ~ N(0, \tau^2) \\ p(\sigma^2) &\propto 1/\sigma^2 \\[3pt] p(\tau) &= \frac{2I(\tau>0)}{\pi(1+\tau^2)} \end{align}

Edit: he incluido $\pi$ p($\tau$). Gracias por señalarlo.

46voto

AdamSane Puntos 1825

Un medio de Cauchy es una de las mitades simétricas de la distribución de Cauchy (si no se especifica, es la mitad derecha que es la intención):

Plot of Cauchy and half-Cauchy densities

Desde el área de la mitad derecha de una de Cauchy es $\frac12$ la densidad, entonces debe ser duplicado. Por lo tanto el 2 en pdf (aunque le falta un $\frac{1}{\pi}$ whuber señaló en los comentarios).

La media de Cauchy tiene muchas propiedades; algunas son útiles propiedades que desea en un caso previo.

Una opción común para una previa en un parámetro de escala es la inversa de la gamma (no menos, porque es conjugada para algunos familiares de los casos). Cuando un débil informativo antes de la deseada, muy pequeños valores de los parámetros utilizados.

La media de Cauchy es bastante pesado de cola y, también, puede ser considerada como bastante débilmente informativo en algunas situaciones. Gelman ([1], por ejemplo) se aboga por la mitad-t priores (incluyendo el medio de Cauchy) a través de la inversa de la gamma debido a que tienen un mejor comportamiento para los pequeños valores de los parámetros, pero sólo se refiere a ella como wealy informativo cuando a gran escala se utiliza el parámetro*. Gelman se ha centrado más en el medio de Cauchy en años más recientes. El papel por Polson y Scott [2] da otras razones para la elección de la mitad de Cauchy en particular.

* Tu post muestra un estándar de la media de Cauchy. Gelman probablemente no elige que de previo. Si usted no tiene ningún sentido en absoluto de la escala, corresponde decir que la escala es de como podría ser superior a 1 como por debajo de 1 (que puede ser lo que usted quiere), pero no encaja con algunas de las cosas Gelman es argumentado.

[1] A. Gelman (2006),
"Antes de las distribuciones de variación de parámetros en los modelos jerárquicos"
El Análisis Bayesiano, Vol. 1, N. 3, pp 515-533
http://www.stat.columbia.edu/~gelman/investigación/publicado/taumain.pdf

[2] N. G. Polson y J. G. Scott (2012),
"En el Medio de Cauchy Antes de un Mundial Parámetro de Escala"
El Análisis Bayesiano, Vol. 7, Nº 4, pp 887-902
https://projecteuclid.org/euclid.ba/1354024466

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X