Приложение 1. Анализ стационарного
демографического состояния

Рассмотрим демографическую характеристику — вероятность Pж для новорождённого прожить не менее T лет:

Pж(T) = F1(T).

Эта зависимость для стационарного состояния имеет универсальный, определяющий характер потому, что зная только её, можно найти большинство других характеристик; это позволяет назвать зависимость F1(T) функцией жизни.

Уточним понятие стационарного демографического состояния. Стационарным будем считать такое состояние, когда не только общая численность людей N0, но и введенная характеристика — функция жизни, остаются неизменными во времени.

Зависимость Pж(T) = F1(T) обладает следующими свойствами:

На рис. П.2 изображены четыре типа таких зависимостей:

1 — вероятность Pж уменьшается равномерно; в этом случае смерть может наступить в любом возрасте от T = 0 до T = Tmax примерно с равной вероятностью;

2 — кривая, характеризующаяся повышенной смертностью в раннем (возможно, детском) возрасте;

3 — кривая, когда существует некоторая наиболее вероятная продолжительность жизни T0. В этом случае вероятность смерти в раннем возрасте мала и существует некоторая, сравнительно небольшая, вероятность дожить почти до максимального возраста Tmax. Наконец,

4 — кривая, при которой продолжительность жизни растёт и приближается к Tmax.

PIC

  Рис. П.2. Четыре типа кривых жизни Pж = F1(T) 

Располагая зависимостью Pж(T), можно найти следующие характеристики:

Рассмотрим связь этих характеристик с функцией Pж(T).

1. Коэффициент рождаемости. Этот коэффициент был определён формулами (1.1) и (1.2) основного текста:

kр = limΔT0 ΔNр N0ΔT = 1 N0 dNр dT ,

где N0 — общее число людей в рассматриваемой группе; ΔNр и dNр — количество родившихся в этой группе за время ΔT и dT. Отсюда:

ΔNр = N0kрΔT, (П1.1) dNр = N0kрdT. (П1.2)

Из dNр (или ΔNр) новорождённых до возраста T доживает только часть, равная dNрF1(T) (или ΔNрF1(T)). В любой момент времени именно столько людей из числа N0 находятся в возрасте от T до T + dT. Поэтому общее число людей в группе равно интегралу:

N0 =0Tmax F1(T)dNр =0Tmax N0kрF1(T)dT = N0kр0Tmax F1(T)dT. (П1.3)

Сокращая это выражение на N0, получим:

kр0Tmax F1(T)dT = 1,илиkр = 1 0TmaxF1(T)dT. (П1.4)

Таким образом, коэффициент рождаемости однозначно определяется функцией F1(T).

2. Относительная плотность распределения людей по возрасту. Гистограммы. Вопрос типа «сколько людей (из группы в N человек) находятся в возрасте 60 лет?», если его понимать буквально (т. е. точно 60 лет), не имеет смысла. Задавая такой вопрос, на самом деле имеют в виду «сколько людей находятся в возрасте от 60 до 61 года?» (ровно 61 год исключается; он не входит в этот диапазон). Говоря о количестве людей определённого возраста, всегда имеют в виду некоторый возрастной диапазон.

Когда всё же подобные характеристики необходимо обрабатывать математически, то вводят понятие о плотности распределения количественной характеристики по непрерывной шкале.

Под плотностью распределения группы из N0 людей по возрасту понимается функция времени (возраста) R(T), которая равна пределу

R(T) = limΔT0N(T;T + ΔT) ΔT , (П1.5)

где ΔT — временной интервал; N(T;T + ΔT) — число людей, возраст которых находится в диапазоне от T до T + ΔT.

Но поскольку группы и число людей в них могут быть разными, то удобнее, разделив R(T) на N0, рассматривать относительную плотность распределения людей по возрасту Rотн. Она равна:

Rотн(T) = R(T) N0 = limΔT0 1 N0 N(T;T + ΔT) ΔT . (П1.6)

При T = 0 получим начальное значение Rотн:

Rотн(0) = R(0) N0 = limΔT0 1 N0 N(0;ΔT) ΔT . (П1.7)

Учтём теперь, что при ΔT 0 N(0;ΔT) есть число рождений ΔNр за время ΔT, а N(T;T + ΔT) равно числу людей из ΔNр, доживших до возраста T:

N(T;T + ΔT) = ΔNрF1(T). (П1.8)

Используя для ΔNр формулу (П1.1) и подставляя это значение в (П1.6) и (П1.7), получим:

Rотн(0) = kр; (П1.9) Rотн(T) = kрF1(T). (П1.10)

