Критерий Акаике
Материал из MachineLearning.
м (→Модификации критерия) |
|||
| (15 промежуточных версий не показаны.) | |||
| Строка 1: | Строка 1: | ||
| - | {{ | + | {{well|Статья написана с использованием LLM '''DeepSeek-V3''' и проверена участником ~~Ilia Vdovin~~}} |
| - | ''' | + | |
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | = | + | = Критерий Акаике (AIC) = |
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | '''Критерий Акаике''' ('''AIC''' — от англ. ''Akaike Information Criterion'') — один из наиболее распространённых [[информационный критерий|информационных критериев]] для [[выбор модели|выбора статистических моделей]] по принципу максимума [[функция правдоподобия|правдоподобия]] с учётом их сложности. AIC оценивает, насколько хорошо модель будет предсказывать новые данные, и штрафует за излишнее количество параметров, предотвращая [[переобучение]]. | |
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | == | + | == Определение и мотивация == |
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | ||
| - | + | При построении статистических моделей возникает фундаментальная проблема: увеличение числа параметров всегда улучшает качество подгонки на обучающей выборке, но может ухудшить предсказания на новых данных. Классические подходы (проверка гипотез) не дают единой меры для сравнения невложенных моделей, а [[кросс-валидация]] вычислительно затратна. | |
| - | | | + | |
| - | + | AIC решает эту проблему, предлагая теоретически обоснованную оценку [[расстояние Кульбака–Лейблера|расстояния Кульбака–Лейблера]] между истинным распределением данных и оцениваемой моделью. Ключевая идея: Акаике показал, что максимизация логарифмического правдоподобия эквивалентна минимизации информационной потери при предсказании новых наблюдений. | |
| - | + | ||
| - | + | == Историческая справка == | |
| - | | | + | |
| - | | | + | Критерий был предложен японским статистиком [[Хиротогу Акаике]] в 1974 году в работе «A new look at the statistical model identification» (IEEE Transactions on Automatic Control) [1]. Акаике исходил из идей [[теория информации|теории информации]] Клода Шеннона и использовал [[энтропия|энтропию]] как меру неопределённости. |
| - | + | ||
| - | + | В своём подходе Акаике соединил две области — теорию информации и математическую статистику, показав, что максимизация логарифмического правдоподобия на обучающей выборке даёт смещённую оценку предсказательной способности модели, и это смещение можно скорректировать. | |
| - | + | ||
| - | + | == Теоретические основы == | |
| - | + | ||
| - | + | === Вывод через KL-дивергенцию === | |
| - | + | ||
| - | + | Пусть <tex>g(x)</tex> — истинное распределение данных (неизвестное), а <tex>f(x \mid \theta)</tex> — кандидатная модель с вектором параметров <tex>\theta</tex> размерности <tex>K</tex>. [[Расстояние Кульбака–Лейблера]] (KL-дивергенция) между ними определяется как: | |
| - | + | ||
| - | }} | + | <tex>D_{KL}(g \parallel f) = \int g(x) \ln \frac{g(x)}{f(x \mid \theta)} \, dx = \mathbb{E}_g[\ln g(x)] - \mathbb{E}_g[\ln f(x \mid \theta)].</tex> |
| - | + | ||
| - | + | '''Важное замечание:''' KL-дивергенция <tex>D_{KL}(g \parallel f)</tex> несимметрична: <tex>D_{KL}(g \parallel f) \neq D_{KL}(f \parallel g)</tex>, поэтому она не является метрическим расстоянием в строгом смысле. Однако она служит естественной мерой того, насколько модель <tex>f</tex> «удалена» от истинного распределения <tex>g</tex>. | |
| - | + | ||
| - | + | Первое слагаемое <tex>\mathbb{E}_g[\ln g(x)]</tex> одинаково для всех моделей, поэтому минимизация <tex>D_{KL}</tex> эквивалентна максимизации <tex>\mathbb{E}_g[\ln f(x \mid \theta)]</tex> — среднего логарифма правдоподобия по истинному распределению. | |
| - | + | ||
| - | + | === Основная идея Акаике === | |
| - | | | + | |
| - | + | Акаике рассуждал следующим образом. Если бы мы знали истинное распределение <tex>g(x)</tex>, то лучшая модель соответствовала бы максимуму <tex>\mathbb{E}_g[\ln f(x \mid \theta)]</tex>. Однако <tex>g(x)</tex> неизвестна, и мы можем лишь оценить параметры <tex>\theta</tex> по обучающей выборке <tex>y = (y_1, \dots, y_n)</tex>, получая оценку <tex>\hat{\theta}</tex>. | |
| - | + | ||
| - | }} | + | Логарифм правдоподобия на обучающей выборке <tex>\ln L(\hat{\theta} \mid y)</tex> является ''смещённой'' оценкой величины <tex>\mathbb{E}_g[\ln f(x \mid \theta)]</tex>: он систематически завышает качество модели, так как параметры подгонялись под те же данные. Акаике показал, что величина смещения асимптотически равна числу параметров <tex>K</tex>: |
| - | + | ||
| - | | | + | <tex>\mathbb{E}_y\left[ \ln L(\hat{\theta} \mid y) - \mathbb{E}_g[\ln f(x \mid \theta)] \right] \approx K.</tex> |
| - | + | ||
| - | + | Отсюда получается несмещённая оценка ожидаемого логарифма правдоподобия на новых данных: | |
| - | + | ||
| - | + | <tex>\widehat{\mathbb{E}_g[\ln f]} = \ln L(\hat{\theta} \mid y) - K.</tex> | |
| - | | | + | |
| - | + | Умножив на <tex>-2</tex> (по историческим причинам, чтобы согласовать с [[хи-квадрат|распределением <tex>\chi^2</tex>]]), получают итоговую формулу: | |
| - | + | ||
| - | # | + | <tex>AIC = -2 \ln L(\hat{\theta} \mid y) + 2K.</tex> |
| - | + | ||
| - | + | === Интерпретация формулы === | |
| - | + | ||
| - | + | * Первое слагаемое <tex>-2 \ln L(\hat{\theta} \mid y)</tex> — характеризует качество подгонки модели на обучающих данных (чем меньше, тем лучше). | |
| + | * Второе слагаемое <tex>2K</tex> — штраф за сложность модели (каждый дополнительный параметр увеличивает AIC на 2). | ||
| + | |||
| + | Таким образом, AIC балансирует между точностью подгонки и сложностью модели. Добавление нового параметра оправдано только в том случае, если оно уменьшает <tex>-2 \ln L</tex> более чем на 2. | ||
| + | |||
| + | == Интерпретация и применение == | ||
| + | |||
| + | AIC является ''относительной'' мерой качества модели. Абсолютное значение AIC не интерпретируется; важны лишь различия между моделями. При сравнении нескольких моделей вычисляют разности: | ||
| + | |||
| + | <tex>\Delta_i = AIC_i - AIC_{min},</tex> | ||
| + | |||
| + | где <tex>AIC_{min}</tex> — минимальное значение среди всех рассматриваемых моделей. Эмпирическое правило: | ||
| + | |||
| + | * <tex>\Delta_i \le 2</tex> — модели практически эквивалентны; | ||
| + | * <tex>4 \le \Delta_i \le 7</tex> — различие заметно; | ||
| + | * <tex>\Delta_i > 10</tex> — модель существенно хуже, её можно исключить из рассмотрения. | ||
| + | |||
| + | Также вычисляют '''веса Акаике''': | ||
| + | |||
| + | <tex>w_i = \frac{\exp(-\Delta_i/2)}{\sum_j \exp(-\Delta_j/2)},</tex> | ||
| + | |||
| + | которые интерпретируются как вероятности того, что модель <tex>i</tex> является наилучшей среди рассматриваемых (в смысле минимизации KL-дивергенции). | ||
| + | |||
| + | === Важные ограничения при применении === | ||
| + | |||
| + | * AIC позволяет сравнивать только модели, построенные на '''одной и той же''' выборке с одинаковым набором наблюдений. | ||
| + | * Зависимая переменная должна быть идентичной (нельзя сравнивать модели с разными преобразованиями <tex>y</tex>). | ||
| + | * Модели должны быть оценены [[метод максимального правдоподобия|методом максимального правдоподобия]]. | ||
| + | |||
| + | === Пример === | ||
| + | |||
| + | Пусть для данных подбираются линейная регрессия с 2 параметрами (<tex>AIC=120</tex>) и полиномиальная модель с 5 параметрами (<tex>AIC=115</tex>). Разность <tex>\Delta=5</tex> указывает на умеренное преимущество полиномиальной модели. Если бы <tex>\Delta>10</tex>, выбор полиномиальной модели был бы однозначным. | ||
| + | |||
| + | == Модификации == | ||
| + | |||
| + | === AICc (исправленный AIC для малых выборок) === | ||
| + | |||
| + | При малых выборках стандартный AIC может давать смещённые оценки. Сугиура (1978) предложил поправку: | ||
| + | |||
| + | <tex>AICc = AIC + \frac{2K(K+1)}{n - K - 1}.</tex> | ||
| + | |||
| + | Рекомендуется использовать при <tex>n/K < 40</tex>. При <tex>n \to \infty</tex> поправка стремится к нулю, и <tex>AICc \to AIC</tex>. | ||
| + | |||
| + | === QAIC (Quasi-AIC) === | ||
| + | |||
| + | Для данных с [[передисперсия|передисперсией]] (например, в экологических моделях) используется квази-правдоподобие. Вводится коэффициент дисперсии <tex>\hat{c}</tex>: | ||
| + | |||
| + | <tex>QAIC = -2 \ln L / \hat{c} + 2K.</tex> | ||
| + | |||
| + | Применяется, когда распределение данных не соответствует стандартным семействам (например, [[биномиальное распределение|биномиальное]] с избыточной вариативностью). | ||
| + | |||
| + | == Ограничения == | ||
| + | |||
| + | * '''Несостоятельность''': AIC не стремится к выбору истинной модели с вероятностью 1 при <tex>n \to \infty</tex> (в отличие от [[байесовский информационный критерий|BIC]]). Если истинная модель конечномерна и входит в множество кандидатов, AIC будет выбирать модели с избыточным числом параметров, дающие сколь угодно малое улучшение правдоподобия. | ||
| + | * Неприменим для сравнения моделей с разными объёмами выборок или разными зависимыми переменными. | ||
| + | * Чувствителен к выбору распределения; при неверной спецификации распределения выводы могут быть ошибочными. | ||
| + | * Требует оценки методом максимального правдоподобия. | ||
| + | |||
| + | == Сравнение с другими критериями == | ||
| + | |||
| + | === BIC (Байесовский информационный критерий) === | ||
| + | |||
| + | Предложен Шварцем (1978): <tex>BIC = -2 \ln L + K \ln n</tex>. Штраф растёт с объёмом выборки, поэтому BIC более консервативен и асимптотически состоятелен. При <tex>n > 8</tex> BIC штрафует за параметры сильнее, чем AIC. BIC предпочтителен, когда цель — идентификация истинной структуры модели, а не прогнозирование. | ||
| + | |||
| + | === DIC и WAIC === | ||
| + | |||
| + | * '''[[DIC]]''' (Deviance Information Criterion) — для байесовских моделей, использует апостериорное среднее отклонений. | ||
| + | * '''[[WAIC]]''' (Widely Applicable Information Criterion) — полностью байесовский критерий, согласованный с [[кросс-валидация|кросс-валидацией]] leave-one-out. | ||
| + | |||
| + | === Кросс-валидация === | ||
| + | |||
| + | Даёт прямую оценку ошибки прогноза, но требует многократного переобучения. При больших данных может быть предпочтительнее AIC. Однако AIC является асимптотически эквивалентным кросс-валидации leave-one-out. | ||
| + | |||
| + | == Практические рекомендации == | ||
| + | |||
| + | * Если цель — прогнозирование и размер выборки мал (<tex>n/K < 40</tex>), используйте AICc. | ||
| + | * Если цель — идентификация истинной структуры модели при большом <tex>n</tex>, предпочтительнее BIC. | ||
| + | * Всегда приводите значения AIC для всех сравниваемых моделей, а также <tex>\Delta_i</tex> и веса Акаике. | ||
| + | * Избегайте «драпировки» моделей (data dredging) — сравнение большого числа моделей без априорных гипотез повышает риск случайных находок. | ||
| + | |||
| + | == Заключение == | ||
| + | |||
| + | Критерий Акаике остаётся одним из самых цитируемых инструментов в статистическом моделировании благодаря своей простоте, теоретической прозрачности и универсальности. Несмотря на разработанные позднее альтернативы, AIC сохраняет актуальность в задачах прогнозирования, особенно при умеренных размерах выборки. Современные тенденции — интеграция AIC с методами регуляризации ([[LASSO]], [[ридж-регрессия]]) и использование информационных критериев в глубоком обучении — расширяют область его применения. | ||
| + | |||
| + | == Литература == | ||
| + | |||
| + | # Akaike, H. (1974). A new look at the statistical model identification. ''IEEE Transactions on Automatic Control'', 19(6), 716–723. | ||
| + | # Burnham, K. P., & Anderson, D. R. (2002). ''Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach''. 2nd ed. Springer. | ||
| + | # McQuarrie, A. D. R., & Tsai, C.-L. (1998). ''Regression and Time Series Model Selection''. World Scientific. | ||
| + | # Konishi, S., & Kitagawa, G. (2008). ''Information Criteria and Statistical Modeling''. Springer. | ||
| + | # Claeskens, G., & Hjort, N. L. (2008). ''Model Selection and Model Averaging''. Cambridge University Press. | ||
| + | |||
| + | Полный промпт, использованный при создании этой статьи, доступен на [[Обсуждение:Критерий Акаике|странице обсуждения]]. | ||
| + | |||
| + | [[Категория:Статистические критерии]] | ||
| + | [[Категория:Выбор модели]] | ||
| + | [[Категория:Теория информации]] | ||
Текущая версия
| | Статья написана с использованием LLM DeepSeek-V3 и проверена участником ~~Ilia Vdovin~~ |
Содержание |
Критерий Акаике (AIC)
Критерий Акаике (AIC — от англ. Akaike Information Criterion) — один из наиболее распространённых информационных критериев для выбора статистических моделей по принципу максимума правдоподобия с учётом их сложности. AIC оценивает, насколько хорошо модель будет предсказывать новые данные, и штрафует за излишнее количество параметров, предотвращая переобучение.
Определение и мотивация
При построении статистических моделей возникает фундаментальная проблема: увеличение числа параметров всегда улучшает качество подгонки на обучающей выборке, но может ухудшить предсказания на новых данных. Классические подходы (проверка гипотез) не дают единой меры для сравнения невложенных моделей, а кросс-валидация вычислительно затратна.
AIC решает эту проблему, предлагая теоретически обоснованную оценку расстояния Кульбака–Лейблера между истинным распределением данных и оцениваемой моделью. Ключевая идея: Акаике показал, что максимизация логарифмического правдоподобия эквивалентна минимизации информационной потери при предсказании новых наблюдений.
Историческая справка
Критерий был предложен японским статистиком Хиротогу Акаике в 1974 году в работе «A new look at the statistical model identification» (IEEE Transactions on Automatic Control) [1]. Акаике исходил из идей теории информации Клода Шеннона и использовал энтропию как меру неопределённости.
В своём подходе Акаике соединил две области — теорию информации и математическую статистику, показав, что максимизация логарифмического правдоподобия на обучающей выборке даёт смещённую оценку предсказательной способности модели, и это смещение можно скорректировать.
Теоретические основы
Вывод через KL-дивергенцию
Пусть — истинное распределение данных (неизвестное), а
— кандидатная модель с вектором параметров
размерности
. Расстояние Кульбака–Лейблера (KL-дивергенция) между ними определяется как:
Важное замечание: KL-дивергенция несимметрична:
, поэтому она не является метрическим расстоянием в строгом смысле. Однако она служит естественной мерой того, насколько модель
«удалена» от истинного распределения
.
Первое слагаемое одинаково для всех моделей, поэтому минимизация
эквивалентна максимизации
— среднего логарифма правдоподобия по истинному распределению.
Основная идея Акаике
Акаике рассуждал следующим образом. Если бы мы знали истинное распределение , то лучшая модель соответствовала бы максимуму
. Однако
неизвестна, и мы можем лишь оценить параметры
по обучающей выборке
, получая оценку
.
Логарифм правдоподобия на обучающей выборке является смещённой оценкой величины
: он систематически завышает качество модели, так как параметры подгонялись под те же данные. Акаике показал, что величина смещения асимптотически равна числу параметров
:
Отсюда получается несмещённая оценка ожидаемого логарифма правдоподобия на новых данных:
Умножив на (по историческим причинам, чтобы согласовать с распределением
), получают итоговую формулу:
Интерпретация формулы
- Первое слагаемое
— характеризует качество подгонки модели на обучающих данных (чем меньше, тем лучше).
- Второе слагаемое
— штраф за сложность модели (каждый дополнительный параметр увеличивает AIC на 2).
Таким образом, AIC балансирует между точностью подгонки и сложностью модели. Добавление нового параметра оправдано только в том случае, если оно уменьшает более чем на 2.
Интерпретация и применение
AIC является относительной мерой качества модели. Абсолютное значение AIC не интерпретируется; важны лишь различия между моделями. При сравнении нескольких моделей вычисляют разности:
где — минимальное значение среди всех рассматриваемых моделей. Эмпирическое правило:
-
— модели практически эквивалентны;
-
— различие заметно;
-
— модель существенно хуже, её можно исключить из рассмотрения.
Также вычисляют веса Акаике:
которые интерпретируются как вероятности того, что модель является наилучшей среди рассматриваемых (в смысле минимизации KL-дивергенции).
Важные ограничения при применении
- AIC позволяет сравнивать только модели, построенные на одной и той же выборке с одинаковым набором наблюдений.
- Зависимая переменная должна быть идентичной (нельзя сравнивать модели с разными преобразованиями
).
- Модели должны быть оценены методом максимального правдоподобия.
Пример
Пусть для данных подбираются линейная регрессия с 2 параметрами () и полиномиальная модель с 5 параметрами (
). Разность
указывает на умеренное преимущество полиномиальной модели. Если бы
, выбор полиномиальной модели был бы однозначным.
Модификации
AICc (исправленный AIC для малых выборок)
При малых выборках стандартный AIC может давать смещённые оценки. Сугиура (1978) предложил поправку:
Рекомендуется использовать при . При
поправка стремится к нулю, и
.
QAIC (Quasi-AIC)
Для данных с передисперсией (например, в экологических моделях) используется квази-правдоподобие. Вводится коэффициент дисперсии :
Применяется, когда распределение данных не соответствует стандартным семействам (например, биномиальное с избыточной вариативностью).
Ограничения
- Несостоятельность: AIC не стремится к выбору истинной модели с вероятностью 1 при
(в отличие от BIC). Если истинная модель конечномерна и входит в множество кандидатов, AIC будет выбирать модели с избыточным числом параметров, дающие сколь угодно малое улучшение правдоподобия.
- Неприменим для сравнения моделей с разными объёмами выборок или разными зависимыми переменными.
- Чувствителен к выбору распределения; при неверной спецификации распределения выводы могут быть ошибочными.
- Требует оценки методом максимального правдоподобия.
Сравнение с другими критериями
BIC (Байесовский информационный критерий)
Предложен Шварцем (1978): . Штраф растёт с объёмом выборки, поэтому BIC более консервативен и асимптотически состоятелен. При
BIC штрафует за параметры сильнее, чем AIC. BIC предпочтителен, когда цель — идентификация истинной структуры модели, а не прогнозирование.
DIC и WAIC
- DIC (Deviance Information Criterion) — для байесовских моделей, использует апостериорное среднее отклонений.
- WAIC (Widely Applicable Information Criterion) — полностью байесовский критерий, согласованный с кросс-валидацией leave-one-out.
Кросс-валидация
Даёт прямую оценку ошибки прогноза, но требует многократного переобучения. При больших данных может быть предпочтительнее AIC. Однако AIC является асимптотически эквивалентным кросс-валидации leave-one-out.
Практические рекомендации
- Если цель — прогнозирование и размер выборки мал (
), используйте AICc.
- Если цель — идентификация истинной структуры модели при большом
, предпочтительнее BIC.
- Всегда приводите значения AIC для всех сравниваемых моделей, а также
и веса Акаике.
- Избегайте «драпировки» моделей (data dredging) — сравнение большого числа моделей без априорных гипотез повышает риск случайных находок.
Заключение
Критерий Акаике остаётся одним из самых цитируемых инструментов в статистическом моделировании благодаря своей простоте, теоретической прозрачности и универсальности. Несмотря на разработанные позднее альтернативы, AIC сохраняет актуальность в задачах прогнозирования, особенно при умеренных размерах выборки. Современные тенденции — интеграция AIC с методами регуляризации (LASSO, ридж-регрессия) и использование информационных критериев в глубоком обучении — расширяют область его применения.
Литература
- Akaike, H. (1974). A new look at the statistical model identification. IEEE Transactions on Automatic Control, 19(6), 716–723.
- Burnham, K. P., & Anderson, D. R. (2002). Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2nd ed. Springer.
- McQuarrie, A. D. R., & Tsai, C.-L. (1998). Regression and Time Series Model Selection. World Scientific.
- Konishi, S., & Kitagawa, G. (2008). Information Criteria and Statistical Modeling. Springer.
- Claeskens, G., & Hjort, N. L. (2008). Model Selection and Model Averaging. Cambridge University Press.
Полный промпт, использованный при создании этой статьи, доступен на странице обсуждения.

