风险模型与风险分析

4217 字
21 分钟
风险模型与风险分析

课程导航#

  1. 课时 1:信念、风险与量化投资流程(Beliefs, Risk, and Process)
  2. 课时 2:投资组合理论(Portfolio Theory)
  3. 课时 3:风险模型与风险分析(Risk Models and Risk Analysis)
  4. 课时 4:Alpha 因子评价(Evaluation of Alpha Factors)
  5. 课时 5:量化因子(Quantitative Factors)
  6. 课时 6:估值技术与价值创造(Valuation Techniques and Value Creation)
  7. 课时 7:多因子 Alpha 模型(Multifactor Alpha Models)
  8. 课时 8:组合换手与最优 Alpha 模型(Portfolio Turnover and Optimal Alpha Model)
  9. 课时 9:高级 Alpha 建模(Advanced Alpha Modeling Techniques)
  10. 课时 10:因子择时模型(Factor Timing Models)
  11. 课时 11:投资约束与信息比率(Portfolio Constraints and Information Ratio)
  12. 课时 12:交易成本与组合实施(Transaction Costs and Portfolio Implementation)

第一节 从 APT 到因子风险模型#

1. 定价命题与风险模型不是同一对象#

套利定价理论(Arbitrage Pricing Theory, APT 是一种多因子资产定价理论:证券收益服从线性因子结构,且同一未来现金流不能长期存在不同价格。在可分散特质风险、无套利等条件下,证券的期望超额收益由其系统性因子暴露及相应风险溢价决定:

E(ri)rf=k=1KXikλk.E(r_i)-r_f=\sum_{k=1}^{K}X_{ik}\lambda_k.
  • E(ri)rfE(r_i)-r_f:股票 ii 相对无风险利率的一期期望超额收益
  • XikX_{ik}:股票 ii 对因子 kk 的敏感度;量纲取决于因子定义与标准化方式。
  • λk\lambda_k:承担一单位因子 kk 暴露所对应的期望风险溢价,单位为收益率除以暴露单位。
  • KK:模型纳入的共同因子数。

多因子风险模型(Multifactor Risk Model) 是用共同因子与特质项估计证券收益协方差的模型。它使用 APT 的线性收益结构,却不要求每个因子都获得非零风险溢价;均值接近零但方差较大的共同冲击仍会影响组合风险。因而“因子能解释协方差”和“因子在均衡中被定价”是两个命题。

将一期因子冲击中心化后,NN 只股票、KK 个因子的收益生成式与协方差结构为:

rμ=Xf+u,E(f)=0,E(u)=0,Cov(f,u)=0,Cov(u)=D,Σ=XFX+D.\begin{aligned} \mathbf r-\boldsymbol\mu&=\mathbf X\mathbf f+\boldsymbol u,\\ E(\mathbf f)&=\mathbf0,\qquad E(\boldsymbol u)=\mathbf0,\\ \operatorname{Cov}(\mathbf f,\boldsymbol u)&=\mathbf0,\qquad \operatorname{Cov}(\boldsymbol u)=\mathbf D,\\ \boldsymbol\Sigma&=\mathbf X\mathbf F\mathbf X^{\prime}+\mathbf D. \end{aligned}
  • r\mathbf rμ\boldsymbol\muN×1N\times1 股票单期收益向量及其条件期望,期限和收益率口径必须一致。
  • X\mathbf XN×KN\times K 股票—因子暴露矩阵;第 ii 行描述股票 ii 的全部因子暴露。
  • f\mathbf fK×1K\times1 零均值因子冲击向量;F=Cov(f)\mathbf F=\operatorname{Cov}(\mathbf f)K×KK\times K 因子协方差矩阵。
  • u\boldsymbol uN×1N\times1 特质收益向量;D\mathbf DN×NN\times N 特质协方差矩阵,常被近似为对角矩阵。
  • Σ\boldsymbol\SigmaN×NN\times N 股票收益协方差矩阵,与 F\mathbf FD\mathbf D 使用同一频率和年化口径。

对总权重 w\mathbf w,组合因子暴露和方差分解为:

bP=Xw,σP2=bPFbP+wDw.\mathbf b_P=\mathbf X^{\prime}\mathbf w, \qquad \sigma_P^2=\mathbf b_P^{\prime}\mathbf F\mathbf b_P +\mathbf w^{\prime}\mathbf D\mathbf w.

w\mathbf wN×1N\times1 总权重列向量,bP\mathbf b_PK×1K\times1 组合因子暴露,σP\sigma_P 是与协方差期限一致的组合波动率。若分析相对基准的主动风险,应把 w\mathbf w 换成主动权重 h=wwb\mathbf h=\mathbf w-\mathbf w_b,此时 σP\sigma_P 的经济口径是跟踪误差,而不是组合总波动率。

“特质”只表示未被当前模型解释。遗漏行业、非线性或状态因子会进入残差;若残差仍有稳定相关结构,对角 D\mathbf D 会低估集中组合的风险。

2. 三类模型的因子与估计方向#

  • 宏观因子模型(Macroeconomic Factor Model):因子是利率、通胀、经济增长、信用利差或油价等可观测时间序列冲击;通常对每只股票做时间序列回归来估计暴露。因子有直接经济含义,但低频发布、数据修订、因子共线性和暴露漂移会削弱稳定性。
  • 基本面因子模型(Fundamental Factor Model):暴露来自决策时已知的行业分类及规模、估值、动量、杠杆等股票特征;每期以横截面回归估计随后实现的因子收益。它便于风险约束和归因,但分类、标准化、回归权重与特质风险设定均会改变结果。
  • 统计因子模型(Statistical Factor Model):因子是从历史收益协方差中提取的潜在共同成分,主成分分析是典型方法。它能压缩高维协方差,却依赖估计窗口;因子的符号、排序和经济含义可能跨期变化。

宏观模型主要从“因子冲击已知”出发估计股票暴露,基本面模型主要从“股票暴露已知”出发估计每期因子收益;两种回归方向不能混写。统计模型同时从收益数据中识别暴露和因子,其样本内解释率也不能证明因子被定价。

3. 主成分截断、动态估计与版本控制#

主成分分析(Principal Component Analysis, PCA 是把协方差矩阵分解为一组彼此正交、按解释方差从大到小排序的线性方向。对样本协方差矩阵 Σ^\widehat{\boldsymbol\Sigma} 做特征分解,并只保留最大的 KK 个主成分,可写成统计因子近似:

Σ^LKΛKLK+DK,RK2=k=1Kλkj=1Nλj.\widehat{\boldsymbol\Sigma} \approx\mathbf L_K\boldsymbol\Lambda_K\mathbf L_K^{\prime}+\mathbf D_K, \qquad R_K^2=\frac{\sum_{k=1}^{K}\lambda_k}{\sum_{j=1}^{N}\lambda_j}.
  • LK\mathbf L_KN×KN\times K 正交特征向量矩阵,行表示股票对保留主成分的载荷。
  • ΛK=diag(λ1,,λK)\boldsymbol\Lambda_K=\operatorname{diag}(\lambda_1,\ldots,\lambda_K):保留特征值构成的 K×KK\times K 对角矩阵,满足 λ1λK0\lambda_1\ge\cdots\ge\lambda_K\ge0
  • DK\mathbf D_K:截断后未由保留主成分解释的特质协方差近似。
  • RK2R_K^2:前 KK 个主成分解释的样本总方差比例,不是样本外预测准确率。

KK 太小会遗漏共同风险,太大会把噪声当成因子;即使 RK2R_K^2 很高,特征向量也可能在新窗口旋转。选择因子数应同时检查样本外协方差误差、特征空间稳定性和最终组合风险,而非只最大化样本内解释率。

基本面模型在时点 tt 用已知暴露 Xt\mathbf X_t 回归下一期横截面收益:

rt+1=at1+Xtft+1+ut+1.\mathbf r_{t+1}=a_t\mathbf1+\mathbf X_t\mathbf f_{t+1}+\boldsymbol u_{t+1}.

ata_t 是截面截距,Xt\mathbf X_t 必须只含时点 tt 可得信息,ft+1\mathbf f_{t+1} 是从 ttt+1t+1 的实现因子收益。多期回归得到的 {ft}\{\mathbf f_t\} 用于估计 F\mathbf F;残差序列及股票特征共同帮助估计 D\mathbf D。市值加权回归提高大盘股拟合度,却不等于对所有股票给出同等准确的风险预测。

权重半衰期(Weight Half-life) 是指数加权中,历史观测权重衰减为当前观测一半所需的期数。若对距当前 \ell 期的观测使用权重 ω\omega^\ell0<ω<10<\omega<1,则:

H=ln2lnω.H=-\frac{\ln 2}{\ln\omega}.

HH 与观测频率使用相同时间单位。较短半衰期对新状态反应更快、参数噪声也更大;较长半衰期更平滑、却可能在相关结构突变后滞后。实务估计还需处理异方差、极端收益、停牌、非同步交易与新股历史不足,并通过收缩或正定化避免协方差矩阵病态。

风险模型必须按时点保存暴露、行业分类、财务数据、股票池、参数和代码版本。事后行业重分类、最终修订财报和幸存股票池会重写历史风险;上线前至少检查矩阵对称性、半正定性、条件数、极端相关和样本外校准。

第二节 边际风险与可加总归因#

1. 边际风险是局部斜率#

  • 绝对风险(Absolute Risk):以总权重 w\mathbf w 计算组合收益 rP=wrr_P=\mathbf w^{\prime}\mathbf r 及波动率 σP\sigma_P
  • 主动风险(Active Risk):以主动权重 h=wwb\mathbf h=\mathbf w-\mathbf w_b 计算主动收益 rA=hrr_A=\mathbf h^{\prime}\mathbf r;其波动率是相对指定基准的跟踪误差。

以下先用绝对风险推导;把 w\mathbf w 换成 h\mathbf h、把 σP\sigma_P 换成跟踪误差即可得到主动风险版本。

边际风险贡献(Marginal Risk Contribution, MRC 是组合波动率对某项权重的局部偏导,表示其他风险资产权重不变时,微小增加该权重造成的波动率变化率:

σP=(wΣw)1/2,MRC=wσP=ΣwσP,MRCi=(Σw)iσP.\sigma_P=(\mathbf w^{\prime}\boldsymbol\Sigma\mathbf w)^{1/2}, \qquad \mathbf{MRC}=\nabla_{\mathbf w}\sigma_P =\frac{\boldsymbol\Sigma\mathbf w}{\sigma_P}, \qquad MRC_i=\frac{(\boldsymbol\Sigma\mathbf w)_i}{\sigma_P}.
  • w\mathbf wN×1N\times1 总权重列向量;若现金被视为无风险且不进入 Σ\boldsymbol\Sigma,单独改变 wiw_i 隐含由现金提供资金。
  • Σ\boldsymbol\SigmaN×NN\times N 证券收益协方差矩阵,单位为每期收益率平方。
  • σP\sigma_P:组合每期波动率,单位为收益率。
  • MRC\mathbf{MRC}N×1N\times1 边际风险向量;MRCiΔwiMRC_i\Delta w_i 是小幅权重变化造成的波动率一阶近似。

MRC 不是证券自身波动率。低波动证券若与当前组合高度同向,MRC 仍可能较高;负 MRC 表示在当前点增加该头寸可局部降低组合风险。

2. 受约束调仓使用方向导数#

完全投资组合不能只改变一项权重。令小额调仓为 Δw=εt\Delta\mathbf w=\varepsilon\mathbf t,则组合波动率的一阶变化为:

ΔσPMRCΔw=εtMRC.\Delta\sigma_P \approx\mathbf{MRC}^{\prime}\Delta\mathbf w =\varepsilon\mathbf t^{\prime}\mathbf{MRC}.
  • t\mathbf tN×1N\times1 调仓方向;完全投资要求 1t=0\mathbf1^{\prime}\mathbf t=0
  • ε\varepsilon:调仓尺度;若 ε=0.01\varepsilon=0.01,方向中一单位对应 1 个百分点权重。
  • tMRC\mathbf t^{\prime}\mathbf{MRC}:沿可行调仓方向的边际风险变化率,而非一个新组合的精确风险差。

等额买入股票 ii、卖出股票 jj 时,tMRC=MRCiMRCj\mathbf t^{\prime}\mathbf{MRC}=MRC_i-MRC_j。若组合还须保持市场 beta 中性,则调仓方向还需满足 βt=0\boldsymbol\beta^{\prime}\mathbf t=0;行业中性、资金中性与 beta 中性必须分别施加。有限调仓存在二阶曲率,幅度较大时应重新计算完整二次型,不能把局部斜率线性外推。

3. 风险贡献与欧拉分解#

  • 风险贡献(Risk Contribution, RC:当前权重乘以对应 MRC,把组合波动率按证券分配。
  • 百分比风险贡献(Percentage Risk Contribution, PCR:单项 RC 占组合波动率的比例;它分配的是当前风险预算,不是预测收益或最大可能损失。
RCi=wiMRCi=wi(Σw)iσP,PCRi=RCiσP=wi(Σw)iσP2.RC_i=w_iMRC_i =\frac{w_i(\boldsymbol\Sigma\mathbf w)_i}{\sigma_P}, \qquad PCR_i=\frac{RC_i}{\sigma_P} =\frac{w_i(\boldsymbol\Sigma\mathbf w)_i}{\sigma_P^2}.

波动率对权重是一阶齐次函数,欧拉定理给出:

i=1NRCi=σP,i=1NPCRi=1.\sum_{i=1}^{N}RC_i=\sigma_P, \qquad \sum_{i=1}^{N}PCR_i=1.

若改按方差归因,证券 ii 的方差贡献是 VCi=wi(Σw)iVC_i=w_i(\boldsymbol\Sigma\mathbf w)_i,贡献和为 σP2\sigma_P^2VCiVC_iRCi=VCi/σPRC_i=VC_i/\sigma_PPCRi=VCi/σP2PCR_i=VC_i/\sigma_P^2 的单位不同,不能在同一图表中混用。空头或对冲头寸可以产生负贡献,表示它降低当前组合风险,不表示该证券本身无风险。

例:从协方差到边际风险、风险贡献与调仓方向

设年化总权重 w=(0.6,0.4)\mathbf w=(0.6,0.4)^{\prime},年化协方差矩阵为:

Σ=[0.040.010.010.09].\boldsymbol\Sigma= \begin{bmatrix} 0.04&0.01\\ 0.01&0.09 \end{bmatrix}.

逐步计算得到:

Σw=(0.028,0.042),σP2=0.0336,σP=18.33%,MRC=(15.28%,22.91%),RC=(9.17%,9.17%),PCR=(50%,50%).\begin{aligned} \boldsymbol\Sigma\mathbf w&=(0.028,0.042)^{\prime},\\ \sigma_P^2&=0.0336,\qquad \sigma_P=18.33\%,\\ \mathbf{MRC}&=(15.28\%,22.91\%)^{\prime},\\ \mathbf{RC}&=(9.17\%,9.17\%)^{\prime},\qquad \mathbf{PCR}=(50\%,50\%)^{\prime}. \end{aligned}

两项方差贡献均为 0.01680.0168,除以 σP\sigma_P 后得到的波动率贡献均为 9.179.17 个百分点。若从资产 2 向资产 1 转移 1 个百分点权重,一阶风险变化为 (15.28%22.91%)×0.01=0.076(15.28\%-22.91\%)\times0.01=-0.076 个百分点;这是当前点附近的调仓近似,不是说资产 1 的 MRC 为负。

线60% 线 16% 20% 24% 28% 0% 20% 40% 60% 80% 100% 60%18.33% 72.7% 1 2 = 1 w 年化波动率
图 3-1:两资产完全投资组合的年化波动率曲线。实线使用示例协方差矩阵重新计算每个权重下的精确风险;虚线是资产 1 权重为 60% 时、沿“买入资产 1 并等额卖出资产 2”方向的局部切线。此时方向斜率为 MRC₁−MRC₂,转移 1 个百分点权重约使波动率下降 0.076 个百分点;远离切点后必须使用曲线而非线性外推。

4. 分组、因子与损失解释#

行业 GG 的风险贡献可以直接加总,RCG=iGRCiRC_G=\sum_{i\in G}RC_i;行业的边际风险则必须先指定一条满足资金和其他约束的调仓方向,不能只把行业内 MRC 相加。

对因子模型,因子与特质方差贡献满足:

k=1KbP,k(FbP)k共同因子方差贡献+i=1Nwi(Dw)i特质方差贡献=σP2.\underbrace{\sum_{k=1}^{K}b_{P,k}(\mathbf F\mathbf b_P)_k}_{\text{共同因子方差贡献}} +\underbrace{\sum_{i=1}^{N}w_i(\mathbf D\mathbf w)_i}_{\text{特质方差贡献}} =\sigma_P^2.

D\mathbf D 为对角矩阵,第 ii 项特质方差贡献简化为 wi2Diiw_i^2D_{ii};将各项方差贡献除以 σP\sigma_P 可得到同一波动率单位下的 RC。因子相关时,单个因子贡献包含协方差项并可能为负;贡献还取决于因子基底,旋转因子会改变归属但不改变总风险。行业与风格暴露不正交时,不能分别归因后再把两套结果相加。

PCR 还有条件性的损失解释。令 \ell 为实现的组合收益,则:

PCRi=Cov(wiri,rP)Var(rP).PCR_i =\frac{\operatorname{Cov}(w_i r_i,r_P)}{\operatorname{Var}(r_P)}.

若证券收益联合正态,条件期望为:

E(wirirP=)=wiμi+PCRi(μP).E(w_i r_i\mid r_P=\ell) =w_i\mu_i+PCR_i(\ell-\mu_P).

μi=E(ri)\mu_i=E(r_i)μP=E(rP)\mu_P=E(r_P)。在短期限内可合理近似各项期望收益为零时,给定组合发生收益 <0\ell<0,证券 ii条件期望损失占比等于 PCRiPCR_i;这不是对某次实际损失的确定性事后归因。对无约束、内部解的均值—方差最优组合,预期收益贡献比例也与 PCR 一致;现实约束或非最优持仓会破坏该等式。

第三节 尾部风险、模型诊断与组合决策#

1. VaR 是损失分位数,不是最坏损失#

定义持有期货币损失为 L=VrPL=-Vr_P风险价值(Value at Risk, VaR 是损失分布在置信水平 qq 的分位数:

VaRq=inf{l:Pr(Ll)q}.VaR_q=\inf\{l:\Pr(L\le l)\ge q\}.
  • LL:指定持有期的货币损失;当组合收益为负时 LL 为正。
  • VV:期初组合市值,单位为货币。
  • qq:置信水平,如 95%95\%99%99\%1q1-q 是超过该损失门槛的尾部概率。
  • VaRqVaR_q:与 LL 同单位的损失阈值,不描述超过阈值后的平均或最大损失。

若线性组合收益近似正态,参数 VaR 为:

VaRq=V(zqσPμP),zq=Φ1(q).VaR_q=V(z_q\sigma_P-\mu_P), \qquad z_q=\Phi^{-1}(q).

μP\mu_PσP\sigma_P 是和持有期匹配的组合期望收益与波动率,Φ1\Phi^{-1} 是标准正态分布的分位数函数。只有在收益独立同分布等额外条件下才可用平方根时间法则缩放波动率;不同持有期或置信水平的 VaR 不能直接比较。

教材级补充:预期尾部损失(Expected Shortfall, ES 是超过 VaR 门槛后的平均损失;对连续损失分布,ESq=E(LLVaRq)ES_q=E(L\mid L\ge VaR_q)ES 描述尾部严重度,而 VaR 只给出进入尾部的门槛,两者都不能替代流动性与情景压力测试。

2. VaR 的边际贡献与欧拉预算#

  • 边际 VaR(Marginal Value at Risk)VaR 对某项权重的局部偏导,记为 MVaRiMVaR_i
  • VaR 贡献(Contribution to Value at Risk):当前权重乘以边际 VaR,记为 RCiVaRRC_i^{VaR}
MVaRi=VaRqwi,RCiVaR=wiMVaRi.MVaR_i=\frac{\partial VaR_q}{\partial w_i}, \qquad RC_i^{VaR}=w_iMVaR_i.

VaR 对头寸可微且是一阶齐次函数时,欧拉分解成立:

i=1NRCiVaR=VaRq.\sum_{i=1}^{N}RC_i^{VaR}=VaR_q.

在线性正态模型中,令 μi\mu_i 为资产 ii 的持有期期望收益,则:

MVaRi=V(zqMRCiμi),RCiVaR=V(zqRCiwiμi).MVaR_i=V(z_qMRC_i-\mu_i), \qquad RC_i^{VaR}=V(z_qRC_i-w_i\mu_i).

忽略短期均值后,RCiVaR=zqVRCiRC_i^{VaR}=z_qVRC_i。沿用第二节两资产示例,取 V=100V=100 万元、q=95%q=95\%zq1.645z_q\approx1.645,一年期参数 VaR 约为 30.1530.15 万元;两项 VaR 贡献各约 15.0815.08 万元。输入协方差只是教学设定,不能视为未来真实参数。

对具有连续联合密度、唯一分位点且 VaR 可微的非正态损失,RCiVaRRC_i^{VaR} 可解释为组合损失恰处于 VaRqVaR_q 边界时,资产 ii 的条件期望损失贡献。该解释只对应这一分位边界;更换 qq 后必须重新计算。历史模拟的离散分位数可能不可微,蒙特卡洛或扰动重估还会引入带宽、抽样误差与数值不稳定。

厚尾、期权型头寸、跳跃、强制平仓和流动性枯竭都会破坏线性正态近似。压力测试应同时改变因子收益、相关性、特质冲击、融资条件和可交易性,不能只把常态协方差乘一个常数。

3. 从预测风险到决策审计#

  • 事前风险(Ex-ante Risk):使用决策时可得的暴露、协方差、持仓与基准版本预测未来风险。
  • 实现风险(Realized Risk):由决策后实际收益路径估计的风险;它混合真实参数变化、抽样波动、成交偏差和持仓漂移。

二者偏差可能源于参数误差、暴露漂移、模型遗漏、约束绑定、成交偏差与状态切换。风险模型不能只靠“长期平均预测接近实现值”通过校验,还应检查:

  • 标准化主动收益残差 zt+1=(rA,t+1α^t)/TEtz_{t+1}=(r_{A,t+1}-\widehat\alpha_t)/TE_t 是否存在持续偏离零、方差偏离 1、自相关或厚尾,其中 α^t\widehat\alpha_tTEtTE_t 是时点 tt 对下一持有期的预期主动收益与事前跟踪误差。
  • 不同置信水平下 VaR 超越频率及超越是否聚集;平均覆盖正确但损失集中在压力期,仍表示条件风险失准。
  • 行业、风格和特质贡献是否稳定,预测贡献较大的来源是否也解释了随后损失。
  • 模型换版本、相关性突变、暴露漂移和流动性压力下,风险排序与组合决策是否发生不可解释跳变。

风险分析服务于组合选择,不要求所有贡献相等。减小最大风险贡献者未必提高组合效用:该头寸可能同时拥有最高预期 alpha,替代交易还会产生成本。可执行比较应同时报告交易前后的预期 alpha、总风险、集中度、流动性、约束余量与净成本,并保存当时可得的模型版本供事后审计。

文章分享

如果这篇文章对你有帮助,欢迎分享给更多人!

评论区

Profile Image of the Author
罗小饭
AI 软件开发工程师,热爱编程、开源与折腾,长期记录 AI、网络、软路由、博客搭建和云服务实践中的经验与思考。
公告
欢迎来到我的博客!这是一则示例公告。
分类
标签
站点统计
文章
47
分类
5
标签
144
总字数
526,682
运行时长
0
最后活动
0 天前

文章目录