回归分析中的岭回归过程

标签:
sas岭回归 |
分类: SAS笔记 |
利用Longley数据做岭回归方程,数据如下所示:
Y
60323
61122
60171
61187
63221
63639
64989
63761
66019
67857
68169
66513
68655
69564
69331
70551
假设数据集为a,则首先对数据进行标准化:
proc standard data=a out=b mean=0 std=1;
run;
然后进行岭回归建模:
proc reg data=b
outest=c;
model y=x1-x6/ridge=0.0 to 1.0 by 0.1;
plot r.*np./ridgeplot;
run;
可以得到如下岭迹图:
http://s8/mw690/0025UTvnzy6JMs14Cure7&690
从上图可以看出当k=0.1时,曲线趋于稳定趋势。因此选取k=0.1,利用SAS建立岭回归方程,如下所示得到各类岭回归系数:
http://s1/mw690/0025UTvnzy6JMsaEE7ec0&690
因此选取k=0.1时的岭回归系数,得到的标准化岭回归方程为:
http://s3/mw690/0025UTvnzy6JMsfCJ5Ed2&690