11

我知道 LIBSVM 在涉及多类 SVM 时只允许一对一分类。但是,我想对其进行一些调整以执行一对一的分类。我试图在下面进行一对一。这是正确的方法吗?

编码:

TrainLabel;TrainVec;TestVec;TestLaBel;
u=unique(TrainLabel);
N=length(u);
if(N>2)
    itr=1;
    classes=0;
    while((classes~=1)&&(itr<=length(u)))
        c1=(TrainLabel==u(itr));
        newClass=c1;
        model = svmtrain(TrainLabel, TrainVec, '-c 1 -g 0.00154'); 
        [predict_label, accuracy, dec_values] = svmpredict(TestLabel, TestVec, model);
        itr=itr+1;
    end
itr=itr-1;
end

我可能犯了一些错误。我想听听一些反馈。谢谢。

第二部分:正如grapeot所说:我需要做Sum-pooling(或投票作为一种简化的解决方案)才能得出最终答案。我不知道该怎么做。我需要一些帮助;我看到了python文件,但仍然不太确定。我需要帮助。

4

3 回答 3

10
%# Fisher Iris dataset
load fisheriris
[~,~,labels] = unique(species);   %# labels: 1/2/3
data = zscore(meas);              %# scale features
numInst = size(data,1);
numLabels = max(labels);

%# split training/testing
idx = randperm(numInst);
numTrain = 100; numTest = numInst - numTrain;
trainData = data(idx(1:numTrain),:);  testData = data(idx(numTrain+1:end),:);
trainLabel = labels(idx(1:numTrain)); testLabel = labels(idx(numTrain+1:end));
%# train one-against-all models
model = cell(numLabels,1);
for k=1:numLabels
    model{k} = svmtrain(double(trainLabel==k), trainData, '-c 1 -g 0.2 -b 1');
end

%# get probability estimates of test instances using each model
prob = zeros(numTest,numLabels);
for k=1:numLabels
    [~,~,p] = svmpredict(double(testLabel==k), testData, model{k}, '-b 1');
    prob(:,k) = p(:,model{k}.Label==1);    %# probability of class==k
end

%# predict the class with the highest probability
[~,pred] = max(prob,[],2);
acc = sum(pred == testLabel) ./ numel(testLabel)    %# accuracy
C = confusionmat(testLabel, pred)                   %# confusion matrix
于 2012-02-04T07:33:05.707 回答
4

从代码中我可以看到您正在尝试首先将标签转换为“某些类”与“非此类”,然后调用 LibSVM 进行训练和测试。一些问题和建议:

  1. 为什么要使用原版TrainingLabel进行培训?在我看来,应该是model = svmtrain(newClass, TrainVec, '-c 1 -g 0.00154');吗?
  2. 使用修改后的训练机制,您还需要调整预测部分,例如使用 sum-pooling 来确定最终标签。在 LibSVM 中使用-bswitch 来启用概率输出也将提高准确性。
于 2012-01-21T13:48:18.847 回答
1

除了概率估计,您还可以使用如下决策值

[~,~,d] = svmpredict(double(testLabel==k), testData, model{k});
prob(:,k) = d * (2 * model{i}.Label(1) - 1);

达到同样的目的。

于 2012-08-24T03:33:27.340 回答