Question

在下面的代码中，我试图创建一个矩阵，列出每个城市的opt.lam。在运行循环时，前两个城市总是工作，然后我得到任何城市的错误。

这是我得到的错误。（coefmatrix运行正常，它只是产生此错误的lambdamatrix）。

[<-中的错误（*tmp* ,, i，值= c（0.577199381062121,0.577199381062121，：下标超出范围

这是我的代码：

lambdamatrix <- matrix(nrow=n,ncol=2)
rownames(lambdamatrix) <- cityIDs
colnames(lambdamatrix) <- c("lambda.min","lambda.1se")
for (i in 1:n) {
  data <- subset(simdata, city==cityIDs[i])
  x <- as.matrix(data.frame(data[,3:24]))
cvfit <- cv.glmnet(x, data$Y, family="poisson", offset=log(data$population))
opt.lam <- c(cvfit$lambda.min, cvfit$lambda.1se)
fit <- glmnet(x, data$Y, family= "poisson", offset=log(data$population))
abline(plot(fit, "lambda", label= TRUE, 
            main = cityIDs[i]), v=log(opt.lam), lty=2, lwd=3, 
                                    col=c("red","dark green"))
coefmatrix[,i] <- coef(fit, s=opt.lam[1])[1:23]
lambdamatrix[,i] <- c(cvfit$lambda.min, cvfit$lambda.1se)[1:n]
}`

Answer 1

在[,i]中，i是列索引器（而[i,]将是行索引器）。

由于您将lambdamatrix定义为matrix(nrow = n, ncol = 2)，一旦超过i=2，您就会要求不存在的列。

`[＆lt; -`（`* tmp *`,,下标超出边界的下标超出范围）出错

1 个答案: