自激马尔可夫调制计数的滤波器和平滑器

2022-5-5 05:47:34

用于自激马尔科夫调制计数过程的滤波器和平滑器。塞缪尔·N·科恩*罗伯特·J·埃利奥特2013年11月26日摘要我们考虑一个自激计数过程，其参数依赖于隐有限状态马尔可夫链。基于对跳跃过程的观察，我们推导出隐藏链的最佳滤波器和平滑器。该过滤器为封闭形式，尺寸有限。我们通过模拟数据和在该框架下分析2010年5月6日的“金融危机”，展示了该过滤器的性能。关键词：非线性滤波器，马尔可夫链，霍克斯过程，高频交易，Flash crashMSC 2010:62M05，60G55，60J28，91G701简介在许多情况下，人们希望从嘈杂的观察中推断隐藏状态过程的性质。在本文中，我们关注一个实际场景，其中观察过程是一个计数过程（即，一个整数值跳跃过程，在任何时间最多增加一个），但一个具有“自激”能力的过程，即跳跃的速率可以依赖于过去的跳跃。我们假设过去跳跃和当前速率之间的精确关系由另一个过程调节，这是不可观察的，但是一个有限状态连续时间马尔可夫链。我们的挑战是从我们的观察中发现这个潜在过程的状态。类似的问题在许多环境中都很常见。近年来，自激过程越来越受到重视，因为模型涉及金融市场的高频交易。这些模型可以在[13,17,1,6,9,7]中找到，并且[2]提供了一个调查。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:47:37

自激过程的其他应用示例包括地震发生[22,24]、神经元环路[8]和犯罪活动[18]。另一方面，隐马尔可夫模型已被用于模拟各种各样的经济、金融和工业现象。这方面的例子构成了本书的重点[11]。对于观测过程是伊藤扩散的情况，使用Zakai的参考概率，这些过程的过滤器在[10]中得到。与我们的情况更接近的是[12]，当观测遵循泊松过程，且速率取决于基础条件时*塞缪尔·科恩（Samuel Cohen）感谢牛津-曼恩定量金融研究所（Oxford–Man Institute for Quantitative Finance）的研究支持，特别是提供了纽约-泛欧交易所（NYSE-Euronext）的数据。状态这些方法的一个关键优势是，由于假设基础过程是一个有限状态马尔可夫链，可以以封闭形式获得基础状态的估计，直至确定性ODE的解。这提供了显著的计算优势，因为这些过滤方程的计算速度非常快。将过滤技术与高频金融数据相结合是一个公认的原则。在[3,15,23]和许多其他文献中可以找到使用各种模型的例子。特别是Frey和Runggaldier[15]的情况与我们的情况非常相似，具有自激马尔可夫调制跳变过程，但他们的滤波器不同（他们提到但没有探讨X是有限状态马尔可夫链时的情况），并且他们不考虑将这些方法应用于任何特定模型或数据集。本文的结构如下。在第2节中，我们介绍了正在考虑的一般模型，以及与霍克斯过程相关的具体案例。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:47:41

在第3节中，扩展了[12]的马尔可夫调制泊松观测情况，我们推导出了滤波和更平滑的方程，以及简单的数值近似。我们还导出了过滤器的一个健壮版本，并指出了为什么在这种情况下，它的效用有限。在第4节中，我们讨论了这种情况下的参数估计，包括说明为什么普遍建议的EM算法无法估计下链的转移矩阵。在第5节中，我们将这些方法应用于模拟数据，在第6节中，我们将这些方法应用于围绕2010年5月6日“火山灰崩塌”的TAQ数据。2自激隐马尔可夫模型我们假设我们有以下情况。设X为N态马氏链，其速率矩阵为。在不失一般性的情况下，我们假设X取RN中标准基向量的值，为了方便起见，我们将其表示为X。然后，正如在[10]中一样，对于RN值的P-鞅，X的表达式为xt=X+Z]0，t]AtXtdt+mtm。设Y是一个计数过程，也就是一个右常数整数值递增过程，跳数最多为1。写D表示Y的路径空间。过滤{Yt}t≥0和{Ft}t≥0，其中yt=σ（Ys；s）≤ t），Ft=σ（Ys，Xs；s≤ t） =Yt∨ σ（Xs；s）≤ t）。设P表示{Yt}t中可预测的σ-代数≥0过滤。我们假设我们有一个概率测度P，使得Y是一个带有P-补偿器的协过程（在{Ft}t中）≥0过滤）λ：X×[0，∞[×D→]0, ∞[根据惯例，这是速率矩阵或速率矩阵的转置。在我们的设置中，元素[At]ij是在时间t从状态j跳到状态i的速率，因此A>是过程的最小生成器。它是B（X）×P可测的，这样s7→ λ（Xs）-, s、 Y（·））是连续的。就是，Yt-Z] 0，t]λ（Xs）-, s、 Y（·））dsa（P，{Ft}t≥0）-鞅。λ（Xs）的值-, s、 Y（·））是Y在时间t的跳跃率。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:47:45

