根据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,代表性消费者的优化问题
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~)(w0−c0)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(a−c)2,则 u ′ ( c ) = a − c u'(c) = a-c u′(c)=a−c,则随机折现因子可以表示为
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~)=δa−c0a−(1+rw~)(w0−c0)=A−Brw~A=δa−c0a−δa−c0w0−c0,B=δa−c0w0−c0
之前在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~]1−1=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} 1⇒1⇒E[rj~]=E[m~(1+rj~)]=E[m~](1+E[rj~])+cov(m~,r~j)=E[m~]1−1−E[m~]cov(m~,r~j)=rf−E[m~]cov(m~,r~j)=rf−E[m~]cov(A−Brw~,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,w≜var(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}]} λw≜E[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
现在我们假设消费者除了拥有初始财富 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+y0−c0)+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~=A−Brw~−Cy1~A=δa−c0a−δa−c0w0+y0−c0,B=δa−c0w0+y0−c0,C=δa−c0a
那么最终会得到
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来消除)
在实践中,以资产定价为目标的投资者往往并不是很关心因子有何经济含义,而只是在乎能否找到对资产回报率有解释力的解释变量来提升自己定价的精度。所以,多因子模型的理论根基其实是 Ross 提出的套利资产定价理论(APT)。
APT的核心思想:所有资产的期望回报率都由一组共同的因子所决定时,基于无套利的思想,不同资产的期望回报率之间会有某种线性关系;
符号记法:
单因子: 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~+(1−w)rj~=w(riˉ+βif~)+(1−w)(rjˉ+βjf~)=[wriˉ+(1−w)rjˉ]+[wβi+(1−w)βj]f~
接下来分两步
假设有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=1∑Kβ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=1∑Nwiri~=i=1∑Nwiriˉ+k=1∑K(i=1∑Nwiβi,k)fk~+i=1∑Nwiϵ~i
进而得到一般情况下的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=1∑Kβi,kλk
构造如下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~0−rf=α0+n=1∑Nβ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
在实践中,多因子模型经常被用来筛选投资标的,最常见的是用来选股。一个因子代表了一个对股票期望回报率有解释力的因素。如果这种解释力很强,那么用因子来给所有股票从好到坏排个序,买入排在前面的股票(卖出排在后面的股票),就应该能获得不错的回报。
统计套利不等于无风险套利,统计套利,则是利用统计分析工具来找出相互联系的资产价格之间长期稳定的数量关系。当观测到现实价格数据大幅偏离这种长期稳定关系时,进行相应的操作来赌这种偏离会消失。所以统计套利不是无风险的。因为我们无法保证过去稳定的数量关系在未来不会破裂。一旦这种数量关系破裂了,那么下注来赌偏离会消失就会亏掉不少钱。
可以计算各只股票过去一段时间的实际回报率。如果有股票实际回报率超过期望回报率,就说明它前段时间股价涨得太快了,有理由预期它接下来一段时间的股价涨幅会慢一些。