「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])