Saturday, January 23, 2010

Simulación de observaciones

Han sido varios los comentarios acerca del método de la grilla descrito en una entrada anterior… Sin embargo, he vuelto a leer lo que escribí y me di cuenta que jamás mencioné cómo simular datos de la distribución horrible que aparecía en esa entrada, cuya forma funcional era más o menos la siguiente:


$latex p(theta)propto sum_{j=1}^p left(frac{1}{sigma^2_j+tau^2}right) prod_{j=1}^p(sigma^2_j+tau^2)^{1/2}expleft( -frac{(y_j-mu)^2}{2sigma^2_j+tau^2}right)$


Lo único que se dijo en ese entonces estaba relacionado con un resumen del método de la grilla (super armamento) y su aplicación a una distribución (munición) exponencial. No sé que me pasó, porque en realidad lo valioso del método es poder simular con distribuciones horribles como la de arriba y no con distribuciones fáciles como la exponencial… En fin, en este documento encontrarán (extracto del nuevo libro) la implementación del método de la grilla en el ambiente computacional de R. Por supuesto el objetivo cumplido es la simulación de datos de esa distribución horrible.

No comments:

Post a Comment