单因子模型到多因子模型

从单因子模型到多因子模型 – 潘登同学的Quant笔记

文章目录

    • 从单因子模型到多因子模型 -- 潘登同学的Quant笔记
  • 单因子模型、多因子模型
    • 拓展到多因子的依据是什么?
      • C-CAPM框架下的单因子
      • C-CAPM框架下的多因子
  • APT
    • 推导多因子模型
      • 单因子两资产
      • 多因子多资产
  • APT的应用
    • α , β \alpha,\beta α,β分离
    • 因子选股
    • 统计套利

单因子模型、多因子模型

根据CAPM:
E [ r j ~ ] = r f + β M , j ( E [ r M ~ ] − r f ) E[\tilde{r_j}] = r_f + \beta_{M,j}(E[\tilde{r_M}]-r_f) E[rj~]=rf+βM,j(E[rM~]rf)
将其写成计量模型,很容易就得出了单因子模型
r j ~ − r f = α i + β i , M ( r M ~ ] − r f ) + ϵ i ~ \tilde{r_j} - r_f = \alpha_i + \beta_{i,M}(\tilde{r_M}]-r_f) + \tilde{\epsilon_i} rj~rf=αi+βi,M(rM~]rf)+ϵi~
Fama和French 1993年指出可以建立一个三因子模型来解释股票回报率
r j ~ − r f = α i + β i , M ( r M ~ ] − r f ) + β i , s S M B i ~ + β i , h H M L i ~ + ϵ i ~ \tilde{r_j} - r_f = \alpha_i + \beta_{i,M}(\tilde{r_M}]-r_f) + \beta_{i,s}\tilde{SMB_i} + \beta_{i,h}\tilde{HML_i} + \tilde{\epsilon_i} rj~rf=αi+βi,M(rM~]rf)+βi,sSMBi~+βi,hHMLi~+ϵi~
其中, S M B i ~ \tilde{SMB_i} SMBi~表示市值因素, H M L i ~ \tilde{HML_i} HMLi~表示账面价值与股票市值之比;

拓展到多因子的依据是什么?

C-CAPM框架下的单因子

回到C-CAPM,代表性消费者的优化问题

max ⁡ c 0 u ( c 0 ) + δ E [ u ( c 1 ~ ) ] s . t . c 1 ~ = ( 1 + r w ~ ) ( w 0 − c 0 ) F O C : 1 = E [ δ u ′ ( c 1 ~ ) u ′ ( c 0 ) ( 1 + r w ~ ) ] \begin{aligned} \max_{c_0} \quad & u(c_0) + \delta E[u(\tilde{c_1})] \\ s.t. \quad & \tilde{c_1} = (1+\tilde{r_w})(w_0-c_0) \\ FOC: & 1 = E[\delta\frac{u'(\tilde{c_1})}{u'(c_0)}(1+\tilde{r_w})] \end{aligned} c0maxs.t.FOC:u(c0)+δE[u(c1~)]c1~=(1+rw~)(w0c0)1=E[δu(c0)u(c1~)(1+rw~)]
其中 r w ~ \tilde{r_w} rw~为市场组合的回报率, w 0 w_0 w0为消费者0时期的财富;

