Estoy buscando un ejemplo de código simple de cómo ejecutar un filtro de partículas en R. El paquete pomp parece admitir el bit matemático del espacio de estado, pero los ejemplos son un poco difíciles de seguir programáticamente para un desarrollador OO simple como yo, particularmente cómo cargar los datos observados en un objeto pomp.
Digamos que tengo un archivo csv con 1 columna de datos ruidosos como entrada, y me gustaría ejecutarlo a través de un filtro de partículas para poder limpiarlo, con la salida como las estimaciones, a otro archivo csv.
y <- read.csv("C:/Dev/VeryCleverStatArb/inputData.csv", header=FALSE)
#CSV to Pomp object ???
#Run Particle Filter
#Write estimates to csv.
La principal dificultad con los ejemplos es cargar datos csv en un objeto pomp.
Un modelo de espacio de estado muy simple debería ser lo suficientemente bueno por ahora.
¿Alguna idea para los R-curiosos?