6

我知道在 MatLab 中这很容易（'-v 10'）。

但我需要在 R 中执行此操作。我确实找到了一条关于添加cross = 10as 参数的评论。但这在帮助文件中没有得到证实，所以我对此持怀疑态度。

svm(Outcome ~. , data= source, cost = 100, gamma =1, cross=10)

任何成功的 R 支持向量机脚本示例也将不胜感激，因为我仍然遇到一些死胡同？

编辑：我忘了在标签之外提到我为此使用 libsvm 包。

4

3 回答 3

7

我也在尝试进行 10 折交叉验证。我认为使用 tune 不是执行它的正确方法，因为此功能用于优化参数，而不是用于训练和测试模型。

我有以下代码来执行 Leave-One-Out 交叉验证。假设该数据集是一个 data.frame，其中存储了您的数据。在每个 LOO 步骤中，添加观察到的与预测的矩阵，因此最后，结果包含全局观察到的与预测的矩阵。

#LOOValidation
for (i in 1:length(dataset)){
    fit = svm(classes ~ ., data=dataset[-i,], type='C-classification', kernel='linear')
    pred = predict(fit, dataset[i,])
    result <- result + table(true=dataset[i,]$classes, pred=pred);
}
classAgreement(result)

所以为了进行 10 折交叉验证，我想我们应该手动划分数据集，并使用折来训练和测试模型。

for (i in 1:10)
    train <- getFoldTrainSet(dataset, i)
    test <- getFoldTestSet(dataset,i)
    fit = svm(classes ~ ., train, type='C-classification', kernel='linear')
    pred = predict(fit, test)
    results <- c(results,table(true=test$classes, pred=pred));

}
# compute mean accuracies and kappas ussing results, which store the result of each fold

我希望这对你有帮助。

于 2012-11-13T16:43:14.370 回答
1

这是一个不使用包创建 10 个测试和训练折叠的简单方法：

#Randomly shuffle the data
yourData<-yourData[sample(nrow(yourData)),]

#Create 10 equally size folds
folds <- cut(seq(1,nrow(yourData)),breaks=10,labels=FALSE)

#Perform 10 fold cross validation
for(i in 1:10){
    #Segement your data by fold using the which() function 
    testIndexes <- which(folds==i,arr.ind=TRUE)
    testData <- yourData[testIndexes, ]
    trainData <- yourData[-testIndexes, ]
    #Use test and train data howeever you desire...
}
于 2014-07-04T20:18:35.353 回答
0

这是我在 cvsegments 的帮助下运行 k 折叠交叉验证以生成索引折叠的通用代码。

# k fold-cross validation
set.seed(1)
k <- 80;
result <- 0;
library('pls');
folds <- cvsegments(nrow(imDF), k);
for (fold in 1:k){
    currentFold <- folds[fold][[1]];
    fit = svm(classes ~ ., data=imDF[-currentFold,], type='C-classification', kernel='linear')
    pred = predict(fit, imDF[currentFold,])
    result <- result + table(true=imDF[currentFold,]$classes, pred=pred);   
}
classAgreement(result)
于 2012-11-14T18:06:04.607 回答