Таким образом, относительная плотность распределения людей по возрасту Rотн(T) с точностью до постоянного множителя (kр) равна функции жизни F1(T). В сущности — это одна и та же функция.

Пользуясь выведенными соотношениями для стационарного состояния, можно определять относительное число людей nотн в любом возрастном диапазоне. Например, для диапазона T1 < T < T2 получим:

nотн(T1;T2) = N(T1;T2) N0 =T1T2 R(T)dT = kр T1T2 F1(T)dT. (П1.11)

Это позволяет строить гистограммы, которые широко применяются в демографических расчётах и литературе. Они отражают относительное количество человек, т. е. процентный состав населения, определённых возрастных групп, например, от 0 до 14, от 15 до 60 и старше 60 лет ([5 и др.]), или с интервалами через 5 лет и т. д. На рис. П.3 показан пример такой гистограммы с интервалами ΔT = 10 лет. Из рисунка видно, что в возрасте от 0 до 10 лет находится 22% всего населения, в возрасте от 10 до 20 лет — 19% и т. д. Видно, что численность людей монотонно уменьшается с возрастом, сначала медленно и равномерно до 50 лет, затем быстрее в возрасте 50–60 лет. Число людей старше 70 лет относительно невелико.

PIC

  Рис. П.3. Иллюстрация гистограммы (ступенчатая линия) и относительной плотности распределения людей по возрасту R(T) (непрерывная кривая)  

Для перехода к абсолютным величинам любую относительную характеристику достаточно умножить на N0 — число людей в рассматриваемой группе.

3. Условная вероятность оставшейся продолжительности жизни. Эта характеристика Pж усл(T|T1) определяет вероятность для человека, прожившего T1 лет, дожить до возраста T > T1. Данная зависимость может быть получена из кривой Pж(T) следующим образом. Разделим все ординаты кривой Pж(T) для значений T T1 на ординату Pж(T1). При этом получится кривая, имеющая смысл для T T1 и равная 1 при T = T1. Её значения для различных T T1 и определяют искомую условную вероятность. На рис. П.4 показана эта процедура.

PIC

  1. Исходная зависимость Pж(T) = F1(T).
  2. Условная вероятность прожить не менее T лет, при условии, что человек прожил T1 < T лет — Pж усл(T|T1).
  Рис. П.4. Построение условной вероятности продолжительности оставшейся жизни  

4. Плотность распределения продолжительности жизни и средняя продолжительность Tср. Продолжительность жизни относится к тем характеристикам, которые должны описываться с помощью плотности вероятности.

Рассмотрим два момента времени T и T + ΔT (см. рис. П4). Вероятность дожить до этих моментов равна соответственно F1(T) и F1(T + ΔT). Следовательно, вероятность смерти в интервале от T до T + ΔT равна:

Pсм(T;T + ΔT) = F1(T) F1(T + ΔT). (П1.12)

Разделив это выражение на ΔT, и переходя к пределу при ΔT 0, получим с одной стороны значение плотности вероятности смерти в возрасте Tρсм(T), а с другой стороны производную от функции F1(T) со знаком минус:

ρсм = limΔT0Pсм(T;T + ΔT) ΔT = limΔT0F1(T) F1(T + ΔT) ΔT = dF1 dT . (П1.13)

Средняя продолжительность жизни равна интегралу:

Tср =0Tmax Tρсм(T)dT = 0Tmax dF1(T) dT TdT = 0Tmax TdF1(T). (П1.14)

Используем теперь формулу интегрирования по частям:

abudv = uv| ab abvdu,

где u = T; v = F1(T); du = dT.

PIC

  Рис. П.5. Иллюстрация к определению средней продолжительности жизни и её стабильности по функции Pж(T) = F1(T).  

При этом:

Tср = TF1(T)|0Tmax +0Tmax F1(T)dT.

Первое слагаемое в этом выражении обращается в нуль, т. к. при T = 0 первый, а при T = Tmax второй сомножитель равен нулю. Поэтому окончательно получим для средней продолжительности жизни:

T1ср =0Tmax F1(T)dT = Sж, (П1.15)

где Sж — площадь, ограниченная функцией F1(T).

На рис П.5 изображен прямоугольник ОАВС, равный по площади Sж, и, следовательно, Tср, и также равный Sж по площади прямоугольник ОDЕF с ординатой, равной 1, и абсциссой Tср: SODEF = 1 Tср = Tср. Общая точка M пересечения сторон AB и EF этих прямоугольников лежит на кривой F1(T)8 и может быть названа точкой средней жизни.

