Регрессионный анализ

Материал из MachineLearning.

(Различия между версиями)
Перейти к: навигация, поиск
(начало)
м (О терминах: исправление очипятки)
 
(30 промежуточных версий не показаны.)
Строка 1: Строка 1:
'''Регрессионный анализ''' — метод моделирования измеряемых данных и исследования их свойств.
'''Регрессионный анализ''' — метод моделирования измеряемых данных и исследования их свойств.
-
Данные состоят из пар значений [[зависимая переменная|зависимой переменной]] (переменной отклика)
+
Данные состоят из пар значений '''зависимой переменной''' (переменной отклика)
-
и [[независимая переменная|независимой переменной]] (объясняющей переменной).
+
и '''независимой переменной''' (объясняющей переменной).
-
Регрессионная модель есть функция независимой переменной и параметров с добавленной [[случайная переменная|случайной переменной]].
+
[[Регрессионная модель]] есть функция независимой переменной и параметров с добавленной [[случайная переменная|случайной переменной]].
Параметры модели настраиваются таким образом, что модель наилучшим образом приближает данные.
Параметры модели настраиваются таким образом, что модель наилучшим образом приближает данные.
Критерием качества приближения (целевой функцией) обычно является [[среднеквадратичная ошибка]]: сумма квадратов разности значений
Критерием качества приближения (целевой функцией) обычно является [[среднеквадратичная ошибка]]: сумма квадратов разности значений
Строка 8: Строка 8:
Регрессионный анализ — раздел [[математическая статистика|математической статистики]] и [[машинное обучение|машинного обучения]].
Регрессионный анализ — раздел [[математическая статистика|математической статистики]] и [[машинное обучение|машинного обучения]].
Предполагается, что зависимая переменная есть сумма значений некоторой модели и [[случайная величина|случайной величины]].
Предполагается, что зависимая переменная есть сумма значений некоторой модели и [[случайная величина|случайной величины]].
-
Относительно характера распределения этой величины делаются распределения, называемые гипотезой порождения данных.
+
Относительно характера распределения этой величины делаются предположения, называемые гипотезой порождения данных.
Для подтверждения или опровержения этой гипотезы выполняются [[статистический тест|статистические тесты]],
Для подтверждения или опровержения этой гипотезы выполняются [[статистический тест|статистические тесты]],
называемые [[анализ остатков|анализом остатков]].
называемые [[анализ остатков|анализом остатков]].
Строка 15: Строка 15:
[[тестирование гипотез|тестирования гипотез]] и выявления скрытых взаимосвязей в данных.
[[тестирование гипотез|тестирования гипотез]] и выявления скрытых взаимосвязей в данных.
-
== Определение ==
+
== Определение регрессионного анализа ==
[[Изображение:Regression_Analysis_Regression.gif|frame|right|Выборка может быть не функцией, а отношением. Например, данные для
[[Изображение:Regression_Analysis_Regression.gif|frame|right|Выборка может быть не функцией, а отношением. Например, данные для
-
построения регрессии могут быть такими: <tex>\{(0,0),(0,1),(0,2),</tex> <tex>
+
построения регрессии могут быть такими: <tex>\{(0,0),(0,1),(0,2),(1,1),(1,2),(1,3)\}</tex>. В такой выборке одному значению
-
(1,1),(1,2),(1,3)\}</tex>. В такой выборке одному значению
+
переменной&nbsp;<tex>x</tex> соответствует несколько значений переменной&nbsp;<tex>y</tex>.]]
переменной&nbsp;<tex>x</tex> соответствует несколько значений переменной&nbsp;<tex>y</tex>.]]
Строка 38: Строка 37:
Требуется найти наиболее вероятные параметры&nbsp;<tex>\bar{\mathbf{w}}</tex>:
Требуется найти наиболее вероятные параметры&nbsp;<tex>\bar{\mathbf{w}}</tex>:
<center><tex>\bar{\mathbf{w}}=\arg\max\limits_{\mathbf{w}\in\mathbb{R}^W}p(y|x,\mathbf{w},f)=p(D|\mathbf{w},f).</tex></center>
<center><tex>\bar{\mathbf{w}}=\arg\max\limits_{\mathbf{w}\in\mathbb{R}^W}p(y|x,\mathbf{w},f)=p(D|\mathbf{w},f).</tex></center>
-
Функция вероятности <tex>p</tex> зависит от гипотезы порождения данных и задается [[Байесовский вывод|Байесовским выводом]] или [[метод наибольшего правдоподобия|методом наибольшего правдоподобия]].
+
Функция вероятности <tex>p</tex> зависит от гипотезы порождения данных и задается [[Байесовский вывод|Байесовским выводом]] или [[метод наибольшего правдоподобия|методом наибольшего правдоподобия]].
== Линейная регрессия ==
== Линейная регрессия ==
 +
{{main|Многомерная линейная регрессия}}
 +
