財(cái)務(wù)的自我修養(yǎng)，微信公眾號(hào):機(jī)智出品(jizhjchupin)
一俄删、系統(tǒng)環(huán)境：
Windows7+Python3.6+Jupyter Notebook
二延柠、目標(biāo)：
以近兩年中國(guó)石化(600028)的收益率作為CAPM模型中的Ri嗤栓，上證指數(shù)(000001)收益率作為Rm佑刷，計(jì)算中國(guó)石化Beta、Alpha值钝计。2016年一年期國(guó)債利率為2.13%恋博。
三、過(guò)程：
step 1：需要的庫(kù)
import pandas as pd
import tushare as ts
import matplotlib.pyplot as plt
import statsmodels.api as sm
step 2：獲取中國(guó)石化和上證指數(shù)數(shù)據(jù)并處理
#tushare獲取市場(chǎng)數(shù)據(jù)
mktcd=ts.get_hist_data('sh',start='2015-10-01',end='2017-10-01')
mktret=pd.Series(mktcd.p_change.values,index=pd.to_datetime(mktcd.index))
mktret.name='mktret'
snp=ts.get_hist_data('600028',start='2015-10-01',end='2017-10-01')
snp.index=pd.to_datetime(snp.index)
snpret=pd.Series(snp.p_change.values,index=pd.to_datetime(snp.index))
snpret.name='snpret'
step 3：合并收益率并計(jì)算風(fēng)險(xiǎn)溢酬
ret=pd.merge(pd.DataFrame(mktret),
pd.DataFrame(snpret),
left_index=True,right_index=True,how='inner')
#計(jì)算無(wú)風(fēng)險(xiǎn)收益率
rf=1.0213**(1/360)-1
#計(jì)算股票超額收益率和市場(chǎng)風(fēng)險(xiǎn)溢酬
Eret=ret-rf
step 4：繪制兩者收益率散點(diǎn)圖
plt.scatter(Eret.values[:,0],Eret.values[:,1])
plt.title('SNP return and market return')
plt.show()
SNP return and market return.png
step 4：擬合CAPM模型
OLS()函數(shù)進(jìn)行擬合
model=sm.OLS(Eret.snpret[1:],sm.add_constant(Eret.mktret[1:]))
result=model.fit()
result.summary()
OLS Regression Results
| Variable: | snpret | squared: | 0.468 |
|---|---|---|---|
| Model: | OLS | Adj. R-squared: | 0.467 |
| Method: | Least Squares | F-statistic: | 427.8 |
| Date: | Mon, 02 Oct 2017 | Prob (F-statistic): | 1.22e-68 |
| Time: | 17:08:22 | Log-Likelihood: | -645.24 |
| No. Observations: | 488 | AIC: | 1294. |
| Df Residuals: | 486 | BIC: | 1303. |
| Df Model: | 1 | ||
| No. Observations: | 488 | ||
| Df Residuals: | 486 |
| coef | std err | t | P>|t| | [0.025 | 0.975] | |
|---|---|---|---|---|---|---|
| const | 0.0295 | 0.041 | 0.715 | 0.475 | -0.051 | 0.110 |
| mktret | 0.7018 | 0.034 | 20.684 | 0.000 | 0.635 | 0.768 |
| Omnibus: | 83.061 | Durbin-Watson: | 1.907 |
|---|---|---|---|
| Prob(Omnibus): | 0.000 | Jarque-Bera (JB): | 186.530 |
| Skew: | 0.897 | Prob(JB): | 3.13e-41 |
| Kurtosis: | 5.440 | Cond. No. | 1.22 |
四私恬、結(jié)果：
根據(jù)OLS擬合結(jié)果债沮，近兩年中國(guó)石化與上證指數(shù)的關(guān)系為：
Ri-Rf=0.0295+0.7018(Rm-Rf)+ε
可根據(jù)該模型預(yù)估中國(guó)石化預(yù)期收益率。
機(jī)智出品.jpg