假设消费者的效用函数是二次型 u ( c ) = − 1 2 ( a − c ) 2 u(c)=-\frac{1}{2}(a-c)^2 u(c)=21(ac)2,则 u ′ ( c ) = a − c u'(c) = a-c u(c)=ac,则随机折现因子可以表示为
m ~ = δ u ′ ( c 1 ~ ) u ′ ( c 0 ) = δ a − ( 1 + r w ~ ) ( w 0 − c 0 ) a − c 0 = A − B r w ~ A = δ a a − c 0 − δ w 0 − c 0 a − c 0 , B = δ w 0 − c 0 a − c 0 \tilde{m} = \delta\frac{u'(\tilde{c_1})}{u'(c_0)} = \delta\frac{a- (1+\tilde{r_w})(w_0-c_0)}{a-c_0} = A - B\tilde{r_w} \\ A=\delta\frac{a}{a-c_0} - \delta\frac{w_0-c_0}{a-c_0},B=\delta \frac{w_0-c_0}{a-c_0} m~=δu(c0)u(c1~)=δac0a(1+rw~)(w0c0)=ABrw~A=δac0aδac0w0c0,B=δac0w0c0
之前在C-CAPM推出来的结论(对所有资产成立,自然对无风险资产成立)
1 = E [ m ~ ( 1 + r f ~ ) ] ⇒ 1 = E [ m ~ ] ( 1 + r f ) + c o v ( m ~ , r f ) ⇒ 1 E [ m ~ ] − 1 = r f 1 = E[\tilde{m}(1+\tilde{r_f})] \\ \Rightarrow 1 = E[\tilde{m}](1+r_f)+cov(\tilde{m},{r}_f) \Rightarrow \frac{1}{E[\tilde{m}]} -1 = r_f 1=E[m~(1+rf~)]1=E[m~](1+rf)+cov(m~,rf)E[m~]11=rf
对于任意资产 i i i:
1 = E [ m ~ ( 1 + r j ~ ) ] ⇒ 1 = E [ m ~ ] ( 1 + E [ r j ~ ] ) + c o v ( m ~ , r ~ j ) ⇒ E [ r j ~ ] = 1 E [ m ~ ] − 1 − c o v ( m ~ , r ~ j ) E [ m ~ ] = r f − c o v ( m ~ , r ~ j ) E [ m ~ ] = r f − c o v ( A − B r w ~ , r ~ j ) E [ m ~ ] = r f + c o v ( r w ~ , r ~ j ) v a r ( r w ~ ) B v a r ( r w ~ ) E [ m ~ ] \begin{aligned} 1 &= E[\tilde{m}(1+\tilde{r_j})]\\ \Rightarrow 1 &=E[\tilde{m}](1+E[\tilde{r_j}])+cov(\tilde{m},\tilde{r}_j) \\ \Rightarrow E[\tilde{r_j}] &= \frac{1}{E[\tilde{m}]} - 1 - \frac{cov(\tilde{m},\tilde{r}_j)}{E[\tilde{m}]} \\ &=r_f - \frac{cov(\tilde{m},\tilde{r}_j)}{E[\tilde{m}]} \\ &=r_f - \frac{cov(A - B\tilde{r_w},\tilde{r}_j)}{E[\tilde{m}]} \\ &=r_f + \frac{cov(\tilde{r_w},\tilde{r}_j)}{var(\tilde{r_w})} \frac{B var(\tilde{r_w})}{E[\tilde{m}]} \\ \end{aligned} 11E[rj~]=E[m~(1+rj~)]=E[m~](1+E[rj~])+cov(m~,r~j)=E[m~]11E[m~]cov(m~,r~j)=rfE[m~]cov(m~,r~j)=rfE[m~]cov(ABrw~,r~j)=rf+var(rw~)cov(rw~,r~j)E[m~]Bvar(rw~)
定义 β \beta β:
β j , w ≜ c o v ( r w ~ , r ~ j ) v a r ( r w ~ ) \beta_{j,w} \triangleq \frac{cov(\tilde{r_w},\tilde{r}_j)}{var(\tilde{r_w})} βj,wvar(rw~)cov(rw~,r~j)
定义常数 λ \lambda λ
λ w ≜ B v a r ( r w ~ ) E [ m ~ ] \lambda_w \triangleq \frac{B var(\tilde{r_w})}{E[\tilde{m}]} λwE[m~]Bvar(rw~)
所以资产j的期望回报率可以写成
E [ r j ~ ] = r f + β j , w λ w E[\tilde{r_j}] = r_f + \beta_{j,w}\lambda_w E[rj~]=rf+βj,wλw

C-CAPM框架下的多因子