Линейная регрессия предполагает, что функция&nbsp;<tex>f</tex> зависит от параметров&nbsp;<tex>\mathbf{w}</tex> линейно.
Линейная регрессия предполагает, что функция&nbsp;<tex>f</tex> зависит от параметров&nbsp;<tex>\mathbf{w}</tex> линейно.
При этом линейная зависимость от свободной переменной&nbsp;<tex>\mathbf{x}</tex> необязательна,
При этом линейная зависимость от свободной переменной&nbsp;<tex>\mathbf{x}</tex> необязательна,
-
<center><tex> y=f(\mathbf{w},\mathbf{x})+\nu=\sum_{j=1}^W w_jg_j(\mathbf{x})+\nu. </tex></center>
+
<center><tex> y=f(\mathbf{w},\mathbf{x})+\nu=\sum_{j=1}^N w_jg_j(\mathbf{x})+\nu. </tex></center>
В случае, когда функция&nbsp;<tex>g\equiv\text{id}</tex> линейная регрессия имеет вид
В случае, когда функция&nbsp;<tex>g\equiv\text{id}</tex> линейная регрессия имеет вид
<center><tex> y=\sum_{j=1}^N w_jx_j+\nu=\langle\mathbf{w},\mathbf{x}\rangle +\nu, </tex></center>
<center><tex> y=\sum_{j=1}^N w_jx_j+\nu=\langle\mathbf{w},\mathbf{x}\rangle +\nu, </tex></center>
Строка 50: Строка 51:
Значения параметров в случае линейной регрессии находят с помощью [[метод наименьших квадратов|метода наименьших квадратов]].
Значения параметров в случае линейной регрессии находят с помощью [[метод наименьших квадратов|метода наименьших квадратов]].
Использование этого метода обосновано предположением о [[гауссовское распределение|гауссовском распределении]] случайной переменной.
Использование этого метода обосновано предположением о [[гауссовское распределение|гауссовском распределении]] случайной переменной.
-
При этом одна из важных оценок критерия качества полученной зависимости называется невязкой:
+
 
 +
Разности <tex>y_i-f(\mathbf{x}_i)</tex> между фактическими значениями зависимой переменной и восстановленными называются '''регрессионными остатками''' (residuals). В литературе используются также синонимы: ''невязки'' и ''ошибки''.
 +
Одной из важных оценок критерия качества полученной зависимости является сумма квадратов остатков:
<center><tex>SSE=\|f(\mathbf{x}_i)-y_i\|_2=\sum_{i=1}^N(y_i-f(\mathbf{w},\mathbf{x}_i))^2.</tex></center>
<center><tex>SSE=\|f(\mathbf{x}_i)-y_i\|_2=\sum_{i=1}^N(y_i-f(\mathbf{w},\mathbf{x}_i))^2.</tex></center>
-
Здесь&nbsp;<tex>SSE</tex>&nbsp;&#151; Sum of Squared Errors, сумма квадратных ошибок.
+
Здесь&nbsp;<tex>SSE</tex>&nbsp;&#151; Sum of Squared Errors.
-
Эти ошибки также называют остатками (residuals).
+
 
