R返回错误：下标超出范围答案

【问题标题】：R giving back error: subscript out of boundsR返回错误：下标超出范围
【发布时间】：2021-07-05 20:17:10
【问题描述】：

我一直在研究 R ISLR College 数据集，我希望在训练集上执行最佳子集选择，并绘制与每种大小的最佳模型相关联的训练集 MSE。

library(ISLR)
library(leaps)
data(College)
head(College)

#splitting the data into 70/30
subset<- sample(nrow(college)*0.7)
collegetrain<- college[subset,]
collegetest<-college[-subset,]

这是我的代码：

regfit.full <- regsubsets(apps ~ ., data = college.train, nvmax = 20)
train.mat <- model.matrix(apps ~ ., data = college.train, nvmax = 20)
val.errors <- rep(NA, 20)
for (i in 1:20) {
coefi <- coef(regfit.full, id = i)
pred <- train.mat[, names(coefi)] %*% coefi
val.errors[i] <- mean((pred - college.train$y)^2)
}
plot(val.errors, xlab = "Number of predictors", ylab = "Training MSE", pch = 19, type = "b")

数据集的结构如下： 777 个观察结果，训练集中有 543 个，测试集中有 234 个。有 18 个变量，其中 17 个是数字，1 个是是和否的因子（这不需要更改）。

我在运行代码时收到的错误消息是： s$which [id, , drop=FALSE] 出错：下标越界

【问题讨论】：

标签： r plot statistics mse

【解决方案1】：

regfit.full <- regsubsets(Apps ~ ., data = collegetrain, nvmax = 20)
train.mat <- model.matrix(Apps ~ ., data = collegetrain, nvmax = 20)

val.errors <- rep(NA, 20)
for (i in 1:17) {
  coefi <- coef(regfit.full, id = i)
  pred <- train.mat[, names(coefi)] %*% coefi
  val.errors[i] <- mean((pred - collegetrain$Apps)^2)
}

plot(val.errors, xlab = "Number of predictors", ylab = "Training MSE", 
     pch = 19, type = "b")

【讨论】：