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

Введение

Модель данных явным образом описывает отношение между переменными predictor и response. Линейная регрессия соответствует модели данных, которая линейна в коэффициентах модели. Наиболее распространенным типом линейной регрессии является least-squares fit, который может соответствовать и линиям и полиномам среди других линейных моделей.

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

MATLAB^® Пользовательский интерфейс Basic Fitting помогает вам соответствовать своим данным, таким образом, можно вычислить коэффициенты модели и построить модель сверху данных. Для примера смотрите Пример: Используя Пользовательский интерфейс Basic Fitting. Также можно использовать MATLAB polyfit и polyval функции, чтобы соответствовать вашим данным к модели, которая линейна в коэффициентах. Для примера см. Программное аппроксимирование.

Если необходимо соответствовать данным нелинейной моделью, преобразуйте переменные, чтобы сделать отношение линейным. В качестве альтернативы попытайтесь соответствовать нелинейной функции непосредственно с помощью любого Statistics and Machine Learning Toolbox™ nlinfit функция, Optimization Toolbox™ lsqcurvefit функция, или путем применения функций в Curve Fitting Toolbox™.

Эта тема объясняет как:

Выполните простую линейную регрессию с помощью \ оператор.
Используйте корреляционный анализ, чтобы определить, связаны ли две величины для обоснования подгонки данных.
Подбор линейной модели к данным.
Оцените качество подгонки путем графического вывода остаточных значений и поиска шаблонов.
Вычислите меры качества подгонки R² и настроенный R²

Простая линейная регрессия

Скрипт Open Live Script

В этом примере показано, как выполнить простую линейную регрессию с помощью accidents набор данных. Пример также показывает вам, как вычислить коэффициент детерминации $R^{2}$ оценивать регрессии. accidents набор данных содержит данные для дорожно-транспортных происшествий со смертельным исходом в штатах США.

Линейная регрессия моделирует отношение между зависимым, или ответ, переменную $y$ и один или несколько независимый, или предиктор, переменные $x_{1}, . . ., x_{n}$ . Простая линейная регрессия рассматривает только одну независимую переменную с помощью отношения

$y = β_{0} + β_{1} x + ϵ,$

где $β_{0}$ свободный член, $β_{1}$ наклон (или коэффициент регрессии), и $ϵ$ остаточный член.

Начните с ряда $n$ наблюдаемые величины $x$ и $y$ данный $(x_{1}, y_{1})$ , $(x_{2}, y_{2})$ , ..., $(x_{n}, y_{n})$ . Используя простое отношение линейной регрессии, эти значения формируют систему линейных уравнений. Представляйте эти уравнения в матричной форме как

$[\begin{array}{c} y_{1} \\ y_{2} \\ ⋮ \\ y_{n} \end{array}] = [\begin{array}{c} 1 & x_{1} \\ 1 & x_{2} \\ ⋮ & ⋮ \\ 1 & x_{n} \end{array}] [\begin{array}{c} β_{0} \\ β_{1} \end{array}] .$

Пусть

$Y = [\begin{array}{c} y_{1} \\ y_{2} \\ ⋮ \\ y_{n} \end{array}], X = [\begin{array}{c} 1 & x_{1} \\ 1 & x_{2} \\ ⋮ & ⋮ \\ 1 & x_{n} \end{array}], B = [\begin{array}{c} β_{0} \\ β_{1} \end{array}] .$

Отношение теперь $Y = X B$ .

В MATLAB можно найти $B$ использование mldivide оператор как B = X\Y.

От набора данных accidents, загрузите данные о несчастном случае в y и данные о населении состояния в x. Найдите отношение линейной регрессии $y = β_{1} x$ между несчастными случаями в состоянии и населении состояния с помощью \ оператор. \ оператор выполняет регрессию наименьших квадратов.

load accidents
x = hwydata(:,14); %Population of states
y = hwydata(:,4); %Accidents per state
format long
b1 = x\y

b1 = 
     1.372716735564871e-04

b1 наклон или коэффициент регрессии. Линейное отношение $y = β_{1} x = 0.0001372 x$ .

Вычислите несчастные случаи на yCalc состояния от x использование отношения. Визуализируйте регрессию путем графического вывода фактических значений y и расчетные значения yCalc.

yCalc1 = b1*x;
scatter(x,y)
hold on
plot(x,yCalc1)
xlabel('Population of state')
ylabel('Fatal traffic accidents per state')
title('Linear Regression Relation Between Accidents & Population')
grid on

Figure contains an axes object. The axes object with title Linear Regression Relation Between Accidents & Population contains 2 objects of type scatter, line.

Улучшите подгонку включением свободного члена $β_{0}$ в вашей модели как $y = β_{0} + β_{1} x$ . Вычислить $β_{0}$ путем дополнения x со столбцом из единиц и использования \ оператор.