Дисперсия остатков вычисляется по формуле
Дисперсия остатков вычисляется по формуле
<center><tex>\bar{\sigma}^2_\nu=\frac{SSE}{N-2}=MSE.</tex></center>
<center><tex>\bar{\sigma}^2_\nu=\frac{SSE}{N-2}=MSE.</tex></center>
Строка 66: Строка 69:
На графиках представлены выборки, обозначенные синими точками, и регрессионные зависимости, обозначенные сплошными линиями.
На графиках представлены выборки, обозначенные синими точками, и регрессионные зависимости, обозначенные сплошными линиями.
По оси абсцисс отложена свободная переменная, а по оси ординат&nbsp;&#151; зависимая.
По оси абсцисс отложена свободная переменная, а по оси ординат&nbsp;&#151; зависимая.
-
Все три зависимости линейны относительно параметров.
+
Все три зависимости линейны относительно параметров.
 +
 
 +
== Нелинейная регрессия ==
 +
{{main|Нелинейная регрессия}}
 +
 
 +
Нелинейные регрессионные модели — [[регрессионная модель|модели]] вида
 +
::<tex>y=f(\mathbf{w},\mathbf{x})+\nu,</tex>
 +
которые не могут быть представлены в виде скалярного произведения
 +
::<tex>f(\mathbf{w},\mathbf{x})=(\mathbf{w},\mathbf{g}(\mathbf{x}))=\sum_{i=1}^n w_i g_i(\mathbf{x}),</tex>
 +
где <tex>\mathbf{w}=[w_1,\ldots, w_n]</tex>&nbsp;— параметры регрессионной модели,
 +
<tex>\mathbf{x}</tex>&nbsp;— свободная переменная из пространства&nbsp;<tex>\mathbb{R}^n</tex>, <tex>y</tex>&nbsp;— зависимая переменная,
 +
<tex>\nu</tex>&nbsp;— случайная величина и&nbsp;<tex>\mathbf{g}=[g_1,\ldots, g_n]</tex>&nbsp;— функция из некоторого
 +
заданного множества.
 +
 
 +
Значения параметров в случае нелинейной регрессии находят с помощью одного из методов градиентного спуска, например [[алгоритм Левенберга-Марквардта|алгоритма Левенберга-Марквардта]].
== О терминах ==
== О терминах ==
-
Термин "регрессия" был введен Фрэнсисом Гальтоном в конце 19-го века.
+
Термин "регрессия" был введён Фрэнсисом Гальтоном в конце 19-го века.
Гальтон обнаружил, что дети родителей с высоким или низким ростом обычно не наследуют выдающийся рост и назвал этот
Гальтон обнаружил, что дети родителей с высоким или низким ростом обычно не наследуют выдающийся рост и назвал этот
феномен "регрессия к посредственности".
феномен "регрессия к посредственности".
Строка 94: Строка 111:
Различают <i>параметрическую</i> и <i>непараметрическую</i> регрессию.
Различают <i>параметрическую</i> и <i>непараметрическую</i> регрессию.
Строгую границу между этими двумя типами регрессий провести сложно.
Строгую границу между этими двумя типами регрессий провести сложно.
-
Сейчас нет существует общепринятого критерия отличия одного типа моделей от другого.
+
Сейчас не существует общепринятого критерия отличия одного типа моделей от другого.
Например, считается, что линейные модели являются параметрическими, а модели, включающие усреднение зависимой переменной по пространству свободной переменной&nbsp;&#151;непараметрическими.
Например, считается, что линейные модели являются параметрическими, а модели, включающие усреднение зависимой переменной по пространству свободной переменной&nbsp;&#151;непараметрическими.
Пример параметрической регресионной модели: линейный предиктор, многослойный персептрон.
Пример параметрической регресионной модели: линейный предиктор, многослойный персептрон.
Строка 129: Строка 146:
Регрессия и классификация тесно связаны друг с другом.
Регрессия и классификация тесно связаны друг с другом.
-
Термин <i>алгоритм</i> в классификации мог бы стать синонимом термина <i>модель</i> в регрессии,
+
Термин <i>[[алгоритм]]</i> в классификации мог бы стать синонимом термина <i>[[регрессионная модель|модель]]</i> в регрессии,
если бы алгоритм не оперировал с дискретным множеством ответов-классов, а модель&nbsp;&#151; с непрерывно-определенной свободной переменной.
если бы алгоритм не оперировал с дискретным множеством ответов-классов, а модель&nbsp;&#151; с непрерывно-определенной свободной переменной.
 +
