Python 投资组合风险管理入门
Dakota Wixom
Quantitative Analyst | QuantCourse.com
$$ R_{P} = $$
$$ RF + \beta_{M}(R_{M}-RF)+b_{SMB} \cdot SMB + b_{HML} \cdot HML + \alpha $$

假设对象 Data 已包含投资组合与市场的超额收益:
import statsmodels.formula.api as smf
model = smf.ols(formula='Port_Excess ~ Mkt_Excess + SMB + HML',
data=Data)
fit = model.fit()
adjusted_r_squared = fit.rsquared_adj
adjusted_r_squared
0.90
要提取 HML 的 p 值,假设工作区中已有回归结果对象 fit:
fit.pvalues["HML"]
0.0063
要检验其是否显著,只需查看它是否小于常用阈值 0.05:
fit.pvalues["HML"] < 0.05
True
要提取 HML 系数,假设工作区中已有回归结果对象 fit:
fit.params["HML"]
0.502
fit.params["SMB"]
-0.243
假设你已将回归结果保存在对象 fit 中:
portfolio_alpha = fit.params["Intercept"]
portfolio_alpha_annualized = ((1 + portfolio_alpha) ** 252) - 1
portfolio_alpha_annualized
0.045
Python 投资组合风险管理入门