请注意，只有“刚才的状态”Xt-出现在速率函数中，但之前观察到的过程Y通常会影响速率（然而，由于可预测性消耗，只有过去的观察才能产生影响）。备注1。当Y出现在强度中时，很难显示λ是否可积分，并且对于概率为1的所有t，Y的值是否为有限（也就是说，我们有一个适当类别的Y密度）。对于特殊形式的λ（t，Y（·））=φZ]-∞,t] h（t）- s）戴斯[4]给出了适当的稳定性条件。我们不必担心这个技术问题，只需假设我们有一个满足这些特性的模型。下面给出了一个形式条件（假设1）。例1。一个重要的例子是λ（Xt-, t、 Y（·））=hα，Xt-i+hβ，Xt-iZ]0，t[e-hγ，Xt-i（t）-s）其中α、β、γ是具有非负项的已知向量。这是前面提到的霍克斯过程的自然变量，参数由马尔可夫链的当前状态决定。注意，超级相似的情况λ=hα，Xt-i+Z]0，t[hβ，Xsie-hγ，Xsi（t-s） DYS不属于我们的模型类别，因为状态X的过去值会影响当前的跳跃率。这一限制是为了确保获得有限尺寸的过滤器。3.过滤和更平滑的方程我们的挑战是，给定Y的观测值，确定状态X。特别是，我们希望能够对任何函数f:X计算E[f（Xt）| YT]→ Rand任意时间t和t。由于X是基向量的空间，很容易看出任何函数f:X→ R可以写成f（Xt）=hf，对于某些向量，可以写成Xti∈ 所以E[f（Xt）| YT]=hf，E[Xt | YT]i。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

kedemingshi

2022-5-5 05:47:47

由于这个原因，我们只要确定E[Xt | YT]就足够了，对于所有时间t，t。对于T<T，因为X是一个马尔可夫链，在时间齐次的情况下，我们有E[Xt | YT]=E[E[Xt | YT]∨ XT]| YT]=E[eA（t）-T）XT | YT]=eA（T-T）E[XT | YT]因此，有必要考虑案例T≥ t、如果A不是时间均匀的，则该方程将略有不同（将使用适当的转移矩阵代替eA（t-T），但同样的简化也是可能的。3.1过滤方程（T=T）我们首先寻求确定E[Xt | Yt]的方程，我们的方法与[12]类似。写λu=λ（u，Xu）-Y（·））表示简单。设Q为测量值，其中Y是标准泊松过程，与X无关，且X具有如上所述的动力学。设Z为方程zt=1+Z]0，t]Zu的解-（λu）- 1）（dYu）- du）=exp-Z] 0，t]（λu）- 1）杜于∈]0，t]λ尤.注意，作为Yt-定义了一个Q-鞅，我们看到Z是一个Z=1的非负Q-局部鞅。假设1。Z是[0，T]上的真鞅。在这个条件下，哪个条件是真的Z] 0，T]（λu）- 1）杜我∞.引理1。如果Z是真鞅，Q是度量，Y是标准泊松过程，X是带有速率矩阵Atat timet的马尔可夫链，那么在度量P定义的bydPdQ下Ft=Zt，Y有补偿器r]0，t]λudu，X和以前一样是带有速率矩阵的马尔可夫链。因此，P=P.Proof。从伊藤的公式中我们可以看到d（ZtYt）Zt-= Yt（λt）- 1）（戴- dt）+λudyt和sod（ZtYt）=Zt-λtdt+（Q-鞅增量），这意味着Y根据需要，在P下具有补偿器r]0，t]λudu。类似地ztxt=X+Z]0，t]Zu-徐-（λu）-1）（dYu）-du）+Z]0，t]Zu-AuXudu+Z]0，t]Zu-dMu。因此X有{Ft}t≥P下的0-compensatorR]0，t]AuXudu，以及具有速率矩阵Au的马尔科夫链。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

2022-5-5 05:47:51

