代码实现如下：
data3.3<-read.csv("C:/Users/Administrator/Desktop/data3.3.csv",head=TRUE)
datas<-data.frame(scale(data3.3[,1:6])) # 对样本数据进行标准化处理并转换为数据框的存储格式
library(MASS)
ridge3.3<-lm.ridge(y~.-1,datas,lambda=seq(0,3,0.1))
# 做岭回归，对于标准化后的数据模型不包含截距项，其中lambda为岭参数k的所有取值
beta<-coef(ridge3.3) # 将所有不同岭参数所对应的回归系数的结果赋给beta
beta
# 绘制岭迹图
k<-ridge3.3$lambda #将所有岭参数赋给k
plot(k,k,type="n",xlab="岭参数k",ylab="岭回归系数",ylim=c(-2.5,2.5))
# 创建没有任何点和线的图形区域
linetype<-c(1:5)
char<-c(18:22)
for(i in 1:5)lines(k,beta[,i],type="o",lty=linetype[i],pch=char[i],cex=0.75)# 画岭迹
legend(locator(1),inset=0.5,legend=c("x1","x2","x3","x4","x5"),cex=0.8,pch=char,lty=linetype)
# 添加图例
输出结果为：
由上图可以看出，变量 x 2 x_{2} x2的岭回归系数 β 2 ^ ( k ) \hat{\beta_{2}}(k) β2^(k)从负值迅速变为正值， ∣ β 1 ^ ( k ) ∣ \left | \hat{\beta_{1}}(k) \right | ∣∣∣β1^(k)∣∣∣和 ∣ β 2 ^ ( k ) ∣ \left | \hat{\beta_{2}}(k) \right | ∣∣∣β2^(k)∣∣∣都迅速减少，两者之和比较稳定。从岭回归的角度来看， x 1 x_{1} x1和 x 2 x_{2} x2只要保留一个就可以了， x 3 , x 4 , x 5 x3,x4,x5 x3,x4,x5的岭回归系数u相对稳定。
通过上面分析，我们剔除 x 1 x_{1} x1后，在进行岭回归。
做出下面修改：
ridge3.3<-lm.ridge(y~.-x1-1,datas,lambda=seq(0,2,0.2))
plot(k,k,type="n",xlab="岭参数k",ylab="岭回归系数",ylim=c(-1,1))
linetype<-c(1:4)
char<-c(18:21)
for(i in 1:4)lines(k,beta[,i],type="o",lty=linetype[i],pch=char[i],cex=0.75)
legend(locator(1),inset=0.5,legend=c("x2","x3","x4","x5"),cex=0.8,pch=char,lty=linetype)
输出结果如下：
剔除 x 1 x_{1} x1后岭回归系数变化幅度减小，从图可以看出，岭参数 k k k大于1.4时，林参数的取值基本稳定，不妨选择 k = 1.4 k=1.4 k=1.4，由此时的回归系数得到样本数据标准化后的岭回归方程。