5. Стабильность продолжительности жизни. Иногда может представлять интерес вопрос о стабильности продолжительности жизни, т. е. о разбросе реальной продолжительности жизни конкретных людей по отношению к Tср. Подобные характеристики чаще всего оцениваются через дисперсию или среднеквадратичное значение разностей между фактическими величинами (в данном случае фактическими возрастами смерти) и средним значением. Однако возможны и другие способы оценки стабильности. В частности, стабильность продолжительности жизни людей может быть оценена величиной площади криволинейных треугольников DEM и MFC (рис. П.5). Из рисунка легко увидеть, что эти две фигуры по площади равны: SDEM = SMFC, поскольку SODEF = SOABC = SODMC. Поэтому для оценки стабильности δж можно использовать одну из этих площадей:

δж = SDEM =0Tср (1 F 1(T))dT = Tср 0TсрF 1(T)dT (П1.16)

На рис. П.5 пунктиром изображены ещё две кривые жизни — 2 и 3, отличающиеся большей, чем кривая 1, стабильностью при одной и той же средней продолжительности жизни Tср. Кривые 2 и 3 всё ближе примыкают к ломаной линии DEFC. По-видимому, человечество заинтересовано в увеличении стабильности продолжительности жизни и приближении Tср к биологическому пределу Tmax.

Безусловно, наряду с относительными представляют интерес абсолютные характеристики: распределение населения по возрасту, рождаемость, общая численность и др. Для решения этих вопросов с помощью зависимости F1(T) достаточно знать какую-либо одну из абсолютных характеристик, например, общее число людей N0 или число людей какой-либо одной возрастной группы и т. д.

Так, абсолютная численность населения в возрасте от T1 до T2 равна относительной численности nотн (формула П1.11), умноженной на N0:

N(T1;T2) = N0nотн = N0kрT1T2 F1(T)dT. (П1.17)

Отсюда, зная N0 и располагая зависимостью F1(T), можно найти N(T1;T2) для любых значений T1 и T2 или построить гистограмму типа рис. П.3 для любых интервалов ΔT. Можно, наоборот, по значению N(T1;T2) для какой-то одной возрастной группы (T1;T2) найти N0 и значения N(Ti;Tj) для любых остальных возрастных групп (Ti;Tj).

Аналогичным способом могут быть найдены и другие абсолютные характеристики.

В связи со сказанным выше следует отметить два момента:

  1. Построение зависимости F1(T) для стационарного или близкого к нему демографического состояния по экспериментальным данным не обязательно требует изучения статистики на длительном временном интервале. Достаточно располагать данными о возрасте умерших в течение сравнительно небольшого промежутка времени (например, 1 года). По этим данным может быть построена зависимость F1(T) — вероятность дожить до возраста T.
  2. Рассмотренные выше соотношения могут оказаться полезными и для обычных демографических расчётов, когда демографический процесс не стационарен и происходит рост или уменьшение общей численности населения вследствие того, что рождаемость не равна смертности или когда изменяется характер кривой F1(T) (например, вследствие уменьшения детской смертности или увеличения величины Tср и т. д.).

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

В основе мнения об экспоненциальном росте в сущности лежат два момента:

  1. Теоретические расчёты показывают, что при неизменных значениях коэффициентов рождаемости kр и смертности kсм изменение численности людей должно происходить по экспоненциальному закону (см. формулу (1.3) основного текста 1-й главы).
  2. Статистические данные о численности людей как будто бы (на первый взгляд и в отдельных случаях) это подтверждают.

Однако можно выдвинуть ряд возражений против такой точки зрения. Они состоят в следующем.

  1. Коэффициенты kр и kсм всё время меняются под воздействием войн и эпидемий, изменений в условиях жизни людей, успехов по борьбе со смертностью, освоением новых источников питания, изменениями в экологической обстановке, в социальной сфере и многих других факторов.
  2. Простые расчёты, отражённые в табл. 1.4, показывают, что длительное существование человечества при неизменных и неравных между собой коэффициентах kр и kсм невозможно.
  3. Существуют биологические пределы для коэффициента рождаемости kр и продолжительности жизни Tmax. А поскольку все люди смертны, то даже если бы все рождающиеся доживали до предельного возраста Tmax, всё равно очень быстро установилось бы равенство между числом рождений и смертей. Вопрос заключается только в том, при какой общей численности людей установилось бы это равенство.
  4. Всякая «гладкая» нелинейная функция на малом участке может быть приближённо представлена в виде экспоненты. Это может создать иллюзию экспоненциального роста, если рассматривать данные на небольшом отрезке времени.
  5. В литературе встречается такая характеристика, как период удвоения — время, в течение которого численность людей возрастает в два раза. Из статистических данных (см. табл. 1.2) видно и отмечается в литературе [7 и др.], что период удвоения всё время и значительно сокращается. Это ни в коем случае не характерно для экспоненты: при экспоненциальном росте период удвоения должен оставаться неизменным. Сокращение периода удвоения говорит о том, что рост населения происходит быстрее, чем по любому экспоненциальному закону.
  6. Это же подтверждается и анализом данных табл. 1.2. Например, логарифм экспоненциальной зависимости должен представлять собой линейную функцию. Однако анализ табл. 1.2 это не подтверждает.

