带着诚实不确定性的回归
普通回归只给出一个斜率,贝叶斯回归给出斜率的整个分布。设股票收益为 y、市场收益为 x,y=βx+e,噪声 e 的方差为 σ2。先给贝塔一个先验,比如以 m0=1 为中心、方差为 v0 的正态分布。数据把斜率推向最小二乘值,先验把它拉回来,这和前几个主题里按精度加权的平均是同一回事。先验精度是 1/v0,数据精度是 ∑x2/σ2,后者在市场波动大时很大。
Pn=v01+σ2∑x2,mn=Pnm0/v0+∑xy/σ2,sd(β)=Pn1(1) 有两点值得注意。第一,以 0 为中心的先验恰好给出岭回归:斜率为 ∑xy/(∑x2+λ),其中 λ=σ2/v0。所以岭回归就是先验不信任大斜率的贝叶斯回归。第二,对新的 x∗ 作预测时,离散度有两个来源:噪声 σ2,以及斜率的不确定性,它贡献 x∗2/Pn。后一项随 x∗ 的平方增长,所以远离数据的预测更不确定。
用数字看贝塔设
∑x2=40、
∑xy=48、
σ2=4。最小二乘斜率为
48/40=1.2。取先验
N(1,0.52),即
v0=0.25:精度为
4+40/4=14,均值
(4+12)/14=1.143,标准差
1/14=0.267。若用以
0 为中心的岭先验,斜率为
48/(40+16)=0.857。在
x∗=2 处的预测均值为
2.29,标准差为
4+4/14=2.07。
- 贝叶斯斜率:精度等于先验精度加数据精度的分布。
- 岭回归是先验以 0 为中心的情形,λ=σ2/v0。
- 预测包含噪声加斜率不确定性,后者随 x∗2 增长。
实验:评分内容
- 后验精度15 分
- 后验斜率15 分
- 斜率的不确定性10 分
- 岭回归斜率15 分
- 预测的离散度10 分
- 非常自信的先验10 分
- 配对术语10 分
- 远离数据时15 分