| Скачать .docx |
Реферат: Парная регрессия
Контрольная работа
по теме: "Парная линейная регрессия"
Данные, характеризующие прибыль торговой компании "Все для себя" за первые 10 месяцев 2004 года (в тыс. руб.), даны в следующей таблице:
| январь | февраль | март | апрель | май | июнь | июль | август | сентябрь | октябрь |
| 367 | 418 | 412 | 470 | 485 | 470 | 525 | 568 | 538 | 558 |
В контрольной работе с использованием табличного процессора Ехсе l необходимо выполнить следующие вычисления и построения:
1. Построить диаграмму рассеяния.
2. Убедится в наличии тенденции (тренда) в заданных значениях прибыли фирмы и возможности принятия гипотезы о линейном тренде.
3. Построить линейную парную регрессию (регрессию вида
). Вычисление коэффициентов b
0
,
b
1
выполнить методом наименьших квадратов.
4. Нанести график регрессии на диаграмму рассеяния.
5. Вычислить значения статистики F и коэффициента детерминации R 2 . Проверить гипотезу о значимости построенного уравнения регрессии.
6. Вычислить выборочный коэффициент корреляции и проверить гипотезу о ненулевом его значении.
7. Вычислить оценку дисперсии случайной составляющей эконометрической модели.
8. Проверить гипотезы о значимости вычисленных коэффициентов b 0 , b 1 .
9. Построить доверительные интервалы для коэффициентов b 0 , b 1 .
10. Построить доверительные интервалы для дисперсии случайной составляющей эконометрической модели.
11. Построить доверительную область для условного математического ожидания М(
)( по оси Х откладывать месяцы январь - декабрь). Нанести границы этой области на диаграмму рассеяния.
12. С помощью линейной парной регрессии сделать прогноз величины прибыли на ноябрь и декабрь месяц и нанести эти значения на диаграмму рассеяния. Сопоставить эти значения с границами доверительной области для условного математического ожидания М(
)
и сделать вывод о точности прогнозирования с помощью построенной регрессионной модели.
Решение.
Используя исходные данные, строим диаграмму рассеяния:

На основе анализа диаграммы рассеяния убеждаемся в наличии тенденции увеличения прибыли фирмы и выдвигаем гипотезу о линейном тренде.
Полагаем, что связь между факторами Х
и У может быть описана линейной функцией
. Решение задачи нахождения коэффициентов b
0
,
b
1
основывается на применении метода наименьших квадратов исводится к решению системы двух линейных уравнений с двумя неизвестными b
0
,
b
1
:
b 0 n + b 1 Уxi = Уyi ,
b 0 Уxi + b 1 Уxi 2 = Уxi yi .
Составляем вспомогательную таблицу:
| № | х | y | x2 | ху | y2 |
| 1 | 1 | 367 | 1 | 367 | 134689 |
| 2 | 2 | 418 | 4 | 836 | 174724 |
| 3 | 3 | 412 | 9 | 1236 | 169744 |
| 4 | 4 | 470 | 16 | 1880 | 220900 |
| 5 | 5 | 485 | 25 | 2425 | 235225 |
| 6 | 6 | 470 | 36 | 2820 | 220900 |
| 7 | 7 | 525 | 49 | 3675 | 275625 |
| 8 | 8 | 568 | 64 | 4544 | 322624 |
| 9 | 9 | 538 | 81 | 4842 | 289444 |
| 10 | 10 | 558 | 100 | 5580 | 311364 |
| сумма | 55 | 4811 | 385 | 28205 | 2355239 |
Для нашей задачи система имеет вид:

Решение этой системы можно получить по правилу Крамера:

Получаем:
,
.
Таким образом, искомое уравнение регрессии имеет вид:
y =364,8 + 21,145x.
4. Нанесем график регрессии на диаграмму рассеяния.

5. Вычислим значения статистики F и коэффициента детерминации R 2 . Коэффициент детерминации рассчитаем по формуле R2 = rxy 2 = 0,9522 = 0,907. Проверим адекватность модели (уравнения регрессии) в целом с помощью F-критерия. Рассчитаем значение статистики F через коэффициент детерминации R2 по формуле:
![]()
Получаем:
. Зададим уровень значимости б =0,01, по таблице находим квантиль распределения Фишера F0,01;1;8
= 11,26, где 1 – число степеней свободы.
Fфакт. > F0,01;1;8 , т.к. 78,098 > 11,26.
Следовательно, делаем вывод о значимости уравнения регрессии при 99% - м уровне значимости.
6. Вычислим выборочный коэффициент корреляции и проверим гипотезу о ненулевом его значении.
Рассчитаем выборочный коэффициент корреляции по формуле:

