Визуализация правил исключения с использованием случайных данных.

Генерировать случайные x и y данные.

xdata = -3 + 6*rand(1,1e4);
ydata = -3 + 6*rand(1,1e4);

В качестве примера исключите данные, которые находятся внутри поля [-1 1 -1 1] или вне домена [-2 2].

outliers1 = ~excludedata(xdata,ydata,'box',[-1 1 -1 1]);
outliers2 = excludedata(xdata,ydata,'domain',[-2 2]);
outliers = outliers1|outliers2;

Постройте график данных, которые не исключаются. Белая область соответствует исключенным областям.

plot(xdata(~outliers),ydata(~outliers),'.')
axis([-3 3 -3 3])
axis square

Загрузите подсчет голосов и названия округов для штата Флорида с президентских выборов 2000 года в США.

load flvote2k

Используйте подсчет голосов для двух основных кандидатов от партии, Буша и Гора, в качестве предикторов подсчета голосов для стороннего кандидата Бьюкенена, и постройте график рассеяния:

plot(bush,buchanan,'rs')
hold on
plot(gore,buchanan,'bo')
legend('Bush data','Gore data')

Предположим, модель, в которой фиксированная доля избирателей Буша или Гора выбирает голосовать за Бьюкенена.

f = fittype({'x'})

f = 
     Linear model:
     f(a,x) = a*x

Исключить данные у открепительных избирателей, которые не воспользовались спорным бюллетенем «бабочка».

nobutterfly = strcmp(counties,'Absentee Ballots');

Выполните жесткую подгонку весов bisquare модели к двум наборам данных, исключая заочных избирателей.

bushfit = fit(bush,buchanan,f,'Exclude',nobutterfly,'Robust','on');
gorefit = fit(gore,buchanan,f,'Exclude',nobutterfly,'Robust','on');

Надежная посадка дает выходам низкий вес, поэтому большие остатки от надежной посадки могут быть использованы для идентификации выходов.

figure
plot(bushfit,bush,buchanan,'rs','residuals')
hold on
plot(gorefit,gore,buchanan,'bo','residuals')

Вычислите остатки.

bushres = buchanan - feval(bushfit,bush);
goreres = buchanan - feval(gorefit,gore);

Определите большие остатки как находящиеся вне диапазона [-500 500].

bushoutliers = excludedata(bush,bushres,'range',[-500 500]);
goreoutliers = excludedata(gore,goreres,'range',[-500 500]);

Отображение графств, соответствующих отклонениям. Округа Майами-Дейд и Броуард соответствуют наибольшим предикторным значениям. Округ Палм-Бич, единственный округ штата, использовавший бюллетень «бабочка», соответствует наибольшим остаточным значениям.

counties(bushoutliers)

ans = 2x1 cell
    {'Miami-Dade'}
    {'Palm Beach'}

counties(goreoutliers)

ans = 3x1 cell
    {'Broward'   }
    {'Miami-Dade'}
    {'Palm Beach'}