Определитель матрицы ▦ Производная функции dydx График 3D Упростить выражение Графический метод решения задач нелинейного программирования ⇲
Примеры решений Коэффициент Спирмена Коэффициент Фехнера Множественная регрессия Нелинейная регрессия Уравнение регрессии Автокорреляция Расчет параметров тренда Ошибка аппроксимации

Пример нахождения коэффициента корреляции

Коэффициент корреляции
где x·y, x, y - средние значения выборок; σ(x), σ(y) - среднеквадратические отклонения.
Кроме того, коэффициент линейной парной корреляции Пирсона может быть определен через коэффициент регрессии b:
коэффициент линейной парной корреляции может быть определен через коэффициент регрессии
где σ(x)=S(x), σ(y)=S(y) - среднеквадратические отклонения, b - коэффициент перед x в уравнении регрессии y=a+bx.

Другие варианты формул:
или

Кxy - корреляционный момент (коэффициент ковариации)
корреляционный момент

Для нахождения линейного коэффициента корреляции Пирсона необходимо найти выборочные средние x и y, и их среднеквадратические отклонения σx = S(x), σy = S(y):

Линейный коэффициент корреляции указывает на наличие связи и принимает значения от –1 до +1 (см. шкалу Чеддока). Например, при анализе тесноты линейной корреляционной связи между двумя переменными получен коэффициент парной линейной корреляции, равный –1. Это означает, что между переменными существует точная обратная линейная зависимость.
Вычислить значение коэффициента корреляции можно по заданным средним выборки, либо непосредственно по исходным табличным данным.
Линейный коэффициент корреляции Пирсона
Геометрический смысл коэффициента корреляции: rxy показывает, насколько различается наклон двух линий регрессии: y(x) и х(у), насколько сильно различаются результаты минимизации отклонений по x и по y. Чем больше угол между линиями, то тем больше rxy.
Знак коэффициента корреляции совпадает со знаком коэффициента регрессии и определяет наклон линии регрессии, т.е. общую направленность зависимости (возрастание или убывание). Абсолютная величина коэффициента корреляции определяется степенью близости точек к линии регрессии.

Свойства коэффициента корреляции

  1. |rxy| ≤ 1;, -1≤x≤1
  2. если X и Y независимы, то rxy=0, обратное не всегда верно;
  3. если |rxy|=1, то Y=aX+b, |rxy(X,aX+b)|=1, где a и b постоянные, а ≠ 0;
  4. |rxy(X,Y)|=|rxy(a1X+b1, a2X+b2)|, где a1, a2, b1, b2 – постоянные.

Поэтому для проверки направления связи выбирается проверка гипотезы при помощи коэффициента корреляции Пирсона с дальнейшей проверкой на достоверность при помощи t-критерия (пример см. ниже).

Количество строк (исходных данных)

Заданы итоговые значения величин (∑x, ∑x2, ∑xy, ∑y, ∑y2)
Здесь будет отображаться решение.
Назначение сервиса. С помощью данного онлайн-калькулятора можно найти параметры уравнения нелинейной регрессии (экспоненциальной, степенной, равносторонней гиперболы, логарифмической, показательной) (см. пример).

Инструкция. Укажите количество исходных данных. Полученное решение сохраняется в файле Word. Также автоматически создается шаблон решения в Excel.

Примечание: если необходимо определить параметры параболической зависимости (y = ax2 + bx + c), то можно воспользоваться сервисом Аналитическое выравнивание.
Ограничить однородную совокупность единиц, устранив аномальные объекты наблюдения можно через метод Ирвина или по правилу трех сигм (устранить те единицы, для которых значение объясняющего фактора отклоняется от среднего более, чем на утроенное среднеквадратичное отклонение).

Пример №1. На основе данных, приведенных в Приложении 1 и соответствующих Вашему варианту (таблица 2), требуется:

  1. Рассчитать коэффициент линейной парной корреляции и построить уравнение линейной парной регрессии одного признака от другого. Один из признаков, соответствующих Вашему варианту, будет играть роль факторного (х), другой – результативного (y). Причинно-следственные связи между признаками установить самим на основе экономического анализа. Пояснить смысл параметров уравнения.
  2. Определить теоретический коэффициент детерминации и остаточную (необъясненную уравнением регрессии) дисперсию. Сделать вывод.
  3. Оценить статистическую значимость уравнения регрессии в целом на пятипроцентном уровне с помощью F-критерия Фишера. Сделать вывод.
  4. Выполнить прогноз ожидаемого значения признака-результата y при прогнозном значении признака-фактора х, составляющим 105% от среднего уровня х. Оценить точность прогноза, рассчитав ошибку прогноза и его доверительный интервал с вероятностью 0,95.
