“For”循环输出问题在R

Question

我在输出分配问题在“for”循环中挣扎：我生成一个空矩阵（“结果”）与适当的（不是???）尺寸在哪里循环应该从功能分配输出（PLS），但我得到一个错误...这里是代码：

library(pls) 
set.seed(10000) 
mat <- replicate(100, rnorm(100)) 
y <- as.matrix(mat[,1], drop=F) 
x <- mat[,2:100] 
min.nc <- 3 
max.nc <- 8 
nc <- seq(min.nc, max.nc) 
results <- matrix(NA, length(nc), 1) 
for(q in nc) { 
    pls <- plsr(y ~ x, ncomp = q) 
    results[q,] <- RMSEP(pls, intercept=F, estimate="train")$val[1,1,q] 
} 

运行后，该错误消息“for”循环说：

Error in `[<-`(`*tmp*`, q, , value = 0.329631604369802) : 
    subscript out of bounds 

......我不明白，因为迭代器“nc”具有“length（nc）”的大小。有任何想法吗？

非常感谢您的时间和帮助！

Chega

Answer 1

你q为3：8，你的矩阵results为1：6。因此，当它试图分配results[q,]时q > 6你得到这个错误。

> results[7,] 
Error in results[7, ] : subscript out of bounds 

的q的值是nc而非seq_along(nc)每个值。

for (q in nc) print(q) 

对战

for (q in seq_along(nc)) print(nc[q])