== Смотри также ==
 +
* [[Регрессионная модель]]
 +
* [[Линейная регрессия (пример)]]
 +
* [[Метод наименьших квадратов]]
 +
* [[Нелинейная регрессия]]
 +
* [[Алгоритм Левенберга-Марквардта]]
 +
* [[Связанный Байесовский вывод]]
 +
* [[Метод группового учета аргументов]]
 +
* [[Анализ регрессионных остатков]]
 +
* [[:Категория:Практика и вычислительные эксперименты|Применение методов регрессионного анализа на практике]]
 +
* [[Регрессионный анализ (рекомендуемые обозначения)]]
== Литература ==
== Литература ==
* Дрейпер&nbsp;Н., Смит&nbsp;Г. Прикладной регрессионный анализ. М.:&nbsp;Издательский дом «Вильямс».&nbsp;2007.
* Дрейпер&nbsp;Н., Смит&nbsp;Г. Прикладной регрессионный анализ. М.:&nbsp;Издательский дом «Вильямс».&nbsp;2007.
 +
* Айвазян&nbsp;С.А. Прикладная статистика и основы эконометрики. М.:&nbsp;Юнити.&nbsp;2001.
 +
* Брандт&nbsp;З. Анализ данных. М.:&nbsp;Мир.&nbsp;2003.
 +
* Стрижов В. В. Методы индуктивного порождения регрессионных моделей. М.: ВЦ РАН. 2008. 55&nbsp;с. [[Media:strijov08ln.pdf|Брошюра, PDF]].
 +
* Стрижов В.В., Крымова Е.А. Методы выбора регрессионных моделей. М.: ВЦ РАН, 2010. 60&nbsp;с. [[Media:Strijov-Krymova10Model-Selection.pdf|Брошюра, PDF]].
-
{{Заготовка}}
+
== Внешние ссылки ==
 +
