Обучите классификатор дерева классификации, а затем выполните перекрестную проверку с помощью пользовательской k-кратной функции потерь.

Загрузите набор данных радужки Фишера.

load fisheriris

Подготовка классификатора дерева классификации.

Mdl = fitctree(meas,species);

Mdl является ClassificationTree модель.

Перекрестная проверка Mdl с использованием 10-кратной перекрестной проверки по умолчанию. Вычислите ошибку классификации (долю неправильно классифицированных наблюдений) для кратных проверке наблюдений.

rng(1); % For reproducibility
CVMdl = crossval(Mdl);
L = kfoldLoss(CVMdl,'LossFun','classiferror')

L = 0.0467

Изучите результат, когда стоимость неправильной классификации цветка как versicolor является 10, и стоимость любой другой ошибочной классификации 1. Создание пользовательской функции noversicolor (показано в конце этого примера). Эта функция присваивает стоимость 10 за неправильную классификацию цветка как versicolor, и стоимость 1 для любой другой неправильной классификации.

Вычислите среднюю ошибку неправильной классификации с помощью noversicolor стоимость.

mean(kfoldfun(CVMdl,@noversicolor))

ans = 0.2267

Этот код создает функцию noversicolor.

function averageCost = noversicolor(CMP,~,~,~,Xtest,Ytest,~)
% noversicolor Example custom cross-validation function
%    Attributes a cost of 10 for misclassifying versicolor irises, and 1 for
%    the other irises.  This example function requires the fisheriris data
%    set.
Ypredict = predict(CMP,Xtest);
misclassified = not(strcmp(Ypredict,Ytest)); % Different result
classifiedAsVersicolor = strcmp(Ypredict,'versicolor'); % Index of bad decisions
cost = sum(misclassified) + ...
    9*sum(misclassified & classifiedAsVersicolor); % Total differences
averageCost = cost/numel(Ytest); % Average error
end