Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Гинзбург Г.Л. Руководство по использованию метода наименьших квадратов в экономико-статистических расчетах учеб. пособие

.pdf
Скачиваний:
19
Добавлен:
23.10.2023
Размер:
4.64 Mб
Скачать

1

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

Соседние файлы в папке книги из ГПНТБ