由于我们的空间是由Y和X的路径定义的，这意味着P=P。我们希望确定[Xt | Yt]=EQ[ZtXt | Yt]EQ[Zt | Yt]的递归关系。使用Zakai参考概率方法，我们将通过找到非标准化密度qt=EQ[ZtXt | Yt]，作为EQ[Zt | Yt]=hqt，1i来实现这一点∈ R、这将允许我们计算我们的期望。定理1。非标准化的“过滤”密度过程满足{Yt}t≥0自适应随机向量ODEqt=q+Z]0，t]（λu- 一）曲-（dYu）- du）+Z]0，t]Auqudu。证据我们知道ztxt=X+Z]0，zi-徐-（λu）-1）（dYu）-du）+Z]0，t]Zu-AuXudu+Z]0，t]Zu-dMu。然后，由于Y在Q下有独立的增量，取条件期望，qt=Q+Z]0，t]EQ[Zu-徐-（λu）- 1） |于|- du）+Z]0，t]Auqudu。Asλ是状态的标量函数-, 我们可以写-λu=λ（u，Xu）-, 徐燕英-= ∧uXu-其中∧是对角线矩阵，其条目（λu）ii=λ（u，ei，Y（·）），因此，至关重要的是，∧uis是可测量的。在我们的方程式中，将其替换为q，我们发现qt=q+Z]0，t]（λu）- 一）曲-（dYu）- du）+Z]0，t]Auqudu。由于Y跳跃的离散性，这个随机系统可以用经典的ODE技术很好地逼近。3.2更平滑的方程（T>T）我们现在希望找到一个rt=EQ[ZTXt | YT]的公式∝ E[Xt | YT]，当T>T时。定理2。每一个t≥ 0存在一个YT可测量向量vt，使得“平滑”密度rtsatiseq[ZTXt | YT]=rt=diag（vt）qt。此外，作为t的函数，v满足方程dvt=-（A>t）- ∧t+I）vt-dt+（λ）-1t- 一） vt-dYt；vT=1。证据首先请注意EQ[ZTXt | YT]=EQ[ZtXtEQ[ZT/ZT | YT∨ [Ft]| YT]。因为X是马尔可夫链，Zt/Zt=exp-Z] t，t]（λu）- 1）杜于∈]t、 t]λ尤不依赖于xs，对于s<t，我们有eq[ZT/ZT | YT∨ Ft]=EQ[ZT/ZT | YT∨ σ（Xt）]=hvt，xtif对于某些YT可测随机变量族{vt}t≤T

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:47:55

