8 votos

¿Por qué no utilizar Beta(1,1) como límite para evitar la prioridad en un parámetro de correlación transformado?

En Análisis bayesiano de datos En el capítulo 13, página 317, segundo párrafo completo, en las aproximaciones modales y distributivas, Gelman et al. escriben

Si el plan es resumir la inferencia por el modo posterior de $\rho$ [el parámetro de correlación en una distribución normal bivariada], nosotros sustituiríamos la distribución a priori U(-1,1) por $p(\rho) \propto (1 - \rho)(1 + \rho)$ , que equivale a una Beta(2,2) sobre el parámetro transformado $\frac{\rho + 1}{2}$ . La densidad previa y la resultante son cero en los límites y, por lo tanto, el modo posterior nunca será -1 o 1. Sin embargo,... la densidad a priori para $\rho$ es lineal cerca de los límites y, por lo tanto, no contradice ninguna probabilidad.

A continuación se muestra un gráfico de la PDF para la distribución Beta(2,2).

PDF plot of Beta(2,2)

Aunque el gráfico se da para el dominio [0,1], la forma es la misma para el dominio [-1,1] obtenido al realizar la inversa de la transformación descrita en la cita anterior. Se trata de una distribución bastante informativa. Da unas siete veces la densidad a $\frac{\rho + 1}{2} = 0.5$ que a $\frac{\rho + 1}{2} = 0.3,0.97$ . Así que de hecho sería contradecir la probabilidad si la probabilidad apuntaba a algo lejos de los límites, pero aún más lejos de $\rho = 0$ . ¿No sería un mejor límite para evitar el previo Beta(1 + $\delta$ ,1 + $\delta$ ), donde $\delta \rightarrow 0$ . Tomemos, por ejemplo, Beta(1,0001, 1,0001), representada a continuación:

PDF plot of Beta(1.0001,1.0001)

El problema con esta prioridad, por supuesto, es que la densidad cae muy bruscamente cerca de cero, lo que puede contradecir la probabilidad de que apunte a un espacio que está muy muy cerca de un límite. Lo que me lleva a mi pregunta:

¿Por qué no establecer la prioridad del parámetro de correlación transformado en Beta(1,1)? Porque la densidad de la distribución beta es cero para $\frac{\rho + 1}{2} = 0,1$ es equivalente a la distribución uniforme sobre el Abrir intervalo (-1,1) en lugar del intervalo cerrado [-1,1], por lo que no es una priorización que evita los límites, y no es preferible a una priorización que cree bastante en la probabilidad de que $\rho = 0$ ¿que sólo es deseable si realmente se tiene esa creencia?

En términos más generales, ¿el uso de la distribución beta no es por definición un límite que evita la prioridad porque su soporte es $0 < \frac{\rho + 1}{2} < 1$ ?

2voto

Christian Hagelid Puntos 121

En los casos en que la función de probabilidad se maximiza en $\rho=-1$ Si se utiliza Beta(1,1) sobre un intervalo abierto, no sirve de nada, ya que entonces no hay un modo posterior bien definido. No hay ningún punto en el intervalo abierto que se pueda decir que es la moda posterior.

El libro quiere que la prioridad sea lineal en los límites, lo que selecciona únicamente Beta(2,2) entre todas las distribuciones Beta. Pero el libro nunca explica con precisión por qué la linealidad es tan importante.

0 votos

¡Y nunca se habla de los efectos secundarios que esbozo en la pregunta! He votado su respuesta, pero aún no la he comprobado como la definitiva. ¡Y me había dado por vencido en esta pregunta!

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