
книги из ГПНТБ / Сивый В.Б. Метод множественной корреляции в анализе и планировании угольных предприятий
.pdfТаблица 7
\ |
X _ ^ о т к |
|
250 |
|
350 |
|
450 |
|
550 |
|
650 |
|
750 |
||
2 0 0 - -300 |
3 0 0 - -400 |
4 0 0 --500 |
5 0 0 - 6 0 0 |
600— 700 |
700 — 800 |
||||||||||
|
у=П |
|
|
то |
т х с р |
т |
т х с р |
т |
т * с р |
т |
т х с р |
то т х ср ТП |
т * 0 р |
||
|
13 |
|
m |
|
|
_ |
|
_ |
|
_ |
|
_ |
|
— |
|
|
1 2 - 1 4 |
т у с р |
— |
|
— |
|
— |
|
— |
|
— |
|
|
|
|
|
15 |
|
то |
_ |
|
_ |
|
_ |
|
9 |
4950 |
5 |
3250 |
3 |
2250 |
|
1 4 - 1 6 |
т у с р |
— |
|
— |
|
— |
|
— |
135 |
— |
75 |
— |
45 |
|
|
|
|
|
|
|
||||||||||
|
17 |
|
т |
_ |
|
_ |
|
7 |
3150 |
15 |
8250 |
7 |
4550 |
2 |
1500 |
|
1 6 - 1 8 |
" ч / Ср |
- |
|
- |
|
— |
113 |
“ |
255 |
— |
119 |
— |
34 |
|
|
19 |
|
т |
_ |
|
5 |
1750 |
20 |
9000 |
12 |
6600 |
2 |
1300 |
1 |
750 |
|
|
|
|
2 2 8 |
|
38 |
|
19 |
|||||||
|
18— 20 |
т |
УС р |
— |
|
|
95 |
— |
380 |
— |
— |
— |
|||
|
21 |
|
т |
1 |
250 |
9 |
3150 |
19 |
8550 |
14 |
7700 |
5 |
3250 |
— |
|
|
2 0 - 2 2 |
т |
Ус р |
- |
21 |
— |
189 |
— |
399 |
— |
294 |
— |
105 |
— |
|
|
|
|
|||||||||||||
|
23 |
|
т |
5 |
1250 |
3 |
1050 |
10 |
4500 |
24 |
13200 |
2 |
1300 |
— |
|
|
2 2 - 2 4 |
т у с р |
|
115 |
|
69 |
|
230 |
|
552 |
|
46 |
|
|
"
25 |
|
т |
6 |
1500 |
7 |
2450 |
25 |
11250 |
9 |
4590 |
_ |
_ |
|
24 — 26 |
т у с р |
— |
150 |
— |
175 |
|
625 |
|
225 |
- |
~ |
- |
|
|
|
|
|||||||||||
27 |
|
т |
15 |
3750 |
7 |
2450 |
17 |
7650 |
1 |
550 |
_ |
— |
|
26— 28 тпуср |
|
405 |
- |
189 |
- |
459 |
— |
27 |
— |
- |
|
||
29 |
|
т |
в |
1500 |
6 |
2100 |
1 |
450 |
— |
|
— |
— |
|
28 —30 |
т у с р |
- |
174 |
- |
174 |
— |
29 |
— |
|
— |
— |
|
|
31 |
|
т |
3 |
750 |
8 |
2800 |
1 |
450 |
— |
|
— |
__ |
|
30— 32 |
т у с р |
_ |
93 |
— |
248 |
' — |
31 |
— |
|
— |
— |
|
|
33 |
|
т |
|
3 |
1050 |
— |
|
— |
|
— |
|
|
|
32— 34 |
т у с р |
|
|
~ |
99 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* т х |
|
|
36 |
48 |
|
|
100 |
|
84 |
21 |
|
6 |
|
|
|
|
|
|
|
|
|
|
|||||
2 ( т х |
у с р ) |
|
958 |
1238 |
|
|
2272 |
|
1716 |
383 |
|
98 |
|
Г т д : х Ср |
|
9000 |
16800 |
|
|
45000 |
|
46200 |
13650 |
4500 |
|||
х Ср 2 ( т х |
Уср) |
239500 |
433300 |
|
1022400 |
943800 |
248950 |
73500 |
|||||
— 2 < т х |
y a v ) |
|
2 6 ,9 |
25,8 |
|
2 2 ,7 |
|
2 0 ,4 |
18,2 |
|
16,3 |
||
* |
2 т х |
|
|
|
|
||||||||
|
|
|
|
/ |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
я —Lqtk |
х= П |
|
13 |
т |
1 2 - 1 4 |
т у ср |
15 |
т |
1 4 - 1 6 |
™Уср |
17 |
т |
1 6 - 1 8 |
т у ср |
19 |
т |
18— 20 |
т Уср |
21 |
т |
20— 22 |
тУср |
23 |
т |
2 2 - 2 4 |
™-УСр |
|
25 |
т |
|
24—26 |
т у ср |
|
'27 |
т |
|
26—28 |
т у ср |
|
29 |
т |
|
28—30 |
т у ср |
|
31 |
т |
|
30—32 |
т у ср |
|
33 |
т |
|
32—34 |
т у ср |
|
Х тх |
|
|
2<тх |
Уср) |
|
Щ тх |
хСр ) |
•■'ср2 »'7"* |
УСр> |
|
_ |
2 (т х Уср) |
|
Ух |
Х тх |
850 |
950 |
80 0 — 900 |
900— 1000 |
1 |
т х ср |
т |
т х ср |
” 1 |
|
|
|
_ |
|
1 |
950 |
—13
2 |
1700 |
1 |
950 |
— |
30 |
|
15 |
1 |
850 |
- |
|
—17
—-
—
—-
—
—-
_ - — -
——
--
--
-
'
- |
|
- |
- |
- |
- |
|
|
|
3 |
|
2 |
47 |
|
28 |
2550 |
|
1900 |
39950 |
|
26600 |
15,7 |
|
14,0 |
Продолжение таб л . 7
|
|
|
|
|
Ху — |
Х т у |
ЖгПу хСр) |
уСр^ |
Уср2 <т ух ср) |
SOriyXgp |
|
|
|
|
|
|
Хту |
1 |
950 |
13 |
, |
12350 |
950 |
20 |
13100 |
300 |
|
196500 |
655 |
32 |
18300 |
544 |
|
311100 |
572 |
40 |
19400 |
760 |
|
368600 |
485 |
48 |
22900 |
1008 |
|
480900 |
477 |
44 |
21300 |
1012 |
|
489900 |
484 |
47 |
20150 |
1175 |
503750 |
429 |
40 |
14400 |
1080 |
388800 |
360 |
13 |
4050 |
377 |
117450 |
312 |
12 |
4000 |
372 |
124000 |
333 |
3 |
1050 |
99 |
34650 |
350 |
300 |
- |
- |
- |
|
- |
- |
6740 |
- |
|
- |
139600 |
- |
- |
|
- |
- |
- |
3028000 |
|
- |
|
- |
|
|
Из расчетной таблицы видно, что каждому значению х соответствует определенное зна
чение ух, т. е. ух = / (х). Для наглядности дан ные функции ух = / (ж) приведены ниже:
X |
250 |
350 |
450 |
550 |
650 |
750 |
850 |
950 |
Ух |
26,6 |
25,8 |
22,7 |
20,4 |
18,2 |
16,3 |
15,7 |
14,0 |
Аналогично может быть представлена ре
грессия Х у по у.
и |
13 |
15 |
17 |
19 |
21 |
23 |
25 |
27 |
29 |
31 |
33 |
Xи |
950 |
655 |
572 |
485 |
477 |
484 |
429 |
360 |
312 |
333 |
350 |
Как указывалось выше, линия, связывающая на корреляционном поле найденные средние
значения ухху в интервалах изменения х (у),
называется эмпирической линией регрессии.
Линия регрессии показывает, как в сред
нем изменяется функция у с увеличением или уменьшением значения аргумента х.
В нашем примере теоретическую линию ре
грессии можно считать прямой и искать ее
уравнение в виде ух = а + Ъх. Для определе ния коэффициентов а и Ь уравнения большей
частью применяется метод наименьших квадра тов. При этом каждая точка (х, ух) учитывается столько раз, сколько данных было использова-
64
йо для ее вычисления. Другими словами, если ух является средней из т х данных, соответ
ствующих значению х, то ухучитывается т х раз. Принцип наименьших квадратов требует,
чтобы сумма квадратов отклонений условных
средних величин от расчетных значений по фор
муле регрессии была наименьшей. Это означает,
что параметры уравнения регрессии прямой ух= = а + Ьх надо выбрать так, чтобы было вы держано условие
2 т х (ух — у)2 = 2 т * (Ух — Ьх — а)2 = минимум. (19)
Суммирование производится по всем значе ниям х.
В наших расчетах величина b выступает в
качестве коэффициента |
регрессии (qx/v) у по |
|||
х. Искомые |
величины |
определяются |
из двух |
|
нормальных |
уравнений |
|
|
|
2 |
У ~ |
Qv/x 2 |
* “Ь |
(20) |
2 ху = |
еУ/* 2 я2 + 2 ха- |
|
Решив эту систему уравнений относительно
Qx/y, ПОЛУЧИМ
2 * у - 2 * - 2 *
2 - 2—( v - ) a
Упростим полученное выражение, разделив
числитель и знаменатель формулы на /г2:
Ci//* = |
(ху)ср — х у |
(21а) |
-т |
5 |
8 |
65 |
Числитель формулы (21а) называется кова
риацией сух. |
Иначе ковариация сух выражается |
|||
как |
частное |
от |
деления суммы |
произведений |
ус р 2 |
(тух ) |
по |
всем интервалам |
изменения у |
(см. табл. 7) на общее число наблюдений п ми
нус произведение средних значений ху.
Аналогично определяется коэффициент ре грессии х по у.
Q x / y — |
(22) |
Геометрический смысл коэффициента регрес
сии Qy/X в случае линейной зависимости заклю
чается в том, что он определяет собой угловой
коэффициент прямой. Реально коэффициент b
уравнения регрессии ух = а + Ъх показывает, на сколько в среднем изменяется у при измене
нии признака х на единицу. Чем больше число
вое значение коэффициента регрессии, тем бо лее значителен прирост функции при том же увеличении аргумента.
Свободный член а определяется путем под
становки известных значений регрессии в си стему уравнений (2 0 ).
Пользуясь данными табл. 7 и формулами (21) и (2 2 ), определяем Qx/y и ду/я:. Значение хуср
определяется делением данных, взятых из стро
ки яСр 2 (яг* г/Ср) или из столбца ycpLmx xcv,
на объем выборки п. Среднее значение аргумен
та х рассчитывается по формуле (3); аналогично
рассчитывается среднее значение функции у.
Вычисленные дисперсии рассмотренных пока зателей равны: Оу = 18,8089 и ох — 16 814,88
66
Определим |
коэффициент |
регрессии х по у*. |
|
3 028 000 |
|
|
|
300 |
• — 22,47 • 465,33 |
— 19,27970. |
|
Qx/V |
18,8089 |
= |
|
|
|
||
Коэффициент регрессии у по х равен |
|||
3 028 000 |
|
|
|
300 |
■ 22,47 • 465,33 |
— 0,02157. |
|
Qy/x — |
16814,88 |
= |
|
|
|
|
|
Определяем |
уравнения |
регрессии: |
|
1 ) регрессия х по у. |
|
|
|
ху — с + |
Qx/уУ — с — 19,27970 • 22,47; |
||
2 ) регрессия у по х: |
|
|
|
ух — а + |
Qy/x = а — 0,02157 • 465,33. |
Для точки с координатами (х, у)
с = 465,33 + 19,27970 • 22.47 = 898,5449; а = 22,47 + 0,02157 • 465,33 = 32,5053.
В окончательном виде уравнения теоретиче ских линий регрессии имеют вид:
Ху = — 19,27970г/ + 898,5449;
ух = — 0,02157х + 32,5053.
Теоретические линии регрессии строятся по
точкам, полученным в результате подстановки
минимального и максимального значений х в приведенные выше уравнения.
Если эмпирические и теоретические линии
регрессии пересекаются примерно в одной точ
ке с координатами (х, у), то это подтверждает правильность выполненных расчетов и графи ческих построений.
5* |
67 |
Внашем примере в результате подстановки
=950 м и oc-i = 250 м соответственно полу
чим: г/i = |
12,0138 |
т/мес и г/2 = |
27,1128 |
т/мес. |
||||||||||
Большой |
наклон |
линии регрессии на |
рис. 6 |
|||||||||||
|
|
|
|
|
|
|
подтверждает |
|
пра |
|||||
|
|
|
|
|
|
|
вильность |
|
прямого |
|||||
|
|
|
|
|
|
|
сопоставления |
|
рас |
|||||
|
|
|
|
|
|
|
сматриваемых |
пока |
||||||
|
|
|
|
|
|
|
зателен. |
Об |
|
этом |
||||
|
|
|
|
|
|
|
свидетельствуют |
так |
||||||
|
|
|
|
|
|
|
же |
корреляционные |
||||||
|
|
|
|
|
|
|
поля и линии регрес |
|||||||
|
|
|
|
|
|
|
сии, |
приведенные на |
||||||
|
|
|
|
|
|
|
рис. |
7. |
|
|
|
|
||
|
|
|
|
|
|
|
Научное |
и прак |
||||||
|
|
|
|
|
|
|
тическое |
|
значение |
|||||
|
|
|
|
|
|
|
корреляционных |
за |
||||||
|
|
|
|
|
|
|
висимостей, |
|
помимо |
|||||
Рис. |
6. |
Корреляционные |
формы |
связи, |
опре |
|||||||||
деляется |
еще другим |
|||||||||||||
поля |
и |
линии |
регрессии |
|||||||||||
свойством |
|
— |
тесно |
|||||||||||
для |
зависимости |
относи |
|
|||||||||||
тельной протяженности |
от |
той |
связи. |
Степень |
||||||||||
каточных |
|
выработок |
на |
тесноты связи харак |
||||||||||
1000 |
т |
месячной |
добычи |
теризует |
|
рассеяние |
||||||||
угля |
LmK |
от |
производи |
|
||||||||||
точек |
корреляцион |
|||||||||||||
тельности труда рабочего П. |
||||||||||||||
ного |
|
поля |
относи |
|||||||||||
|
|
|
|
|
|
|
|
|||||||
тельно |
линии |
регрессии |
(рис. |
8 .) |
Чем |
ближе |
или плотнее располагаются точки к линии ре
грессии, тем меньшее влияние на исследуемую
величину |
оказывают прочие |
факторы |
и |
|
тем более |
точна |
и определенна |
зависимость |
|
у по х. |
|
|
|
за |
При тесной и близкой к функциональной |
||||
висимости |
между |
показателями каждая отдель |
||
на я точка корреляционного поля |
находится не |
посредственно около линии регрессии, а у прак тически полностью зависит от х.
Если точки корреляционного поля распола
гаются на значительном расстоянии от линии
регрессии, т. е. связь слабая, то действие како
го-либо аргумента на функцию можно пред сказать лишь ориентировочно.
Таким образом, степень тесноты связи яв ляется весьма важным свойством корреляцион ной зависимости.
В экономических исследованиях с помощью
корреляции, как правило, исследуются зави
симости между отдельными переменными, имею
щими неодинаковые единицы измерения. По этому для оценки тесноты связи между призна ками разной природы пользуются коэффициен том корреляции (гху), который представляет
собой коэффициент регрессии у по г в стандар
тизованной системе единиц. Стандартизован
ный коэффициент регрессии (коэффициент кор реляции) показывает, на какую часть своего квадратического отклонения ау изменяется у при изменении х на величину своего квадрати ческого отклонения ах, когда остальные факто
ры, не связанные с х, остаются неизменными, а
связанные с х — изменяются неопределенно.
Коэффициент корреляции может иметь как по ложительное, так и отрицательное значение в
пределах (—1 ) ^ Гух < 1 .
Теснота связи в случае линейной корреля
ционной зависимости определяется по формуле
rvx — Qy/xQx/y (23)
Сопряженные коэффициенты регрессии бе рутся при этом со своими знаками ( + или —).
69
Рис. 7. Корреляционные поля и линии регрессии для зависимости между производи тельностью труда рабочего П и скоростью проведения подготовительных выработок Уподг (“)> относительной протяженностью основных выработок на 1000 т месячной до
бычи угля Lnp0B (б), скоростью подвигания очистной линии забоев Роч (в), относи тельной протяженностью ремонтируемых выработок на 1000 т месячной добычи угля
£рем (®)*
а знак коэффициента корреляции совпадает со
знаком коэффициента регрессии и показывает
характер зависимости между показателями — возрастающий или убывающий.
У |
У |
. |
У |
Рис. 8. Рассеивание опытных точек (наблюде ний) при различной тесноте связи между слу чайными величинами у и х:
а — тесная связь; б — слабая связь; в — отсутствие связи.
Коэффициент линейной корреляции можно
также определить, пользуясь формулой
г у х ---- |
ух |
(” У)ср — |
(23а) |
ЯД, |
«Д , |
|
|
|
|
т. е. как частное от деления ковариации на произведение среднеквадратичных отклонений (сгжи ау) . Эта формула может быть представлена в ином виде:
п 2, |
х у — 2 У • 2 ж________ |
|
Гу* |
— |
. (236) |
V '[» 2 * 2 |
(2 *)2] [П S V* -Т 2 Э Д |
При выполнении расчетов вручную более
удобна формула (23а), а на счетных машинах — формула (236).
П