我也是新手。。不明白为什么你的不对，我刚才对着你的一遍，没有找出来，但是我的可以运行，给你看看我的代码吧。


shrinkage<-function(fit,k=10){
require(bootstrap)
theta.fit<-function(x,y){lsfit(x,y)}
theta.predict<-function(fit,x){cbind(1,x)%*%fit$coef}
x<-fit$model[,2:ncol(fit$model)]
y<-fit$model[,1]
results<-crossval(x,y,theta.fit,theta.predict,ngroup=k)
r2<-cor(y,fit$fitted.values)^2
r2cv<-cor(y,results$cv.fit)^2
cat("Original R-square=", r2,"\n")
cat(k,"Fold Cross-Validated R-square =",r2cv,"\n")
cat("Change =",r2-r2cv,"\n")
}
fit<-lm(Murder ~ Population + Income + Illiteracy +Frost,
data=states)
shrinkage(fit)

谢谢谢谢！！