因此，我们可以写我们可以写我们可以写我们可以写的EQ[EQ[ZT[ZT[ZT[ZtXtX>TTTTTTTTTTTTTTTTTTTTT[TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT[TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT[TTTTTTTTTTTTTTTT[TTTTTTTTTTTTTTTTTTTTTTTTTT[TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTqti=dhvt，qti=hdvt，qti+hvt，dqti+h及物动词，qti=hHtvtdt+KtvtdYt，qti+hvt，Atqtdt+（λt- 一） qt（dYt）- dt）i+hKtvt，（λt）- 一） qtidYt=Dvt，H> t+At- （λt）- （一）qtEdtDvt，K> t+（λt）- 一） +K>t（λt）- （一）qtEdYt。因此，我们看到vt必须具有dynamicsvt=-（A>- ∧t+I）vt-dt+（λ）-1t- 一） vt-戴特。备注2。可以看到，在Y跳跃的时间τ处，我们有vτ=vτ-+ (Λ-1τ- 一） vτ-= Λ-1τvτ-所以vτ-= ∧τvτ。备注3。在滤波和平滑两种情况下，我们都推导出了非标准化密度q和r的方程。然而，实际上，通常很方便“在平面上”归一化，以增加数值稳定性。由于我们所有的方程在q和v中都是线性的（r只是q和v的乘积），通过标量过程进行乘法归一化不会引起任何困难，并且结果仍然正确，直到乘以常数。3.3具体示例：霍克斯过程我们现在将考虑一个具体示例，我们可以推导出一个明确的过滤器。在这种情况下，我们也将推导出更高阶的近似值。设λ由λ（Xt）给出-, t、 Y（·））=hα，Xt-i+hβ，Xt-iZ[0，t[e]-hγ，Xt-i（t）-s）对于已知的非负向量α、β、γ。我们考虑两个关键设置，一个是在连续时间内观察到脉冲，另一个是在离散时间网格上观察到计数数据。3.3.1连续观测当在连续时间内观测到跳跃时，原则上可以完全实施滤波。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

大多数88

2022-5-5 05:47:59

这涉及到求解跳跃之间的以下等式系统：d∧t=-diag（γ）（λt）- diag（α））dtdqt=（A- （λt）- 一））qtdtdvt=-（A>- （λt）- 一））vtdtAt跳跃：∧t+=∧t+diag（β）qt=∧t-qt-及物动词-= λt-vT边界值∧=diag（α），q=X，vT=1（注意∧是左连续的，而q和v是cadlag。）解决这个系统的唯一困难是我们有q和v的非自治ODE，并且需要soa数值积分技术。例如，使用指数积分器方法可以给出Y在ti处跳跃之间的解-1和ti，∧ti=diag（α）+exp(-diag（γ）（ti）- 钛-1））（λti-1+diag（β- α））qti≈ ∧tiexp（A）- （λti）-1.- 一））（ti- 钛-1)qti-1vti-1.≈ 经验（A>- （λti）-1.- 一））（ti- 钛-1)根据问题的参数，在跳跃之间加入更多的点，并在常微分方程的解中使用这些点也将提高性能。3.3.2离散观测当观测仅离散观测时，即我们有固定的时间网格{ti}并观测Yti值，则另一种方法是合适的。首先注意∧ti=diag（α）+exp(-diag（γ）（ti）- 钛-1））（λti-1.- diag（α））+diag（β）Z[ti-1，ti[exp(-diag（γ）（s）- 钛-1））未提供的γ（ti）-钛-1）不是很大，Y近似为泊松过程，在区间[ti]上具有恒定速率-1，ti]，所以如果我们观察间隔[ti]内的跳跃-1，ti]，那么这些跳跃在区间内大致均匀分布。因此∧ti≈ diag（α）+exp(-diag（γ）（ti- 钛-1））（λti-1.- diag（α））+diagβ1 - 经验(-γ（ti）- 钛-1））γ（ti- 钛-1)δYti（1），其中第二行括号内的操作是可读的。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:48:02

（采用这种近似可以略微减少偏差vsa naive方案，其中最后一项仅为+diag（β）δYti。）由于q和v的方程是指数形式，因此QTI给出了有效近似值≈ 经验（A）- （λti）-1.- 一））（ti- 钛-1） +δYtilog（λti）-1)· qti-1vti-1.≈ 经验（A>- （λti）-1.- 一））（ti- 钛-1） +δYtilog（λti）-1)· vti。（2）结合这些方程，可以近似计算连续时间滤波器。3.4稳健滤波器为了完整性，我们现在概述了一个“稳健”滤波器，在这里我们可以避免对滤波方程中的Y进行积分。定理3。定义矩阵值过程Γt=exp-Z] 0，t]（λs）- 一） ds·Y∈]0，t]（λs）是的。然后“qt:=Γ”-1Tqt和“vt=Γtvtsatisfy”qt=q+Z]0，t]-1 sΓsΓqsdsΓvt=\'v-Z] 0，t]ΓsA>sΓ-1s“vsds=ΓT1+Z]t，t]ΓsA>sΓ-注意，这些方程不直接涉及观测过程Y（asit仅通过Γ出现）。证据首先要注意的是，Γ和∧都是对角矩阵，因此是可交换的。从定义上，我们可以看出Γ-有活力吗-1t=（λt）- 一） Γ-1tdt+（λ）-1t- 一） Γ-1t-戴斯。现在根据伊藤法则，d（Γ）-1tqt）=Γ-1t-λt- 一） qt-（戴- dt）+Aqtdt+（λt）- 一） Γ-1tdt+（λ）-1t- 1)Γ-1t-戴斯qt-+ (Λ-1t- 1)Γ-1t-（λt）- 一） qt-dYs=Γ-1tAqtdtso，写“qt=Γ”-1tqt，我们有\'qt=\'q+Z]0，t]Γ-1sAΓsΓqsds。也可以通过写入qt，vti=h\'qt，\'vti=hΓ导出相应的平滑器-1tqt，从中我们可以看到。应用伊藤法则，这满足了“vt=”v-Z] 0，t]ΓsA>-当vT=1时，vT=T1+Z]t，t]sA>-1s’VSD。备注4。从数值上看，这些方程具有显著的影响。当Γ与指数有关时，除非∧sis接近I，否则通常至少会有Γ和Γ中的一个-1快速成长。由于不可交换性，这意味着Γ的某些成分-1沙Γ会很大。随着“q”像exp（Γ）一样增长-1sAΓs），这会导致双指数，迅速变得数值不稳定。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

kedemingshi

2022-5-5 05:48:05

这限制了这些稳健过滤器的实际适用性。4参数估计为了实现这些过滤器，通常需要使用训练数据集估计过程的参数。我们将我们的估计分为两部分——确定λ的参数的估计和链的传递矩阵A的估计。这种区别很有用，因为相对于基础链中的跳数，我们通常会在Y中观察到大量跳数——事实上，正是这种“多尺度”行为允许过滤器工作。由于A的估计取决于X的跳跃次数，这意味着我们通常能够很好地估计λ的参数，但是我们对A的估计需要大量的数据（相应地，成本也要高得多）。4.1观测参数我们继续寻找给定状态路径X的最大似然估计器，然后我们可以将其与EM算法相结合。下面的结果是标准引理2（例如参见[7]）。给定X的对数似然函数和区间[0，T]上的观测值，加上常数Xu：Yu6=0logλ（u，Xu，Y（·））-ZTλ（u，Xu，Y（·））然后，假设概率^ru=E[Xu | Yu]，这就得到了部分信息似然函数∝ 鲁。我们用λ（u，ei，Y（·））为向量写λuf。引理3。在区间[0，T]上观察到的对数似然函数，加上常数Xu：Yu6=0h^ru，对数（λu）i-ZTh^ru，λuidut根据λ的性质，可以对该函数进行数值优化或解析优化。然后，EM算法允许我们在给定λr的参数（可以使用数值优化技术实现）和给定λ参数（可以使用滤波方程实现）之间进行迭代。备注5。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:48:09

在某些情况下，可能很难为λ的参数指定良好的初始值，然而，绘制观测数据可能会通过检查揭示自然近邻聚类。使用这种特殊聚类作为r的初始值，我们可以首先计算λ的参数，然后通过过滤器实现这些参数。下面使用这种技术。备注6。对于实际数据，这种方法可能会有问题，因为它对理论模型的偏差非常敏感。确定能够产生稳定、稳健方法的修改是一个持续研究的领域。考虑到EM算法中的反馈效应，这个问题可能非常重要，并且很容易出现在实际数据中。4.1.1示例：离散观测霍克斯过程在离散时间内，这种对数似然的各种近似是可能的。给定（1）和（2）中λ和r的先前计算值，xid^rti给出了一个简单的近似值，δYtilog（λti）-1)- λti-1（ti）- 钛-1）可以快速计算。然后，可以使用各种数值方法使该函数最大化。4.2基础链动态为了估计基础链转移矩阵，通常成本非常高，因为与观测数量相比，链的跳跃很少。这意味着，即使对基础链进行了完美的观察，过渡矩阵的估计也会很差。一种建议的方法（例如参见[10]）是使用EM算法，其中隐藏变量是状态之间的转换次数以及每个状态的占用时间。然而，这种方法在这种情况下会退化，如下所示。设J是带有条目Jij的随机矩阵，其中jiji是观察期内从状态i到状态J的转换次数，Jii=-Pjjij对于所有i.让K是每个状态的占用时间的向量。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:48:11

然后，从X的结构，如^rt=E[Xt | YT]我们知道E[K | YT]=EhZ]0，T]Xt-dtYTi=Z]0，T]rtdtE[J|YT]=EhZ]0，T]Xt-dX>tYTi=EhZ]0，T]Xt-十> t型-A> dt+Z]0，T]Xt-dM>tYTi=Z]0，T]E[Xt-十> t-|YT]A>dt=diagZ] 0，T]rtdtA> 。因此，给出A的初始估计值，我们有[J | YT]≈ 诊断Z] 0，T]^rtdt^A>=diag（E[K|YT]）^A>。给定跃迁次数（J）和占用时间（K）的最大似然估计为J>diag（K）-因此EM算法给出了估计值[J | YT]>diag（E[K | YT]）-1=^Adiag（E[K | YT]）diag（E[K | YT]）-1=^A因此，从一些初始估计开始，迭代EMalgorithm，我们得到估计序列^A=^A=。。。由此可知，以这种方式应用的EM算法确实会产生速率矩阵的一致性估计。我们在对真实数据建模时采用的另一种观点是，将模型中的速率矩阵视为过滤器的“调整参数”，应在校准期间使用专家判断来确定。从这个角度来看，我们可以选择一个参数化的速率矩阵族，例如nA=-1 11 -1.; > 0o然后在校准期间选择这使得过滤器在确定基本状态时具有可接受的性能。当如果规模较大，则允许频繁切换状态，而如果速度较慢，则在过滤器检测到状态变化之前需要strongerdata。然后，应通过考虑样本外性能来确定该参数选择的有效性。5模拟数值结果为了评估这些方法的准确性，方程式（1）和（2）给出的算法在R中实现。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

