Tengo esta función:
$$\pi(\theta|\vec{x}) = \frac{2^{-n/2-2}\theta^{-n/2-3}\left(\sum_\limits{i = 1}^{n} x_i^2+4\right)^{n/2+2}\exp\left({-{\left(\sum_\limits{i = 1}^{n} x_i^2 + 4\right)}/{2\theta}}\right)}{\Gamma(n/2 +2)} $$
y estoy tratando de tomar la media de la misma con respecto a $\theta$ .
Esto es lo que hice:
$$E(\theta|\vec{x}) = \int_{0}^{\infty} \theta \, \frac{2^{-n/2-2}\theta^{-n/2-3}\left(\sum_\limits{i = 1}^{n} x_i^2+4\right)^{n/2+2}\exp\left({-{\left(\sum_\limits{i = 1}^{n} x_i^2 + 4\right)}/{2\theta}}\right)}{\Gamma(n/2 +2)} d\theta$$ Utilizando wolfram alpha, la respuesta resultó ser $$=\frac{n (\sum_{i = 1}^{n} x_i^2 + 4)\Gamma(n/2)}{4 \Gamma(n/2+2)}$$
Cuando simplifiqué esto a mano, obtuve:
$$ = \frac{n (\sum_{i = 1}^{n} x_i^2 + 4)\Gamma(n/2)}{4 (n/2+2)(n/2+1)\Gamma(n/2)} = \frac{n (\sum_{i = 1}^{n} x_i^2 + 4)}{(n+4)(n+2)}$$
Sin embargo, wolfram demostró que se trataba de una simplificación: $$= \frac{(\sum_{i = 1}^{n} x_i^2 + 4)}{(n+2)}$$ que no es equivalente a la respuesta original de la integral al introducir los números.
Agradecería cualquier sugerencia sobre dónde pueden estar fallando las cosas. La respuesta que obtuve de la integral y la simplificación manual aparentemente no son correctas.
0 votos
Llame a $Y = \sum_{i=1}^nX_i^2 + 4$ e intenta ajustarlo a una distribución conocida para obtener la respuesta correcta (probablemente chi-cuadrado o gamma)