R语言实战回归部分交叉验证问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

学习R语言实战，在回归的交叉验证部分遇到了问题。
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")
}

states <- as.data.frame(state.x77[,c("Murder", "Population","Illiteracy", "Income", "Frost")])
fit <- lm(Murder ~ Population + Income + Illiteracy + Frost, data=states)
shrinkage(fit)

----------------------------------------------------------------
这里面主要是搞不懂在shrinkage函数中，x,y还没被定义出来呢就已经在theta.fit、theta.predict中用了，这样也可以吗？
此外cbind(1,x)%*%fit$coef这个语句应该是想计算预测的值，但是我不用函数的方式就是先回归再计算
states <- as.data.frame(state.x77[,c("Murder", "Population","Illiteracy", "Income", "Frost")])
fit <- lm(Murder ~ Population + Income + Illiteracy + Frost, data=states)
x <- fit$model[,2:ncol(fit$model)]
y <- fit$model[,1]
cbind(1,x)%*%fit$coef

就显示
> cbind(1,x)%*%fit$coef
Error in cbind(1, x) %*% fit$coef :
  requires numeric/complex matrix/vector arguments
好像cbind(1,x)莫名奇妙的就变成list了，为什按照原来的语句就能顺畅的执行了呢？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：R语言实战 交叉验证 R语言 shrinkage Age

诶好像dollar号都没打出来

states <- as.data.frame(state.x77[,c("Murder", "Population","Illiteracy", "Income", "Frost")])
fit <- lm(Murder ~ Population + Income + Illiteracy + Frost, data=states)
x <- fit$model[,2:ncol(fit$model)]
y <- fit$model[,1]
cbind(1,x)%*%fit$coef

代码贴上来有点问题，原代码在R in action 2的194页

同问啊。cbind(1,x)这个应该是个data.frame，不能够用%*%，  但是为什么源代码就能够顺利执行。  是因为function有什么特殊功能吗