kedemingshi

2022-5-5 05:48:15

首先模拟马尔科夫链（在连续时间内），然后使用绪方[21]的方法和分支方法（如[19]的算法1所建议，忽略边缘效应）的组合来模拟马尔科夫链跳跃之间的霍克斯过程，从而模拟数据。选择参数值以给出以下行为。我们有两个状态链，在模拟的时间尺度上有一些跳跃。在FirstState中，Y过程具有明显的自激性，但基频率较低，即α较小，β和γ较大。在第二种状态下，该过程是可忽略的自激过程，但具有较高的基速率，即α为中等大，β为小，γ为中等小。然后选择这些参数，这样，如果基础状态保持不变，那么每个状态下的长期平均跳变率是相同的。对于模拟视界ofT=1000，选择的值为过渡矩阵=-0.01 0.010.01 -0.01观测参数αβγ状态1 10/7状态2 18 0.01 0.1在没有状态变化的情况下，在长期内，每个状态的平均跳跃率为20=α1-单位时间的β/γ。记录每次间隔长度为0.1的跳跃次数。下面显示了一个典型的示例路径（图1）。图1：马尔可夫调制霍克斯过程的模拟观察路径对于上述样本路径，粗略的初始聚类是指在大约50、200、250、300、400、650和950的时间段存在状态变化，这些变化可以通过肉眼进行。然后，可以使用最大似然估计器，使用该初始分配来估计参数。转移矩阵被认为是已知的。使用optim命令中的内置Nelder–Mead方法对可能性进行优化。根据这些估计，可以应用过滤器来确定隐藏状态，从而重新校准和重复。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:48:18