На рис. 1.3 (основной текст, гл. 1) показан общий рост населения. Видно, что в период с 400 г. до н. э. по 1600 г. рост не был монотонным; наблюдаются периоды спада, например, в 400 и 1400 годах. Быстрый рост, начавшийся с этого времени, привёл к тому, что численность людей достигла 1 миллиарда в 1800 году, 1,5 млрд. в 1900 г., несколько превысила 5 млрд. в 1990 и, по прогнозам ООН, обещает перевалить за 6 млрд. в 2000 году.

На рис. П.6 представлен рост населения Земли в целом и по основным регионам за тот же период в логарифмическом масштабе. По оси абсцисс отложена функция r = lgN106 = lgN 6. Из этих графиков можно увидеть следующее.

  1. Наблюдаются колебания численности во всех регионах, хотя из- за логарифмического масштаба они выглядят сглаженными. Так, например, уменьшению функции r в Америке с 1400 по 1600 годы от значения 1,59 до 1,08 соответствует трёхкратное уменьшение населения с 39 до 13 миллионов человек.
  2. Наблюдается падение численности в Азии с 0-го по 400 и с 800 по 1200 годы, в Европе с 400 по 800 годы, на территории России (СССР) с 1200 по 1400 годы, в Америке с 1400 по 1600 годы, в Австралии и Океании с 1750 по 1850 годы. Как показывает исторический анализ, это было вызвано нашествиями иноземных завоевателей, войнами, колонизацией, эпидемиями, голодом.
  3. Примерно с 1700–1800 годов во всех регионах начался бурный рост населения; в развивающихся странах он продолжается и поныне.
  4. Общий рост населения планеты с 1700–1800 годов происходит быстрее, чем по экспоненциальному закону. Это видно хотя бы из того, что кривая 1 на рис. П.6 не является линейной, а всё более круто устремляется вверх. Дополнительная иллюстрация к этому представлена на рис. П.7. Здесь построена экспонента 1, примерно отражающая рост населения в 1700–1800 годах — на отрезке времени, где логарифмическая кривая 1 на рис. П.6 имеет приблизительно линейный участок. Эта экспонента достаточно хорошо выражает (аппроксимирует) численность в 1700, 1750, 1800, 1850 годах. Но начиная с 1850 года население растёт всё быстрее и быстрее и кривая 2, отражающая реальную численность, всё более уклоняется вверх от экспоненты.

На рис. П.8 , построенном по данным табл. 1.3 и [5], показан рост населения на всей планете во второй половине 20 века. На этом рисунке по оси абсцисс отложена численность в миллиардах человек с шагом в 10 лет.

Всего за 40 лет с 1950 по 1990 годы население планеты возросло в 2,1 раза; в том числе в Азии (без СССР) — в 2,2 раза, в Африке — в 2,9 раза, в Америке — в 2,2 раза, в Европе (без СССР) — в 1,27 раза, в СССР — в 1,6 раза. В Австралии и Океании, которые на рис. 1.2 не отражены, также наблюдалось 2-кратное увеличение.9

Вместе с тем даже на этом графике видно снижение темпов роста в Европе и СССР. Это явление характерно для развитых стран; наблюдается оно также в США и Канаде (если их рассматривать отдельно), хотя и не в такой степени, как в Европе.

PIC

Функция r = lg N 106 = lgN 6.

Число людей N выражено в миллионах.

  Рис. П.6. Рост численности людей во всём мире и по регионам с 400 г. до н. э. по 1980 г. Логарифмический масштаб.  

PIC

  1. Экспонента, аппроксимирующая рост населения в 1700 — 1800 годах.
  2. Фактический рост населения.
  Рис. П.7. Отклонение фактического роста населения мира от экспоненциальной зависимости.  

PIC

  1. Экспонента, аппроксимирующая рост населения в 1700 — 1800 годах.
  2. Фактический рост населения.
  Рис. П.8. Рост населения во всём мире и по регионам с 1950 по 2000 годы.