книги из ГПНТБ / Гинзбург Г.Л. Руководство по использованию метода наименьших квадратов в экономико-статистических расчетах учеб. пособие
.pdf1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
1922 |
|
|
|
|
|
|
|
|
13 |
III |
5 |
25 |
125 |
284 |
2,45 |
12,25 |
61,25 |
306,25 |
14 |
VI |
7 |
49 |
341- |
317 |
2,50 |
17,50 |
122,50 |
857,50 |
15 |
IX |
9 |
81 |
721 |
1458 |
3,16 |
28,44 |
255,96 |
2303,64 |
16 |
XII |
11 |
12il |
1331 |
7593 |
3,88 |
42,68 |
469,48 |
5164,28 |
|
1923 |
|
|
|
|
|
|
|
|
17 |
III |
13 |
169 |
2197 21186 |
4,33 |
56,29 |
731,77 |
9513,01 |
|
18 |
VI |
15 |
225 |
3375 110000 |
5,04 |
75,60 |
1134,00 |
17010,00 |
|
19 |
IX |
17 |
289 |
4913 |
11107 |
7,04 |
136,68 |
2323,56 |
39500,52 |
20 |
XII |
19 |
361 |
6859 421011 |
12,62 |
239,78 |
4555,82 |
86560,78 |
|
|
|
|
|
|
|
61,75 |
472,87 |
11475,21 |
135833,51 |
Теперь параметры вычислим по формулам (§ 2):
10= 6,97759 — 5,50236 = 1,47523;
Х1== 1,14752 + 1,13110 = — 0,01641;
ла = — 0,02961 + 0,04085 = 0,01124;
Х3 = — 0,0039948 + 0,0048099 = 0,0008151
Уравнение (3, 25) после подстановки полученных парамет ров примет вид
щ __ 1Q1,4в—0 , 0 1 1<J+0,0008P |
(3 26) |
Ввиду малости коэффициента Хз, есть основание пола гать, что выравнивание по суперпозиции 2-го порядка тоже даст хороший результат.
Для получения |
выровненных значений yt подставляем |
в формулу (3, 25) |
значения t, которые мы берем из табли |
цы 24. |
|
Чтобы упростить расчет yt и построить график получен ной функции, прологарифмируем обе части выражения
(3, 26)
\gyx - 1,48—0,016* +0,011** + 0.0008*3
Найденные значения lgyt см. в табл. 25.
101
Т а б л и ц а 24
t |
|
—19 —17 —15 |
— 13 |
—11 |
—9 |
—7 |
—5 |
—3 |
||
lg yt |
|
0,993 |
1,527 |
1,863 |
2,036 |
2,078 2,029 |
1,906 |
1,758 |
1,615 |
|
|
|
|
|
|
|
|
Продолжение таблицы 25 |
|||
—1 |
1 |
3 |
5 |
7 |
9 |
11 |
13 |
15 |
17 |
19 |
t
1,484311,475 1,545 1,754 2,134 2,715 3,546 4,644 6,049 7,793 9,911
После получения значений lgyt наносим на график вы ровненную кривую. На рис. 20 изображен график изменения денежной массы во время инфляции в Германии (1919— 1923 гг.) '.
1 Проф. А.Я. Боярский в работе [2] пишет, что... «в экономической области трудно найти более яркий пример преимуществ, которые имеет логарифмическая шкала, чем при изображении изменений в области финан сов и цен в период крупной инфляции».
102
По оси абсцисс откладываем / в линейном масштабе, а по
оси ординат lgyi н lgyt
§9. Приспособление методологии прогноза
кручному счету
Одним из нерешенных в настоящее время вопросов приме нения математико-статистических методов в практике прогно зирования является приспособление существующей методоло гии к «ручному счету». Это приводит к тому, что расчеты ве дутся короткими цепочками, прерываемыми неформализуемыми в настоящее время операциями: планирования «от базы» и «экспертными оценками». Перевод таких расчетов на вычис лительную технику весьма сложен: требуются надежные тео ретические модели оценки параметров различных по форме зависимостей данных наблюдений.
Весьма распространенным аналитическим методом прогно за является экстраполяция динамического ряда, выровненно го при помощи полиномов различных порядков, а также при помощи различных суперпозиций: полинома и экспоненты, по линома и .ригоер'болы .и т. д. При использовании суперпозиций высоких степеней применяются различные методы, сводящие процесс отыскания параметров к серии простых операций на счетно-клавишных машинах и облегчающие труд экономистаисследователя. Такие методы часто применяются при прове дении обширных экономических расчетов. •■
Наиболее распространенными методами выбора степени полинома являются подбор кривых различных видов по мини мальной величине остаточной дисперсии (критерий F), приме нение полиномов Чебышева и др. в разумном сочетании с эко номическими соображениями, «экспертными оценками» и т. д.
Однако повышение степени полинома ведет к усложнению вида искомой функции и увеличению числа параметров, под лежащих определению. Эти и другие соображения, подтверж денные расчетами, приводят к мысли, что задача экстраполя ции на незначительный период оптимально решается с по мощью полиномов и суперпозиций степени не выше третьей, а задачу прогноза на продолжительный период времени луч ше всего решать с помощью полиномов и суперпозиций степе ни не выше второй, так как они дают более «устойчивый» тренд.
Пример. Рассмотрим данные о численности населения .
СССР за период с 1950 по 1969 годы в оценке на 1 января, приведенные в статистическом ежегоднике «Народное хозяй ство СССР в 1969 году».
103
В таблице 26 выписан ряд динамики за п= 20 лет (число точек эксперимента — четное п—2 к, откуда &=10). Произве
дем выравнивание этого |
ряда по кубической |
параболе у = |
|
= а^+а\Х-\-а2Х2-\-а-;,хг, для |
чего |
осуществим |
линейное пре |
образование х' — 2(х—к)— 1, т. |
е. перенесем |
начало коорди |
|
нат в середину динамического ряда. Параметры системы нор мальных уравнений (3, 11). определим по формулам (3, 12).
Значения й0, Ь\, Ь2, Ь3 |
находим в нижней строке табл. 26, |
|||||||
а значения 7^ -—•в приложении 5. Получаем: |
|
|
||||||
|
«о—210,79645; |
|
|
а, = |
1,82772; |
|
|
|
|
а2= 0 ,00467; |
|
|
а3= |
- 0,00647. |
|
|
|
|
|
|
|
|
|
Т а б л и ц а 25 |
||
п |
Годы |
У |
|
t |
ty |
r-y |
|
t3y |
1 |
2 |
3 |
|
4 |
5 |
6 |
|
7 |
1 |
1950 |
178,5 |
—19 |
—3391,5 |
64438,5 |
—1224331,5 |
||
2 |
1951 |
181,6 |
—17 |
—3087,2 |
52482,4 |
— 892200,8 |
||
3 |
1952 |
184,8 |
—15 |
—2772,0 |
'41580,0 |
— 623700,0 |
||
4 |
1953 |
188,0 |
—13 |
—2444,0 |
31772,0 |
— 413036,0 |
||
5 |
1954 |
191,0 |
—11 |
—2Ш1.0 |
23Ы1,0 |
— 254221,0 |
||
6 |
4955 |
194,4 |
— 9 |
—1749,6 |
15746,1 |
— |
141717,6 |
|
7 |
1956 |
197,9 |
— 7 |
—1385,3 |
9697,1 |
— |
67879,7 |
|
8 |
1957 |
201,4 |
— 5 |
—1007,0 |
5035,0 |
— |
25175,0 |
|
9 |
1958 |
204,9 |
— 3 |
— 614,7 |
1744,1 |
— |
3532,3 |
|
10 |
1959 |
208,8 |
— |
1 |
— 208,8 |
208,8 |
— |
208,8 |
11 |
1960 |
212,4 |
|
1 |
212,3 |
212,3 |
|
212,3 |
12 |
1961 |
216,3 |
|
3 |
648,6 |
1945,8 |
|
5837,4 |
13 |
1962 |
220,0 |
|
5 |
1099,0 |
5495,0 |
|
27475,0 |
14 |
1963 |
223,5 |
|
7 |
1562,4 |
10936,8 |
|
76557,6 |
15 |
1964 |
226,7 |
|
9 |
2037,6 |
18338,4 |
|
165045,6 |
16 |
1965 |
229,6 |
|
11 |
2522,3 |
27745,3 |
|
305198,3 |
•17 |
1966 |
232,2 |
|
13 |
3013,4 |
39174,2 |
|
509264,6 |
18 |
1967 |
234,8 |
|
15 |
3516,0 |
52740,0 |
|
791100,0 |
19 |
1968 |
237,2 |
17 |
4023,9 |
68406,3 |
|
И62907Л |
|
20 |
1969 |
239,5 |
|
19 |
4539,1 |
86242,9 |
|
1638615,1 |
|
&„ =4203,5 |
|
|
5, =4451,1 |
62 =557752,3 6з = 1043759,7 |
|||
Уравнение кубической параболы ийеет вид |
|
|
||||||
|
yt= |
210,8 + |
1,83^ — 0,005^ — 0,0006^ |
|
|
|||
•Остаточная дисперсия |
|
|
|
|
|
|
||
|
аг = Ш —1)г = |
= о, 061. |
|
|
||||
|
|
п |
|
|
20 |
|
|
|
104
Произведем экстраполяцию на 1970 год, для чего в уравне ние кубической параболы подставим t=2\.
Получим, что у t= 241,1. Фактически по переписи 1970 го да численность населения СССР составила 241,6 млн. чел. (на
15 января), |
т. е. абсолютная погрешность Д=241,7—241,1 = |
= 0,6 млн. |
чел. |
Относительная погрешность 6=0,25%,.
Такая точность в большинстве экономических исследова ний считается достаточной, так как аппроксимируемые вели чины точно известны не всегда и невозможно учесть все дей ствующие факторы.
Выровняем тот же ряд по суперпозиции полинома и экспо ненты третьего порядка вида
|
у |
л o Q о ~\~ |
1 |
“Ь ЛоЛГ2 - f #а-£2 |
/ |
/- |
гч г*\ |
||
|
Ю 0 г |
1 |
^ |
3 |
(см. табл. |
26). |
|||
Параметры |
суперпозиции |
находятся |
по формулам |
(3, |
12). |
||||
Значения Ь0,Ьх,Ь^Ьг |
находим в нижней строке таблицы 26, |
||||||||
а значения |
yi — |
в приложении 5. |
|
|
|
|
|||
Получаем: |
|
|
|
|
|
|
|
|
|
|
<2о = 2,32373; |
|
|
щ = 0,00380; |
|
|
|||
|
а2= — 0,00002; |
|
|
=0,000001. |
|
|
|||
Уравнение |
суперпозиции будет |
|
|
|
|
|
|
||
=1Q2,32+0,0038;—0,00002f’ —0,0000 Ог:3
|
|
|
|
|
|
|
Таблица 26 |
|
п |
Годы |
У |
t |
lg У |
Hg у |
Pig у |
Pig У |
|
Г |
2 |
3 |
4 |
5 |
6 |
7 |
|
8 |
1 |
I960 |
178,5 |
—19 |
2,2516 |
—42,7804 |
812,8276 |
—.15443,7244 |
|
2 |
1951 |
181,6 |
—12 |
2,2591 |
—38,4047 |
652,8799 |
—11098,9583 |
|
3 |
1952 |
184,8 |
—15 |
2,2667 |
—34,0005 |
510,0075 |
— 7650,1125 |
|
4 |
195? |
188,0 |
—13 |
2,2742 |
—29,5646 |
■ 384,3398 |
— 4996,4174 |
|
5 |
1954 |
191,0 |
—111 |
2,2810 |
—25,0910 |
276,0010 |
— 3036,0110 |
|
6 |
1956 |
194,4 |
— 9 |
2,2887 |
—20,5983 |
185.3847 |
— |
1668,4623 |
7 |
1956 |
197,9 |
— 2 |
2,2964 |
—16,0748 |
112,5236 |
— |
787,6652 |
8 |
1957 |
201,4 |
— 5 |
2,3040 |
—11,5200 |
57,6000 |
— |
288,0000 |
9 |
1958 |
204,9 |
— 3 |
2,3115 |
— 6,9345 |
20,8035 |
— |
62,4105 |
10 |
1956 |
208,8 |
— 1 |
2,3198 |
— 2,3198 |
2,3198 |
— |
2,3198 |
11 |
196С |
212,4 |
1 |
2,3269 |
2.3269 |
2,3269 |
|
2,3269 |
12 |
1961 |
216,3 |
3 |
2,3349 |
7,0047 |
21,0141 |
|
63,0423 |
13 |
1962 |
220,0 |
5 |
2.3420 |
11,7100 |
58,5500 |
|
292,7500 |
14 |
1963 |
223,5 |
7 |
2,5487 |
16,4409 |
115,0863 |
|
805,6041 |
105
1 |
2 |
3 |
4 |
5 |
6 |
|
7 |
8 |
15 |
1964 |
226,7 |
9 |
2,3549 |
21,1941 |
|
190,7469 |
1746,7221 |
16 |
1965 |
229,6 |
11 |
2,3604 |
25,9644 |
|
285,6084 |
3141 6924 |
17 |
1966 |
232,2 |
13 |
2,3651 |
30,7463 |
|
399,7019 |
5196,1247 |
18 |
1967 |
234,8 |
15 |
2,3699 |
35,5485 |
|
533,2275. |
7998,4125 |
19 |
1968 |
237,2 |
17 |
2,3742 |
40,3614 |
|
686,1438 |
11664,4446 |
20 |
1969 |
239,5 |
19 |
2,3783 |
45,1877 |
|
858,5663 |
16312,7597 |
|
6. =46,4144 |
|
|
bt =9,2720 |
b2= |
b3= 2175,2730 |
||
|
|
|
=6166,7584 |
|||||
Остаточная дисперсия с2 |
2.45 |
0,122. |
|
|||||
20 |
|
|||||||
|
|
|
|
|
|
|
|
|
Подставим в полученное уравнение t= 21.
Тогда
lg yt = 2,3840, |
yt ~ 242,1; |
Д= 242,1—241,7=0,4 млн. чел.;
6=0,15%,.
В данном случае экстраполяция привела к более точному результату.
Если учесть, что оценка численности населения в статисти ческом ежегоднике дается на 1 января, а перепись населения произведена 15 января 1970 г., то результат экстраполяции окажется практически достоверным.
Произведем прогноз численности населения СССР на
1980 год:
а) по параболе 2-го порядка вида
у = а0-j- CL\t -f- |
. |
Параметры параболы определены по формулам:
а0 = тп h — Г is Ьл=0 ,11297 -4203,5-0,00047 -557752,3 =
=210,79645;
04= 7226! = 0,000376-4451,1 = 1,67335;
а2= — 7з1 Ь0+ т13 Ь2= — 0,00047 •4203,5+0,00000356 ■ ■557752,3 = 0,00467.
Уравнение параболы будет
yt = 210,8+ 1,67^ — 0,0047^
106
Подставив в это уравнение / = 41, получим прогнозируемую величину =271,5 млн. чел.;
б) по суперпозиции полинома и экспоненты вида
у = 1 + +
Параметры суперпозици по формулам (3, 12) соответст венно будут:
сг0= 5 ,24359; «1= 0,00349; а2 = — 0,00002.
Уравнение суперпозиции получит вид
105'24+0’00з;~0’00002''2
При t = 41, г/t = 102'428 lg г/t = 2,428.
Прогнозируемая величина г/t =267,8 млн. чел.
§ 10. Экспериментальная проверка модифицированного метода
В мае—июне 1968 г. с целью испытания методов выравни вания статистических рядов и составления корреляционных
уравнений, а также проверки эффективности |
предлагаемого |
в данной работе модифицированного метода, |
в МИНХ им. |
Г. В. Плеханова был поставлен педагогический эксперимент.
Две группы студентов II курса финансового факультета (№№ 921 и 922) в течение 6 часов изучали метод последова тельного исключения и модифицированный метод. Столько же времени они обучались выполнению расчетов на клавишных арифмометрах. После этого дважды была проведена экспе риментальная контрольная на составление корреляционного уравнения и выравнивание динамического ряда по параболе второго порядка.
В аудитории имелись малые вычислительные машины раз личных типов, на каждой машине работали по три студента.
Кроме того, 922 группа была снабжена специальными таб лицами коэффициентов для выполнения расчетов модифици рованным методом. (Количество таблиц было недостаточным).
Все исходные данные были выписаны на доске.
107
Задание 1 предусматривало решение задачи на составле ние корреляционного уравнения. Требовалось определить за висимость урожайности от количества удобрений.
Внесено минер, удоб- |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
рений, и,!га |
|
|
|
|
|
|
|
|
|
|
Урожайность, ц/га |
10 |
13 |
15 |
16 |
19 |
18 |
21 |
23 |
22 |
23 |
Группа 921 решала |
задачу методом |
последовательного |
||||||||
исключения неизвестных. Первый правильный результат по ступил из 922 группы через 30 мин., из 921 группы через 1 час. 30 мин. Таким образом применение модифицированного мето да при решении этой, довольно простой, задачи дало эконо мию времени в три раза.
В задании 2 была предложена задача на выравнивание
динамического ряда. |
Имеются |
данные ассоциации |
воздуш |
||||
ного транспорта США. |
|
|
|
|
|
|
|
|
|
|
|
Годы |
|
|
|
|
1946 |
1948 |
1950 |
1952 |
1954 |
1956 |
1958 |
Полный ежегодный |
|
|
|
|
|
|
|
объем пассажирских пе |
|
|
|
|
|
|
|
ревозок по воздуху, млн. |
5910 |
5910 |
7955 |
12461 |
16696 |
22276 |
28250 |
пассажиромиль. |
|||||||
Первый правильный результат поступил из 922 группы че рез 1 час, из 921 группы через 4 часа. Таким образом приме нение модифицированного метода дало экономию во времени в 4 раза '.
§ 11. Практические рекомендации и выводы
При выравнивании статистических рядов целесообразно переносить начало координат в середину рассматриваемого периода. В этом случае формулы для определения парамет ров уравнения связи еще более упрощаются, а объем вычис лительной работы сводится до минимума.
Таблицы, применявшиеся для общего, нечетного и четного случаев, с успехом могут быть использованы для нахождения параметров экспоненциальной формы связи.
1 Несколько студентов из обеих групп с задачей вообще не справились,
ГОВ
Решение конкретной задачи обычным способом с помощью ЭВМ при больших значениях т и п требует ввода в машину 10— 15 значных чисел и ЭВМ «решает» систему «округлен ных» нормальных уравнений. При вводе в машину наших фор мул получается точный ответ через число k в таком виде, ко торый пригоден для окончательных расчетов.
§ 12. Замечания к пользованию таблицами
Все числа в таблицах |
записаны |
в нормальной форме |
(с плавающей запятой), |
т. е. состоят |
только из значащих |
цифр. |
|
|
При пользовании таблицами необходимо иметь в виду, что перед каждой мантиссой нуль опущен. Округление всюду про изведено с недостатком. Запись чисел в плавающей форме с одинаковым количеством знаков в мантиссе обеопечивает для всех чисел одинаковую относительную погрешность. Пе ред тем как возникает необходимость воспользоваться табли цами, рекомендуется просмотреть решение одного из приме ров, приведенных в работе.
Лример. Найти 7з2 при га= 17, пользуясь приложением 1.
Так как пь—3, то |
в строке, соответствующей га=17, в пятом |
столбце (732 = Т2з) |
находим 2321981410~02 , что означает |
в естественной форме (с фиксированной запятой)
Тз2 = 0,0023219814
ЗАКЛЮЧЕНИЕ
Результаты, полученные в,работе, уже опубликованы в не которых статьях и успешно внедряются в практику в таких организациях как Отдел внедрения математических методов и ЭВМ в планирование народного хозяйства Госплана РСФСР, Отдел экономико-статистического анализа ЦСУ -РСФСР, Отдел выборочных методов НИИ ЦСУ СССР, Все союзный научно-исследовательский институт торгового маши ностроения, ЦЫИЛ при .Главном управлении государствен ных материальных резервов, МИНХ им. Г. В. Плеханова, Мос ковский экоиомикощтатистический институт, Высшая школа профдвижения, Всесоюзный .науяно-исследовательский инсти тут по изучению спроса на товары широкого потребления и
•конъюнктуры торговли.
Полученные .результаты дают возможность сэкономить -время и труд -экономиотов-практиков и исследователей при решении широкого жруга экономико-статистических задач, а так же высвободить время работников ВЦ для выполнения 'массовых 'трудоемких 'расчетов.
.109
Использованию предлагаемого метода должен предшест вовать экономический анализ, позволяющий предположить наличие математического закона изменения данного процесса или явления. Метод прост, удобен и надежен. Одним из его достоинств является то, что для его применения нет необходи мости изучать всю предпосланную ему теорию. Вычислитель ные процедуры легко формализуются на счетно-клавишных машинах.
Решать задачи регрессивного анализа и выравнивания ди намических рядов смогут специалисты среднего звена, не имеющие высшего или специального образования. Метод дает экономию в расчетах и во времени в 3—4 раза и сокращает количество записей в 6—8 раз. Точность расчетов по данному методу весьма высока, относительные погрешности не пре высили 0,07%,.
Еще П. Л. Чебышев говорил: «Практика ищет всегда самого лучшего, самого выгодного» [15].
Для использования приведенных в работе вычислитель ных методов желательно применять счетно-клавишные маши ны, обладающие запасом не менее десяти знаков. Применение этих машин сокращает затраты труда при вычислениях в 2,4 раза ’ .
Большое количество экономических примеров иллюстри рует применение предлагаемого метода в экономических ис следованиях. Прежде чем решать данным методом специаль ную задачу экономист может выбрать для ориентации наибо лее подходящий для него пример.
В настоящее время в сфере материального производства, во всех звеньях социалистического строительства и особенно в области экономики, резко возрастает значение труда эконо мистов, инженерно-технических работников, работников науч ных учреждений, которые затрачивают огромное количество времени на выполнение различных расчетов.
Внедрение прогрессивных систем и методов расчетов с ис пользованием ЭВМ даст огромную экономию затрат общест венного труда.
Автор будет считать свою задачу выполненной и цель до стигнутой, если предлагаемый модифицированный метод вы равнивания статистических рядов пробудит интерес у эконо- миста-исследователя и будет использован при решении ак туальных экономических проблем.
Чрезвычайная признательность заранее выражается за всякие указания и, в особенности, за постановку новых эконо мико-статистических задач, в которых с успехом могут быть применены методы, изложенные в работе.
1 В. С. К р а сн и ц к и й. Применение ЭЦВМ в математической стати стике, Куйбышев, 1967 г.
110
