Zarinelli, Treccani, Farmer & Lillo (2014) — Beyond the Square Root

作者: Elia Zarinelli(LIST S.p.A., Trieste)、Michele Treccani(Mediobanca;QUANTLab Pisa)、J. Doyne Farmer(Oxford Martin School 数学研究所;Santa Fe Institute)、Fabrizio Lillo(Scuola Normale Superiore;QUANTLab;Santa Fe Institute)。

发表: arXiv:1412.2152(2014-12-05);正刊 Market Microstructure and Liquidity 1(2), 2015,DOI 10.1142/S2382626615500045。PDF 38 页。

⚠️ 标题更正:本库早期记录误记为「Overcoming the Pro-Rata Curse」,错误。那是同组作者的早期工作论文标题;本文正刊标题是「Beyond the Square Root: Evidence for Logarithmic Dependence of Market Impact on Size and Participation Rate」。backlog 里 7.4 与 7.5 实为同一篇的不同版本。

一句话总结

平方根公式在订单规模的两个数量级上拟合得很好,但换成对数函数能覆盖近五个数量级——而且冲击其实依赖「执行时长」和「参与率」两个变量,平方根的「只依赖规模」是数据聚合造成的假象。

数据

  • Ancerno 数据库的机构 metaorder(母单)执行记录,美股 2007–2009
  • 过滤后 6,944,883 笔 metaorder;冲击曲面分析用其中约 500 万笔(限定 \eta > 10^{-3} 以避免强噪声区)
  • Kibot 历史数据补充:Russell 3000 指数 3500 只股票的 1 分钟 OHLCV
  • 日波动代理 \sigma_D = (S_h - S_l)/S_o

三个刻画 metaorder 的变量:

记号 含义
\pi 日均相对份额 = 成交量 / 日成交量
\eta 参与率 = 交易速度 / 市场速度
F 持续时间(占总交易日的比例)

\pi = \eta \cdot F——两个变量就足以导出第三个。

三个参数的分布(论文的发现之一):

参数 分布 幂律指数
参与率 \eta 截断幂律(10^{-4} \le \eta \le 0.1 区间) a = -0.864 \pm 0.001
持续时间 F 截断幂律(0.01 \le F \le 0.5) a = -0.932 \pm 0.003
日均份额 \pi 明显不是幂律 —

论文特别指出 \pi 不是幂律这件事很重要:

这可能是一个重要结果,因为某些市场冲击理论的预测依赖于这一点,而且一般假设是幂律行为。

\eta 与 F 的线性相关只有 -0.022(几乎无关);贡献主要来自极端区域——\eta 很大意味着 F 很小,反之亦然(激进的 metaorder 通常很短)。

一个数据质量问题(论文的处理值得学):有些 metaorder 的起止时间明显不可靠(持续时间恰好是 410 分钟、从 9:30 到 16:20、参与率系统性偏低)。论文判断这些订单报给 Ancerno 的时间不准确,因此直接剔除——「代价是显著缩小了样本」。

核心结果一:对数函数拟合好得多

评估指标:加权均方根误差 E_{RMS}(按每个观测的标准误加权),越小越好。

两种函数形式:

f(\pi|Y,\delta) = Y\,\pi^{\delta}, \qquad g(\pi|a,b) = a\log_{10}(1+b\pi)

全样本拟合结果:

函数 最优参数 E_{RMS}
幂律(平方根) Y = 0.15 \pm 0.01,\delta = 0.47 \pm 0.02 6.70
对数 a = 0.028 \pm 0.001,b = 465 \pm 33 2.80

\delta = 0.47 与既有文献一致(Potters-Bouchaud、Lillo-Farmer-Mantegna、Gabaix 等)。

论文对双对数图上「看起来是直线」的辨析(这段是全文最关键的读法):

\pi 在 10^{-3} 到 10^{-1} 范围内,数据点第一近似落在一条直线上。然而明显的凹性是存在的——对很大的和很小的 \pi,冲击曲线都向下弯。

分年份与分市值的结果(表 4):

分组 幂律 \delta 幂律 E_{RMS} 对数 a 对数 b 对数 E_{RMS}
2007 0.41 5.47 0.029 491 1.36
2008 0.41 4.96 0.025 547 2.11
2009 0.46 3.18 0.032 316 2.23
大市值 0.51 1.38 0.030 441 0.40
中市值 0.46 1.32 0.030 400 0.69
小市值 0.42 0.69 0.027 428 1.03

除小市值分组外,对数函数总是更好。

一个重要的稳健性结论:临时冲击曲线的形状几乎不依赖交易时段(2007/08/09)也不依赖市值——这印证了平方根律「鲁棒」的说法,也说明换成对数形式后同样鲁棒。

核心结果二:平方根的「只依赖规模」是聚合假象

这是全文最有价值的发现。

平方根公式只依赖日均份额 \pi——这意味着临时冲击应当同时独立于 F 和 \eta。

