11 votos

¿Cuáles son la media y la varianza de la distribución Gamma?

Existen dos formas para la distribución Gamma, cada una con diferentes definiciones para los parámetros de forma y escala. En lugar de preguntar cuál es la forma que se utiliza para la gsl_ran_gamma probablemente sea más fácil pedir las definiciones asociadas para la media y la desviación estándar en términos de los parámetros de forma y escala.

Se agradecerá cualquier indicación sobre las definiciones.

2 votos

La documentación ambigua (o la ausencia de ella) es una señal de alarma, porque sugiere que los implementadores son demasiado inexpertos como para ser conscientes de que hay diferentes convenciones y que la suya necesita una documentación detallada. Por lo tanto, además de averiguar qué convención se utiliza, sería conveniente realizar pruebas exhaustivas de la implementación.

1 votos

@whuber: la documentación proporciona la forma del pdf explícitamente. En este caso, la forma dada es la misma que se utiliza en, por ejemplo, Wikipedia.

0 votos

@whuber: La documentación de GSL es clara y sin ambigüedades; error del usuario. La versión HTML que recibí era lo suficientemente clara como para no sospechar que hubiera un problema.

16voto

Niall Puntos 51

Si el parámetro de forma es $k>0$ y la escala es $\theta>0$ una parametrización tiene la función de densidad

$$p(x) = x^{k-1} \frac{ e^{-x/\theta} }{\theta^{k} \Gamma(k)}$$

donde el argumento, $x$ es no negativo. Una variable aleatoria con esta densidad tiene la media $k \theta$ y la varianza $k \theta^{2}$ (esta parametrización es la que se utiliza en la página de la wikipedia sobre la distribución gamma).

Una parametrización alternativa utiliza $\vartheta = 1/\theta$ como parámetro de velocidad (parámetro de escala inverso) y tiene la densidad

$$p(x) = x^{k-1} \frac{ \vartheta^{k} e^{-x \vartheta} }{\Gamma(k)}$$

Con esta elección, la media es $k/\vartheta$ y la varianza es $k/\vartheta^{2}$ .

0 votos

Gracias por la rápida respuesta; ¿alguna idea de qué forma se utiliza en gsl_ran_gamma?

0 votos

No, pero podrías simular algunas para valores dados de la forma y la escala y ver si la media muestral se acerca más a $k \theta$ o $k/\vartheta$ .

1 votos

@Aengus: Según la documentación, GSL utiliza la parametrización con media $k\theta$ . En la notación de la documentación la media sería $ab$ .

3voto

Tjad Clark Puntos 1

En realidad, además de lo que dijo Macro, hay una tercera forma para la distribución gamma Con un parámetro de forma $v$ y un parámetro medio $\mu$

$ p(x\mid \mu,v) = constant \times x^{\frac{v-2}{2}} e^{-\frac{xv}{2\mu}} $

si $x \sim G(\mu,v)$ entonces $ E(x) = \mu$ y $var(x) =\dfrac{2\mu^2}{v}$

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