X = [ones(length(x),1) x];
b = X\y

Этот результат представляет отношение $y = β_{0} + β_{1} x = 142.7120 + 0.0001256 x$ .

Визуализируйте отношение путем графического вывода его на той же фигуре.

yCalc2 = X*b;
plot(x,yCalc2,'--')
legend('Data','Slope','Slope & Intercept','Location','best');

Figure contains an axes object. The axes object with title Linear Regression Relation Between Accidents & Population contains 3 objects of type scatter, line. These objects represent Data, Slope, Slope & Intercept.

От фигуры две подгонки выглядят подобными. Один метод, чтобы найти лучшую подгонку должен вычислить коэффициент детерминации, $R^{2}$ . $R^{2}$ одна мера того, как хорошо модель может предсказать данные и падения между $0$ и $1$ . Выше значение $R^{2}$ , лучше модель при предсказании данных.

Где $y_{}^{ˆ}$ представляет расчетные значения $y$ и $y_{}^{‾}$ среднее значение $y$ , $R^{2}$ задан как

$R^{2} = 1 - \frac{\sum_{i = 1}^{n} {(y_{i} - {y_{}^{ˆ}}_{i})}^{2}}{\sum_{i = 1}^{n} {(y_{i} - y_{}^{‾})}^{2}} .$

Найдите лучший припадок двух подгонок путем сравнения значений $R^{2}$ . Как $R^{2}$ значения показывают, вторая подгонка, которая включает свободный член, лучше.

Rsq1 = 1 - sum((y - yCalc1).^2)/sum((y - mean(y)).^2)

Rsq1 = 
   0.822235650485566

Rsq2 = 1 - sum((y - yCalc2).^2)/sum((y - mean(y)).^2)

Rsq2 = 
   0.838210531103428

Невязки и качество подгонки

Остаточные значения являются различием между наблюдаемыми величинами ответа (зависимая) переменная и значения, которые предсказывает модель. Когда вы подбираете модель, которая подходит для ваших данных, остаточные значения аппроксимируют независимые случайные ошибки. Таким образом, распределение остаточных значений не должно показывать заметный шаблон.

Создание подгонки с помощью линейной модели требует минимизации суммы квадратов остаточных значений. Эта минимизация дает к тому, что называется методом наименьших квадратов. Можно получить сведения о “качестве” соответствия путем визуального исследования графика невязок. Если остаточный график имеет шаблон (то есть, остаточные точки данных, кажется, не имеют случайное рассеяние), случайность указывает, что модель правильно не соответствует данным.

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

Одной мерой качества подгонки является coefficient of determination или R² (объявленный r-квадратом). Эта статистическая величина указывает, как тесно значения, которые вы получаете из подбора кривой модели, совпадают с зависимой переменной, которую модель предназначается, чтобы предсказать. Статистики часто задают R² использование остаточного отклонения от подобранной модели:

R² = 1 – _SSresid / _SStotal

_SSresid является суммой квадратов остатков от регрессии. _SStotal является суммой различий в квадрате от среднего значения зависимой переменной (полная сумма квадратов). Оба - положительные скалярные величины.

Изучить, как вычислить R² когда вы будете использовать инструмент Basic Fitting, см. R2, Коэффициент детерминации. Узнать больше о вычислении R² статистическая величина и ее многомерное обобщение, продолжите читать здесь.

Пример: вычисление R² от аппроксимаций полиномом

Можно вывести R² от коэффициентов параболической регрессии, чтобы определить, сколько отклонения в y линейная модель объясняет, как следующий пример описывает:

Создайте две переменные, x и y, из первых двух столбцов count переменная в файле данных count.dat:
```
load count.dat
x = count(:,1);
y = count(:,2);
```
Использование polyfit вычислить линейную регрессию, которая предсказывает y от x:
```
p = polyfit(x,y,1)

p =
    1.5229   -2.1911
```
p(1) наклон и p(2) прерывание линейного предиктора. Можно также получить коэффициенты регрессии с помощью Пользовательского интерфейса Basic Fitting.
Вызвать polyval использовать p предсказать y, вызов результата yfit:
```
yfit = polyval(p,x);
```
Используя polyval избавляет вас от ввода уравнения подгонки самостоятельно, который в этом случае похож:
```
yfit =  p(1) * x + p(2);
```
Вычислите остаточные значения как вектор из чисел со знаком:
```
yresid = y - yfit;
```
Придайте остаточным значениям квадратную форму и составьте их, чтобы получить остаточную сумму квадратов:
```
SSresid = sum(yresid.^2);
```
Вычислите полную сумму квадратов y путем умножения отклонения y количеством наблюдений минус 1:
```
SStotal = (length(y)-1) * var(y);
```
Компьютер² использование формулы, данной во введении этой темы:
```
rsq = 1 - SSresid/SStotal

rsq =
    0.8707
```
Это демонстрирует что линейное уравнение 1.5229 * x -2.1911 предсказывает 87% отклонения в переменной y.

