我正在处理糖尿病数据集。我正在努力将从我之前的任务中获得的"最佳"模型拟合到训练集上,并报告该模型的测试错误。代码低于
library(nclSLR)
library(leaps)
library(glmnet)
data(diabetes)
set.seed(1)
(n = nrow(diabetes))
(p = ncol(diabetes) - 1)
train = diabetes [1:350,]
test = diabetes [351:442,]
grid=10^seq(10,-2,length=100)
best.subset=c("sex", "bmi", "map", "tc", "ldl", "ltg")
diabetes = diabetes[,best.subset]
#Train lm model
lm.mod=lm(sex~., diabetes[train,])
我最后一行写错了。显示
Error in xj[i] : invalid subscript type 'list'.
我试过取消登录((,但仍然没有成功。
问题是我们使用了i
,即行索引作为数据帧。这里,我们需要已经创建的数据帧的子集,即"训练">
lm.mod=lm(sex~., train)