1. 项目概述当收益率曲线遇见宏观经济“变脸”在固定收益领域收益率曲线建模是每个从业者绕不开的核心课题。无论是为债券定价、管理利率风险还是解读央行的政策信号一个靠谱的模型都是我们手中的“导航仪”。从业十几年我见过太多模型在样本内拟合得天花乱坠一到现实检验就“水土不服”。问题的核心往往在于传统的线性模型比如经典的动态Nelson-Siegel模型假设市场环境是稳定不变的但现实中经济总是在扩张、收缩、通胀、通缩等不同“状态”间切换。这就好比用一套固定的驾驶模式去应对城市拥堵、高速公路和乡村泥泞路显然会力不从心。宏观工具化机制转换动态Nelson-Siegel模型正是为了解决这个痛点而生。它不再将收益率曲线的动态视为一成不变而是承认其演化规律会随着宏观经济状态如高利率环境、低通胀时期的改变而“变脸”。这个项目的核心创新在于它没有采用传统隐马尔可夫模型HMM那种“黑箱”式的状态划分而是直接用我们看得见、摸得着的宏观经济变量如联邦基金利率、通货膨胀率作为“标尺”通过一个贝叶斯决策树结构清晰地将历史时期划分成几个具有明确经济含义的“机制”。然后模型允许收益率曲线的核心参数如因子的持续性、波动性以及宏观因子对收益率因子的影响在这些不同机制下取不同的值。这套方法的技术价值极高。首先它极大地提升了模型对现实尤其是对金融危机、政策转向等极端时期的刻画能力。其次它为我们检验“宏观跨越”假说——即当前收益率曲线是否已经包含了所有相关的宏观经济信息——提供了一个更精细的透镜也许在“风平浪静”的常态下宏观变量确实没有额外的预测力但在“波涛汹涌”的非常时期宏观冲击对曲线的影响可能至关重要。对于从事利率衍生品定价、资产负债管理或宏观经济研究的同行来说理解这种非线性特征意味着能更精准地评估风险、捕捉机会。2. 核心思路与模型架构拆解2.1 从静态到动态Nelson-Siegel模型的演进要理解这个项目的精妙之处得从基础说起。经典的Nelson-Siegel模型是一个静态的收益率曲线拟合公式。在任意一个时间点t对于期限为τ的利率y_t(τ)模型将其分解为三个因子的线性组合y_t(τ) L_t S_t * [(1 - e^{-λτ})/(λτ)] C_t * [(1 - e^{-λτ})/(λτ) - e^{-λτ}]这里的L_t,S_t,C_t就是大名鼎鼎的水平、斜率、曲率因子。L_t代表长期利率水平S_t是短长期利差通常S_t y_t(短期) - y_t(长期)为负时曲线倒挂C_t则捕捉曲线中间部分的凸度。λ是一个衰减参数控制曲率因子载荷函数的形态。这个模型的优美之处在于其极强的经济学直觉和惊人的简洁性仅用三个因子就能灵活拟合出从向上倾斜到向下倾斜从陡峭到平坦的各种曲线形态。然而静态模型无法描述这些因子如何随时间变化。Diebold和Li在2006年的开创性工作将其动态化提出了动态Nelson-Siegel模型。他们让L_t,S_t,C_t这三个因子遵循一个向量自回归过程比如一阶VAR(1)F_t A * F_{t-1} η_t其中F_t [L_t, S_t, C_t]。这样一来模型就变成了一个状态空间模型测量方程是上面的NS公式将可观测的各类期限利率与不可观测的潜在因子联系起来状态转移方程则描述了这些因子自身的动态演化。这个DNS模型成为了学术和业界的一个基准。注意在早期实践中衰减参数λ常常被预先校准为一个固定值如0.0609以简化估计。但在更严谨的建模中尤其是在引入机制转换后将λ作为待估参数可以增加模型的灵活性使其能更好地适应不同期限结构的数据特征。2.2 引入宏观经济与机制转换模型的两次关键升级单纯的DNS模型仍是一个“纯收益率”模型。Diebold, Rudebusch和Aruoba在2006年进一步将其扩展为**“收益率-宏观”模型**。思路很直接既然利率是宏观经济的重要变量那么将宏观因子与收益率因子放在同一个状态空间里建模是顺理成章的。通常我们会选取几个核心宏观变量例如产能利用率代表实体经济活动的“量”。联邦基金利率代表中央银行货币政策立场的“价”。通货膨胀率代表名义经济活动的“价格水平”。将这三个宏观变量m_t [CU_t, FFR_t, INFL_t]与三个收益率因子f_t [L_t, S_t, C_t]拼接形成一个六维的状态向量FF_t [f_t, m_t]。其状态转移方程变为一个六维的VAR(1)过程。这个扩展使得我们可以直接研究宏观冲击如何影响收益率曲线以及反过来收益率曲线是否蕴含了未来的宏观信息。但上述模型仍是“单机制”的即假设整个样本期内状态转移矩阵A和扰动项协方差矩阵H是恒定不变的。这与经济周期和金融市场的结构性变化事实相悖。因此本项目引入了机制转换。其核心思想是允许模型的关键参数包括状态转移矩阵A、因子均值µ、因子创新协方差矩阵H在不同的“机制”或“状态”下取不同的值。如果用z_t ∈ {1, 2, ..., G}表示t时期所处的机制那么状态转移方程就变为FF_t A_{z_{t-1}} * FF_{t-1} η_t且η_t ~ N(0, H_{z_t})。这意味着在高通胀机制和低增长机制下收益率曲线的演化规律、波动率以及宏观与收益率的互动关系可能是完全不同的。2.3 核心创新基于宏观变量的树结构机制检测传统的机制转换模型如马尔可夫机制转换模型存在一个重大缺陷机制z_t是隐状态通过统计方法估计出来但缺乏清晰的经济解释。我们可能得到“机制A”和“机制B”但很难说清A到底对应着什么样的经济环境。本项目的最大亮点就是彻底解决了这个“黑箱”问题。它借鉴了机器学习中分类与回归树的思想但目标不是预测而是基于模型拟合优度进行机制划分。具体步骤如下选择分割变量与阈值候选分割变量是标准化后的宏观经济变量如FFR, INFL候选阈值是其历史分位数如0.2, 0.4, 0.6, 0.8。例如一个候选分割规则是“FFR_t ≤ 0.6”即联邦基金利率是否低于其历史60%分位数。评估分割效果对于每一个候选分割规则它将当前样本比如根节点即所有时期划分成两个子集对应两个潜在的机制。然后在这两个子集分别估计“收益率-宏观”模型并计算整个样本的边际似然。选择最优分割选择那个能使整体边际似然最大的分割规则。边际似然是贝叶斯模型比较的黄金标准因为它综合考虑了模型拟合度和复杂度通过积分掉所有参数实现。递归分割对分割后产生的每个子节点即每个新机制重复步骤1-3继续寻找最优分割直到满足停止条件如机制数达到上限3或某个机制内的样本数少于24个月。这个过程最终会生成一个决策树。例如实证结果发现的第一层分割是FFR ≤ 0.6第二层在FFR 0.6的分支上对INFL ≤ 0.6进行分割。这样就得到了三个具有明确经济含义的机制机制1FFR ≥ 0.6即高利率环境。机制2FFR 0.6且INFL 0.6即低利率且低通胀环境。机制3FFR 0.6且INFL ≥ 0.6即低利率但高通胀环境。这种方法的优势是革命性的机制是可观测、可解释的直接由宏观经济状况定义。它完美地将经济直觉与统计建模结合了起来。3. 模型估计的实战细节与贝叶斯实现3.1 状态空间模型与贝叶斯估计框架将上述思路转化为可估计的模型我们得到一个包含机制转换的“收益率-宏观”状态空间模型。测量方程将不同期限的观测收益率y_t与潜在因子包括收益率因子和宏观因子联系起来转移方程则描述了因子在特定机制下的VAR(1)动态。模型的参数集Θ非常庞大包括载荷矩阵Λ由衰减参数λ决定。各机制下的状态转移矩阵A_g、因子均值µ_g、因子创新协方差矩阵H_g。测量误差的协方差矩阵Q通常假设为对角阵。机制序列{z_t}由决策树确定。面对如此复杂的模型频率学派的极大似然估计会异常困难尤其是在处理潜在因子和机制转换时。因此贝叶斯方法搭配马尔可夫链蒙特卡洛估计成为了自然且强大的选择。3.2 先验分布的选择与设定技巧贝叶斯估计始于先验。合理的先验能稳定估计特别是在样本量有限或模型复杂时。状态转移矩阵A_g的先验这是关键。我们希望对A_g进行稀疏化估计即认为不是所有变量之间都存在强烈的动态关联。为此对A_g的非对角元素使用了钉螺先验。具体来说对于A_g的第(j,k)个非对角元素a_{jk}^g我们设定一个混合先验π(a_{jk}^g | γ_{jk}^g) ∼ (1 - γ_{jk}^g) * N(0, ξ_0^2) γ_{jk}^g * N(0, ξ_1^2)其中γ_{jk}^g ∼ Bernoulli(0.5)是一个潜变量ξ_0^2是一个很小的值“钉”ξ_1^2是一个较大的值“螺”。如果数据支持a_{jk}^g不为零MCMC采样会使γ_{jk}^g倾向于1从而让a_{jk}^g从一个方差较大的正态分布中采样反之则倾向于从方差极小的正态分布中采样 effectively将其收缩至0。这相当于一种贝叶斯变量选择能自动识别出重要的跨因子影响。协方差矩阵的先验对于因子创新协方差矩阵H_g我们使用共轭的逆Wishart先验。对于对角化的测量误差协方差矩阵Q其每个对角元素σ_i^2使用逆Gamma先验。这些都是状态空间模型中的标准选择。其他参数因子均值µ_g使用正态先验衰减参数λ使用均匀先验。3.3 吉布斯采样与卡尔曼滤波的协同在给定机制序列{z_t}由决策树确定后模型的估计通过吉布斯采样器进行。这是一个迭代算法在每一步我们依条件分布从其他参数和潜在因子的当前值中抽样出某一个参数或因子。核心难点在于潜在因子F_t的抽样。由于模型是线性和高斯的我们可以利用卡尔曼滤波和平滑算法来高效地从F_t的全条件后验分布中抽样。具体步骤是给定模型参数A_g,H_g,Q等运行前向的卡尔曼滤波得到每一步的预测分布。运行后向的平滑算法从整个序列F_{1:T}的联合后验分布中抽取一条路径。这个过程在吉布斯采样的每次迭代中都要执行以更新潜在因子的值。其他参数如A_g,H_g在给定因子F_t和机制z_t后其全条件后验分布由于共轭先验的存在通常是标准分布如多元正态、逆Wishart易于抽样。实操心得在实现这个采样器时对卡尔曼滤波的数值稳定性要格外小心。对于长期数据协方差矩阵可能变得非正定。一个实用的技巧是在滤波过程中使用平方根滤波或添加一个极小的扰动项来确保数值稳定。此外由于机制转换我们需要为每个机制存储一套参数并在时间序列上根据z_t切换使用哪套参数在编程实现时索引管理要清晰。3.4 边际似然的计算机制选择的关键如前所述决策树选择分割规则的核心标准是边际似然p(Y, M | 树结构)。计算状态空间模型的边际似然没有解析解。本项目采用了一个巧妙的数值方法首先在给定参数下通过对潜在因子F_t积分得到“部分边际似然”一个多元正态密度这可以通过卡尔曼滤波的一次前向传递计算出来。然后对于剩下的参数A_g,H_g,Q等我们利用吉布斯采样得到的后验样本{Θ^{(s)}}进行数值积分。具体来说计算每个后验样本s对应的部分边际似然值p(Y, M | Θ^{(s)}, 树结构)然后对这些值取平均。这相当于用后验样本来近似对参数空间的积分。这个计算量很大因为每评估一个候选分割规则都需要运行一次完整的MCMC来获取后验样本并计算平均边际似然。但这是确保机制划分基于整体模型拟合度而非过拟合某个局部特征的可靠方法。4. 实证分析美国国债收益率曲线的三机制解读4.1 数据准备与描述性分析我们使用Liu和Wu2021构建的美国国债零息收益率月度数据期限从3个月到120个月共13个期限样本从1971年8月到2022年12月。宏观变量采用制造业产能利用率、联邦基金利率和年化通货膨胀率。从描述性统计中我们可以验证几个经典事实收益率曲线通常向上倾斜长期利率高于短期短期利率的波动性远大于长期利率长期利率的自相关性更强。这些都为后续建模提供了背景。4.2 机制识别结果与经济解释运行前述的树结构算法后模型识别出了三个机制其划分清晰而富有经济意义机制1高利率环境FFR ≥ 0.6。对应美联储采取紧缩货币政策、主动提高利率以对抗通胀或经济过热的时期。样本中包含191个月。机制2低利率低通胀环境FFR 0.6且INFL 0.6。这是常见的“金发姑娘”经济状态增长温和、通胀可控、货币政策宽松。这是样本中最主要的机制包含310个月。机制3低利率高通胀环境FFR 0.6且INFL ≥ 0.6。这是一种相对棘手的状态通胀已经抬头但政策利率仍处于低位。这可能发生在通胀初期或央行出于对经济增长的担忧而延迟加息时。包含116个月。这个划分完美捕捉了二战后美国宏观经济政策的几个典型阶段例如上世纪70-80年代的“大通胀”与沃尔克加息时期机制1为主2008年金融危机后的长期低通胀、零利率下限时期机制2为主以及2021-2022年疫情后供给冲击带来的高通胀与利率快速上升期从机制3向机制1转换。4.3 参数估计与“宏观跨越”假说的再检验模型估计出的参数揭示了不同机制下截然不同的动态。下表摘要展示了三个机制下状态转移矩阵A中部分关键系数的后验均值括号内为标准差机制参数组水平(L)斜率(S)曲率(C)通胀(INFL)对曲率(C)的影响机制1因子持续性0.98*0.92*0.78*0.13*(高利率)(A矩阵对角元)(0.02)(0.03)(0.09)(0.04)机制2因子持续性0.99*0.99*0.93*0.00(低利率低通胀)(A矩阵对角元)(0.01)(0.01)(0.03)(0.01)机制3因子持续性0.96*0.95*0.92*0.00(低利率高通胀)(A矩阵对角元)(0.02)(0.02)(0.03)(0.01)*注*表示95%后验可信区间不包含0。最重要的发现关乎“宏观跨越”假说在机制2和机制3中宏观变量通胀、产能利用率、联邦基金利率对三个收益率因子的预测系数几乎都不显著后验区间包含0。这意味着在这两种状态下当期的收益率曲线已经基本包含了所有关于未来宏观经济的相关信息。宏观变量没有提供额外的预测能力“宏观跨越”假说大致成立。这符合传统仿射模型的理论预期。然而在机制1高利率环境下情况完全不同。通货膨胀对曲率因子有显著的正向影响系数0.13。这意味着在高利率机制下通胀的上升会显著影响收益率曲线的中期形态曲率。此时宏观变量具有独立的预测内容“宏观跨越”假说被拒绝。这个发现具有深刻的政策含义。它表明在货币政策紧缩期通胀动态会以一种非线性的方式“渗入”收益率曲线的定价中这是单机制或线性模型无法捕捉的。对于交易员而言在高利率环境下需要格外关注通胀数据对中期利率如5-7年期的潜在冲击。4.4 模型拟合效果与机制对比从拟合优度看无论是全样本还是分机制模型拟合的收益率曲线与真实曲线都高度吻合尤其是在中长期限。这证明了模型在捕捉不同形态曲线方面的能力。通过比较不同机制下参数的后验密度分布如图7所示我们可以更直观地看到差异因子持续性水平因子A_{1,1}在三个机制下都很高且接近但斜率因子A_{2,2}和曲率因子A_{3,3}的持续性在机制1中明显更低说明在高利率环境下曲线的斜率和形态更易发生改变。波动性因子创新方差H的对角元素在后验密度上重叠很少。例如水平因子的波动性H_{1,1}在机制1高利率中最大在机制2低利率低通胀中最小。这符合直觉货币政策剧烈变动期整体利率水平的不确定性也最高。5. 常见问题、挑战与实战建议5.1 模型选择与过拟合问题Q为什么选择最多三个机制为什么机制内最少需要24个月的数据A这是模型复杂性与样本信息量之间的权衡。我们的月度数据跨度约50年600个月。过多的机制如4个或5个会导致每个机制内的样本量过少参数估计不精确模型稳定性差。限制每个机制至少24个月2年是确保参数估计有一定可靠性的经验法则。这三个机制已经能够清晰对应主要的经济状态紧缩、温和、滞胀/类滞胀。在实际应用中可以根据数据频率和长度调整这些阈值。Q如何避免决策树过拟合A我们的方法本质上是基于模型比较的预设定而非纯粹的数据挖掘。分割标准是边际似然它天然包含了对模型复杂度的惩罚贝叶斯奥卡姆剃刀。此外我们只考虑由宏观经济分位数定义的、有经济意义的简单分割规则如“利率是否高于历史中位数”而不是在连续空间上搜索最优分割点这大大减少了过拟合的风险。当然在样本外预测中检验模型的稳健性始终是必要的。5.2 计算效率与实现难点Q贝叶斯估计计算量巨大有何优化技巧A确实这是本项目的主要挑战。以下几点可以提升效率并行计算评估不同的候选分割规则是相互独立的可以完美并行化。简化模型在构建决策树的搜索阶段可以使用简化版的模型例如固定λ或对H_g使用更简单的先验进行快速评估以筛选出最有希望的几个分割规则再用完整模型精算边际似然。利用共轭性质在吉布斯采样中确保先验选择是共轭的可以保证从标准分布中抽样速度远快于Metropolis-Hastings等通用算法。代码优化使用C、Rcpp或Julia重写卡尔曼滤波和采样循环等核心部分能带来数量级的速度提升。Q如何处理参数的可识别性问题A在状态空间模型中潜在因子的尺度和旋转存在不确定性。在Nelson-Siegel框架下通常通过固定载荷矩阵Λ的结构由λ决定来识别水平、斜率、曲率因子。此外对因子均值µ_g施加适当的先验如将其中心化有助于稳定估计。在机制转换模型中确保不同机制的参数有足够差异也是避免标签切换问题的关键。5.3 模型的应用与扩展Q这个模型如何用于预测A预测需要分两步机制预测首先需要预测未来时段th处于哪个机制。由于机制由可观测宏观变量FFR, INFL决定我们需要对宏观变量进行预测可以使用简单的AR模型或专业的宏观预测模型。将宏观变量的预测值代入训练好的决策树即可得到未来的机制标签z_{th}。收益率预测在给定预测机制z_{th}后使用该机制下估计的参数A_{z}, H_{z}, µ_z和当前的状态向量FF_t通过状态转移方程迭代向前预测FF_{th}再通过测量方程得到收益率曲线y_{th}(τ)的预测分布。Q模型可以如何扩展A有几个有前景的方向纳入更多宏观或金融变量如信用利差、股市波动率VIX、期限溢价等作为机制划分的候选变量或直接纳入状态向量。时变衰减参数λ允许λ随机制或时间变化可以增加模型对短期和超长期利率的拟合灵活性。应用于其他市场将模型应用于欧元区、英国、日本等国债市场检验机制划分的普适性。连接至无套利定价在DNS模型基础上引入无套利约束使其更适用于衍生品定价同时保留机制转换特征。在我个人的实践中将这种具有清晰经济解释的机制转换模型应用于投资组合的久期管理和凸性调整效果显著优于单机制模型。它帮助我们在2022年利率快速上行初期更早地识别出机制从“低利率低通胀”向“低利率高通胀”再向“高利率”转换的风险从而及时降低了组合的利率风险暴露。模型的价值不在于其复杂性而在于它为我们理解市场所处的“状态”提供了一个量化的、基于数据的视角。