Задача 1.

1.1 Выяснить,  существует ли связь между потреблением дизельного топлива (y) и объёмом валовой продукции (). (Для этого построить поле рассеяния. На основе его визуального анализа выдвинуть гипотезу о виде статистической зависимости y от ). Найти точечные оценки неизвестных параметров модели. Выяснить, существует ли связь между потреблением дизельного топлива (y) и объёмом капитальных вложений (). Найти оценки неизвестных параметров модели.

Годы


Валовая продукция промышленности

( млрд. р.)

Объём капитальных вложений

( млрд. р.)

Объём потребления дизельного топлива

(y млн. т.)

1985

7.1

4.9

1.6

1986

6.9

4.8

1.5

1987

6.8

4.7

1.5

1988

6.8

4.6

1.5

1989

7

4.7

1.6

1990

7.1

4.9

1.6

1991

7.2

5.1

1.7

1992

7.5

5.4

1.7

1993

7.7

5.7

1.7

1994

7.8

5.8

1.8

1995

7.8

5.6

1.8

1996

7.7

5.5

1.6

1997

7.9

5.9

1.8

1998

8

5.8

1.9

1999

8

5.8

1.8

2000

7.9

6

1.8


Решение: По исходным данным построим поля рассеяния переменной у в зависимости от  и , нанесем линии тренда и эллипсы рассеяния

Вид полей и эллипсов рассеивания позволяют выдвинуть гипотезу о том, что зависимость потребления дизельного топлива (y) от объёмов продукции () описывается линейной моделью вида:

где  и  - неизвестные постоянные коэффициенты, а e – случайное отклонение, вызванное влиянием неучтённых факторов и погрешностями измерений.

Аналогично, между y и   зависимость описывается моделью:

Найдем уравнения линейной регрессии

 и

неизвестные коэффициенты находятся по формулам (используя метод наименьших квадратов (МНК)):

      .

        .

Найдет сначала коэффициенты первого уравнения. Вычисления поясним с помощью таблицы:

 

y

1

7.1

4.9

1.6

50.41

24.01

11.36

7.84

34.79

2.56

2

6.9

4.8

1.5

47.61

23.04

10.35

7.2

33.12

2.25

3

6.8

4.7

1.5

46.24

22.09

10.2

7.05

31.96

2.25

4

6.8

4.6

1.5

46.24

21.16

10.2

6.9

31.28

2.25

5

7

4.7

1.6

49

22.09

11.2

7.52

32.9

2.56

6

7.1

4.9

1.6

50.41

24.01

11.36

7.84

34.79

2.56

7

7.2

5.1

1.7

51.84

26.01

12.24

8.67

36.72

2.89

8

7.5

5.4

1.7

56.25

29.16

12.75

9.18

40.5

2.89

9

7.7

5.7

1.7

59.29

32.49

13.09

9.69

43.89

2.89

10

7.8

5.8

1.8

60.84

33.64

14.04

10.44

45.24

3.24

11

7.8

5.6

1.8

60.84

31.36

14.04

10.08

43.68

3.24

12

7.7

5.5

1.6

59.29

30.25

12.32

8.8

42.35

2.56

13

7.9

5.9

1.8

62.41

34.81

14.22

10.62

46.61

3.24

14

8

5.8

1.9

64

33.64

15.2

11.02

46.4

3.61

15

8

5.8

1.8

64

33.64

14.4

10.44

46.4

3.24

16

7.9

6

1.8

62.41

36

14.22

10.8

47.4

3.24

сумма

119.2

85.2

26.9

891.08

457.4

201.19

144.09

638.03

45.47

Средние

7.45

5.325

1.68125



 

 

 



n=16,

  

  

  

 

Тогда  

    

Таким образом,

Аналогично находятся оценки коэффициентов модели

 а именно,

1.2. По найденным в п. 1.1. уравнениям регрессии построить доверительные интервалы потребления дизельного топлива, соответствующие вероятности 0,9 при следующих значениях независимой переменной: Построить доверительную полосу для уравнения регрессии. Изобразить на графике поля рассеяния, прямые регрессии и доверительные полосы.

 Решение.

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

