Intento ajustar un modelo de serie temporal de población en stan/rstan(2.7.0) donde la tasa de mortalidad depende de la generación anterior (n-1) pero la reproducción depende de una generación desconocida (n-x). No he encontrado una manera de estimar x ya que stan no tiene opciones para distribuciones de enteros y la función floor devuelve reales que no pueden ser utilizados como índices de arrays. ¿Alguien conoce una buena solución para estimar los desfases enteros en stan?
Respuesta
¿Demasiados anuncios?
La estimación de los parámetros enteros en STAN es imposible. HMC depende de los cálculos de gradiente de la densidad posterior. Los parámetros enteros no son adecuados para esos cálculos. Redondear o truncar un número real a un entero no ayudará porque la superficie posterior será plana en una región a lo largo de ese eje, por lo que no habrá información para el muestreador sobre dónde ir después.
Sin embargo, existen varias alternativas.
- Especifique el número de rezagos como datos para un rango de rezagos plausibles y compare las probabilidades posteriores de estos modelos. Tenga en cuenta que tendrá que utilizar las funciones que calculan la probabilidad posterior completa, en lugar de la probabilidad posterior hasta una constante de proporcionalidad. Consulte el manual de STAN para obtener más detalles.
- Utilice una rutina de muestreo estándar de Metropolis/Gibbs.