参数估计的最初几次迭代如下所示。迭代αβγ1 7.35 20.23 0.813 8.6×10-081.34 0.2072 7.54 19.85 0.869 0.0032 1.485 0.1783 7.64 19.59 0.891 0.0050 1.538 0.1644 7.60 19.45 0.893 0.0064 1.534 0.169这些估计值有一些误差，但是下面给出了使用这些最终参数的过滤和平滑路径（图2），以及用于比较的真实（未知）状态。从这些路径中，我们可以看到滤波器在这些噪声参数值下具有良好的性能。我们注意到有一些状态变化（在t=700和t=800之间）是平滑器无法检测到的。这并不奇怪，因为这些变化只是一个短暂的持续时间，所以越平稳的人就越容易确定跳跃率的任何明显变化更有可能是由于随机变化。0 200 400 600 800 10000.0 0.4 0.8时间过滤状态0 200 400 600 800 10000.0 0.4 0.8时间平滑状态0 200 400 600 800 10000.0 0 0.4 0.8时间真实状态图2：图1中观察路径的过滤、平滑和真实隐藏状态（尤其是Xt=e的概率）。我们发现，对于这些参数值，过滤器在检测从状态2到状态1的变化时表现良好（变化非常剧烈），但在检测从状态1到状态2的变化时效果较差。图3显示了每个ei的λ（ei，t，Y（·））值。当nx=e时（图中的情况），每个状态下的电位率差异非常明显。跳跃的低基本速率意味着λ（e，t，Y（·））通常较低，偶尔出现大的尖峰。相反，当X=e时，跳跃的速度要稳定得多。因此λ（e，t，Y（·））和λ（e，t，Y（·））相对接近，并且λ（e，t，Y（·））的挥发性比nx=e时小。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

何人来此

2022-5-5 05:48:21

由于在确定滤波器时，速率的相对差异非常重要，因此我们发现，从状态2到状态1的真实状态变化会导致λ的可能值变得更加明显，因此滤波器会很快检测到这种变化。图3:timest在图1中观察路径的每个状态下的速率∈ [150, 300]. 函数λ（e，t，Y（·））对应于X为1时的速率，以蓝色显示，而λ（e，t，Y（·））为绿色，λ（Xt，t，Y（·））为黑色。6现实世界的应用我们现在概述这些方法对现实世界数据的可能应用。考虑到为现实世界建立良好模型的复杂性，本节应被视为这些方法实用性的经验支持，而不是建模的实际处方。特别是，我们将使用现实但特别的参数选择，以证明过滤器的有效性。我们将在上午9:30后给出时间，单位为“t=秒”。2010年5月6日，市场经历了众所周知的“泡沫崩盘”，道琼斯指数在5分钟内下跌600点（5-6%），然后迅速恢复。其中一个主要因素是在芝加哥商品交易所Globex平台上的E-mini S&P期货合约上进行交易。关于这一事件的报告[20]详细描述了市场行为，进一步的研究可以在[16]中找到。在一天的波动交易之后，一名基本面交易员在下午2:32（t=18120）对电子迷你合约发出了一份慷慨的卖出指令，该指令在20分钟内执行，不考虑价格或时间。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:48:24