где  соответственно верхняя и нижняя границы доверительного интервала; значение независимой переменной  для которого определяется доверительный интервал, квантиль распределения Стьюдента, доверительная вероятность, (n-2) – число степеней свободы;

,

Сначала рассмотрим уравнение . По условию задачи  число степеней свободы 16 тогда, по таблице распределения Стьюдента находим t0.90 = 1.76. 

Поясним вычисления с помощью таблицы:

i

1

1.5909

1.0228

0.1706

2

1.5392

0.0015

-0.5500

3

1.5134

0.0002

-0.6500

4

1.5134

0.0002

-0.6500

5

1.5650

0.0012

-0.4500

6

1.5909

0.0001

-0.3500

7

1.6167

0.0069

-0.2500

8

1.6942

0.0000

0.0500

9

1.7458

0.0021

0.2500

10

1.7716

0.0008

0.3500

11

1.7716

0.0008

0.3500

12

1.7458

0.0213

0.2500

13

1.7975

0.0000

0.4500

14

1.8233

0.0059

0.5500

15

1.8233

0.0005

0.5500

16

1.7975

0.0000

0.4500

сумма

 

1.0644

 

Получим

 и заполним таблицу:





6.800

1.513

0.125

1.293

1.734


7.450

1.681

0.069

1.560

1.803

8.000

1.823

0.112

1.626

2.021


График уравнения регрессии у от х1, доверительные интервалы и доверительная полоса приведены на следующем рисунке.

Таким же образом находим доверительные интервалы для уравнения . Используя вспомогательную таблицу:

I

1

1.5842

0.0003

-0.4250

2

1.5613

0.0038

-0.5250

3

1.5385

0.0015

-0.6250

4

1.5156

0.0002

-0.7250

5

1.5385

0.0038

-0.6250

6

1.5842

0.0003

-0.4250

7

1.6299

0.0049

-0.2250

8

1.6984

0.0000

0.0750

9

1.7669

0.0045

0.3750

10

1.7898

0.0001

0.4750

11

1.7441

0.0031

0.2750

12

1.7212

0.0147

0.1750

13

1.8126

0.0002

0.5750

14

1.7898

0.0122

0.4750

15

1.7898

0.0001

0.4750

16

1.8354

0.0013

0.6750

сумма

0.0508 


График уравнения регрессии, и доверительные интервалы и доверительная полоса приведены на следующем рисунке

Задача 2.

2.1. Найти все коэффициенты парной корреляции, проверить их значимость и проанализировать тесноту линейной связи между всеми парами переменных.

Решение.

Коэффициент парной корреляции находится по формуле:

Подставляя соответствующие значения, получим 0,84. Так как

>1.76=t1-a/2, то существенно отличается от 0 и существует сильная линейная положительная связь между y и .

Аналогично проверим неравенство для  :

8,987>1,761, значит, также существенно отличается от 0 и существует сильная линейная положительная связь между y и. .

 Для :

7,813>1,761, значит,  существенно отличается от 0 и существует сильная линейная положительная связь между и .

2.2. Найти по методу наименьших квадратов оценки коэффициентов линейной регрессионной модели

Решение.

Уравнение регрессии ищем в виде:

Обозначения

тогда

Вектор  находится по формуле:

В нашем случае

,

Тогда

Таким образом, .


2.3. Найти коэффициенты множественной корреляции и детерминации.

Решение.

Коэффициент R множественной корреляции определяется по формуле:

Воспользуемся вспомогательной таблицей:

i

1

2,335

-1,1025

0,565

2

2,889

-0,5485

-0,289

3

2,765

-0,6725

-0,165

4

2,395

-1,0425

0,305

5

2,457

-0,9805

0,543

6

2,951

-0,4865

0,249

7

3,036

-0,4015

0,194

8

3,723

0,2855

-0,323

9

4,063

0,6255

-0,363

10

4,202

0,7645

-0,202

11

3,462

0,0245

0,338

12

3,631

0,1935

-0,431

13

4,341

0,9035

-0,341

14

3,971

0,5335

0,329

15

3,971

0,5335

0,129

16

4,480

1,0425

-0,480

Сумма квадратов


7,9623

1,966


Тогда .

