Построение модели множественной регрессии

 

Цель работы: овладеть способами построения модели множественной линейной регрессии и выработать умения и навыки нахождения параметров уравнения, оценки надежности уравнения регрессии и его параметров, проведения экономической интерпретации полученных результатов.

 

Содержание работы: на основании опытных данных требуется

1. Рассмотреть результативные и факторные признаки и определить форму связи между ними. Построить все корреляционные поля (попарно). И на основании этого записать уравнение множественной регрессии.

2. Исключить лишние факторы, обосновав математически данный выбор.

3. Измерить тесноту связи между факторами.

4. Используя метод наименьших квадратов оценить уравнение регрессии.

5. Проверить на адекватность уравнение регрессии используя следующие способы:

– используя критерий Фишера;

– с помощью коэффициента детерминации;

– используя среднюю ошибку аппроксимации.

6. Дать экономическую интерпретацию найденных оценок уравнения регрессии.

Задача. Исходные данные для признаков _{, , ,} – для различных нефтегазодобывающих управлений — приведены в табл. 2. 21.

 

Таблица 2.21

Признаки	Значение признаков на различных НГДУ
	0,92	0,93	0,89	0,90	0,90	0,89	0,92	0,91	0,93	0,89

 

В таблице обозначено: – коэффициент эксплуатации скважин; – дебит скважин (тн/сут.); – уровень автоматизации труда (%); – производительность труда (тн/чел.).

 

 

Выполнение работы

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

, .

, .

, .

, .

.

Таблица 2.22

Коэффициенты парной корреляции равны соответственно:

,

.

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

,

.

Рис. 2.9.

Рис. 2.11.

Рис. 2. 13.

Рис. 2.10

Рис. 2.12.

Рис. 2.14.

X1

X1

X1

X 2

X 2

X3

X 2

X 3

X 3

Y

Y

Y

 

 

Так как _,то между признаками и связь сильнее, чем между и . Этот факт подтверждается путем вычисления коэффициентов частной корреляции и :

 

,

.

 

Поэтому из модели множественной линейной регрессии исключаем фактор . Тогда в модель будут включены факторы и и уравнение регрессии запишется в виде:

 

.

 

Включение фактора в модель обосновано значимостью коэффициента парной корреляции :

 

.

 

Для выяснения вопроса о силе линейной связи между факторами, включенными в модель, вычисляем множественный коэффициент корреляции R:

.

 

Так как в нашем примере объем выборки небольшой (), то произведем корректировку R:

 

.

 

Проверяем значимость по критерию Стьюдента. Вычисляем среднеквадратическую ошибку :

 

,

 

Вычисляем статистику (ф.1.102):

 

 

По таблице критических точек распределения Стьюдента при уровне значимости с числом степеней свободы находим (приложение 5). Так как , то делаем вывод, что значим.

Для нахождения оценок , , уравнения регрессии решаем систему нормальных уравнений:

 

 

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

Проверяем адекватность уравнения регрессии. Используем коэффициент детерминации , полагая . Для полученной модели . Это означает, что полученная модель приблизительно на 66% объясняет изменение производительности труда в зависимости от изменения включенных в модель факторов и , что является не плохим показателем.

Проведем проверку модели на адекватность по критерию Фишера – Снедекора. Найдем статистику по формуле (ф.1.78), полагая в ней :

 

.

 

По таблице критических точек распределения Фишера — Снедекора при уровне значимости и числах степеней свободы , (p — число факторов , включенных в модель, n — объем выборки) находим . Так как , то найденное уравнение регрессии, устанавливающее зависимость производительности труда на десяти нефтегазодобывающих управлениях (НГДУ) от коэффициента эксплуатации скважин и дебита скважин , значимо описывает опытные данные и может быть принято для руководства.

Оценим адекватность уравнения регрессии по средней ошибке аппроксимации , которую вычислим по формуле (ф.1.63):

 

.

 

Для нахождения суммы составляем расчетную табл. 2.23.

Таблица 2.23
Y
	35,2	0,2	0,000114
	36,8	0,8	0,017778
	31,5	0,5	0,005952
	33,8	0,8	0,019394
	32,7	1,3	0,049706
	32,6	0,6	0,008
	36,3	1,7	0,076053
	34,7	0,7	0,014412
	37,1	0,1	0,00027
	32,2	0,8	0,019394
			0,211073

 

По данной таблице находим:

.

Среднеквадратическая ошибка небольшая, что дает основание считать, что построенная модель адекватно описывает опытные данные.

Итак, все три метода проверки модели на адекватность подтвердили гипотезу о том, что уравнение регрессии в целом статистически значимо и хорошо соответствует данным наблюдений.

Дадим экономическую интерпретацию найденных коэффициентов уравнения регрессии. Значение свободного члена характеризует влияние неучтенных в модели факторов, в частности фактора (уровень автоматизации труда). Знак минус говорит о том, что отсутствие этого фактора в модели отрицательно сказывается на повышении производительности труда. Величина коэффициента показывает, что при увеличении коэффициента эксплуатации на 0,01 производительность труда увеличивается в среднем на 86,3271 тн/чел. Коэффициент показывает, что при увеличении дебита скважин на одну тонну производительность труда увеличивается в среднем на 0,360611 тн/чел.

Лабораторная работа № 6.