Получаем:
![]()
Проверка существенности отличия коэффициента корреляции от нуля проводится по схеме:если
, то гипотеза о существенном отличии коэффициента корреляции от нуля принимается, в противном случае отвергается.
Здесь t1-б/2, n -2 – квантиль распределения Стьюдента, б - уровень значимости или уровень доверия, n – число наблюдений, (n-2) – число степеней свободы. Значение б задается. Примем б = 0,05, тогда t1-б/2, n -2 = t0,975,8 = 2,37. Получаем:
.
Следовательно, коэффициент корреляции существенно отличается от нуля и существует сильная линейная связь между х и у.
С использованием табличного процессора Ехсеl проведем регрессионную статистику:
Вывод итогов:
| Регрессионная статистика | |
| Множественный R | 0,952409 |
| R-квадрат | 0,907083 |
| Нормированный R-квадрат | 0,895468 |
| Стандартная ошибка | 21,7332 |
| Наблюдения | 10 |
| Дисперсионный анализ | |||||
| df | SS | MS | F | Значимость F | |
| Регрессия | 1 | 36888,245 | 36888,25 | 78,09816 | 2,119E-05 |
| Остаток | 8 | 3778,6545 | 472,3318 | ||
| Итого | 9 | 40666,9 | |||
| Коэфф. | Станд. ошибка | t-статистика | P-Значение | Нижние 95% | Верхние 95% | |
| Y-пересечение | 364,8 | 14,846599 | 24,57128 | 8,04E-09 | 330,56368 | 399,0363 |
| Переменная X 1 | 21,14545 | 2,3927462 | 8,837316 | 2,12E-05 | 15,627772 | 26,66314 |
Вычисленные значения коэффициентов b 0 , b 1 , значения статистики F , коэффициента детерминации R 2 выборочного коэффициента корреляции rxy совпадают с выделенными в таблице.
7. Оценка дисперсии случайной составляющей эконометрической модели вычисляется по формуле
.
Используя результаты регрессионной статистики, получаем:
.
8. Проверим значимость вычисленных коэффициентов b 0 , b 1 по t-критерию Стьюдента. Для этого проверяем выполнение неравенств:
и
,
где
,
,
,
.
Используем результаты регрессионной статистики:
| Коэффициенты | Стандартная ошибка | t-статистика | P-Значение | Нижние 95% | Верхние 95% | |
| Y-пересечение | 364,8 | 14,846599 | 24,57128 | 8,04E-09 | 330,56368 | 399,0363 |
| Переменная X 1 | 21,14545 | 2,3927462 | 8,837316 | 2,12E-05 | 15,627772 | 26,66314 |
Получаем:
;
Примем б = 0,05, тогда t1-б/2,
n
-2
= t0,975,8
= 2,37.
Так как
и
, делаем вывод о значимости коэффициентов линейного уравнения регрессии.
9. Доверительные интервалы для коэффициентов b 0 , b 1 получаем с помощью результатов регрессионной статистики.
Доверительный интервал для коэффициента b 0 уравнения регрессии:
![]()
Доверительный интервал для коэффициента b 1 уравнения регрессии:
![]()
10. Построим доверительный интервал для дисперсии случайной составляющей эконометрической модели по формуле:
.
Примем б = 0,05, тогда по таблице для 10-элементной выборки q = 0,65.
Получаем:
,
.
11. Построим доверительную область для условного математического ожидания М(
).
Доверительные интервалы для уравнения линейной регрессии:
находятся по формуле:
![]()
где
соответственно верхняя и нижняя границы доверительного интервала;
значение независимой переменной
для которого определяется доверительный интервал,
квантиль распределения Стьюдента,
доверительная вероятность, (n-2) – число степеней свободы;

![]()
Рассмотрим уравнение: y =364,8 + 21,145x. Пусть
тогда
. Зная
и
, заполним таблицу:
| 1 | 385,95 | 20,25 | 4,634 | 377,327 | 394,564 |
| 2 | 407,09 | 12,25 | 5,215 | 397,391 | 416,791 |
| 3 | 428,24 | 6,25 | 5,738 | 417,564 | 438,908 |
| 4 | 449,38 | 2,25 | 6,217 | 437,819 | 460,945 |
| 5 | 470,53 | 0,25 | 6,661 | 458,138 | 482,917 |
| 6 | 491,67 | 0,25 | 7,078 | 478,508 | 504,838 |
| 7 | 512,82 | 2,25 | 7,471 | 498,921 | 526,715 |
| 8 | 533,96 | 6,25 | 7,845 | 519,372 | 548,556 |
| 9 | 555,11 | 12,25 | 8,202 | 539,854 | 570,365 |
| 10 | 576,25 | 20,25 | 8,544 | 560,363 | 592,146 |
| сумма | 82,5 | ||||
| 11 | 597,4 | 30,25 | 8,873 | 580,897 | 613,903 |
| 12 | 618,55 | 42,25 | 9,190 | 601,453 | 635,638 |
График уравнения регрессии, доверительная полоса, диаграмма рассеяния:

12. С помощью линейной парной регрессии сделаем прогноз величины прибыли на ноябрь и декабрь месяц:
597,4,
618,55.
Нанесем эти значения на диаграмму рассеяния.

Эти значения сопоставимы с границами доверительной области для условного математического ожидания М(
).
Точность прогнозирования: с вероятностью 0,95 прибыль в ноябре находится в интервале (487,292; 515,508); прибыль в декабре находится в интервале (497,152; 526,376).