Решение. Уравнение имеет вид y = ax + b
Средние значения



Дисперсия

Среднеквадратическое отклонение


Коэффициент корреляции
Связь между признаком Y фактором X сильная и прямая (определяется по шкале Чеддока).
Уравнение регрессии
Коэффициент регрессии: k = a = 4.01
Коэффициент детерминации
R 2 = 0.99 2 = 0.97, т.е. в 97% случаев изменения х приводят к изменению y. Другими словами - точность подбора уравнения регрессии - высокая. Остаточная дисперсия: 3%.
xyx2y2x·yy(x)(yi-y)2(y-y(x))2(x-xp)2
1107111449107103.19333.0614.530.25
2109411881218107.2264.063.2320.25
3110912100330111.21232.561.4712.25
41131612769452115.22150.064.956.25
51202514400600119.2327.560.592.25
61223614884732123.2410.561.550.25
71234915129861127.265.0618.110.25
812864163841024131.277.5610.672.25
913681184961224135.28115.560.526.25
10140100196001400139.29217.560.5112.25
11145121210251595143.3390.062.920.25
12150144225001800147.31612.567.2530.25
7815036501906171034315032366.2566.23143

Примечание: значения y(x) находятся из полученного уравнения регрессии:
y(1) = 4.01*1 + 99.18 = 103.19
y(2) = 4.01*2 + 99.18 = 107.2
... ... ...

Значимость коэффициента корреляции

Выдвигаем гипотезы:
H0: rxy = 0, нет линейной взаимосвязи между переменными;
H1: rxy ≠ 0, есть линейная взаимосвязь между переменными;
Для того чтобы при уровне значимости α проверить нулевую гипотезу о равенстве нулю генерального коэффициента корреляции нормальной двумерной случайной величины при конкурирующей гипотезе H1 ≠ 0, надо вычислить наблюдаемое значение критерия (величина случайной ошибки):
По таблице Стьюдента находим tтабл (n-m-1;α/2) = (10;0.025) = 2.228
Поскольку Tнабл > tтабл, то отклоняем гипотезу о равенстве 0 коэффициента корреляции. Другими словами, коэффициент корреляции статистически - значим.
Интервальная оценка для коэффициента корреляции (доверительный интервал)
Интервальная оценка для коэффициента корреляции

r - Δr ≤ r ≤ r + Δr
Δr = ±tтаблmr = ±2.228 • 0.0529 = 0.118
0.986 - 0.118 ≤ r ≤ 0.986 + 0.118
Доверительный интервал для коэффициента корреляции: 0.868 ≤ r ≤ 1

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





Sa=0.2152

Доверительные интервалы для зависимой переменной

Рассчитаем границы интервала, в котором будет сосредоточено 95% возможных значений Y при неограниченно большом числе наблюдений и X = 7
(122.4;132.11)
Проверка гипотез относительно коэффициентов линейного уравнения регрессии

1) t-статистика


Статистическая значимость коэффициента регрессии подтверждается (18.63>2.228).

Статистическая значимость коэффициента регрессии подтверждается (62.62>2.228).
Доверительный интервал для коэффициентов уравнения регрессии
Определим доверительные интервалы коэффициентов регрессии, которые с надежность 95%  будут следующими (tтабл=2.228):
(a - tтабл·Sa; a + tтабл·S a)
(3.6205;4.4005)
(b - tтабл·Sb; b + tтабл·Sb)
(96.3117;102.0519)

2) F-статистики


Fkp = 4.96. Поскольку F > Fkp, то коэффициент детерминации статистически значим (см. критерий Фишера).

см. также Корреляционный анализ. Примеры решения задач.