论文按 F 与 \eta 分组检验,结论是否定的:

条件变量 幂律指数的变化
按持续时间 F 分组 从 0.58 降到 0.42
按参与率 \eta 分组 从 0.74 降到 0.47

(指数误差均为 0.02。)

机制(论文的解释,这段值得记住):

由于 \pi = \eta F,按 \pi 条件化意味着对 \eta-F 平面上从左下到右上的对角条带取平均。这个平均包含了正负残差、它们部分抵消,从而产生观察到的数据坍缩(collapse)。

所以「平方根律的普适性」在很大程度上是聚合的产物,而不是真的与执行参数无关。

核心结果三:冲击曲面与局部指数

「冲击曲面」概念:把 I_{tmp}(\Omega=\{\eta,F\}) 看作 (\eta, F) 平面上的曲面。

两种拟合的对比:

f(\eta,F|Y,\delta,\gamma_1) = Y\,\eta^{\delta}F^{\gamma_1}, \qquad g(\eta,F|a,b,c) = a\log_{10}(1+b\eta)\cdot\log_{10}(1+cF)

残差图(表 8):幂律拟合的残差有宏观结构(正负残差在平面上系统分布);对数拟合的残差这个模式被强烈削弱。

局部指数的估计结果(图 10):

指数 行为
\hat\delta(\eta,F) 小 \eta 时一致更大;大 \eta 时一致小于 0.5;变化范围显著,从约 0.1 到 0.9
\hat\gamma_1(\eta,F) 小 F 时接近 1;中等 F 时接近 0.3;大 F 时行为更复杂(低参与率时高指数,反之亦然)

论文的结论:

平方根预测的值 \gamma = \delta = 0.5 只在 \eta-F 平面的中央区域有效。

核心结果四:从冲击反推「潜在订单簿」

这是本文的方法论贡献。

Tóth 等 (2011) 的理论:真实订单簿不反映实际供需(参与者不透露真实意图),潜在订单簿只有在价格移动时才显形,因此可以从冲击形状反推。

方程(20):V(x|b,n) 是潜在订单簿在对数价格 x 处的量,尺寸为 \pi 的 metaorder 产生满足

\pi = \int_{x_0}^{x_0 + I(\pi)} V(x|b,n)\,dx

的冲击 I(\pi)。

关键的理论对照(论文明确列出):

潜在订单簿形状 冲击函数
对数价格线性 V \propto x 平方根(Tóth 等的主张)
常数 $V = $ const 线性

本文的参数化形式(式 21):

V(x|b,n) := \frac{1}{Y}\cdot\frac{x^n\exp(bx)}{\int_0^1 dy\, y^n\exp(by)}, \qquad x\in[0,1]

x \ll 1/b 时按多项式增长,x \gg 1/b 时按指数增长。

两个特例:

  • n = 0:恢复对数函数 I(\pi|Y,b,0) = Y\frac{\log(1+c\pi)}{\log(1+c)},其中 c = e^b - 1。\pi \ll 1/c 时线性增长,\pi \gg 1/c 时对数增长
  • n = 1:小 x 时按平方根增长,大 x 时对数增长

拟合结果:

模型 结论
n=0 比 n=1 好
线性/对数的分界 1/c \pi^\star \approx 2\times 10^{-3}
三参数模型 I(\pi\mid Y,b,n) 明显改善拟合;\hat n = 0.22

\hat n = 0.22 接近零,意味着最优买卖价附近的潜在订单簿几乎是平的(论文说左侧噪声较大)。

论文的结论:数据与大成交量下潜在订单簿呈指数形状一致。

核心结果五:执行不是 VWAP,且价格在执行结束前就开始回撤

论文问的问题(很直白):

给定两个参与率相同、持续时间不同的 metaorder(F_1 < F_2),第一个在 t = F_1 时达到的冲击,应该与第二个相同吗?

先验理由:此时可获得的公开信息相同,所以从市场冲击的角度它们应该无法区分。

实测(图 11):明确的否定答案。价格轨迹明显偏离由圆圈标出的临时冲击面。

两个发现:

  1. 小参与率时轨迹明显高于冲击面(更极端);参与率越大轨迹越接近圆圈
  2. 有些情形下价格在 metaorder 结束前就开始回撤

模型能解释多少(论文对比了两个模型):

模型 能解释 不能解释
Almgren-Chriss 前置加载的执行轨迹高于冲击面 不能重现冲击的凹性(它预测线性冲击)
传播子模型(Bouchaud 等、Lillo-Farmer) 凹性 ✓、平方根律 ✓、前置加载轨迹高于曲面 ✓ —
价格可能在执行结束前开始回撤 ✗

传播子模型怎么解释回撤:论文说——存在衰减核,使冲击有可能在执行结束前就开始回撤;如果交易压力弱于市场的恢复力量,冲击就开始回撤。

论文的推断:

这些结果强烈提示,我们数据库中经纪商的典型交易剖面不是广泛使用的 VWAP 交易剖面,而是前置加载的执行方案。这可能是出于风险厌恶,或为了避免错失盈利机会。

诚实的边界:论文明确说——「不幸的是,这些替代方案无法在 Ancerno 数据库包含的信息内检验。」

还有一个方法论说明:Ancerno 数据库不提供每笔子订单的数量、时间、成交量、价格,所以论文用 1 分钟分辨率的公开价格来跟踪价格动力学。论文说这个方法低估了执行期间回撤的幅度(因为只跟到 v = F_i)。

核心结果六:永久冲击约为峰值的 2/3,但精确值依赖 \eta 与 F

背景:Farmer 等 (2004) 在「metaorder 尺寸分布有指数 1.5 的幂律尾」的假设下,预测冲击衰减到峰值的大约 2/3。

论文立刻指出这个假设在数据里不成立:

注意这个假设在我们的数据中似乎没有被满足——如第 2 节所见,\pi 明显不是幂律分布。这有严重后果,因为在这种情况下他们的理论不预测平方根冲击函数。

做法:把时间重标度为 z = v/F(这样 metaorder 完成时统一为 z=1),冲击也除以完成时的冲击值,得到归一化的衰减轨迹 I_{ren}(z)。

结论:

平均而言冲击松弛到峰值冲击的约 2/3,但精确的渐近价格取决于 metaorder 的参与率与持续时间。

论文给出的解释:metaorder 之间的羊群效应(herding)。

羊群效应的观察证据(第 2 节图 1):在大多数交易日存在一种「情绪(mood)」——给定一天内大多数 metaorder 方向相同。

对既有理论的批评(第 5 节)

论文对三个既有理论框架做了评价:

框架 论文的评价
Tóth 等的潜在订单簿方法 本文直接沿用并做参数化扩展(引入 n);结论 \hat n = 0.22(接近平坦)
Farmer 等的公允定价方法 论文指出:其 2/3 永久冲击预测依赖 \pi 幂律假设,而该假设在数据中不成立
O'Hara 等(平方根律的理论来源之一) 批评要点:平方根律的普适性表象部分来自对不同 \eta、F 的订单取聚合平均

论文自陈的边界

  • Ancerno 不提供子订单级别信息——无法检验执行剖面的具体形状(这是第 4.3 节结论的直接限制)
  • 价格跟踪用 1 分钟公开数据,而非逐笔——低估执行期间的回撤幅度
  • 剔除了时间不可靠的 metaorder,代价是样本显著缩小
  • 只覆盖美股 2007–2009(含金融危机)
  • 冲击曲面分析限定 \eta > 10^{-3} 以避免噪声区

对本 wiki 的位置

与另外两篇「拒绝平方根」的对照

今天收录的三篇里,有两篇拒绝平方根——但对象不同:

论文 拒绝什么 数据
Almgren 等 2005 拒绝临时冲击 h(v) 的指数 =1/2(实测 \beta=3/5) 花旗美股大盘股经纪执行
本篇 拒绝总冲击对规模 Q 的指数恒为 1/2(对数拟合 E_{RMS} 2.80 vs 幂律 6.70) Ancerno 的机构 metaorder(含中小盘)
Sato-Kanazawa 2025 不拒绝——确认 \langle\delta\rangle = 0.489 东交所全账户

三者的差别可能来自数据源:Almgren 与本篇都没有全账户数据(前者只见花旗客户、后者只见 Ancerno 覆盖的机构);而 Sato-Kanazawa 用的是全账户。这是一个仍未解决的争议,详见 concepts/microstructure/冲击指数的三个口径 (Three Measures of the Impact Exponent)。

数值速查

量 值
样本 Ancerno metaorder,美股 2007–2009,6,944,883 笔(曲面分析约 500 万笔)
\eta 的幂律指数 -0.864 ± 0.001(10^{-4}–0.1 区间)
F 的幂律指数 -0.932 ± 0.003(0.01–0.5 区间)
\pi 的分布 明显不是幂律
\eta 与 F 的相关 -0.022
幂律全样本拟合 Y=0.15,\delta=\mathbf{0.47},E_{RMS}=\mathbf{6.70}
对数全样本拟合 a=0.028,b=465,E_{RMS}=\mathbf{2.80}
按 F 分组的指数 0.58 → 0.42
按 \eta 分组的指数 0.74 → 0.47
局部 \hat\delta 的范围 约 0.1 到 0.9
局部 \hat\gamma_1 小 F 时≈1,中 F 时≈0.3
潜在订单簿参数 \hat n = 0.22(接近平坦)
线性/对数的分界 \pi^\star \approx 2\times 10^{-3}
永久冲击/峰值 约 2/3(但精确值依赖 \eta、F)
冲击回撤 价格在执行结束前就开始(非 VWAP 执行)