##example 7.4
dat<-read.csv(file="D:/xt10.3.csv",sep=",",header=TRUE)
dat
M<-dat[,2:ncol(dat)]
M<-as.data.frame(as.matrix(M))
colnames(M)<-c("y",paste("x",1:2,sep=""))
traceback()
library(MASS)
lmr<lm.ridge(y~x1+x2,X,lambda=seq(0,10,0.5),model=T)************************************
traceback()
matplot(lmr$lambda,t(lmr$coef),xlab="k",
ylab="estimated cofficients beta_hat(k)",type="l",
lty=l,col=1:5,main="Ridge Trace")
abline(0,0,lty=2)
legend(x=8,y=-0.05,legend=c("x1","x2"),
col=1:5,lty=1,cex=0.75)
#ridge regression in a smaller window
lmr<-lm.ridge(y~x1+x2+x3+x4+x5,X,lambda=seq(3,6,0.2),model=T)
matplot(lmr$lambda,t(lmr$coef),xlab="k",
ylab="estimated cofficients beta_hat(k)",type="p",pch=1:5,lty=0,
col=1:5,main="Ridge Trace")
lmr$coef[,'4.0']/p
你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答
本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。
因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。