Пример №2
1. Расчет средних значений x, y:
x = ∑xin = 660.611 = 60.05
y = ∑yin = 333.9411 = 30.36
x·y=∑xi·yin=19952.0711= 1813.82
2. Расчет дисперсий:
S2(x) =xi2n-x2=40337.211- 60.052= 60.47
S2(y) =yi2n-y2=10329.5211- 30.362= 17.43
3. Расчет среднеквадратических отклонений:
S(x) =S2(x)=60.47= 7.78
S(y) =S2(y)=17.43= 4.17
4. Расчет линейного коэффициента корреляции Пирсона:
rxy=x·y-x·yS(x)·S(y)=1813.82-60.05·30.367.78·4.17= -0.2872 Линейный коэффициент корреляции принимает значения от –1 до +1.
Связи между признаками могут быть слабыми и сильными (тесными). Их критерии оцениваются по шкале Чеддока:
0.1 < rxy < 0.3: слабая;
0.3 < rxy < 0.5: умеренная;
0.5 < rxy < 0.7: заметная;
0.7 < rxy < 0.9: высокая;
0.9 < rxy < 1: весьма высокая;
В нашем примере связь между признаком Y фактором X слабая и обратная.

x yx2y2x·yy(x)(yi-y)2(y-y(x))2
68.5 22.39 4692.25 501.31 1533.72 29.06 63.49 44.44
75.7 29.24 5730.49 854.98 2213.47 27.95 1.25 1.67
52.7 32.92 2777.29 1083.73 1734.88 31.49 6.56 2.04
60.2 33.52 3624.04 1123.59 2017.9 30.34 10 10.14
62.3 30.98 3881.29 959.76 1930.05 30.01 0.39 0.94
48.3 37.17 2332.89 1381.61 1795.31 32.17 46.4 25
56.5 32.12 3192.25 1031.69 1814.78 30.91 3.1 1.47
65.9 31.76 4342.81 1008.7 2092.98 29.46 1.97 5.3
56.2 28.48 3158.44 811.11 1600.58 30.95 3.53 6.11
51.1 23.17 2611.21 536.85 1183.99 31.74 51.67 73.42
63.2 32.19 3994.24 1036.2 2034.41 29.87 3.36 5.37
660.6333.9440337.210329.5219952.07333.94191.71175.9

Значимость линейного коэффициента корреляции Пирсона.
tнабл= rxy·n-21-rxy2= 0.2872·91-0.28722= 0.9
По таблице Стьюдента с уровнем значимости α=0.05 и степенями свободы k=n-m-1=11-1-1=9 находим tкрит: tкрит(n-m-1;α/2) = tкрит(9;0.025) = 2.262, где m=1 - количество объясняющих переменных.
Если tнабл > tкритич, то полученное значение коэффициента корреляции Пирсона признается значимым (нулевая гипотеза, утверждающая равенство нулю коэффициента корреляции, отвергается).
Поскольку tнабл < tкрит, то принимаем гипотезу о равенстве 0 коэффициента корреляции. Другими словами, коэффициент корреляции статистически - не значим
В парной линейной регрессии t2r = t2b и тогда проверка гипотез о значимости коэффициентов регрессии и корреляции равносильна проверке гипотезы о существенности линейного уравнения регрессии.

Интервальная оценка для линейного коэффициента корреляции Пирсона
(rxy- tкрит·1-rxy2n;rxy+ tкрит·1-rxy2n)
Доверительный интервал для коэффициента корреляции
(0.29 - 2.262·1-0.29211;0.29 + 2.262·1-0.29211)
Доверительный интервал для линейного коэффициента корреляции Пирсона: r(-0.9129;0.3386)

Перейти к онлайн решению своей задачи

Пример нахождения доверительных интервалов коэффициентов регрессии

В таблице 5 представлена динамика российского нефтяного экспорта и цены нефти марки 11га1з в 1997-2003 гг.

Задание:
1. Постройте поле корреляции и сформулируйте гипотезу о форме связи.
2. Постройте уравнение зависимости экспорта нефти от цены на нефть.
3. Рассчитайте среднюю ошибку аппроксимации и коэффициент детерминации. Оценить статистическую значимость параметров регрессии и уравнения в целом.
4. Оцените полученные результаты, выводы оформите в аналитической записке.
Таблица 5

ГодыЦена нефти марки Urals (Россия), долл/барр.Экспорт нефти и нефтепродуктов, млн.т.
199718,3360,6
199811,8353,8
199917,3056,9
200026,6361,9
200122,9770,8
200223,7375,0
200327,0476,4

Решение:

Уравнение имеет вид y = ax + b
1. Параметры уравнения регрессии.
Средние значения


Дисперсия


Среднеквадратическое отклонение
Среднеквадратическое отклонение

