Бустинг, AdaBoost — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
(Пример на языке Scala)
(Пример на языке Scala)
Строка 114: Строка 114:
 
   libraryDependencies '''+=''' "com.github.haifengl" '''%%''' "smile-scala" '''%''' "1.5.2"
 
   libraryDependencies '''+=''' "com.github.haifengl" '''%%''' "smile-scala" '''%''' "1.5.2"
 
Пример классификации датасета и вычисления F1 меры<ref>[https://en.wikipedia.org/wiki/F1_score F1 мера]</ref> используя smile.classification.adaboost<ref>[https://haifengl.github.io/smile/classification.html#adaboost Smile, AdaBoost]</ref>:
 
Пример классификации датасета и вычисления F1 меры<ref>[https://en.wikipedia.org/wiki/F1_score F1 мера]</ref> используя smile.classification.adaboost<ref>[https://haifengl.github.io/smile/classification.html#adaboost Smile, AdaBoost]</ref>:
   import smile.classification._
+
   '''import '''smile.classification._
   import smile.data._
+
   '''import '''smile.data._
   import smile.plot._
+
   '''import '''smile.plot._
   import smile.read
+
   '''import '''smile.read
   import smile.validation.FMeasure
+
   '''import '''smile.validation.FMeasure
  
   val iris: AttributeDataset = read.table("iris.csv", delimiter = ",", response = Some((new NumericAttribute("class"), 2)))
+
   '''val '''iris: AttributeDataset = read.table("iris.csv", delimiter = ",", response = Some(('''new '''NumericAttribute("class"), 2)))
   val x: Array[Array[Double]] = iris.x()
+
   '''val '''x: Array[Array['''Double''']] = iris.x()
   val y: Array[Int] = iris.y().map(_.toInt)
+
   '''val '''y: Array['''Int'''] = iris.y().map(_.toInt)
   val ada: AdaBoost = adaboost(x, y, ntrees = 500, maxNodes = 2)
+
   '''val '''ada: AdaBoost = adaboost(x, y, ntrees = 500, maxNodes = 2)
   val predictions: Array[Int] = x.map(ada.predict)
+
   '''val '''predictions: Array['''Int'''] = x.map(ada.predict)
   val f1Score = new FMeasure().measure(predictions, y)
+
   '''val '''f1Score = '''new '''FMeasure().measure(predictions, y)
 
   plot(x, y, ada)
 
   plot(x, y, ada)
  

Версия 11:17, 22 января 2019

Описание

Бустинг (англ. boosting) — это композиционный мета-алгоритм обучения машин[на 18.01.19 не создан], не использующий параллельное обучение базовых классификаторов как бэггинг (англ. bagging). Бустинг также схож со стэкингом (англ. stacking), но стэкинг комбинирует результаты различных алгоритмов, получая тем самым более точный ответ. Основной идеей бустинга является комбинирование слабых функций, которые строятся в ходе итеративного процесса, где на каждом шаге новая модель обучается с использованием данных об ошибках предыдущих. В отличие от слабого алгоритма, сильный обучающий алгоритм является классификатором, хорошо коррелирующим с верной классификацией.

Одним из недостатков бустинга является то, что он может приводить к построению громоздких композиций, состоящих из сотен алгоритмов. Такие композиции исключают возможность содержательной интерпретации, требуют больших объёмов памяти для хранения базовых алгоритмов и существенных затрат времени на вычисление классификаций.

Алгоритмы бустинга

Большинство алгоритмов бустинга состоит из итеративного обучения слабых классификаторов с целью сборки их в сильный классификатор. Когда они добавляются, им обычно приписываются некоторым образом веса, которые, обычно, связаны с точностью обучения[на 18.01.19 не создан]. После того, как слабый классификатор добавлен, веса пересчитываются, что известно как «пересчёт весовых коэффициентов». Неверно классифицированные входные данные получают больший вес, а правильно классифицированные экземпляры теряют вес. Тем самым последующее слабое обучение фокусируется больше на примерах, где предыдущие слабые обучения дали ошибочную классификацию.

Основное расхождение между многими алгоритмами бустинга заключается в методах определения весовых коэффициентов точек тренировочных данных[на 18.01.19 не создан] и гипотез. Первым алгоритмом, который смог адаптироваться к слабому обучению был AdaBoost (сокр. Adaptive Boosting), предложенный Шапире и Фройндом.

Алгоритмы бустинга могут основываться на выпуклых или невыпуклых алгоритмах оптимизации. Выпуклые алгоритмы, такие как AdaBoost и LogitBoost[1], могут «потерпеть крушение» из-за случайного шума, так как не могут обучить базовым и поддающимся научению комбинациям слабых гипотез. Алгоритмы бустинга, основанные на невыпуклой оптимизации, такие как BrownBoost[2], могут быть обучены из данных с шумами и лежащий в основе классификатор Лонг–Серведио[3] для набора данных может быть обучен.

Прикладное использование алгоритмов бустинга

Задача классификации объектов

Если даны изображения, содержащие различные известные в мире объекты, классификатор может быть обучен на основе них для автоматической классификации объектов в будущих неизвестных изображениях. Простые классификаторы, построенные на основе некоторых признаков изображения объекта, обычно оказываются малоэффективными в классификации. Использование методов бустинга для классификации объектов — это путь объединения слабых классификаторов специальным образом для улучшения общей возможности классификации.

Классификация признаков является типичной задачей компьютерного зрения, где определяется, содержит ли изображение некоторую категорию объектов или нет. Идея тесно связана с распознаванием, идентификацией и обнаружением. Классификация по обнаружению объекта обычно содержит выделение признаков[на 18.01.19 не создан], обучение классификатора и применение классификатора к новым данным. Есть много способов представления категории объектов, например по анализу формы, с помощью модели «мешок слов», с помощью локальных описателей, таких как SIFT[4], и так далее. Примерами классификаторов с учителем служат наивные байесовские классификаторы[на 18.01.19 не создан], методы опорных векторов[на 18.01.19 не создан], смесь гауссиан и нейронные сети. Однако исследования показали, что категории объектов и их положение в изображениях могут быть обнаружены также с помощью обучения без учителя.

Задача ранжирования выдачи поисковых систем

Благодаря AdaBoost в мире появился градиентный бустинг (англ. gradient boosting) или GBM. Задачу ранжирования выдачи поисковых запросов рассмотрели с точки зрения функции потерь, которая штрафует за ошибки в порядке выдачи, поэтому было удобно внедрить GBM в ранжирование.

AdaBoost

Описание

Алгоритм может использоваться в сочетании с несколькими алгоритмами классификации для улучшения их эффективности. Алгоритм усиливает классификаторы, объединяя их в «комитет». AdaBoost является адаптивным в том смысле, что каждый следующий комитет классификаторов строится по объектам, неверно классифицированным предыдущими комитетами. AdaBoost чувствителен к шуму в данных и выбросам. Однако он менее подвержен переобучению по сравнению с другими алгоритмами машинного обучения.

AdaBoost вызывает слабые классификаторы в цикле [math]t = 1,...,T[/math]. После каждого вызова обновляется распределение весов [math]D_t[/math], которые отвечают важности каждого из объектов обучающего множества для классификации. На каждой итерации веса каждого неверно классифицированного объекта возрастают, таким образом новый комитет классификаторов «фокусирует своё внимание» на этих объектах.

Достоинства и недостатки

Достоинства:

  1. Простота реализации
  2. Хорошая обобщающая способность. В реальных задачах удаётся строить композиции, превосходящие по качеству базовые алгоритмы. Обобщающая способность может улучшаться по мере увеличения числа базовых алгоритмов.
  3. Время построения композиции практически полностью определяется временем обучения базовых алгоритмов.
  4. Возможность идентифицировать выбросы. Это наиболее «трудные» объекты $x_i$, для которых в процессе наращивания композиции веса $w_i$ принимают наибольшие значения.

Недостатки:

  1. Склонен к переобучению при наличии значительного уровня шума в данных.
  2. Требует достаточно длинных обучающих выборок. Другие методы линейной коррекции, в частности, бэггинг, способны строить алгоритмы сопоставимого качества по меньшим выборкам данных.


Алгоритм для задачи построения двоичного классификатора

Пакет AdaBoost может быть использован для распознавания лиц как пример двоичной классификации. Две категории — это лица и фон. Общий алгоритм выглядит следующим образом:

Дано: [math](x_1,y_1),...,(x_m,y_m)[/math], где [math]x_i \in X, y_i \in Y = \{-1,+1\}[/math]
Инициализируем [math]D_1(i) = \frac{1}{m},i=1,...,m[/math].
Для каждого [math]t=1,...,T[/math] пока не выполнен критерий останова:
  1. Находим классификатор [math]h_t:X\to \{-1,+1\}[/math] который минимизирует взвешенную ошибку классификации: [math]h_t = \arg \min_{h_j \in \mathcal{H}} \epsilon_j[/math], где [math]\epsilon_j = 
      \sum\limits_{i=1}^{m} D_t(i) [y_i\neq h_j(x_i)][/math]
  2. Если величина [math]\epsilon_t \geqslant 0.5[/math], то останавливаемся.
  3. Выбираем [math]\alpha_t \in \mathbf{R}[/math], обычно [math]\alpha_t = \frac{1}{2}\ln\frac{1-\epsilon_t}{\epsilon_t}[/math], где [math]\epsilon_t[/math] взвешенная ошибка классификатора 
     [math]h_t[/math]
  4. Обновляем: [math]D_{t+1}(i) = \frac{D_t(i)\exp^{-\alpha_t y_i h_t(x_i)}}{Z_t}[/math], где [math]Z_t[/math] является нормализующим параметром (выбранным так, чтобы [math]D_{t+1}[/math] являлось распределением вероятностей, то есть [math]\sum\limits_{i-1}^{m} D_{t+1}(i) = 1[/math]).
Строим результирующий классификатор:
 [math]H(x) = \textrm{sign}(\sum\limits_{t=1}^{T} \alpha_t h_t(x))[/math]
Выражение для обновления распредления [math]D_t[/math] должно быть сконструировано таким образом, чтобы выполнялось условие:
 [math]\exp^{\alpha_t y_i h_t(x_i)} \begin{cases}\lt 1,\ y(i) = h_t(x_i) \\ \gt 1,\ y(i) \neq h_t(x_i)\end{cases}[/math]

Таким образом, после выбора оптимального классификатора [math]h_t[/math] для распределения [math]D_t[/math], объекты [math]x_i[/math], которые классификатор [math]h_t[/math] идентифицирует корректно, имеют веса меньшие, чем те, которые идентифицируются некорректно. Следовательно, когда алгоритм тестирует классификаторы на распределении [math]D_{t+1}[/math], он будет выбирать классификатор, который лучше идентифицирует объекты неверно распознаваемые предыдущим классификатором.

Построение линейной комбинации классификаторов

Пусть [math]N(b,U^{l})[/math] — суммарный вес ошибочных классификаций, где $b$ — алгоритм классификации, [math]U^l = (u_1,...,u_l)[/math] — нормированный вектор весов объектов обучающей выборки.

Дано: [math]X^l, Y^l[/math] — обучающая выборка, $T$ — максимальное число базовых алгоритмов

1. инициализация весов объектов: [math]w_i = 1/l,\ i = 1,...,l;[/math]
2. для всех [math]t = 1,...,T[/math], пока не выполнен критерий останова:
  1. [math]b_t = \arg \min_{b} N(b; W^l)[/math];
  2. [math]\alpha_t = \frac{1}{2} \ln\frac{1 - N(b_t;W^l)}{N(b_t;W^l)}[/math];
  3. пересчет весов объектов: [math]w_i = w_i\exp{-\alpha_t y_i b_t(x_i)},\ i = 1,...,l[/math];
  4. нормировка весов объектов: [math]w_0 = \sum\limits_{j=1}^{l} w_j; w_i = \frac{w_i}{w_0},\ i = 1,...,l[/math];

Пример кода на python для scikit-learn

Классификатор sklearn.ensemble.AdaBoostClassifier имеет 5 параметров: base_estimator, n_estimators, learning_rate, algorithm, random_state. Наиболее важными являются:

  1. base_estimator — базовый алгоритм. По умолчанию используется DecisionTreeClassifier(max_depth=1)
  2. n_estimators — максимальное количество оценок, после которого бустинг прекращается. Если произойдет полное совпадение, то закончится раньше.
  3. learning_rate — вклад каждой модели в весовые коэффициенты и значение по умолчанию равно $1$. Снижение этого параметра будет означать, что весовые коэффициенты буду увеличиваться или уменьшаться в небольшой степени, вынуждая модель дольше обучаться (но иногда повышается производительность).
from sklearn.ensemble import AdaBoostClassifier
from sklearn import datasets
from sklearn.model_selection import train_test_split
from sklearn import metrics

iris = datasets.load_iris()

X = iris.data
y = iris.target

X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.3)
abc = AdaBoostClassifier(n_estimators=50, learning_rate=1)

model = abc.fit(X_train, y_train)

y_pred = model.predict(X_test)

print("Accuracy:",metrics.accuracy_score(y_test, y_pred))
Accuracy: 0.8888888888888888

Теперь рассмотрим алгоритм с SVC в качестве базы:

from sklearn.svm import SVC

svc=SVC(probability=True, kernel='linear')

abc = AdaBoostClassifier(base_estimator=svc, n_estimators=50, learning_rate=1)

model = abc.fit(X_train, y_train)

y_pred = model.predict(X_test)

print("Accuracy:",metrics.accuracy_score(y_test, y_pred))
Accuracy: 0.9555555555555556

Пример на языке Scala

SBT зависимость:

 libraryDependencies += "com.github.haifengl" %% "smile-scala" % "1.5.2"

Пример классификации датасета и вычисления F1 меры[5] используя smile.classification.adaboost[6]:

 import smile.classification._
 import smile.data._
 import smile.plot._
 import smile.read
 import smile.validation.FMeasure
 val iris: AttributeDataset = read.table("iris.csv", delimiter = ",", response = Some((new NumericAttribute("class"), 2)))
 val x: Array[Array[Double]] = iris.x()
 val y: Array[Int] = iris.y().map(_.toInt)
 val ada: AdaBoost = adaboost(x, y, ntrees = 500, maxNodes = 2)
 val predictions: Array[Int] = x.map(ada.predict)
 val f1Score = new FMeasure().measure(predictions, y)
 plot(x, y, ada)

См. также

Примечания

Источники информации

  1. AdaBoostClassifier — документация scikit-learn
  2. AdaBoost — статья на machinelearning.ru
  3. AdaBoost — презентация по AdaBoost