Вычисление настроенный R² для параболических регрессий

Обычно можно уменьшить остатки в модели, установив полином более высокой степени. Когда вы добавляете больше терминов, вы увеличиваете коэффициент детерминации, R². Вы получаете более близкую подгонку к данным, но за счет более сложной модели, для который R² не может считать. Однако улучшение этой статистической величины, настроенного R², действительно включает штраф за количество терминов в модели. Настроенный R², поэтому, более подходит для сравнения, как различные модели соответствуют к тем же данным. Настроенный R² задан как:

R²_{настроенный} = 1 - (_SSresid / _SStotal) * ((n-1) / (n-d-1))

где n является количеством наблюдений в ваших данных, и d является степенью полинома. (Линейная подгонка имеет степень 1, квадратичные подходящие 2, кубическое соответствие 3, и так далее.)

Следующий пример повторяет шаги предыдущего примера, Примера: Вычисление R2 от Аппроксимаций полиномом, но выполняет кубическое (степень 3) подгонка вместо линейного (степень 1) подгонка. От кубического соответствия вы вычисляете и простой и настроенный R² значения, чтобы оценить, улучшают ли дополнительные термины предсказательную силу:

Создайте две переменные, x и y, из первых двух столбцов count переменная в файле данных count.dat:
```
load count.dat
x = count(:,1);
y = count(:,2);
```
Вызвать polyfit сгенерировать кубическое соответствие, чтобы предсказать y от x:
```
p = polyfit(x,y,3)

p =
   -0.0003    0.0390    0.2233    6.2779
```
p(4) точка пересечения кубического предиктора. Можно также получить коэффициенты регрессии с помощью Пользовательского интерфейса Basic Fitting.
Вызвать polyval использовать коэффициенты в p предсказать y, именование результата yfit:
```
yfit = polyval(p,x);
```
polyval оценивает явное уравнение, которое вы могли вручную ввести как:
```
yfit =  p(1) * x.^3 + p(2) * x.^2 + p(3) * x + p(4);
```
Вычислите остаточные значения как вектор из чисел со знаком:
```
yresid = y - yfit;
```
Придайте остаточным значениям квадратную форму и составьте их, чтобы получить остаточную сумму квадратов:
```
SSresid = sum(yresid.^2);
```
Вычислите полную сумму квадратов y путем умножения отклонения y количеством наблюдений минус 1:
```
SStotal = (length(y)-1) * var(y);
```
Вычислите простой R² для кубического соответствия с помощью формулы, данной во введении этой темы:
```
rsq = 1 - SSresid/SStotal

rsq =
    0.9083
```
Наконец, вычислите настроенный R² с учетом степеней свободы:
```
rsq_adj = 1 - SSresid/SStotal * (length(y)-1)/(length(y)-length(p))

rsq_adj =
    0.8945
```
Настроенный R², 0.8945, меньше, чем простой R². 9083. Это обеспечивает более надежную оценку степени вашей полиномиальной модели предсказать.

Во многих моделях параболической регрессии, добавляя термины в уравнение увеличивает обоих R² и настроенный R²В предыдущем примере использование кубической подгонки увеличило обе статистики по сравнению с линейной подгонкой. (Можно вычислить настроенный R² для линейного пригодного для себя, чтобы продемонстрировать, что это имеет нижнее значение.) Однако это не всегда верно, что линейная подгонка хуже, чем подгонка высшего порядка: более сложная подгонка может иметь ниже настроенный R² чем более простая подгонка, указывая, что увеличенная сложность не выравнивается по ширине. Кроме того, в то время как R² всегда находится в интервале между 0 и 1 для моделей параболической регрессии, которые генерирует инструмент Basic Fitting, настроил R² поскольку некоторые модели могут быть отрицательными, указав, что модель, которая имеет слишком много терминов.

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

Подбор кривой данным с функциями Curve Fitting Toolbox

Программное обеспечение Curve Fitting Toolbox расширяет базовую функциональность MATLAB путем включения следующих соответствующих данным возможностей:

Линейный и нелинейный параметрический подбор кривой, включая стандартный линейный метод наименьших квадратов, нелинейный метод наименьших квадратов, метод взвешенных наименьших квадратов, метод наименьших квадратов с ограничениями и робастные методы подбора
Непараметрический подбор кривой
Статистика для определения качества подгонки
Экстраполяция, дифференцирование и интегрирование
Диалоговое окно, которое упрощает разделение данных и сглаживание
Сохранение результатов подгонки в различных форматах, включая файлы кода MATLAB, MAT-файлы и переменные рабочей области

Для получения дополнительной информации см. документацию Curve Fitting Toolbox.

Документация