- 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$coef和fit$model是什么意思(在help中查出coef()是求模型系数的函数,至于model不知道是什么意思)?
先谢谢各位了


雷达卡






京公网安备 11010802022788号







