我正在寻找一个如何在R中运行粒子过滤器的简单代码示例。pomp程序包似乎支持状态空间数学位,但是对于像我这样的简单OO开发人员,这些示例很难以编程方式进行遵循如何将观察到的数据加载到pomp对象中。
可以说我有一个带有1列嘈杂数据的csv文件作为输入,并且我希望通过“粒子过滤器”运行它,以便希望将其清理为输出,并将其输出为另一个csv文件。
y <- read.csv("C:/Dev/VeryCleverStatArb/inputData.csv", header=FALSE)
#CSV to Pomp object ???
#Run Particle Filter
#Write estimates to csv.
这些示例的主要困难是将csv数据加载到pomp对象中。
现在,一个非常简单的状态空间模型应该已经足够了。
对R好奇有什么想法吗?
这对于任何从事成对交易或一般算法交易的人都是有用的,在这些交易中,经济关系充满了噪音。
恕我直言,您最好自己对过滤器进行编码...
—
G博士2012年
太棒了!样本/提示/指示器,供所有人查看?替代解决方案总比没有解决方案好。
—
user1180428'2
@ user1180428:我已经编辑了答案,现在可以为您提供一个可能的选择。
—
韦恩