Минимизация эмпирического риска
Материал из MachineLearning.
| (7 промежуточных версий не показаны.) | |||
| Строка 1: | Строка 1: | ||
{{well|Статья написана с использованием LLM '''Gemini 3.1 Pro Preview''' и проверена участником Polina Khadralinova}} | {{well|Статья написана с использованием LLM '''Gemini 3.1 Pro Preview''' и проверена участником Polina Khadralinova}} | ||
| + | {{TOCright}} | ||
| - | ''' | + | '''Минимизация эмпирического риска''' (англ. ''Empirical Risk Minimization, ERM'') — это фундаментальный принцип машинного обучения. Он заключается в том, что из заданного семейства моделей мы выбираем ту, которая показывает наименьшую ошибку на доступных тренировочных данных. |
| - | + | Поскольку мы никогда не знаем всех возможных данных в природе, мы не можем измерить «истинную» ошибку модели. Принцип ERM постулирует, что вместо недоступной истинной ошибки мы можем опираться на её практическую оценку — эмпирический риск. | |
| - | == | + | == Введение == |
| - | + | Задача машинного обучения с учителем формулируется так. У нас есть множество объектов <tex>X</tex> (например, фотографии) и множество допустимых ответов <tex>Y</tex> (например, метки «кошка» или «собака»). Существует скрытая закономерность — идеальное правило <tex>y^*: X \to Y</tex>. | |
| - | + | Мы не знаем эту закономерность целиком. Нам доступна лишь конечная обучающая выборка <tex>X^{\ell} = \{ (x_i, y_i) \}_{i=1}^{\ell}</tex>. | |
| - | = | + | Наша цель — найти такой алгоритм <tex>a: X \to Y</tex> из заданного семейства моделей <tex>A = \{ a(x, w) \mid w \in W \}</tex>, который не просто запомнит правильные ответы для обучающей выборки, но и сможет эффективно (с минимальным количеством ошибок) работать на новых данных. Ошибку модели на конкретном объекте мы измеряем с помощью функции потерь <tex>\mathcal{L}(y, a(x))</tex>. |
| - | + | == Историческая справка == | |
| - | + | ||
| - | + | ||
| - | + | Истоки принципа минимизации эмпирического риска восходят к философской концепции [[Эмпирическая индукция|эмпирической индукции]] Фрэнсиса Бэкона (1620 г.), утверждавшего, что законы природы необходимо выводить путём обобщения фактов опыта, а не постулировать умозрительно. | |
| - | + | ||
| - | + | В строгой математической форме этот принцип впервые был применён Карлом Фридрихом Гауссом в 1795 году при разработке [[Метод наименьших квадратов|метода наименьших квадратов]] для расчёта орбит небесных тел. Гаусс предложил искать параметры модели, минимизируя сумму квадратов отклонений предсказанных значений от наблюдаемых. Позже, в 1936 году, Рональд Фишер применил схожий принцип для задачи классификации ([[Линейный дискриминантный анализ|линейный дискриминантный анализ]]). | |
| - | + | Своё современное теоретико-вероятностное обоснование принцип ERM получил в конце 1960-х — начале 1970-х годов в трудах советских математиков В. Н. Вапника и А. Я. Червоненкиса. В рамках созданной ими статистической теории обучения ([[Теория Вапника-Червоненкиса|теории Вапника-Червоненкиса]]) были строго сформулированы условия состоятельности принципа минимизации эмпирического риска и получены оценки скорости сходимости ошибки. | |
| - | + | == Ожидаемый и эмпирический риск == | |
| - | + | В теории машинного обучения мы считаем, что данные не появляются из ниоткуда. Существует неизвестное вероятностное распределение <tex>P(x, y)</tex> на пространстве <tex>X \times Y</tex>. Все наши объекты и ответы генерируются независимо из этого распределения. | |
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | В идеальном мире мы хотели бы создать модель, которая ошибается как можно реже на любых возможных данных. Эта идеальная мера ошибки называется '''ожидаемым риском''' (или истинным риском). Математически это математическое ожидание функции потерь по всему распределению: | |
| + | ::<tex>R(a) = \mathbb{E}_{x,y \sim P} [\mathcal{L}(y, a(x))] = \int_{X \times Y} \mathcal{L}(y, a(x)) dP(x, y)</tex> | ||
| + | '''Интуитивная аналогия:''' Представьте, что вы готовитесь к сложному экзамену. Ожидаемый риск — это ваша реальная оценка на самом экзамене, где вам может попасться абсолютно любой билет по предмету. Мы хотим сдать экзамен на «отлично», то есть свести ожидаемый риск <tex>R(a)</tex> к минимуму. | ||
| + | |||
| + | Проблема в том, что заранее узнать все вопросы экзамена (распределение <tex>P(x, y)</tex>) невозможно. Поэтому вычислить ожидаемый риск <tex>R(a)</tex> напрямую нельзя. | ||
| + | |||
| + | '''Зачем нам нужен показатель, который невозможно вычислить?''' Ожидаемый риск — это наш главный теоретический ориентир. Он показывает, насколько хорошо модель работает глобально. Хотя мы не можем рассчитать его напрямую, вся математическая теория машинного обучения построена на поиске путей приблизиться к нему. | ||
| + | |||
| + | Поскольку у нас нет всех данных мира, мы используем билеты прошлых лет — нашу обучающую выборку <tex>X^{\ell}</tex>. Закон больших чисел говорит нам, что теоретическое математическое ожидание можно оценить на практике, просто усреднив ошибку на всех доступных данных. Так мы получаем '''эмпирический риск''': | ||
| + | ::<tex>Q(a, X^{\ell}) = \frac{1}{\ell} \sum_{i=1}^{\ell} \mathcal{L}(y_i, a(x_i))</tex> | ||
| + | |||
| + | В нашей аналогии эмпирический риск <tex>Q(a, X^{\ell})</tex> — это доля ошибок, которые вы делаете, решая старые билеты дома. | ||
| + | |||
| + | Метод минимизации эмпирического риска делает логичный шаг. Раз мы не можем минимизировать ошибку на всех данных мира, давайте найдём алгоритм <tex>a^*</tex> (или вектор весов <tex>w^*</tex>), который доставляет минимум этому функционалу на обучающей выборке: | ||
| + | ::<tex>a^* = \arg\min_{a \in A} Q(a, X^{\ell})$</tex> | ||
| + | |||
| + | == Условия состоятельности и переобучение == | ||
| + | |||
| + | Главной проблемой принципа ERM является явление '''[[Переобучение|переобучения]]''' (overfitting). Возвращаясь к нашей аналогии: если вы просто зазубрите наизусть ответы к билетам прошлых лет, ваш эмпирический риск дома будет равен нулю. Но на реальном экзамене при малейшем изменении формулировки вопроса вы провалитесь (истинный риск <tex>R(a)</tex> окажется высоким). Модель теряет способность к обобщению. | ||
| + | |||
| + | Чтобы гарантировать, что наш метод обучения на компьютере вообще работает и не является делом чистой случайности, математики ввели понятие состоятельности. | ||
| + | |||
| + | === Зачем нужна состоятельность и равномерная сходимость? === | ||
| + | Метод ERM называется '''строго состоятельным''', если при увеличении объёма данных до бесконечности (<tex>\ell \to \infty</tex>) ошибка на обучении гарантированно приближается к реальной ошибке в жизни. | ||
| + | |||
| + | Однако простой сходимости для одной фиксированной модели недостаточно. Поскольку во время обучения мы ищем лучшую модель среди огромного семейства <tex>A</tex>, всегда есть риск найти «везучую» модель, которая случайно показала нулевую ошибку на тренировке, но будет бесполезна в реальности. | ||
| + | |||
| + | Чтобы исключить фактор случайного везения, сходимость должна быть '''равномерной по всему семейству моделей'''. Математически это означает, что даже для самой «везучей» модели разница между тренировочной и реальной ошибкой должна стремиться к нулю. Это записывается через предел супремума: | ||
| + | ::<tex>\lim_{\ell \to \infty} \mathbb{P} \left( \sup_{a \in A} |R(a) - Q(a, X^\ell)| > \epsilon \right) = 0</tex> | ||
| + | |||
| + | === Оценка Вапника — Червоненкиса === | ||
| + | Связующим мостом между теорией бесконечных данных и реальной практикой является знаменитая теорема Вапника — Червоненкиса. Она не просто постулирует предел, а математически доказывает его и показывает, с какой скоростью мы к нему приближаемся. | ||
| + | |||
| + | Согласно теореме, для конечной выборки объёмом <tex>\ell</tex> с вероятностью не менее <tex>1 - \delta</tex> для любого алгоритма <tex>a \in A</tex> справедливо неравенство: | ||
| + | ::<tex>R(a) \le Q(a, X^{\ell}) + \sqrt{\frac{h \left( \ln\left(\frac{2\ell}{h}\right) + 1 \right) - \ln\left(\frac{\delta}{4}\right)}{\ell}}</tex>, | ||
| + | где <tex>h</tex> — [[Размерность Вапника-Червоненкиса|VC-размерность]] (мера сложности) семейства моделей. | ||
| + | |||
| + | Правое слагаемое с квадратным корнем выступает в роли '''«штрафа за сложность»''' модели. Из него легко увидеть два фундаментальных правила: | ||
| + | * '''Влияние объёма данных:''' Чем больше у нас данных <tex>\ell</tex>, тем меньше штрафной корень. При <tex>\ell \to \infty</tex> этот корень стремится к нулю, математически доказывая вышеописанный предел равномерной сходимости. | ||
| + | * '''Влияние сложности модели:''' Чем сложнее семейство моделей <tex>h</tex>, тем больше штрафной корень и тем выше риск переобучения при фиксированном объёме данных. | ||
| + | == Регуляризация == | ||
| + | |||
| + | Для борьбы с переобучением метод ERM модифицируется. Мы не просто минимизируем ошибку на данных, но и штрафуем модель за излишнюю сложность. Этот подход известен как принцип структурной минимизации риска. | ||
| + | |||
| + | Оптимизируемый функционал принимает вид: | ||
| + | ::<tex>Q_{\text{reg}}(w, X^{\ell}) = \sum_{i=1}^{\ell} \mathcal{L}(y_i, a(x_i, w)) + \tau \mathcal{R}(w) \to \min_{w}</tex>, | ||
| + | где <tex>\mathcal{R}(w)</tex> — функция-регуляризатор, запрещающая параметрам модели принимать слишком большие или сложные значения, а <tex>\tau > 0</tex> — коэффициент регуляризации. | ||
| + | |||
| + | Наиболее часто используемые регуляризаторы: | ||
| + | * '''<tex>L_2</tex>-регуляризация''' (гребневая, Ridge): <tex>\mathcal{R}(w) = \|w\|_2^2 = \sum_{j=1}^n w_j^2</tex>. Не даёт весам модели неограниченно расти. | ||
| + | * '''<tex>L_1</tex>-регуляризация''' (Lasso): <tex>\mathcal{R}(w) = \|w\|_1 = \sum_{j=1}^n |w_j|</tex>. Способна занулять малозначимые признаки, упрощая итоговую модель. | ||
| + | |||
| + | == Основные типы функций потерь == | ||
| + | |||
| + | Выбор функции потерь <tex>\mathcal{L}(y, a(x))</tex> строго зависит от типа решаемой задачи и природы данных. | ||
| + | |||
| + | === В задачах регрессии === | ||
| + | В задачах [[Регрессионный анализ|регрессии]] мы предсказываем действительные числа (<tex>Y = \mathbb{R}</tex>). Типичные функции потерь: | ||
| + | * '''Квадратичная функция потерь''' (MSE, метод наименьших квадратов): | ||
| + | :<tex>\mathcal{L}(y, a(x)) = (a(x) - y)^2</tex> | ||
| + | :* ''Особенности:'' Очень чувствительна к аномалиям в данных (выбросам). Поскольку ошибка возводится в квадрат, алгоритм будет изо всех сил пытаться подстроиться под редкие опечатки или шумы в данных, жертвуя качеством на нормальных объектах. | ||
| + | * '''Абсолютная ошибка''' (MAE, модуль отклонения): | ||
| + | :<tex>\mathcal{L}(y, a(x)) = |a(x) - y|</tex> | ||
| + | :* ''Особенности:'' Является робастной (устойчивой) функцией потерь. Она одинаково штрафует за малые и большие отклонения, поэтому просто игнорирует редкие аномальные выбросы, выстраивая наиболее устойчивую зависимость. | ||
| + | |||
| + | === В задачах классификации === | ||
| + | В бинарной [[Классификация|классификации]] ответом является класс (<tex>Y = \{-1, +1\}</tex>). Самый естественный выбор — пороговая функция потерь, которая просто штрафует на единицу за несовпадение ответов и даёт ноль при правильном предсказании: | ||
| + | ::<tex>\mathcal{L}(y, a(x)) = [a(x) \neq y]</tex> | ||
| + | |||
| + | Однако минимизация такой функции на конечной выборке является дискретной NP-полной задачей комбинаторной оптимизации, так как функция потерь имеет вид «ступеньки» (у неё нулевой градиент почти везде, и её невозможно продифференцировать). | ||
| + | |||
| + | Для линейного классификатора <tex>a(x) = \text{sign}(\langle w, x \rangle)</tex> вводится понятие отступа (margin): <tex>M = y \langle w, x \rangle</tex>. Пороговая функция потерь выражается через отступ как <tex>[M < 0]</tex>. Чтобы применить методы градиентной оптимизации, пороговую функцию заменяют её гладкими или кусочно-гладкими верхними оценками <tex>\mathcal{L}(M)</tex>: | ||
| + | * '''Логистическая функция потерь''' (применяется в [[Логистическая регрессия|логистической регрессии]]): | ||
| + | :<tex>\mathcal{L}(M) = \log_2(1 + e^{-M})</tex> | ||
| + | :* ''Особенности:'' Гладкая функция. Она плавно штрафует модель даже за правильные, но неуверенные предсказания (когда отступ <tex>M</tex> положителен, но мал). Это заставляет классификатор не просто угадывать класс, а делать это с максимальной уверенностью. | ||
| + | * '''Кусочно-линейная функция потерь''' (Hinge loss, применяется в [[Метод опорных векторов|методе опорных векторов]], SVM): | ||
| + | :<tex>\mathcal{L}(M) = \max(0, 1 - M)</tex> | ||
| + | :* ''Особенности:'' Даёт ровно нулевую ошибку, если объект классифицирован правильно и с хорошим запасом (когда отступ <tex>M \ge 1</tex>). Это свойство приводит к тому, что оптимальное решение зависит только от небольшой доли объектов выборки (опорных векторов), делая решение разреженным. | ||
| + | * '''Экспоненциальная функция потерь''' (применяется в алгоритме AdaBoost): | ||
| + | :<tex>\mathcal{L}(M) = e^{-M}</tex> | ||
| + | :* ''Особенности:'' Очень жестко (экспоненциально) наказывает за ошибки (когда отступ <tex>M < 0</tex>). Позволяет быстро настраивать композиции алгоритмов, но крайне неустойчива к сильному шуму в данных. | ||
| + | |||
| + | == Методы оптимизации == | ||
| + | |||
| + | Когда функция потерь непрерывна и имеет производные по параметрам <tex>w</tex> (или хотя бы субградиенты), задачу минимизации эмпирического риска решают методами численной оптимизации. | ||
| + | |||
| + | === Классический градиентный спуск === | ||
| + | Исторически базовым алгоритмом является '''градиентный спуск''' (Gradient Descent, GD): | ||
| + | ::<tex>w_t = w_{t-1} - h \nabla_w Q(w_{t-1}, X^{\ell})</tex> | ||
| + | |||
| + | Распишем эту формулу подробнее, чтобы понять её физический смысл. Поскольку градиент суммы равен сумме градиентов, мы можем подставить формулу эмпирического риска прямо внутрь градиентного шага: | ||
| + | ::<tex>w_t = w_{t-1} - \frac{h}{\ell} \sum_{i=1}^{\ell} \nabla_w \mathcal{L}(y_i, a(x_i, w_{t-1}))</tex> | ||
| + | |||
| + | Разберём, из чего состоит это выражение: | ||
| + | * <tex>w_{t-1}</tex> — вектор весов на предыдущем шаге оптимизации; | ||
| + | * <tex>h</tex> — размер шага (темп обучения, learning rate); | ||
| + | * <tex>\nabla_w \mathcal{L}(y_i, a(x_i, w_{t-1}))</tex> — вектор частных производных (градиент) функции потерь на конкретном <tex>i</tex>-м объекте. Он показывает направление наискорейшего роста ошибки на этом объекте. Знак «минус» перед ним заставляет нас двигаться в противоположную сторону — в сторону уменьшения ошибки. | ||
| + | |||
| + | === Почему классический градиентный спуск работает медленно? === | ||
| + | Обрати внимание на знак суммы <tex>\sum_{i=1}^{\ell}</tex> в формуле выше. Чтобы сделать '''всего один маленький шаг''' обновления весов <tex>w_t</tex>, компьютер вынужден совершить так называемый «полный проход по всей выборке» (full batch pass). То есть процессор должен: | ||
| + | # Вычислить градиенты функции потерь отдельно для каждого из <tex>\ell</tex> объектов нашей обучающей выборки. | ||
| + | # Сложить все эти векторы градиентов между собой. | ||
| + | # Разделить сумму на <tex>\ell</tex>, чтобы найти среднее направление спуска. | ||
| + | |||
| + | Если наш датасет огромен (например, в эпохе Big Data выборка составляет <tex>\ell = 1\,000\ 000</tex> объектов), то для совершения одной крошечной итерации компьютеру придётся проделать один миллион тяжелых математических операций вычисления градиентов. Классический градиентный спуск в таких условиях работает катастрофически медленно. | ||
| + | |||
| + | === Почему стохастический метод (SGD) работает быстрее? === | ||
| + | В индустрии машинного и глубокого обучения повсеместно применяется '''[[Стохастический градиентный спуск|стохастический градиентный спуск]]''' (Stochastic Gradient Descent, SGD). | ||
| + | |||
| + | Идея метода заключается в том, что вместо вычисления точного среднего градиента по миллиону объектов, мы оцениваем его по '''одному случайно выбранному''' объекту <tex>x_i</tex>. На каждом шаге веса обновляются мгновенно: | ||
| + | ::<tex>w_t = w_{t-1} - h \nabla_w \mathcal{L}(y_i, a(x_i, w_{t-1}))</tex> | ||
| + | |||
| + | Вместо вычисления 1 000 000 градиентов мы вычисляем всего один. Пока классический алгоритм будет делать один-единственный тяжелый шаг по всей выборке, метод SGD успеет сделать миллион быстрых шагов и уйти далеко вперёд в сторону оптимума. Хотя траектория движения SGD получается шумной и «скачущей» (стохастической), в среднем он сходится к решению в сотни раз быстрее. | ||
| + | |||
| + | На практике чаще всего используется промежуточный пакетный вариант (Mini-batch SGD). При нём шаг делается не по одному объекту, а по небольшой группе объектов (батчу, обычно от 32 до 256 штук). Это снижает случайный шум стохастического градиента и позволяет эффективно распараллеливать вычисления на видеокартах (GPU). | ||
== См. также == | == См. также == | ||
| - | * [[ | + | * [[Теория Вапника-Червоненкиса]] |
| - | + | ||
* [[Переобучение]] | * [[Переобучение]] | ||
| - | * [[ | + | * [[Регуляризация]] |
== Примечания == | == Примечания == | ||
| Строка 43: | Строка 151: | ||
== Литература == | == Литература == | ||
| - | * | + | * {{книга | автор = Вапник В. Н., Червоненкис А. Я. | заглавие = Теория распознавания образов | место = М. | издательство = Наука | год = 1974 | страниц = 416 }} |
| - | * | + | * {{книга | автор = Вапник В. Н. | заглавие = Восстановление зависимостей по эмпирическим данным | место = М. | издательство = Наука | год = 1979 | страниц = 448 }} |
| + | * {{книга | автор = Hastie T., Tibshirani R., Friedman J. | заглавие = The Elements of Statistical Learning, 2nd edition | издательство = Springer | год = 2009 | страниц = 745 }} | ||
| + | * {{книга | автор = Воронцов К. В. | заглавие = Математические методы обучения по прецедентам | место = М. | издательство = МФТИ | год = 2012 }} | ||
[[Категория:Машинное обучение]] | [[Категория:Машинное обучение]] | ||
| - | [[Категория: | + | [[Категория:Математические методы]] |
| + | [[Категория:Теория вычислительного обучения]] | ||
Текущая версия
| | Статья написана с использованием LLM Gemini 3.1 Pro Preview и проверена участником Polina Khadralinova |
|
Минимизация эмпирического риска (англ. Empirical Risk Minimization, ERM) — это фундаментальный принцип машинного обучения. Он заключается в том, что из заданного семейства моделей мы выбираем ту, которая показывает наименьшую ошибку на доступных тренировочных данных.
Поскольку мы никогда не знаем всех возможных данных в природе, мы не можем измерить «истинную» ошибку модели. Принцип ERM постулирует, что вместо недоступной истинной ошибки мы можем опираться на её практическую оценку — эмпирический риск.
Введение
Задача машинного обучения с учителем формулируется так. У нас есть множество объектов (например, фотографии) и множество допустимых ответов
(например, метки «кошка» или «собака»). Существует скрытая закономерность — идеальное правило
.
Мы не знаем эту закономерность целиком. Нам доступна лишь конечная обучающая выборка .
Наша цель — найти такой алгоритм из заданного семейства моделей
, который не просто запомнит правильные ответы для обучающей выборки, но и сможет эффективно (с минимальным количеством ошибок) работать на новых данных. Ошибку модели на конкретном объекте мы измеряем с помощью функции потерь
.
Историческая справка
Истоки принципа минимизации эмпирического риска восходят к философской концепции эмпирической индукции Фрэнсиса Бэкона (1620 г.), утверждавшего, что законы природы необходимо выводить путём обобщения фактов опыта, а не постулировать умозрительно.
В строгой математической форме этот принцип впервые был применён Карлом Фридрихом Гауссом в 1795 году при разработке метода наименьших квадратов для расчёта орбит небесных тел. Гаусс предложил искать параметры модели, минимизируя сумму квадратов отклонений предсказанных значений от наблюдаемых. Позже, в 1936 году, Рональд Фишер применил схожий принцип для задачи классификации (линейный дискриминантный анализ).
Своё современное теоретико-вероятностное обоснование принцип ERM получил в конце 1960-х — начале 1970-х годов в трудах советских математиков В. Н. Вапника и А. Я. Червоненкиса. В рамках созданной ими статистической теории обучения (теории Вапника-Червоненкиса) были строго сформулированы условия состоятельности принципа минимизации эмпирического риска и получены оценки скорости сходимости ошибки.
Ожидаемый и эмпирический риск
В теории машинного обучения мы считаем, что данные не появляются из ниоткуда. Существует неизвестное вероятностное распределение на пространстве
. Все наши объекты и ответы генерируются независимо из этого распределения.
В идеальном мире мы хотели бы создать модель, которая ошибается как можно реже на любых возможных данных. Эта идеальная мера ошибки называется ожидаемым риском (или истинным риском). Математически это математическое ожидание функции потерь по всему распределению:
Интуитивная аналогия: Представьте, что вы готовитесь к сложному экзамену. Ожидаемый риск — это ваша реальная оценка на самом экзамене, где вам может попасться абсолютно любой билет по предмету. Мы хотим сдать экзамен на «отлично», то есть свести ожидаемый риск к минимуму.
Проблема в том, что заранее узнать все вопросы экзамена (распределение ) невозможно. Поэтому вычислить ожидаемый риск
напрямую нельзя.
Зачем нам нужен показатель, который невозможно вычислить? Ожидаемый риск — это наш главный теоретический ориентир. Он показывает, насколько хорошо модель работает глобально. Хотя мы не можем рассчитать его напрямую, вся математическая теория машинного обучения построена на поиске путей приблизиться к нему.
Поскольку у нас нет всех данных мира, мы используем билеты прошлых лет — нашу обучающую выборку . Закон больших чисел говорит нам, что теоретическое математическое ожидание можно оценить на практике, просто усреднив ошибку на всех доступных данных. Так мы получаем эмпирический риск:
В нашей аналогии эмпирический риск — это доля ошибок, которые вы делаете, решая старые билеты дома.
Метод минимизации эмпирического риска делает логичный шаг. Раз мы не можем минимизировать ошибку на всех данных мира, давайте найдём алгоритм (или вектор весов
), который доставляет минимум этому функционалу на обучающей выборке:
Условия состоятельности и переобучение
Главной проблемой принципа ERM является явление переобучения (overfitting). Возвращаясь к нашей аналогии: если вы просто зазубрите наизусть ответы к билетам прошлых лет, ваш эмпирический риск дома будет равен нулю. Но на реальном экзамене при малейшем изменении формулировки вопроса вы провалитесь (истинный риск окажется высоким). Модель теряет способность к обобщению.
Чтобы гарантировать, что наш метод обучения на компьютере вообще работает и не является делом чистой случайности, математики ввели понятие состоятельности.
Зачем нужна состоятельность и равномерная сходимость?
Метод ERM называется строго состоятельным, если при увеличении объёма данных до бесконечности () ошибка на обучении гарантированно приближается к реальной ошибке в жизни.
Однако простой сходимости для одной фиксированной модели недостаточно. Поскольку во время обучения мы ищем лучшую модель среди огромного семейства , всегда есть риск найти «везучую» модель, которая случайно показала нулевую ошибку на тренировке, но будет бесполезна в реальности.
Чтобы исключить фактор случайного везения, сходимость должна быть равномерной по всему семейству моделей. Математически это означает, что даже для самой «везучей» модели разница между тренировочной и реальной ошибкой должна стремиться к нулю. Это записывается через предел супремума:
Оценка Вапника — Червоненкиса
Связующим мостом между теорией бесконечных данных и реальной практикой является знаменитая теорема Вапника — Червоненкиса. Она не просто постулирует предел, а математически доказывает его и показывает, с какой скоростью мы к нему приближаемся.
Согласно теореме, для конечной выборки объёмом с вероятностью не менее
для любого алгоритма
справедливо неравенство:
,
где — VC-размерность (мера сложности) семейства моделей.
Правое слагаемое с квадратным корнем выступает в роли «штрафа за сложность» модели. Из него легко увидеть два фундаментальных правила:
- Влияние объёма данных: Чем больше у нас данных
, тем меньше штрафной корень. При
этот корень стремится к нулю, математически доказывая вышеописанный предел равномерной сходимости.
- Влияние сложности модели: Чем сложнее семейство моделей
, тем больше штрафной корень и тем выше риск переобучения при фиксированном объёме данных.
Регуляризация
Для борьбы с переобучением метод ERM модифицируется. Мы не просто минимизируем ошибку на данных, но и штрафуем модель за излишнюю сложность. Этот подход известен как принцип структурной минимизации риска.
Оптимизируемый функционал принимает вид:
,
где — функция-регуляризатор, запрещающая параметрам модели принимать слишком большие или сложные значения, а
— коэффициент регуляризации.
Наиболее часто используемые регуляризаторы:
-
-регуляризация (гребневая, Ridge):
. Не даёт весам модели неограниченно расти.
-
-регуляризация (Lasso):
. Способна занулять малозначимые признаки, упрощая итоговую модель.
Основные типы функций потерь
Выбор функции потерь строго зависит от типа решаемой задачи и природы данных.
В задачах регрессии
В задачах регрессии мы предсказываем действительные числа (). Типичные функции потерь:
- Квадратичная функция потерь (MSE, метод наименьших квадратов):
- Особенности: Очень чувствительна к аномалиям в данных (выбросам). Поскольку ошибка возводится в квадрат, алгоритм будет изо всех сил пытаться подстроиться под редкие опечатки или шумы в данных, жертвуя качеством на нормальных объектах.
- Абсолютная ошибка (MAE, модуль отклонения):
- Особенности: Является робастной (устойчивой) функцией потерь. Она одинаково штрафует за малые и большие отклонения, поэтому просто игнорирует редкие аномальные выбросы, выстраивая наиболее устойчивую зависимость.
В задачах классификации
В бинарной классификации ответом является класс (). Самый естественный выбор — пороговая функция потерь, которая просто штрафует на единицу за несовпадение ответов и даёт ноль при правильном предсказании:
Однако минимизация такой функции на конечной выборке является дискретной NP-полной задачей комбинаторной оптимизации, так как функция потерь имеет вид «ступеньки» (у неё нулевой градиент почти везде, и её невозможно продифференцировать).
Для линейного классификатора вводится понятие отступа (margin):
. Пороговая функция потерь выражается через отступ как
. Чтобы применить методы градиентной оптимизации, пороговую функцию заменяют её гладкими или кусочно-гладкими верхними оценками
:
- Логистическая функция потерь (применяется в логистической регрессии):
- Особенности: Гладкая функция. Она плавно штрафует модель даже за правильные, но неуверенные предсказания (когда отступ
положителен, но мал). Это заставляет классификатор не просто угадывать класс, а делать это с максимальной уверенностью.
- Особенности: Гладкая функция. Она плавно штрафует модель даже за правильные, но неуверенные предсказания (когда отступ
- Кусочно-линейная функция потерь (Hinge loss, применяется в методе опорных векторов, SVM):
- Особенности: Даёт ровно нулевую ошибку, если объект классифицирован правильно и с хорошим запасом (когда отступ
). Это свойство приводит к тому, что оптимальное решение зависит только от небольшой доли объектов выборки (опорных векторов), делая решение разреженным.
- Особенности: Даёт ровно нулевую ошибку, если объект классифицирован правильно и с хорошим запасом (когда отступ
- Экспоненциальная функция потерь (применяется в алгоритме AdaBoost):
- Особенности: Очень жестко (экспоненциально) наказывает за ошибки (когда отступ
). Позволяет быстро настраивать композиции алгоритмов, но крайне неустойчива к сильному шуму в данных.
- Особенности: Очень жестко (экспоненциально) наказывает за ошибки (когда отступ
Методы оптимизации
Когда функция потерь непрерывна и имеет производные по параметрам (или хотя бы субградиенты), задачу минимизации эмпирического риска решают методами численной оптимизации.
Классический градиентный спуск
Исторически базовым алгоритмом является градиентный спуск (Gradient Descent, GD):
Распишем эту формулу подробнее, чтобы понять её физический смысл. Поскольку градиент суммы равен сумме градиентов, мы можем подставить формулу эмпирического риска прямо внутрь градиентного шага:
Разберём, из чего состоит это выражение:
-
— вектор весов на предыдущем шаге оптимизации;
-
— размер шага (темп обучения, learning rate);
-
— вектор частных производных (градиент) функции потерь на конкретном
-м объекте. Он показывает направление наискорейшего роста ошибки на этом объекте. Знак «минус» перед ним заставляет нас двигаться в противоположную сторону — в сторону уменьшения ошибки.
Почему классический градиентный спуск работает медленно?
Обрати внимание на знак суммы в формуле выше. Чтобы сделать всего один маленький шаг обновления весов
, компьютер вынужден совершить так называемый «полный проход по всей выборке» (full batch pass). То есть процессор должен:
- Вычислить градиенты функции потерь отдельно для каждого из
объектов нашей обучающей выборки.
- Сложить все эти векторы градиентов между собой.
- Разделить сумму на
, чтобы найти среднее направление спуска.
Если наш датасет огромен (например, в эпохе Big Data выборка составляет объектов), то для совершения одной крошечной итерации компьютеру придётся проделать один миллион тяжелых математических операций вычисления градиентов. Классический градиентный спуск в таких условиях работает катастрофически медленно.
Почему стохастический метод (SGD) работает быстрее?
В индустрии машинного и глубокого обучения повсеместно применяется стохастический градиентный спуск (Stochastic Gradient Descent, SGD).
Идея метода заключается в том, что вместо вычисления точного среднего градиента по миллиону объектов, мы оцениваем его по одному случайно выбранному объекту . На каждом шаге веса обновляются мгновенно:
Вместо вычисления 1 000 000 градиентов мы вычисляем всего один. Пока классический алгоритм будет делать один-единственный тяжелый шаг по всей выборке, метод SGD успеет сделать миллион быстрых шагов и уйти далеко вперёд в сторону оптимума. Хотя траектория движения SGD получается шумной и «скачущей» (стохастической), в среднем он сходится к решению в сотни раз быстрее.
На практике чаще всего используется промежуточный пакетный вариант (Mini-batch SGD). При нём шаг делается не по одному объекту, а по небольшой группе объектов (батчу, обычно от 32 до 256 штук). Это снижает случайный шум стохастического градиента и позволяет эффективно распараллеливать вычисления на видеокартах (GPU).
См. также
Примечания
Литература
- Вапник В. Н., Червоненкис А. Я. Теория распознавания образов. — М.: Наука, 1974. — 416 с.
- Вапник В. Н. Восстановление зависимостей по эмпирическим данным. — М.: Наука, 1979. — 448 с.
- Hastie T., Tibshirani R., Friedman J. The Elements of Statistical Learning, 2nd edition. — Springer, 2009. — 745 с.
- Воронцов К. В. Математические методы обучения по прецедентам. — М.: МФТИ, 2012.