这种卖出指令很可能主要被高频交易员、期货市场的基本面买家和跨市场套利者（尤其是通过与标准普尔500指数中的基本面证券进行配对交易）吸收。随着交易员调整仓位，在下午2:41至2:44的三分钟内（t=18660至18840），E-mini和SPY的价格因综合抛售压力下降了3%。今年秋天产生了一个烫手山芋效应，导致进一步下跌SPDR S&P500 ETF（SPY）是一只交易所交易基金，复制了在纽约证券交易所交易的S&P500。直到下午2:45:28（t=18928），迷你合约暂停五秒钟的交易。从下午2点41分到交易暂停，E-mini的价格下跌了5%，而间谍的价格下跌了6%。这些合约随后以波动的方式恢复，直到下午3点左右（t=19800），价格稳定在接近原始水平。这一事件导致标准普尔500指数波动指数（VIX）上涨31.7%。我们将考虑纽约-泛欧交易所TAQ数据作为间谍合同。我们将尝试使用我们的过滤器，仅使用SPY上交易的频率，自动识别金融崩溃。我们将考虑的数据源只包括带有时间戳的间谍交易，时间仅为1秒精度。我们将把注意力集中在坠机当天，以及相关数据上（而不是引用）。我们将考虑交易数量，而不是交易量，因为这将更容易显示高频交易者的自我激励行为。我们在图4中绘制了当天的交易数据（价格和交易数量）。Kirilenko等人[16]将14:32–15:08（t=18120到20280）这段时间确定为坠机时间，如图4中的红色所示，图5给出了这段时间的更高分辨率图。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:48:28

从记录的交易价格中，我们可以清楚地看到崩盘的时期，以及今天下午交易量的增加。我们还可以看到，数据似乎包含明显的记录错误，尤其是在危机期间和之后，一些交易的记录价格与大部分市场价格相差甚远（或可能在错误的时间）。这在高频数据[5,14]中很常见，但由于我们的方法只使用交易数量，因此可以合理地避免这些错误，并且我们不进行进一步的数据清理。图4:2010年3月6日纽约-泛欧交易所间谍交易（9:30–16:00）、记录价格（蓝色）和交易数量（黑色）。14:32–15:08期间的交易以红色显示。图5:2010年3月6日纽约-泛欧交易所间谍交易（14:32–15:08），记录价格（蓝色）和交易数量（黑色）。我们将把市场建模为两种状态之一，一种是正常状态X=e，另一种与高交易期X=e有关，我们希望这将反映崩溃。对于粗略的近似值，我们将使用13:43:20到16:00的时间段作为崩盘影响的时间段，而09:30到13:43:20的时间段将用于校准通常的交易条件。（这段时间是上午9:30后的15000秒。）由于自激过程的朴素EMalgorithm对模型偏差非常敏感，因此我们将使用单通道对该分组进行校准，而不是应用EM算法。与泊松分布相比，每秒的交易数量似乎明显过度分散，达到了无法通过在1秒时间尺度上使用自激模型（这是我们模型的最小分辨率）进行校正的程度，这一事实突显了这一点。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2022-5-5 05:48:32