Коэффициент множественной детерминации  равен квадрату коэффициента множественной корреляции, то есть  =0,76.

 

2.4. В 2005г. планируется увеличение объёма валовой продукции на 1 млрд. р. по сравнению с 2000г., а объём капитальных вложений – на 0,3 млрд. р. Дать точечный и интервальный прогноз среднего потребления дизельного топлива в 2005г. при уровне доверия 0,9. (Считая, что объёмы валовой продукции и капитальных вложений в 2000г. будут равны запланированным).

Решение.

Для =9,5; = 8,2 получаем точечный прогноз: = 3,74.

Для нахождения интервального прогноза вычислим значения всех параметров, входящих в формулу: , где ,  - соответственно верхняя и нижняя границы доверительного интервала, - вектор независимых переменных, для которого определяется интервал, =1,78 - квантиль распределения Стьюдента, - доверительная вероятность, n – количество наблюдений, (n-3)- число степеней свободы,

=30,37;    =1,966/13=0,1535;    S =0,392;    =2,16;

И тогда = 7,5848; = -0,1048.

2.5. На основе полученных в задачах 1, 2 статистических характеристик  провести содержательный экономический анализ зависимости потребления дизельного топлива от объёмов валовой продукции  сельского хозяйства и капитальных вложений.

Решение.

На основании проведённых расчётов и полученных статистических характеристик можно сделать определённые выводы относительно взаимосвязей между исследуемыми экономическими показателями.

Так как 0,84 и проверка значимости показала его существенное отличие от 0, то есть основания утверждать, что между y и  существует достаточно тесная положительная линейная зависимость, которая может быть отражена с помощью найденного уравнения регрессии . Коэффициент = 0,6858  в данном случае не имеет  экономического смысла. Коэффициент =0,7606 характеризует размер прироста потребления дизельного топлива, обусловленного приростом объёма валового выпуска продукции сельского хозяйства на единицу.

Значение  0,9216 свидетельствует о тесной линейной связи между y и :   Коэффициент b1 = 1,0017 в уравнении показывает, какого прироста потребления дизельного топлива следует ожидать при увеличении объема капитальных вложений на единицу.

Коэффициент = -0,77  в уравнении    показывает, что при росте валового выпуска продукции промышленности на 1 млрд. р. и неизменном объёме капитальных  вложений следует ожидать уменьшение потребления  дизельного топлива на 0,77 млн. тонн. Коэффициент  = 2,16 показывает, что при увеличении объёма капитальных вложений на 1 млрд. р.  и неизменном объеме валовой продукции промышленности следует ожидать увеличение потребления дизельного топлива на 2,16 млн. тонн.


Задача 3.

3.1. Построить ломанную кривую изменения потребления дизельного топлива во времени. Выдвинуть гипотезу о виде зависимости объёма потребления дизельного топлива от времени. Записать трендовую модель, отражающую изменение потребления дизельного топлива во времени. Оценить неизвестные параметры модели методом наименьших квадратов.

Решение.

           На основании визуального наблюдения ломанной кривой, отражающей характер изменения по годам объема потребления дизельного топлива, выдвигаем гипотезу о линейном тренде. Следовательно, трендовая модель, отражающая изменение потребления дизельного топлива, запишется в виде

 

где - неизвестные параметры, -случайное отклонение.

Коэффициенты регрессионного уравнения тренда  находятся по методу наименьших квадратов и равны:

Воспользуемся вспомогательной таблицей:


t

y

ty

1

1

2,9

2,9

2,637

0,263

2

2

2,6

5,2

2,744

-0,144

3

3

2,6

7,8

2,851

-0,251

4

4

2,7

10,8

2,958

-0,258

5

5

3

15

3,065

-0,065

6

6

3,2

19,2

3,172

0,0028

7

7

3,5

24,5

3,279

0,221

8

8

3,4

27,2

3,386

0,014

9

9

3,7

33,3

3,493

0,207

10

10

4

40

3,6

0,4

11

11

3,8

41,8

3,707

0,093

12

12

3,2

38,4

3,814

-0,614

13

13

4

52

3,921

0,079

14

14

4,3

60,2

4,028

0,272

15

15

4,1

61,5

4,135

-0,035

16

16

4

64

4,242

