R中的for循环错误“尝试选择少于一个元素”

Question

我已经为我的向量创建了正确数量的索引，并且我正在尝试从 for 循环中输入第 i 个元素作为索引来保存分类错误值。但我得到了错误：

indeces.gen_error[[i]] 中的错误 尝试选择少于一个元素

uspscl.txt

uspsdata.txt

我的代码：

library(e1071)
library(caret)

set.seed(733)
uspscldf = read.table('uspscl.txt', header=F, sep=',')
uspsdatadf = read.table('uspsdata.txt', header=F, sep='\t')

trainIndex <- createDataPartition(uspscldf$V1,list=FALSE, p = .80,times=1)
dataTrain <- uspsdatadf[ trainIndex,]
dataTest  <- uspsdatadf[-trainIndex,]

classTrain <- uspscldf[ trainIndex,]
classTest  <- uspscldf[-trainIndex,]

indeces = seq(0.00001, 1, by=0.001)

indeces.gen_error = NULL
indeces.softmargin = NULL
for (i in seq(0.00001, 1, by=0.001)){
    # For svm(): soft margin is "cost" 
    # Gaussian kernel bandwidth (sigma) = is implicitly defined by "gamma"
    # kernal=radial is non-linear while kernal=linear is linear
    svm.model <- svm(classTrain ~ ., data = dataTrain, cost = i,type="C-classification",kernal = "linear")

    svm.pred <- predict(svm.model, dataTrain)

    # confusion matrix
    tab <- table(pred = svm.pred, true = classTrain)

    classification_error <- 1- sum(svm.pred == classTrain)/length(svm.pred)

    indeces.gen_error[[i]] <- paste(classification_error)
    indeces.softmargin[[i]]<-i
}

我在第一次迭代中打印了第一个 i，它给出了 1e-5，这是正确的，所以我不知道为什么它说我选择的元素少于一个。 任何帮助，将不胜感激。谢谢

回答::: 在我自己解决答案之前，我没有看到皮埃尔对此的回答，但他的解释更好，所以我接受他的回答。我现在的新代码是：

indeces = seq(0.00001, 1, by=0.001)

indeces.gen_error = NULL
indeces.softmargin = NULL
count=0
for (i in indeces){
  count=count+1
    # For svm(): soft margin is "cost" 
    # Gaussian kernel bandwidth (sigma) = is implicitly defined by "gamma"
    # kernal=radial is non-linear while kernal=linear is linear
    svm.model <- svm(classTrain ~ ., data = dataTrain, cost = i,type="C-classification",kernal = "linear")

    svm.pred <- predict(svm.model, dataTrain)

    # confusion matrix
    tab <- table(pred = svm.pred, true = classTrain)

    classification_error <- 1- sum(svm.pred == classTrain)/length(svm.pred)

    indeces.gen_error[[count]] <- paste(classification_error)
    indeces.softmargin[[count]]<-i
}

Answer 1

#Example
x <- NULL
for( i in seq(0.01, 1, .01)) {
  a <- 10 * i
  x[[i]] <- paste("b", a)
}
# Error in x[[i]] <- paste("b", a) : 
#   attempt to select less than one element

#The right way
x <- NULL
myseq <- seq(0.01, 1, 0.01)
for( i in 1:length(myseq)) {
  a <- 10 * myseq[i]
  x[i] <- paste("b", a)
}

为什么第一种方法会失败for( i in seq(0.01, 1, .01)) 将使用序列作为“i”。每当循环失败时，第一种排除故障的方法是逐个尝试每个循环。因此，每个循环都会获取序列的一个值，并将其输入到有i 的任何地方。第一个循环看起来像：

for (i in seq(0.00001, 1, by=0.001)){

    svm.model <- svm(classTrain ~ ., data = dataTrain, cost = 0.00001,type="C-classification",kernal = "linear")

    svm.pred <- predict(svm.model, dataTrain)

    # confusion matrix
    tab <- table(pred = svm.pred, true = classTrain)

    classification_error <- 1- sum(svm.pred == classTrain)/length(svm.pred)

    indeces.gen_error[[0.00001]] <- paste(classification_error)
    indeces.softmargin[[0.00001]]<- 0.00001
}

你看到问题了吗？使用indeces.gen_error[[0.00001]]，请注意这里发生的事情。你不是故意的。您的意思是让indeces.gen_error[[1]] 成为第一个条目。

您正在使用小数进行子集化。如果我们有：

x <- 1:10

您认为x[2.5] 会发生什么？我们要求 R 获取位置 2.5 处的元素。那没有意义。没有半位置。有第2个或第3个。试试看返回什么。

在您的循环中，您向 R 询问 indeces.gen_error[[0.00001]]。因此，您请求的是第 1/100,000 个位置。那没有意义。评估器将强制子集为整数。它转到0。我们得到一个错误。