MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   资料存档 (https://www.labfans.com/bbs/forumdisplay.php?f=72)
-   -   MATLAB中10倍SVM分类的示例 (https://www.labfans.com/bbs/showthread.php?t=23543)

poster 2019-12-10 20:41

MATLAB中10倍SVM分类的示例
 
我需要一个说明性的示例,展示如何对两类数据集进行10倍SVM分类。 MATLAB文档中只有一个示例,但并非十折。有人能帮我吗?



[B]回答:[/B]

这是一个完整的示例,使用了生物信息学工具箱中的以下功能: [URL="http://www.mathworks.com/help/stats/svmtrain.html"]SVMTRAIN[/URL] , [URL="http://www.mathworks.com/help/stats/svmclassify.html"]SVMCLASSIFY[/URL] , [URL="http://www.mathworks.com/help/bioinfo/ref/classperf.html"]CLASSPERF[/URL]和[URL="http://www.mathworks.com/help/bioinfo/ref/crossvalind.html"]CROSSVALIND[/URL] 。

load fisheriris %# load iris dataset groups = ismember(species,'setosa'); %# create a two-class problem %# number of cross-validation folds: %# If you have 50 samples, divide them into 10 groups of 5 samples each, %# then train with 9 groups (45 samples) and test with 1 group (5 samples). %# This is repeated ten times, with each group used exactly once as a test set. %# Finally the 10 results from the folds are averaged to produce a single %# performance estimation. k=10; cvFolds = crossvalind('Kfold', groups, k); %# get indices of 10-fold CV cp = classperf(groups); %# init performance tracker for i = 1:k %# for each fold testIdx = (cvFolds == i); %# get indices of test instances trainIdx = ~testIdx; %# get indices training instances %# train an SVM model over training instances svmModel = svmtrain(meas(trainIdx,:), groups(trainIdx), ... 'Autoscale',true, 'Showplot',false, 'Method','QP', ... 'BoxConstraint',2e-1, 'Kernel_Function','rbf', 'RBF_Sigma',1); %# test using test instances pred = svmclassify(svmModel, meas(testIdx,:), 'Showplot',false); %# evaluate and update performance object cp = classperf(cp, pred, testIdx); end %# get accuracy cp.CorrectRate %# get confusion matrix %# columns:actual, rows:predicted, last-row: unclassified instances cp.CountingMatrix 输出:

ans = 0.99333 ans = 100 1 0 49 0 0 我们获得了99.33%准确度,其中只有一个被错误分类为“非setosa”的“ setosa”实例

[B]更新[/B] :SVM功能已移至R2013a中的“统计”工具箱



[url=https://stackoverflow.com/questions/3070789]更多&回答...[/url]


所有时间均为北京时间。现在的时间是 01:04

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.