现在我们假设消费者除了拥有初始财富 w 0 w_0 w0 外,还在 0 期与 1 期分别拥有工资性收入 y 0 y_0 y0 y 1 ~ \tilde{y_1} y1~,则最优化问题为
max ⁡ c 0 u ( c 0 ) + δ E [ u ( c 1 ~ ) ] s . t . c 1 ~ = ( 1 + r w ~ ) ( w 0 + y 0 − c 0 ) + y 1 ~ F O C : 1 = E [ δ u ′ ( c 1 ~ ) u ′ ( c 0 ) ( 1 + r w ~ ) ] \begin{aligned} \max_{c_0} \quad & u(c_0) + \delta E[u(\tilde{c_1})] \\ s.t. \quad & \tilde{c_1} = (1+\tilde{r_w})(w_0+y_0-c_0) + \tilde{y_1} \\ FOC: & 1 = E[\delta\frac{u'(\tilde{c_1})}{u'(c_0)}(1+\tilde{r_w})] \end{aligned} c0maxs.t.FOC:u(c0)+δE[u(c1~)]c1~=(1+rw~)(w0+y0c0)+y1~1=E[δu(c0)u(c1~)(1+rw~)]
则随机折现因子表示为
m ~ = A − B r w ~ − C y 1 ~ A = δ a a − c 0 − δ w 0 + y 0 − c 0 a − c 0 , B = δ w 0 + y 0 − c 0 a − c 0 , C = δ a a − c 0 \tilde{m} = A - B\tilde{r_w} - C\tilde{y_1} \\ A=\delta\frac{a}{a-c_0} - \delta\frac{w_0+y_0-c_0}{a-c_0},B=\delta \frac{w_0+y_0-c_0}{a-c_0},C = \delta \frac{a}{a-c_0} m~=ABrw~Cy1~A=δac0aδac0w0+y0c0,B=δac0w0+y0c0,C=δac0a
那么最终会得到
E [ r j ~ ] = r f + β j , w ′ λ w ′ + β j , y ′ λ y ′ E[\tilde{r_j}] = r_f + \beta_{j,w}'\lambda_w' + \beta_{j,y}'\lambda_y' E[rj~]=rf+βj,wλw+βj,yλy
这样,我们就把资产的期望回报率用一个两因子的模型给表示了出来。现在,决定资产期望回报率的,既有资产回报率与市场总回报率之间的相关性,也有资产回报率与工资收入之间的相关性。

从以上的推演,我们可以看到多因子模型背后的直觉。所谓因子,实际上是会影响随机折现因子的不确定性来源。

因子一定会贡献系统性风险(否则可以通过diversification来消除)

APT

在实践中,以资产定价为目标的投资者往往并不是很关心因子有何经济含义,而只是在乎能否找到对资产回报率有解释力的解释变量来提升自己定价的精度。所以,多因子模型的理论根基其实是 Ross 提出的套利资产定价理论(APT)。

APT的核心思想:所有资产的期望回报率都由一组共同的因子所决定时,基于无套利的思想,不同资产的期望回报率之间会有某种线性关系;

推导多因子模型

符号记法:

  • factor risk: 因子的系统性风险
  • loading: 因子前的系数
  • idiosyncratic risk: 与因子风险无关的剩余风险, ϵ i ~ \tilde{\epsilon_i} ϵi~

单因子两资产