因此，作为临时修正，我们还将简单地将每秒的交易数除以25，因为这将产生一个工作模型。这个重新缩放的过程我们称为Y。我们假设一个马尔可夫调制的霍克斯过程，其类型更为复杂，带有一个额外的多项式项，这使得跳跃的效果更加灵活。也就是说，我们认为Y也近似于一个带速率的纯跳跃过程hα，Xti+hβ，XtiZ[0，t[e]-hγ，Xti（t-s）戴斯hζ，Xti。在1秒尺度上离散积分，并使用初始分组，我们得到最大似然估计状态αβγζX=e1。0014241 1.0416288 0.9996515 0.8987939X=e0。5222101 1.7255265 0.5095281 0.6821643底层马尔可夫链的速率矩阵A可以作为我们模型的调整参数。采取一种形式=-1 11 -1.,我们发现 = 10-7产生良好的表现，每115天对应一次自我激励事件。这些参数的选择显然是临时的，是为了证明过滤器的潜在性能。图6和图7显示了过滤和平滑的结果。图6:2010年5月6日间谍交易和价格，以及经过过滤和平滑的州估计。可以看出，在碰撞发生的当天，除碰撞本身外，平滑度越高，检测不到任何状态变化。在过滤器表明可能发生状态变化的地方，观察到了一些异常点，但图7：间谍交易和价格2010年5月6日14:32–15:08，以及经过过滤和平滑的状态估计。所有问题都会在接下来的几秒钟内迅速解决（除了一天结束时的问题）。更仔细地考虑图7，我们可以看到，正如预期的那样，过滤器在更平稳的情况下检测到崩溃，但仍然在价格急剧下跌之前。在价格恢复的中途，过滤器恢复到正常状态。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

大多数88

2022-5-5 05:48:35

从价格下跌开始，越平稳的人成功地检测到崩溃，并恢复到价格复苏开始时的正常状态。回想一下，这种计算是在不观察价格的情况下进行的，因为过滤和平滑完全基于交易数量。7结论我们发现，这些过滤方法具有检测异常交易行为的巨大潜力，允许在正常和异常状态下存在自激。这种方法的一个关键优势是，它在自激结构的选择上非常灵活，并且，由于过滤和平滑的方程以闭合形式可用，因此可以以非常高的速度实现。鉴于自动高频交易的增加，以及由此可能产生的潜在问题，这一点尤为重要。这项工作的未来发展将需要为贸易数据构建更现实、更灵活的模型。特别是，我们对过度分散的特别修正（交易数量除以25）并不令人满意，而开发多级自激模型，可以直接管理数据集的“尖峰”性质，将是更好的选择。可靠统计方法的发展也是一个标志，特别是考虑到THEM算法的使用和高频数据常见的误差。参考文献[1]E.Bacry、S.Delattre、M.Homann和J.F.Muzy。用相互激励的点过程模拟微观结构噪声。http://arxiv.org/pdf/1101.3422.pdf.[2] Luc Bauwens和Nikolaus Hautsch。使用点过程对金融高频数据进行建模。《金融时间序列手册》，第953-979页，2009年。[3] 彼得·J·博兰德和杰罗姆·T·康纳。用于高频金融数据价格估计的约束神经网络卡尔曼滤波器。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

何人来此

2022-5-5 05:48:39

《国际神经系统杂志》，8:399–415，1997年。[4] Pierre Br’emaud和Laurent Massouli’e.非线性Hawkes过程的稳定性。《概率年鉴》，24:1563-15881996。[5] C.T.布朗利和G.M.加洛。超高频金融经济计量分析：数据处理问题。《计算统计与数据分析》，51:2232–22452006。[6] 阿尔瓦罗·卡特亚、塞巴斯蒂安·贾蒙格尔和杰森·里奇。低买高卖：高频交易视角。http://ssrn.com/abstract=1964781.[7] 陈峰和彼得·霍尔。非平稳自激发点过程的推断及其在超高频金融数据建模中的应用。应用概率杂志，即将出版。[8] E.S.乔诺博伊、L.P.施拉姆和A.F.卡尔。神经点过程系统的最大似然辨识。生物控制论，59:265–275，1988。[9] Rama Cont.高频金融数据的统计建模：事实、模型和挑战。IEEE信号处理杂志（金融应用信号处理特刊），28:16–252011。http://ssrn.com/abstract=1748022.[10] 罗伯特·J·埃利奥特。噪声观测马尔可夫链的新有限维滤波器和平滑器。IEEE信息论学报，39:265–271，1993年。[11] 罗伯特·J·埃利奥特、L·阿贡和J·B·摩尔。隐马尔可夫模型：估计与控制。Springer Verlag，柏林海德堡纽约，1994年。[12] 罗伯特·J·埃利奥特和W·P·马尔科姆。马尔科夫调制泊松观测的一般平滑公式。IEEE自动控制交易，50:1123–11342005。[13] 罗伯特·F·恩格尔和杰弗里·R·拉塞尔。自回归条件持续时间：不规则间隔事务数据的新模型。《计量经济学》，66:1127-11621998。[14] 托马斯·尼尔·法肯贝里。高频数据过滤。Tick Data Inc.白皮书。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