эконометрика. На уровне значимости а 0,05 проверить гипотезу о согласии линейной регрессии с результатами наблюдений
Скачать 160.15 Kb.
|
Задача 1 Найти оценки параметров линейной регрессии у на х. Построить диаграмму рассеяния и нанести прямую регрессии на диаграмму рассеяния. На уровне значимости а = 0,05 проверить гипотезу о согласии линейной регрессии с результатами наблюдений. С надежностью p = 0,95 найти доверительные интервалы для параметров линейной регрессии. Таблица
Решение: Метод наименьших квадратов дает наилучшие (состоятельные, эффективные и несмещенные) оценки параметров уравнения регрессии. Но только в том случае, если выполняются определенные предпосылки относительно случайного члена (ε) и независимой переменной (x). Формально критерий МНК можно записать так: = ∑(yi - y*i)2 → min Система нормальных уравнений. a•n + b∑x = ∑y∑x + b∑x2 = ∑y•x Для наших данных система уравнений имеет вид a + 724 b = 492 a + 40134 b = 26907 Домножим уравнение (1) системы на (-51.71), получим систему, которую решим методом алгебраического сложения. -724a -37438.04 b = -25441.32 a + 40134 b = 26907 Получаем: .96 b = 1465.68 Откуда b = 0.5435 Теперь найдем коэффициент «a» из уравнения (1): a + 724 b = 492 a + 724 • 0.5435 = 492 a = 98.51= 7.0361 Получаем эмпирические коэффициенты регрессии: b = 0.5435, a = 7.0361 Уравнение регрессии (эмпирическое уравнение регрессии):= 0.5435 x + 7.0361 Таблица. Для расчета параметров регрессии построим расчетную таблицу
Линейный коэффициент корреляции принимает значения от -1 до +1. Связи между признаками могут быть слабыми и сильными (тесными). Их критерии оцениваются по шкале Чеддока: .1 .3 .5 .7 .9 Линейное уравнение регрессии имеет вид y = 0.54 x + 7.04 Коэффициентам уравнения линейной регрессии можно придать экономический смысл. Коэффициент регрессии b = 0.54 показывает среднее изменение результативного показателя (в единицах измерения у) с повышением или понижением величины фактора х на единицу его измерения. В данном примере с увеличением на 1 единицу y повышается в среднем на 0.54. Коэффициент a = 7.04 формально показывает прогнозируемый уровень у, но только в том случае, если х=0 находится близко с выборочными значениями. Но если х=0 находится далеко от выборочных значений х, то буквальная интерпретация может привести к неверным результатам, и даже если линия регрессии довольно точно описывает значения наблюдаемой выборки, нет гарантий, что также будет при экстраполяции влево или вправо. Подставив в уравнение регрессии соответствующие значения х, можно определить выровненные (предсказанные) значения результативного показателя y (x) для каждого наблюдения. Связь между у и х определяет знак коэффициента регрессии b (если > 0 - прямая связь, иначе - обратная). В нашем примере связь прямая. Коэффициент детерминации. Квадрат (множественного) коэффициента корреляции называется коэффициентом детерминации, который показывает долю вариации результативного признака, объясненную вариацией факторного признака. Для оценки качества параметров регрессии построим расчетную таблицу
Значимость коэффициента корреляции. Выдвигаем гипотезы:0: rxy = 0, нет линейной взаимосвязи между переменными;1: rxy ≠ 0, есть линейная взаимосвязь между переменными; По таблице Стьюдента с уровнем значимости α=0.05 и степенями свободы k=12 находим tкрит: крит (n-m-1;α/2) = (12;0.025) = 2.179 где m = 1 - количество объясняющих переменных. Если |tнабл| >tкритич, то полученное значение коэффициента корреляции признается значимым (нулевая гипотеза, утверждающая равенство нулю коэффициента корреляции, отвергается). Поскольку |tнабл| >tкрит, то отклоняем гипотезу о равенстве 0 коэффициента корреляции. Другими словами, коэффициент корреляции статистически - значим В парной линейной регрессии t2r = t2b и тогда проверка гипотез о значимости коэффициентов регрессии и корреляции равносильна проверке гипотезы о существенности линейного уравнения регрессии. Интервальная оценка для коэффициента корреляции (доверительный интервал). Доверительный интервал для коэффициента корреляции. (0.813;1) Индивидуальные доверительные интервалы для Y при данном значении X. (a + bxi ± ε) где крит (n-m-1;α/2) = (12;0.025) = 2.179 Таблица
С вероятностью 95% можно гарантировать, что значения Y при неограниченно большом числе наблюдений не выйдет за пределы найденных интервалов. Проверка гипотез относительно коэффициентов линейного уравнения регрессии. ) t-статистика. Критерий Стьюдента. крит (n-m-1;α/2) = (12;0.025) = 2.179 Поскольку 13.51 > 2.179, то статистическая значимость коэффициента регрессии b подтверждается (отвергаем гипотезу о равенстве нулю этого коэффициента). Поскольку 3.27 > 2.179, то статистическая значимость коэффициента регрессии a подтверждается (отвергаем гипотезу о равенстве нулю этого коэффициента). Доверительный интервал для коэффициентов уравнения регрессии. Определим доверительные интервалы коэффициентов регрессии, которые с надежность 95% будут следующими: (b - tкритSb; b + tкритSb) (0.54 - 2.179 • 0.0402; 0.54 + 2.179 • 0.0402) (0.456;0.631) С вероятностью 95% можно утверждать, что значение данного параметра будут лежать в найденном интервале. (a - tкритSa; a + tкрит Sa) (7.036 - 2.179 • 2.15; 7.036 + 2.179 • 2.15) (2.344;11.728) С вероятностью 95% можно утверждать, что значение данного параметра будут лежать в найденном интервале. F-статистика. Критерий Фишера. где m - число факторов в модели. Оценка статистической значимости парной линейной регрессии производится по следующему алгоритму: . Выдвигается нулевая гипотеза о том, что уравнение в целом статистически незначимо: H0: R2=0 на уровне значимости α. 2. Далее определяют фактическое значение F-критерия: где m=1 для парной регрессии. Задача 2 Исследуется зависимость производительности труда (Yi) от уровня механизации работ (Xi %) и среднего возраста работников (Xi лет) по данным 14 промышленных предприятий (i - порядковый номер предприятия). Статистические данные приведены в таблице . Требуется: Вычислить ковариации и составить ковариационную матрицу. Найти оценки параметров множественной линейной регрессии и составить уравнение плоскости регрессии у = b0+b1x +b2x На уровне значимости а = 0,05 проверить гипотезу о согласии линейной множественной регрессии с результатом наблюдений. С надежностью p = 0,95 найти доверительные интервалы для параметров множественной линейной регрессии. 2.1. Таблица
Решение: Определим вектор оценок коэффициентов регрессии. Согласно методу наименьших квадратов, вектор s получается из выражения: s = (XTX)-1XTY Таблица
Матрица Y
Таблица. Матрица XT
В матрице, (XTX) число 14, лежащее на пересечении 1-й строки и 1-го столбца, получено как сумма произведений элементов 1-й строки матрицы XT и 1-го столбца матрицы X Умножаем матрицы, (XTY) Находим обратную матрицу (XTX)-1
Вектор оценок коэффициентов регрессии равен Уравнение регрессии (оценка уравнения регрессии) = 1.74 + 0.53X1 + 0.16X2 Матрица парных коэффициентов корреляции R. Число наблюдений n = 14. Число независимых переменных в модели равно 2, а число регрессоров с учетом единичного вектора равно числу неизвестных коэффициентов. С учетом признака Y, размерность матрицы становится равным 4. Матрица, независимых переменных Х имеет размерность (14 х 4). Матрица, составленная из Y и X
Транспонированная матрица.
Матрица ATA.
Полученная матрица имеет следующее соответствие:
Найдем парные коэффициенты корреляции.
Матрица парных коэффициентов корреляции R:
Оценка дисперсии равна: e2 = (Y - X*Y(X))T(Y - X*Y(X)) = 46.76 Коэффициент детерминации R2 = 0.945 Коэффициент детерминации. 0.9722 = 0.945 Проверка гипотез относительно коэффициентов уравнения регрессии (проверка значимости параметров множественного уравнения регрессии). Число v = n - m - 1 называется числом степеней свободы. Считается, что при оценивании множественной линейной регрессии для обеспечения статистической надежности требуется, чтобы число наблюдений, по крайней мере, в 3 раза превосходило число оцениваемых параметров. t-статистика табл (n-m-1;α/2) = (11;0.025) = 2.201 (bi - tiSbi; bi + tiSbi) b0: (1.74 - 2.201 • 5.12 ; 1.74 + 2.201 • 5.12) = (-9.53;13.01)1: (0.53 - 2.201 • 0.0426 ; 0.53 + 2.201 • 0.0426) = (0.43;0.62)2: (0.16 - 2.201 • 0.14 ; 0.16 + 2.201 • 0.14) = (-0.15;0.47) Задача 3 Исследуется зависимость себестоимости единицы продукции (у тыс. р.) от объема произведенной продукции (х тыс. шт.) по данным 15 предприятий (г - порядковый номер предприятия). Статистические данные приведены в таблице. Требуется: Построить диаграмму рассеяния. Убедиться, что между себестоимостью и объемом произведенной продукции существует нелинейная связь. Считая, что регрессия у по х представляется многочленом второй степени, найти оценки параметров параболической регрессии и составить уравнение линии регрессии. Построить кривую регрессии и нанести ее на диаграмму рассеяния.
С помощью средств MS Excel нанесем точки рассеивания на координатную плоскость. Анализируя, расположение точек на диаграмме, можем утверждать наличие нелинейной связи между факторами. Составим уравнение регрессии |