再用R软件来进行。需要调用库MASS中的函数lm.ridge
且在进行之前需要对数据进行标准化处理。
d2=scale(d1)
d2=as.data.frame(d2)
> library(MASS)
> res3=lm.ridge(y~-1+.,data=d2,lambda=seq(0,1,0.05))
> res3
x1 x2 x3 x4 x5
0.00 2.4473858 -2.48509745 -0.08314028 0.5305377 0.5635372
0.05 0.2132936 -0.23113004 -0.11628955 0.6448621 0.4322741
0.10 0.1641667 -0.08529681 -0.11045247 0.5875122 0.3871112
0.15 0.1629117 -0.01491928 -0.10567036 0.5441832 0.3568321
0.20 0.1695221 0.02964903 -0.10192920 0.5113297 0.3349291
0.25 0.1772741 0.06116693 -0.09893708 0.4856791 0.3183414
0.30 0.1845911 0.08486816 -0.09647920 0.4651031 0.3053549
0.35 0.1911164 0.10342168 -0.09441044 0.4482176 0.2949227
0.40 0.1968339 0.11837013 -0.09263208 0.4340951 0.2863673
0.45 0.2018226 0.13068103 -0.09107542 0.4220934 0.2792302
0.50 0.2061802 0.14099738 -0.08969147 0.4117552 0.2731903
0.55 0.2099992 0.14976558 -0.08844448 0.4027457 0.2680156
0.60 0.2133599 0.15730628 -0.08730784 0.3948148 0.2635348
0.65 0.2163301 0.16385638 -0.08626137 0.3877714 0.2596184
0.70 0.2189665 0.16959492 -0.08528953 0.3814676 0.2561670
0.75 0.2213158 0.17465995 -0.08438014 0.3757866 0.2531029
0.80 0.2234174 0.17915966 -0.08352357 0.3706352 0.2503646
0.85 0.2253036 0.18318016 -0.08271210 0.3659381 0.2479027
0.90 0.2270020 0.18679080 -0.08193947 0.3616335 0.2456772
0.95 0.2285356 0.19004815 -0.08120056 0.3576707 0.2436554
1.00 0.2299238 0.19299879 -0.08049116 0.3540073 0.2418102
从以上结果可以看出，其计算结果与用SPSS计算不同，从第二行就开始不同，不知是何因。
打算手动用矩阵来进行计算，看结果有何变化没有？
待续！
扫码加好友，拉您进群
全部版块
我的主页
收藏