* [http://en.wikipedia.org/wiki/Category:Regression_analysis Wikipedia.org, Категория «Regression analysis»]
 +
* [http://www.nlreg.com/ Non linear regression and curve fitting, NLREG software]
 +
* [http://www.ncrg.aston.ac.uk/netlab/ Netlab neural network software by Yan Nabney]
 +
* [http://curvefit.com/ Curvefit.com, A complete guide to nonlinear regression]
[[Категория:Регрессионный анализ]]
[[Категория:Регрессионный анализ]]
 +
[[Категория:Энциклопедия анализа данных]]
 +
[[Категория:Популярные и обзорные статьи]]

Текущая версия

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

Содержание

Определение регрессионного анализа

Выборка может быть не функцией, а отношением. Например, данные для построения регрессии могут быть такими: . В такой выборке одному значению переменной  соответствует несколько значений переменной .
Выборка может быть не функцией, а отношением. Например, данные для построения регрессии могут быть такими: \{(0,0),(0,1),(0,2),(1,1),(1,2),(1,3)\}. В такой выборке одному значению переменной x соответствует несколько значений переменной y.

Регрессия — зависимость математического ожидания (например, среднего значения) случайной величины от одной или нескольких других случайных величин (свободных переменных), то есть E(y|\mathbf{x})=f(\mathbf{x}). Регрессионным анализом называется поиск такой функции f, которая описывает эту зависимость. Регрессия может быть представлена в виде суммы неслучайной и случайной составляющих.

 y=f(\mathbf{x})+\nu,

где f — функция регрессионной зависимости, а \nu — аддитивная случайная величина с нулевым матожиданием. Предположение о характере распределения этой величины называется гипотезой порождения данных. Обычно предполагается, что величина \nu имеет гауссово распределение с нулевым средним и дисперсией \sigma^2_\nu.

Задача нахождения регрессионной модели нескольких свободных переменных ставится следующим образом. Задана выборка — множество \{\mathbf{x}_1,...,\mathbf{x}_N|\mathbf{x}\in\mathbb{R}^M\} значений свободных переменных и множество \{y_1,...,y_N| y\in\mathbb{R}\} соответствующих им значений зависимой переменной. Эти множества обозначаются как D, множество исходных данных \{(\mathbf{x},y)_i\}. Задана регрессионная модель — параметрическое семейство функций f(\mathbf{w},\mathbf{x}) зависящая от параметров \mathbf{w}\in\mathbb{R} и свободных переменных \mathbf{x}. Требуется найти наиболее вероятные параметры \bar{\mathbf{w}}:

\bar{\mathbf{w}}=\arg\max\limits_{\mathbf{w}\in\mathbb{R}^W}p(y|x,\mathbf{w},f)=p(D|\mathbf{w},f).

Функция вероятности p зависит от гипотезы порождения данных и задается Байесовским выводом или методом наибольшего правдоподобия.

Линейная регрессия

Линейная регрессия предполагает, что функция f зависит от параметров \mathbf{w} линейно. При этом линейная зависимость от свободной переменной \mathbf{x} необязательна,

 y=f(\mathbf{w},\mathbf{x})+\nu=\sum_{j=1}^N w_jg_j(\mathbf{x})+\nu.

В случае, когда функция g\equiv\text{id} линейная регрессия имеет вид

 y=\sum_{j=1}^N w_jx_j+\nu=\langle\mathbf{w},\mathbf{x}\rangle +\nu,

здесь x_j — компоненты вектора \mathbf{x}.

Значения параметров в случае линейной регрессии находят с помощью метода наименьших квадратов. Использование этого метода обосновано предположением о гауссовском распределении случайной переменной.

Разности y_i-f(\mathbf{x}_i) между фактическими значениями зависимой переменной и восстановленными называются регрессионными остатками (residuals). В литературе используются также синонимы: невязки и ошибки. Одной из важных оценок критерия качества полученной зависимости является сумма квадратов остатков:

SSE=\|f(\mathbf{x}_i)-y_i\|_2=\sum_{i=1}^N(y_i-f(\mathbf{w},\mathbf{x}_i))^2.

Здесь SSE — Sum of Squared Errors.

Дисперсия остатков вычисляется по формуле

\bar{\sigma}^2_\nu=\frac{SSE}{N-2}=MSE.

Здесь MSE — Mean Square Error, среднеквадратичная ошибка.

Изображение:Regression_Analysis_Linear.gif

Изображение:Regression_Analysis_Quadratic.gif

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

Нелинейная регрессия

Основная статья: Нелинейная регрессия

Нелинейные регрессионные модели — модели вида

y=f(\mathbf{w},\mathbf{x})+\nu,

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

f(\mathbf{w},\mathbf{x})=(\mathbf{w},\mathbf{g}(\mathbf{x}))=\sum_{i=1}^n w_i g_i(\mathbf{x}),

где \mathbf{w}=[w_1,\ldots, w_n] — параметры регрессионной модели, \mathbf{x} — свободная переменная из пространства \mathbb{R}^n, y — зависимая переменная, \nu — случайная величина и \mathbf{g}=[g_1,\ldots, g_n] — функция из некоторого заданного множества.

Значения параметров в случае нелинейной регрессии находят с помощью одного из методов градиентного спуска, например алгоритма Левенберга-Марквардта.

О терминах

Термин "регрессия" был введён Фрэнсисом Гальтоном в конце 19-го века. Гальтон обнаружил, что дети родителей с высоким или низким ростом обычно не наследуют выдающийся рост и назвал этот феномен "регрессия к посредственности". Сначала этот термин использовался исключительно в биологическом смысле. После работ Карла Пирсона этот термин стали использовать и в статистике.

Аппроксимация функций: непрерывная функция  приближает непрерывную или дискретную функцию 
Аппроксимация функций: непрерывная функция f приближает непрерывную или дискретную функцию u

В статистической литературе различают регрессию с участием одной свободной переменной и с несколькими свободными переменными — одномерную и многомерную регрессию. Предполагается, что мы используем несколько свободных переменных, то есть, свободная переменная — вектор \mathbf{x}\in\R^N. В частных случаях, когда свободная переменная является скаляром, она будет обозначаться x. Различают линейную и нелинейную регрессию. Если регрессионную модель не является линейной комбинацией функций от параметров, то говорят о нелинейной регрессии. При этом модель может быть произвольной суперпозицией функций g из некоторого набора. Нелинейными моделями являются, экспоненциальные, тригонометрические и другие (например, радиальные базисные функции или персептрон Розенблатта), полагающие зависимость между параметрами и зависимой переменной нелинейной.

Различают параметрическую и непараметрическую регрессию. Строгую границу между этими двумя типами регрессий провести сложно. Сейчас не существует общепринятого критерия отличия одного типа моделей от другого. Например, считается, что линейные модели являются параметрическими, а модели, включающие усреднение зависимой переменной по пространству свободной переменной —непараметрическими. Пример параметрической регресионной модели: линейный предиктор, многослойный персептрон. Примеры смешанной регрессионной модели: функции радиального базиса. Непараметрическая модель — скользящее усреднение в окне некоторой ширины. В целом, непараметрическая регрессия отличается от параметрической тем, что зависимая переменная зависит не от одного значения свободной переменной, а от некоторой заданной окрестности этого значения.

Интерполяция: функция  задана значениями узловых точек
Интерполяция: функция f задана значениями узловых точек

Есть различие между терминами: "приближение функций", "аппроксимация", "интерполяция", и "регрессия". Оно заключается в следующем.

Приближение функций. Дана функция u дискретного или непрерывного аргумента. Требуется найти функцию f из некоторого параметрическую семейства, например, среди алгебраических полиномов заданной степени. Параметры функции  f должны доставлять минимум некоторому функционалу, например,

\rho(u,f)=\left(\frac{1}{b-a}\int_a^b|f(x)-g(x)|^2dx\right)^{\frac{1}{2}}.

Термин аппроксимация — синоним термина "приближение функций". Чаще используется тогда, когда речь идет о заданной функции, как о функции дискретного аргумента. Здесь также требуется отыскать такую функцию f, которая проходит наиболее близко ко всем точкам заданной функции. При этом вводится понятие невязки — расстояния между точками непрерывной функции f и соответствующими точками функции u дискретного аргумента.

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

Регрессия и классификация тесно связаны друг с другом. Термин алгоритм в классификации мог бы стать синонимом термина модель в регрессии, если бы алгоритм не оперировал с дискретным множеством ответов-классов, а модель — с непрерывно-определенной свободной переменной.

Смотри также

Литература

  • Дрейпер Н., Смит Г. Прикладной регрессионный анализ. М.: Издательский дом «Вильямс». 2007.
  • Айвазян С.А. Прикладная статистика и основы эконометрики. М.: Юнити. 2001.
  • Брандт З. Анализ данных. М.: Мир. 2003.
  • Стрижов В. В. Методы индуктивного порождения регрессионных моделей. М.: ВЦ РАН. 2008. 55 с. Брошюра, PDF.
  • Стрижов В.В., Крымова Е.А. Методы выбора регрессионных моделей. М.: ВЦ РАН, 2010. 60 с. Брошюра, PDF.

Внешние ссылки

Личные инструменты