Коэффициент корреляции.
Коэффициент корреляции
Связь между признаком Y фактором X сильная и прямая
Уравнение регрессии

Уравнение регрессии

Коэффициент детерминации
R2= 0.922 = 0.84, т.е. в 84% случаев изменения х приводят к изменению y. Другими словами - точность подбора уравнения регрессии - высокая
xyx2y2x ∙ yy(x)(y-y)2(y-y(x))2(x-xp)2
119298.121416188875.5335476.28219.63232120.86160.5624362.01
203481.0341209231389.8697649.09521.1689328.761610.265196.01
281539.1278961290650.37151492.72801.1557979.4268658.5135.01
305653.5793025427153.74199338.85887.315961.5954628.94895.01
381987.66145161975472.28376298.461160.1143160.4129738.5711218.34
3631252.851317691569633.12454784.551095.5223673.0324760.357729.34
3891276.881513211630422.53496706.321188.83246980.017753.5712977.01
3871396.701497691950770.89540522.91181.65380430.9346248.0412525.34
315952.0399225906361.12299889.45923.1929625.58831.491593.34
217619.9647089384350.4134531.32571.4125583.742356.853373.67
149384.4022201147763.3657275.6327.32156427.53258.2315897.01
192516.5936864266865.2399185.28481.6769336.981219.246902.84
33019358.9110107558869708.452943150.829358.911570608.75247224.62102704.92

2. Оценка параметров уравнения регрессии
Значимость коэффициента корреляции

По таблице Стьюдента находим Tтабл (n-m-1;a) = (10;0.05) = 1.812
Поскольку Tнабл > Tтабл , то отклоняем гипотезу о равенстве 0 коэффициента корреляции. Другими словами, коэффициента корреляции статистически - значим.

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





S a = 0.4906
Доверительные интервалы для зависимой переменной

Рассчитаем границы интервала, в котором будет сосредоточено 95% возможных значений Y при неограниченно большом числе наблюдений и X = 1
(-587.75;179.86)
Проверка гипотез относительно коэффициентов линейного уравнения регрессии
1) t-статистика


Статистическая значимость коэффициента регрессии a подтверждается (7.32>1.812)

Статистическая значимость коэффициента регрессии b не подтверждается (1.46<1.812)
Доверительный интервал для коэффициентов уравнения регрессии
Определим доверительные интервалы коэффициентов регрессии, которые с надежность 95%  будут следующими (tтабл=1.812):
(a - tтабл·S a; a + tтабл·Sa)
(2.7006;4.4786)
(b - tтабл·S b; b + tтабл·Sb)
(-465.5454;50.4796)
2) F-статистики


Fkp = 4.96
Поскольку F > Fkp, то коэффициент детерминации статистически значим.

Перейти к онлайн решению своей задачи

Доверительные интервалы для зависимой переменной

Уравнение тренда имеет вид y = at2 + bt + c
1. Находим параметры уравнения методом наименьших квадратов.
Система уравнений

Для наших данных система уравнений имеет вид (см. таблицу).

Получаем a0 = -11.37, a1 = 88.47, a2 = 2151.09
Уравнение тренда: y = -11.37t2+88.47t+2151.09
Оценим качество уравнения тренда с помощью ошибки абсолютной аппроксимации.


Поскольку ошибка больше 15%, то данное уравнение не желательно использовать в качестве тренда
Средние значения



Дисперсия

Среднеквадратическое отклонение

Индекс детерминации


т.е. в 87.35 % случаев влияет на изменение данных. Другими словами - точность подбора уравнения тренда - высокая

tyt2y2x ∙ yy(t)(y-ycp)2(y-y(t))2(t-tp)2(y-y(t)) : yt3t4t2y
12225.314951960.092225.32228.1965.60998.352166431.117112225.3
22254.945084574.014509.82282.55462.25764.5225962347.9858169019.6
32332.395439623.296996.92314.179781.21328.6969442284.599278120990.7
42365.8165597009.649463.22323.0517529.761827.56251101137.956425637852.8
52295.4255268861.16114772309.193844190.1641031653.56612562557385
62303.9365307955.2113823.42272.594970.25980.3161172135.109216129682940.4
72166.7494694588.8915166.92213.254448.892166.90254100859.8853432401106168.3
82080.4644328064.1616643.22131.17234092577.59299105621.9085124096133145.6
92075.9814309360.8118683.12026.3524806.252455.202516102860.8457296561168147.9
4520100.628544981997.2698988.820100.5189317.219911299.31260625332.9644050306661235751.2

