Фільтр частинок на прикладі R - тривіального коду


9

Я шукаю простий приклад коду, як запустити фільтр частинок в Р. Здається, що пакет пам’яті підтримує біт математики простору стану, але приклади трохи складні, щоб програматично наслідувати простого розробника OO, такого як я, зокрема як завантажити спостережувані дані в помпезний об’єкт.

Скажімо, у мене є вхідний файл CSV з 1 стовпцем галасливих даних, і я хотів би запустити його через фільтр частинок для того, щоб сподіватися очистити його, з вихідним підрахунком, до іншого файлу csv.

 y <- read.csv("C:/Dev/VeryCleverStatArb/inputData.csv", header=FALSE)
 #CSV to Pomp object ???
 #Run Particle Filter
 #Write estimates to csv.

Основна складність прикладів - завантаження даних csv в помповий об'єкт.

На сьогодні дуже проста модель простору стану повинна бути достатньо хорошою.

Якісь ідеї для R-цікавих?


Це може бути корисним для тих, хто шукає торгівлю парами або алгоритмічну торгівлю взагалі, де це економічні відносини, оточені шумом.

ІМХО, вам краще самостійно кодувати фільтр ...
Д-р G

Дивовижно! Зразки / поради / покажчики для всіх, хто переглядає це? Альтернативне рішення краще, ніж рішення.
користувач1180428

@ user1180428: я відредагував свою відповідь, яка тепер може стати можливою альтернативою для вас.
Вейн

Відповіді:


1

EDIT: Схоже, зараз більшість пакетів фільтрів для частинок зникла. Однак я грав у LaplacesDemon (байєсівський пакет MCMC), і він має функцію PMC (Населення Монте-Карло), яка реалізує PMC, який є типом фільтра для частинок. Можливо, занадто багато техніки для швидкого фільтрування частинок, але речі, які варто вивчити.

Ви можете знайти пакет та підручники у CRAN .

ОРИГІНАЛ: Якщо чесно, у найпростішому випадку pompце важко використовувати. Це дуже гнучко для всього, що ви можете зробити, але це як використовувати космічний корабель, щоб перейти до продуктового магазину.

Ви спробували переглянути фільтри Kalman (якщо ваші дані можуть задовольнити припущення щодо фільтра Калмана), включаючи базові функції tsSmoothта StructTS(лише універсаріат) та пакет dlm? Я також роздивився б loessі інші плавніші.

Я сподіваюся, що я помиляюся, і хтось швидко скакає сюди, "Ось як це зробити для простих універсальних даних, таких як у вас, зі скромними припущеннями". Я хотів би мати можливість використовувати пакет сам.


1
Був там зламаний. На жаль, простий ковзний середній показник здає, що в цьому є корисний сигнал краще, ніж фільтр Калмана, та багато інших прикладів - Kalman: link , SMA: link Дані стаціонарні, настільки, що Dickey Fuller p <0,01. Можливо, я роблю це неправильно. Сценарій для запуску фільтра частинок на цих даних та інших пар, що торгують кандидатами, був би дивним (я думаю).
користувач1180428
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.