单因子: f f f,将其标准化 E [ f ~ ] = 0 E[\tilde{f}]=0 E[f~]=0,构造两个资产
{ r i ~ = r i ˉ + β i f ~ r j ~ = r j ˉ + β j f ~ \begin{cases} \tilde{r_i} = \bar{r_i} + \beta_i \tilde{f} \\ \tilde{r_j} = \bar{r_j} + \beta_j \tilde{f} \\ \end{cases} {ri~=riˉ+βif~rj~=rjˉ+βjf~

构造资产组合p:
r p ~ = w r i ~ + ( 1 − w ) r j ~ = w ( r i ˉ + β i f ~ ) + ( 1 − w ) ( r j ˉ + β j f ~ ) = [ w r i ˉ + ( 1 − w ) r j ˉ ] + [ w β i + ( 1 − w ) β j ] f ~ \begin{aligned} \tilde{r_p} &= w \tilde{r_i} + (1-w) \tilde{r_j} \\ &= w(\bar{r_i} + \beta_i \tilde{f}) + (1-w) (\bar{r_j} + \beta_j \tilde{f}) \\ &= [w\bar{r_i} + (1-w) \bar{r_j}] + [w\beta_i + (1-w) \beta_j] \tilde{f} \\ \end{aligned} rp~=wri~+(1w)rj~=w(riˉ+βif~)+(1w)(rjˉ+βjf~)=[wriˉ+(1w)rjˉ]+[wβi+(1w)βj]f~
接下来分两步

  1. 找到 w 0 w_0 w0使得 w β i + ( 1 − w ) β j = 0 ⇒ w 0 = β j β j − β i w\beta_i + (1-w) \beta_j=0\Rightarrow w_0=\frac{\beta_j}{\beta_j-\beta_i} wβi+(1w)βj=0w0=βjβiβj
    r ~ p 0 = β j r i ˉ − β i r j ˉ β j − β i = r f ⇒ r i ˉ − r f β i = r j ˉ − r f β j \tilde{r}_{p0} = \frac{\beta_j \bar{r_i} - \beta_i \bar{r_j}}{\beta_j - \beta_i} = r_f\\ \Rightarrow \frac{\bar{r_i}-r_f}{\beta_i} = \frac{\bar{r_j}-r_f}{\beta_j} \\ r~p0=βjβiβjriˉβirjˉ=rfβiriˉrf=βjrjˉrf
    则对任意一个资产,
    λ ≜ r i ˉ − r f β i ∀ i r i ˉ = r f + β i λ \lambda \triangleq \frac{\bar{r_i}-r_f}{\beta_i} \qquad \forall i \\ \bar{r_i} = r_f + \beta_i \lambda λβiriˉrfiriˉ=rf+βiλ
    对于任意一个资产 i , j i,j i,j,
    λ = r i ˉ − r j ˉ β i − β j ( a b = c d = a + c b + d ) \lambda = \frac{\bar{r_i} - \bar{r_j}}{\beta_i - \beta_j} \qquad (\frac{a}{b} = \frac{c}{d} = \frac{a+c}{b+d}) λ=βiβjriˉrjˉ(ba=dc=b+da+c)
  2. 找到 w 1 w_1 w1,使得 w 1 β i + ( 1 − w 1 ) β j = 1 ⇒ w 1 = 1 − β j β i − β j w_1\beta_i + (1-w_1)\beta_j = 1\Rightarrow w_1=\frac{1-\beta_j}{\beta_i-\beta_j} w1βi+(1w1)βj=1w1=βiβj1βj
    r ~ p 1 = [ 1 − β j β i − β j r i ˉ + ( 1 − 1 − β j β i − β j ) r j ˉ ] + f ~ = β i r j ˉ − β j r i ˉ β i − β j + r i ˉ − r j ˉ β i − β j + f ~ = r f + λ + f ~ \begin{aligned} \tilde{r}_{p1} &= [\frac{1-\beta_j}{\beta_i-\beta_j}\bar{r_i}+(1-\frac{1-\beta_j}{\beta_i-\beta_j})\bar{r_j}] + \tilde{f}\\ &= \frac{\beta_i\bar{r_j}-\beta_j\bar{r_i}}{\beta_i-\beta_j} + \frac{\bar{r_i}-\bar{r_j}}{\beta_i-\beta_j} + \tilde{f}\\ &=r_f + \lambda + \tilde{f}\\ \end{aligned} r~p1=[βiβj1βjriˉ+(1βiβj1βj)rjˉ]+f~=βiβjβirjˉβjriˉ+βiβjriˉrjˉ+f~=rf+λ+f~
    两边同时取期望
    r ˉ p 1 = r f + λ ⇒ λ = r ˉ p 1 − r f ( 表示 f a c t o r p r e m i u m ) λ = r i ˉ − r f β i = r ˉ p 1 − r f ⇒ r i ˉ − r f = β i ( r ˉ p 1 − r f ) \begin{aligned} \bar{r}_{p1} &= r_f + \lambda \\ \Rightarrow \lambda &= \bar{r}_{p1} - r_f \quad (表示factor \quad premium) \\ \lambda &= \frac{\bar{r_i}-r_f}{\beta_i} = \bar{r}_{p1} - r_f \\ \Rightarrow \bar{r_i}-r_f &=\beta_i(\bar{r}_{p1} - r_f) \end{aligned} rˉp1λλriˉrf=rf+λ=rˉp1rf(表示factorpremium)=βiriˉrf=rˉp1rf=βi(rˉp1rf)

多因子多资产

假设有K个因子,N个资产
r i = r i ˉ + ∑ k = 1 K β i , k f k ~ + ϵ ~ i , i = 1 , … , N , N > K r_i = \bar{r_i} + \sum_{k=1}^K \beta_{i,k} \tilde{f_k} + \tilde{\epsilon}_i,i=1,\dots,N, N>K ri=riˉ+k=1Kβi,kfk~+ϵ~i,i=1,,N,N>K
假设 E [ f k ~ ] = 0 , E [ ϵ ~ i ] = 0 , E [ f k ~ 2 ] = 1 , E [ ϵ ~ i 2 ] = σ ϵ 2 < ∞ , E [ f ~ k f ~ k ′ ] = E [ ϵ ~ k ϵ ~ k ′ ] = E [ f ~ k ϵ ~ k ] = 0 E[\tilde{f_k}]=0,E[\tilde{\epsilon}_i]=0,E[\tilde{f_k}^2]=1,E[\tilde{\epsilon}_i^2]=\sigma_{\epsilon}^2<\infty,E[\tilde{f}_k\tilde{f}_{k'}]=E[\tilde{\epsilon}_k\tilde{\epsilon}_{k'}]=E[\tilde{f}_k\tilde{\epsilon}_{k}]=0 E[fk~]=0,E[ϵ~i]=0,E[fk~2]=1,E[ϵ~i2]=σϵ2<,E[f~kf~k]=E[ϵ~kϵ~k]=E[f~kϵ~k]=0,
r p ~ = ∑ i = 1 N w i r i ~ = ∑ i = 1 N w i r i ˉ + ∑ k = 1 K ( ∑ i = 1 N w i β i , k ) f k ~ + ∑ i = 1 N w i ϵ ~ i \tilde{r_p} = \sum_{i=1}^N w_i\tilde{r_i} = \sum_{i=1}^N w_i\bar{r_i} + \sum_{k=1}^K (\sum_{i=1}^Nw_i\beta_{i,k})\tilde{f_k} + \sum_{i=1}^N w_i \tilde{\epsilon}_i rp~=i=1Nwiri~=i=1Nwiriˉ+k=1K(i=1Nwiβi,k)fk~+i=1Nwiϵ~i

  1. 所有loading=0
    ( ∑ i = 1 N w i β i , k ) = 0 ∀ k , k = 1 , … , K (\sum_{i=1}^Nw_i\beta_{i,k}) = 0 \quad \forall k,k=1,\dots,K (i=1Nwiβi,k)=0k,k=1,,K
    [ β 1 , 1 … β N , 1 ⋮ ⋱ ⋮ β 1 , K … β N , K ] K × N [ w 1 ⋮ w N ] N × 1 = [ 0 ⋮ 0 ] K × 1 \begin{bmatrix} \beta_{1,1} & \dots & \beta_{N,1} \\ \vdots & \ddots & \vdots \\ \beta_{1,K} & \dots & \beta_{N,K} \\ \end{bmatrix}_{K\times N } \begin{bmatrix} w_1\\ \vdots\\ w_N \end{bmatrix}_{N\times 1} = \begin{bmatrix} 0\\ \vdots\\ 0 \end{bmatrix}_{K\times 1} β1,1β1,KβN,1βN,K K×N w1wN N×1= 00 K×1
    因为 N > K N>K N>K,且 E [ f ~ k f ~ k ′ ] = 0 E[\tilde{f}_k\tilde{f}_{k'}]=0 E[f~kf~k]=0,因子间相互独立,所以系数矩阵满秩,所以方程组有解; 接下来继续分析 r p 0 ~ \tilde{r_{p0}} rp0~
    r p 0 ~ = ∑ i = 1 N w i r i ˉ + ∑ i = 1 N w i ϵ ~ i ( w i 是各个 β 的函数 ) \tilde{r_{p0}} = \sum_{i=1}^N w_i\bar{r_i} + \sum_{i=1}^N w_i \tilde{\epsilon}_i \quad (w_i是各个\beta的函数) rp0~=i=1Nwiriˉ+i=1Nwiϵ~i(wi是各个β的函数)
    两边同时取方差
    σ 2 ( r p 0 ~ ) = ( ∑ i = 1 N w i 2 ) σ ϵ 2 \sigma^2(\tilde{r_{p0}}) = (\sum_{i=1}^N w_i^2)\sigma_{\epsilon}^2 σ2(rp0~)=(i=1Nwi2)σϵ2
    当资产的数量N很大的时候,每个权重就大概为 1 N \frac{1}{N} N1,因此, σ 2 ( r p ~ ) \sigma^2(\tilde{r_p}) σ2(rp~)的数量级就为
    ( 1 N ) 2 × N × σ ϵ 2 = σ ϵ 2 N (\frac{1}{N})^2 \times N \times \sigma_{\epsilon}^2 = \frac{\sigma_{\epsilon}^2}{N} (N1)2×N×σϵ2=Nσϵ2
    N → ∞ , σ 2 ( r p 0 ~ ) → 0 , r p 0 N\to \infty,\sigma^2(\tilde{r_{p0}}) \to 0,r_{p0} N,σ2(rp0~)0,rp0就是无风险组合
    r p 0 ~ ≈ ∑ i = 1 N w i r i ˉ = r f \tilde{r_{p0}} \approx \sum_{i=1}^N w_i\bar{r_i} = r_f rp0~i=1Nwiriˉ=rf
  2. 对K个因子,各个因子的loading=1进行构建,就能得到 λ k \lambda_k λk(表示第k个因子的因子溢价)

进而得到一般情况下的APT
r i ˉ = r f + ∑ k = 1 K β i , k λ k \bar{r_i} = r_f + \sum_{k=1}^K \beta_{i,k}\lambda_k riˉ=rf+k=1Kβi,kλk

APT的应用

α , β \alpha,\beta α,β分离

构造如下APT模型
r ~ 0 − r f = α 0 + ∑ n = 1 N β 0 , n r ~ n + ϵ ^ 0 \tilde{r}_0 - r_f = \alpha_0 + \sum_{n=1}^N \beta_{0,n}\tilde{r}_n + \hat{\epsilon}_0 r~0rf=α0+n=1Nβ0,nr~n+ϵ^0
r ~ 0 \tilde{r}_0 r~0表示0这个资产组合,可以用另一个资产组合的收益 ∑ n = 1 N β 0 , n r ~ n \sum_{n=1}^N \beta_{0,n}\tilde{r}_n n=1Nβ0,nr~n来表示;

那么只需要买入0,卖出 ∑ n = 1 N β 0 , n r ~ n \sum_{n=1}^N \beta_{0,n}\tilde{r}_n n=1Nβ0,nr~n这个资产组合就能对冲掉因子所形成的系统性风险,得到了一个 α 0 + ϵ ^ 0 \alpha_0 + \hat{\epsilon}_0 α0+ϵ^0

因子选股

在实践中,多因子模型经常被用来筛选投资标的,最常见的是用来选股。一个因子代表了一个对股票期望回报率有解释力的因素。如果这种解释力很强,那么用因子来给所有股票从好到坏排个序,买入排在前面的股票(卖出排在后面的股票),就应该能获得不错的回报。

统计套利

统计套利不等于无风险套利,统计套利,则是利用统计分析工具来找出相互联系的资产价格之间长期稳定的数量关系。当观测到现实价格数据大幅偏离这种长期稳定关系时,进行相应的操作来赌这种偏离会消失。所以统计套利不是无风险的。因为我们无法保证过去稳定的数量关系在未来不会破裂。一旦这种数量关系破裂了,那么下注来赌偏离会消失就会亏掉不少钱。

可以计算各只股票过去一段时间的实际回报率。如果有股票实际回报率超过期望回报率,就说明它前段时间股价涨得太快了,有理由预期它接下来一段时间的股价涨幅会慢一些。

你可能感兴趣的:(计量经济学,算法)