2. Анализ точности определения оценок параметров уравнения тренда.

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




S a = 4.8518
Доверительные интервалы для зависимой переменной
Доверительные интервалы для зависимой переменной

По таблице Стьюдента находим Tтабл (n-m-1;a) = (7;0.05) = 1.895
Рассчитаем границы интервала, в котором будет сосредоточено 95% возможных значений Y при неограниченно большом числе наблюдений и t = 6
2151.09 + 88.47*6 + -11.37*62 - 1.895*39.911 ; 2151.09 + 88.47*6 + -11.37*62 - 1.895*39.911
(-55.3814;95.8814)
Интервальный прогноз.
Определим среднеквадратическую ошибку прогнозируемого показателя.


где L - период упреждения; уn+L - точечный прогноз по модели на (n + L)-й момент времени; n - количество наблюдений во временном ряду; Sy - стандартная ошибка прогнозируемого показателя; Tтабл - табличное значение критерия Стьюдента для уровня значимости а и для числа степеней свободы, равного n — 2.
Точечный прогноз, t = 10: y(10) = -11.37*102 + 88.47* + 2151.09 = 1898.79
K1 = 247.4924
1898.79 - 247.4924 = 1651.2976 ; 1898.79 + 247.4924 = 2146.2824
t = 10: (1651.2976;2146.2824)
Точечный прогноз, t = 11: y(11) = -11.37*112 + 88.47* + 2151.09 = 1748.49
K2 = 261.9213
1748.49 - 261.9213 = 1486.5687 ; 1748.49 + 261.9213 = 2010.4113
t = 11: (1486.5687;2010.4113)
Точечный прогноз, t = 12: y(12) = -11.37*122 + 88.47* + 2151.09 = 1575.45
K3 = 278.0099
1575.45 - 278.0099 = 1297.4401 ; 1575.45 + 278.0099 = 1853.4599
t = 12: (1297.4401;1853.4599)
Точечный прогноз, t = 13: y(13) = -11.37*132 + 88.47* + 2151.09 = 1379.67
K4 = 295.4871
1379.67 - 295.4871 = 1084.1829 ; 1379.67 + 295.4871 = 1675.1571
t = 13: (1084.1829;1675.1571)
Точечный прогноз, t = 14: y(14) = -11.37*142 + 88.47* + 2151.09 = 1161.15
K5 = 314.1213
1161.15 - 314.1213 = 847.0287 ; 1161.15 + 314.1213 = 1475.2713
t = 14: (847.0287;1475.2713)
3. Проверка гипотез относительно коэффициентов линейного уравнения тренда.
1) t-статистика. Критерий Стьюдента.


Статистическая значимость коэффициента уравнения подтверждается

Статистическая значимость коэффициента тренда подтверждается
Доверительный интервал для коэффициентов уравнения тренда
Определим доверительные интервалы коэффициентов тренда, которые с надежность 95% будут следующими (tтабл=1.895):
(a - tтабл·Sa; a + tтабл·Sa)
(-20.5642;-2.1758)
(b - t табл·Sb; b + tтаблb)
(36.7313;140.2087)
2) F-статистика. Критерий Фишера.


Fkp = 5.32
Поскольку F > Fkp, то коэффициент детерминации статистически значим
4. Тест Дарбина-Уотсона на наличие автокорреляции остатков для временного ряда.
yy(x)ei= y-y(x)e2(ei- ei-1)2
2225.32228.19-2.898.35210
2254.92282.55-27.65764.5225613.0576
2332.32314.1718.13328.69692095.8084
2365.82323.0542.751827.5625606.1444
2295.42309.19-13.79190.16413196.7716
2303.92272.5931.31980.31612034.01
2166.72213.25-46.552166.90256062.1796
2080.42131.17-50.772577.592917.8084
2075.92026.3549.552455.202510064.1024
11299.312124689.8824



Критические значения d1 и d2 определяются на основе специальных таблиц для требуемого уровня значимости a, числа наблюдений n и количества объясняющих переменных m.
Не обращаясь к таблицам, можно пользоваться приблизительным правилом и считать, что автокорреляция остатков отсутствует, если 1.5 < DW < 2.5. Для более надежного вывода целесообразно обращаться к табличным значениям.
d1 < DW и d2 < DW < 4 - d2.