-0,242

cумма

136


503,8

55,032

-0,0572

сумма квадратов

1496






И получим  =1,49;  = 0,0.225. Следовательно, уравнение регрессии будет иметь вид:

3.2. Для найденного уравнения регрессии построить доверительную полосу при уровне доверия 0,9. Нарисовать её на графике вместе с линией регрессии.

 Решение.

Доверительный интервал для линейного тренда  находится по формуле:

где

В нашем случае:

*  = 8; ; S = 0,267.

Результат запишем в виде таблицы:

год

t

1985

1

2,637

0,1215

0,4232

2,8508

1992

8

3,386

0,0667

3,2686

3,5034

2000

16

4,242

0,1331

4,0077

4,4763


На рисунке. изображены график тренда, доверительные интервалы (для t=1,8,16), и доверительная полоса .

3.3. По линейному уравнению тренда найти точечный и интервальный прогнозы среднего потребления дизельного топлива в 2005. и 2007г. (доверительную вероятность принять равной 0,9). Изобразить на графике точечный и интервальный прогноз.

Решение.

 Тогда =4,777 ; = 4,991. Аналогично пункту 3.2 решение запишем в виде таблицы:

годы

t




2005

21

4,777

0,1987

4,4273

5,1267

2007

23

4,991

0,2260

4,5932

5,3888


На рисунке изображены графики линейного тренда , точечные прогнозы  (обозначены точками на прямой) и доверительные интервалы прогноза потребления дизельного топлива в 2005, 2007гг.


Задача 4.

4.2. Для уравнений регрессии:

проверить наличие или отсутствие автокорреляции, используя критерий Дарбина-Уотсона и уровне значимости .

Решение.

Критерий Дарбина-Уотсона имеет вид:

где - отклонения от линии регрессии, i=1,..n.

Для регрессионной модели :

Используя таблицу:



1

-0,565


2

0,289

0,854

3

0,165

-0,124

4

-0,305

-0,47

5

-0,543

-0,238

6

-0,249

0,294

7

-0,194

0,055

8

0,323

0,517

9

0,363

0,04

10

0,202

-0,161

11

-0,338

-0,54

12

0,431

0,769

13

0,341

-0,09

14

-0,329

-0,67

15

-0,129

0,2

16

0,480

0,609



Посчитаем d=3,1573/1,966=1,5818.

У нас n=16, m=2, , следовательно, условие 0,98=>d>=1,54 не выполняется, значит, автокорреляция присутствует.  

Для  используя таблицу: 


I

-0,263


2

0,144

0,407

3

0,251

0,107

4

0,258

0,007

5

0,065

-0,193

6

-0,0028

-0,0678

7

-0,221

-0,2182

8

-0,014

0,081

9

-0,207

-0,193

10

-0,4

-0,167

11

-0,093

0,89

12

0,614

1,5314

13

-0,079

-0,693

14

-0,272

-0,193

15

0,035

0,307

16

0,242

0,207

сумма квадратов

1,0012

4,1302


Посчитаем d=4,1302/1,0012=4,1252.

У нас n=16, m=2, , следовательно, условие 0,98=>d>=1,54   не выполняется, значит, автокорреляция присутствует.

Для сначала найдём коэффициенты:


ei

1

-0,3382


2

-0,1389

0,1993

3

0,1604

0,2993

4

0,3597

0,1993

5

0,359

-0,0007

6

0,1583

-0,2007

7

0,0576

-0,1007

8

-0,1431

-0,2007

9

-0,3438

-0,2007

10

-0,3445

-0,0007

11

-0,0452

0,2993

12

0,1541

0,1993

13

-0,1466

-0,3007

14

0,0527

0,1993

15

0,152

0,0993

16

0,0513

-0,1007

сумм квадратов

0,7792

0,5794


Тогда d=0,5794/0,7792=0,7436. У нас n=16, m=1, , следовательно, условие 1,1=<d<=1,37 не выполняется, значит, автокорреляция присутствует.

4.2. Для уравнения проверить наличие или отсутствие мультиколлинеарности.

Решение.

В пункте 2.1 нашли, что 0,902 – достаточно близок к 1, значит, имеет место мультиколлинеарность.