留一法这个R语言函数哪里有问题么？ 和我手动算的一直不一样。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

函数如下：

cv.lm <- function(obj) {
data <- data.frame(y = obj$model[,1], x = obj$model[, -1])
N <- nrow(data)
{
   ytest <- data[, 1]
   newrsq <- numeric(N)
  ypred <- numeric(N)
   for(i in 1:N){
     newtrain <- data[-i, ]
     xtest <- data[i, -1]
   newfm <- lm(y ~., data = newtrain)
    newrsq <- summary(newfm)$r.squared
   ypred <- predict(newfm, xtest)
      }
  }
q.squared <- 1 - (sum((ytest - ypred)^2) / sum((ytest - mean(ytest))^2))
SDEP <- sqrt(sum((ytest - ypred)^2) / N)
return(list(q.squared = q.squared, SDEP = SDEP, newrsq = newrsq))
}




q.squared我计算出来的一直是负数……

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：留一法 交叉验证