УДК 519.24, 316.4
СТОХАСТИЧЕСКАЯ МОДЕЛЬ СТРУКТУРЫ СОЦИАЛЬНОЙ
ГРУППЫ1
Савельев Л.Я., Гончарова Г.С.
Новосибирский государственный университет, Новосибирск, Тувинский государственный университет, НИЛ «Теории вероятностей и ее приложений»,
Кызыл,
Институт философии и права СО РАН, Новосибирск THE STOCHASTIC MODELS OF SOCIAL GROUPS STRUCTURE
Savelyev L.J., Goncharova G.S. Novosibirsk State University, Novosibirsk, Tuvan State University, Laboratory "Theory of Probability and its Applications", Kyzyl. Institute of Philosophy and Law, Russian
Academy of Sciences, Novosibirsk
В статье рассматриваются простая марковская стохастическая модель социальной группы и примеры ее применения к исследованию динамики структуры конкретных групп и ее прогноза. Описываются группы учащейся молодежи некоторых регионов Сибири. Исследуются условия достижения выбранной стационарной структуры рассматриваемой социальной группы.
Ключевые слова: вероятность, марковское свойство, социокультурные типы, коллективисты, индивидуалисты.
The article discusses a simple stochastic Markov model of a social group, and examples of its application to the study of the dynamics of the structure of particular groups and its prognosis. Describes a group of students in some regions of Siberia. We study the conditions for achieving the selected stationary structure considered a social group.
Key words: probability, Markov property, socio-cultural type, collectivists, individualists.
Введение. Рассматриваются формализованные классификации индивидов данного социума по выбранным личным, общественным и культурным признакам. В простейшей двоичной классификации по каждому выбранному признаку индивид относится к одному из двух уровней: высокому или низкому. В более подробной троичной классификации по каждому выбранному признаку индивид относится к одному из трех уровней: высокому, среднему или низкому. Содержание выбранных признаков и их уровней определяется рассматриваемым социумом и поставленной задачей. Предлагается стохастическая модель, позволяющая описывать состояние социума в данный момент и его динамику в рассматриваемом индивидуально-социально-культурном пространстве, а также делать краткосрочный, среднесрочный и долгосрочный прогнозы состояния
1 Материалы статьи подготовлены по гранту РГНФ № 12-03-00546 а. «Социокультурные типы молодежи современной России: этническое и региональное измерение»
Тувинский государственный университет
социума. Термин классификация употребляется в математическом смысле: установление отношения эквивалентности или разбиение множества на части.
Среди стохастических моделей, описывающих различные естественные и социальные процессы, выделяются марковские. Они отличаются сравнительной простотой и часто с успехом применяются. Для описания социальных явлений эффективно используются марковские процессы с дискретным временем и конечным множеством состояний, называемые конечными марковскими цепями [1; 3]. Определения основных понятий элементарной теории вероятностей есть в книге [4]. Важность и эффективность применения статистических методов в социологии подробно описаны в книге [5].
Использование предлагаемых простых математических моделей предполагает существенное упрощение выбранных характеристик сложных представлений о личности, обществе и культуре. Но простые модели позволяют достаточно адекватно описывать некоторые фундаментальные закономерности структуры рассматриваемых общественных групп, связанные с личнымикачествами составляющих их индивидов, социальными и культурными особенностями группы. Деление на два уровня дает возможность по имеющимся неполным данным достаточно уверенно относить индивида в одну из двух групп. Если двух уровней мало для решения поставленной задачи и нужно более подробная характеристика составляющих группу индивидов, их личных, социальных и культурных качеств, то можно использовать три уровня. Деление на три уровня часто тоже дает возможность по имеющимся неполным данным достаточно уверенно относить индивида в одну из трех групп.
В теории вероятностей и математической статистике разработаны модели с любым числом признаков и любым числом уровней. Если есть надежный статистический материал, то можно использовать более общие модели, чем предложенная двоичная и троичная. Использование модели с большим числом признаков и большим числом уровней связано с большим объемом вычислений. Применение современных компьютерных программ позволяет с ним справиться. Главными недостатками применения моделей с большим числом признаков и большим числом уровней является трудность учета взаимной зависимости выбранных признаков и размытость границ выбранных уровней. Кроме того, адекватность модели существенно снижается использованием неточных описательных определений рассматриваемых признаков и уровней в имеющихся материалах. Эти недостатки существенно снижают адекватность моделей с большим числом признаков и большим числом уровней. Наконец, большое число признаков и уровней затрудняет содержательную интерпретацию полученных результатов. Часто для надежной интерпретации приходится группировать признаки и уровни, получая те же результаты, что и при изначальном использовании модели с небольшим числом признаков и небольшим числом уровней.
Особенно хорошо разработаны марковские стохастические модели. Их применение связано с теоретическим или статистическим определением
вероятностей перехода рассматриваемой системы из одного выделенного состояния в другое. Определение таких вероятностей обычно требует большого объема статистических данных. Часто нужных данных нет, а их получение связано со значительными трудностями. Для моделей с большим числом признаков и большим числом уровней эти трудности нередко бывают непреодолимы.
Оценить такие трудности можно на примере классификации по трем признакам (личные, общественные и культурные качества) и пяти уровням (очень высокий, высокий, средний, низкий, очень низкий). Точное описание составленной из получающихся 125 групп структуры и ее динамики представляется сложной проблемой. Два уровня (высокий, низкий) дают 8 групп, а три уровня (высокий, средний, низкий) - 27 групп. В этих случаях можно рассчитывать на успех.
Серьезным ограничением для применения марковских стохастических моделей является необходимость выполнения марковского свойства: распределение состояний системы в следующий момент зависит только от ее состояния в данный момент и не зависит от состояний в предыдущие моменты. Марковское свойство часто поясняется фразой: при известном настоящем прошлое и будущее стохастическинезависимы. Термины состояние и момент имеют очень общий смысл, и поэтому модель может описывать широкий круг процессов любой природы. Этим объясняется широкое применение марковских стохастических моделей в самых разных областях. Исключительно важным свойством этих моделей является возможность прогнозировать поведение рассматриваемой системы. Марковские прогнозы часто хорошо согласуются с реальностью.
1. Марковские случайные последовательности. Содержание понятий событие, случайность, вероятность хорошо знакомо каждому. Опишем кратко простейшую математическую модель, используемую в статье.
Основным понятием простейшей стохастической модели является вероятностное пространство составленное из конечного множества
и = {и} исходов и и семейства р = {р[и\ 1 элементарных вероятностей р[и] > 0, = 1
. Исходами могут быть элементы любой природы. Часто
их удобно нумеровать, отождествлять исходы с их номерами и считать I/= {X 2,..., т} тдет ~ число исходов. Части Л £ Е/ множества называются событиями. Сумма Р Ш элементарных вероятностей РМ, составляющих событие ^ исходов и , пазываетсявероятностъю этого
события: цел Ясно, что 0 < Р\А\ ^ 1- В частности, вероятность
пустого (невозможного) события ® равна нулю, а вероятность полного {достоверного) события V равна единице: ^М = 0, Р[ЕЛ=1 . Вместо г иногда пишут
Рг^. Числа? М и
Р\Л\
служат мерами реализуемости исхода
и и события ^ .
Для вероятности, ^ верно, правило сложения: + В] = Р[,4] + для непересекающихся (не имеющих общих исходов) событий , В . Равенство Р[АВ]
Р[А\В] =
Р[В] выражает условную вероятность события ^ при условии, что произошло событие & . Верно правило умно.»сепия: Р\АВ] = Я[АЦ?]Р|В] Здесь сумма А + В обозначает объединение событий -4. В а произведение - их пересечение (общую часть). Предполагается, что Р\Щ > О
(событие
В возможно). Если = р[Д|Р[В] т0 события ^ и В называются
стохастически независимыми или просто независимыми. Независимость событий А и В выражается также равенством P\A\S\ = P\Ä\, Понятие независимости является фундаментальным для теории вероятностей.
Любая функция на множестве исходов называется случайной переменной. Она может принимать абстрактные значения. Выделяются случайные переменные со шкалированными значениями, в частности числовыми. Числовые случайные переменные называются случайными величинами. Они являются главным объектом теории вероятностей. Активно развивается теория возможностей для шкалированных случайных переменных [6]. Применения теории шкал в социологии подробно описаны в книге [7].
Случайная величина / на вероятностном пространстве^'Р^порождает вероятностное пространство составленное из множества = Ы ее
значений х~ /М и семейства Я ~ flflli элементарных вероятностей дМ = ^
Л"]=х (складываются элементарные вероятности? всех исходов и ,
М = 1
для которых /М = I ) Ясно, что t/Dt] > О _ Числоравно в
пространстве р! вероятности P\fkf\u\ = x 1 того, что случайная величина / примет значение х . Семейство Ч ~ fyMJ называете ^распределением случайной величины /Важнейшей характеристикой случайной величины / является ее
я/ = /мрм
среднеезначение u&i . По связанной с азартными играми традиции
число Ef называют также математическим ожиданием. Верно удобное для
Ef = У хцЫ
вычислений равенство хех . Для оценки отклонений значений случайной величины ) от среднего значения u ~ Ef часто применяется ее дисперсия, равная среднему квадратичному отклонению: ~ — ß)2] Верно удобное
Df = - d)zqlx\
для вычислений равенство хех . Используется также стандартное
отклонение а ~ Vö/. С ним связано знаменитое правило трех сигм: значения
/И. для которых I* а| > 3а _ считаются невозможными. Это правило
часто применяют необоснованно.
Для лучшего восприятия вероятностной модели полезна ее механическая интерпретация. Семейство Р ~ описывает распределение единицы массы
на множестве ^ : в точке и сосредоточена масса РМ. а на множестве ^ сосредоточена масса Р\Л\, Случайная величина / переносит массу на множество ^ вещественной прямой: в точке х сосредотачивается масса Я М. Среднее значение случайной величины / является центром масс, а дисперсия - моментом инерции такой механической системы. Применение вероятностных моделей в статистике обосновывает закон больших чисел, утверждающий, что при определенных условиях вероятность события в определенном смысле близка частоте его реализаций.
Семейства случайных переменных исследует теория случайных процессов. Элементы этой теории вместе с основами теории вероятностей излагаются в книге [8]. Опишем кратко используемое в статье достаточно простое семейство случайных переменных, не уточняя вероятностное пространство, на котором они определены. Рассмотрим последовательность случайных переменных
1. = ■.■ 1 2 ... с конечным множеством значений. Эти значения могут быть любыми, но удобно их перенумеровать и заменить номерами. Множеством значений тогда становится конечное множество номеров С — {1,2,...,) Последовательность с, называется марковской цепью, если она обладает марковским свойством', при известном настоящем прошлое я <
будущее
{■ТОО- стохастически независимы.
В рассматриваемом конечном случае марковское свойство выражается в том, что совместные распределения конечных семейств {?(0. О < < п) задаются произведениями матриц, составленных из начальных и переходных вероятностей. Строка ^ ~ ^ С = 1, — начальных вероятностей состоит из вероятностей й; = Рт\£(О) = ¿] Того, что в начальный момент £ = О значением случайной величины
т будет номер ; . Матрица у переходных
вероятностей состоит из условных вероятностей Яч ~ ^^ + У — ./К(£) ~~ ' ]
того, что в следующий момент
значением случайной величины
будет номер / , если в данный момент £ значением случайной величины является номер 1 . Подчеркнем, что здесь термин момент служит синонимом термина номер или термина индекс и может не иметь никакого отношения ко времени.
Распределение случайной величины выражает строка V' ~
составленная из вероятностей РГ = (0 = г'] того, что в момент £ значением случайной величины £(£) является номер ; . В частности, распределение случайной величины <Г(0) выражается равенствами
Р =~ Р' ~ 1 ■■■>Замена состояний их номерами упрощает
обозначения. Состояния могут описываться целыми фразами или сложными графическими изображениями. Матрицы умножаются по правилу строка на
столбец:
Переход через ^ моментов, от £ к £ + fc описывается матричным равенством Pt+ =PtQ . В частности, распределение случайной величины выражается равенством р11 — AQU При та =2,3 эти распределения легко вычисляются в общем виде [4]. Получены также общие формулы для их важных характеристик, связанных со структурой серий в реализациях таких марковских последовательностях.
2. Двоичная и троичная формальные классификации. Опишем подробнее предлагаемые классификации индивидов данного социума по личным , социальным
и культурным (Ю признакам. Двоичная классификация. По каждому признаку К социум делится на типы А и В . В результате получаются группы РА, РВ, К А, К В Тип
Л описывает высокий уровень индивида (личный, социальный, культурный), тип В . низкий. Конкретное содержание этих типов определяется в зависимости от рассматриваемого социума и поставленной задачи.
Из полученных групп образуются 8 классов
А А А, ААВ, ABA, В А А, АВВ, ВАВ, В В А, ВВВ g класс ААА относятся индивиды с
высоким уровнем по всем трем признакам. В класс ААВ относятся индивиды с высоким уровнем по первым двум признакам и низким по третьему. Остальные классы тоже легко интерпретируются. В последний класс ВВВ относятся индивиды с низким уровнем по всем трем признакам.
Троичная классификация. По каждому признаку Р> К социум делится на типы А, С, В В результате получаются группы
РА, РС.РВ, £„4, SC.SB, К А, К С, К В тип А описывает высокий уровень
индивида (личный, социальный, культурный), тип ^ - средний, тип ^ -низкий. Конкретное содержание этих типов определяется в зависимости от рассматриваемого социума и поставленной задачи. Уровни качеств можно определять, не уточняя общие определения самих качеств. Это существенно облегчает задачу.
Из полученных групп образуются 27 классов:
шмм^шм^шщ^тттдщштшщмшщмшм _
В класс относятся индивиды с высоким уровнем по всем трем признакам. В класс '4,4(7 относятся индивиды с высоким уровнем по первым двум признакам и средним по третьему. В класс ААВ относятся индивиды с высоким уровнем по первым двум признакам и низким по третьему. Остальные классы тоже легко
интерпретируются. В последний класс относятся индивиды с низким уровнем по всем трем признакам. Отступление от алфавитного порядка объясняется соответствием обозначений предложенных классификаций: всюду ^ - высокий уровень, а ® - низкий.
От троичной классификации можно перейти к двоичной, если разделить множество индивидов уровня С на два множества и
(1 — Х)С, 0 < А < 1 ^ присоединив ^ - долю множества к множеству ^ , а оставшуюся (1 — Я) - долю - к множеству В. Объединенные группы
определят новые высокий и низкий уровни в
двоичной классификации.
Выбранный порядок классов в определенной степени отражает их значимость: от высоких уровней к более низким. Это позволяет при анализе результатов переходить от абсолютной шкалы к порядковой. При одинаковом конечном числе градаций эти шкалы асимптотически эквивалентны [9]. Применение марковской модели дает возможность исследовать динамику психологической, социальной и культурной структур выделенного общественного слоя. Можно выделить периоды стабильности и делать краткосрочный, среднесрочный и долгосрочный прогнозы.
Общая классификация. В общем случае выбранное множество делится по определенным признакам на п попарно непересекающихся (не имеющих общих элементов) частей, которым присваиваются номера 1,2, ..., п . Назовем эти части большими. Каждая большая часть в свою очередь по определенным правилам делится на т попарно непересекающихся мелких частей, которым присваиваются внутренние номера 1, 2, ..., т . Кроме того, мелкой части с внутренним номером 1 из большой части с номером / присваивается двойной номер У . В итоге множество разбивается на т" частей. В частности, при п = 2,т = 2 получается части; при п = 3,т = 2 получается частей; при и = 3, т — 3 получается частей. Содержание такой модели может быть самым разным. Деление больших частей на различные числа маленьких технически усложняет модель.
В последнее время стали иногда вместо множеств рассматривать функции со значениями в отрезке [0,1], описывающие степень принадлежности элемента множеству [6]. Такие функции называются нечеткими множествами. Нередко они более адекватно описывают содержание задачи, но их корректное применение связано с математическими трудностями.
3. Марковская модель с 2 состояниями. Рассмотрим двоичную марковскую последовательность \ случайных переменных С > 0 с
множеством значении ^ — {1,2} (номера состояний), начальной и переходной матрицами
!и (¡12\ IV
Л» Ч12Л / V А = 1а1з а^] = {я, 1 - а! * \д21 Чгг' и-1
Здесь а1 — Рг{№ ~ вероятность появления в начальный
момент
номера
а
цч Рг{?,\1 I И Ц
вероятность появления в момент
номера
)
при условии, что в момент
появился номер
I (1 < г',/ < 2) Будем использовать более простые обозначения °> Р- Ч и
, _1 - ч
а = р + ц - 1 1-Й'
Нетрудно доказать, что верно следующее выражение для £ - ой степени О матрицы С? :
Из матричного равенства Р" = следуют формулы
р! = Рг[<(0 = 1] = Ъ + (а - ьур,
VI = М<(0 = 2] = 1 - Ы- (а -
для вероятностей появления в момент
номеров
верны приближенные
И<1 , то при достаточно больших значениях Е равенства Р1 к , Рз и 1 ~ & .
Постепенно последовательность (О стабилизируется. Распределение случайной величины все ближе к ГЬД - Ь] и все меньше зависит от
начального распределения 1 ~ я).
Распределение
{Ь,1~Ь1ста щонарно: 1 - ¿0- С = 1ЬД-Ь1 Если в некоторый моменте распределение случайной величины £,00 оказывается равным \Ъ,1-Ъ} х0 эх0 распределение имеют и все следующие случайные величины Процесс переходит в стационарный режим. Чем меньше
тем быстрее убывает ИГ и быстрее достигается стационарный режим. Выделяется случай Р + Ч = 1 . В этом случае <¿ = 0 и р^ = Ь = р
t> О
и
1-й
при
и любом начальном распределении, 1 ~~ а }. Марковская последовательность превращается в последовательность независимых случайных переменных. Стационарный режим достигается уже в первый момент.
В статье [10, с. 82] описываются 6 социокультурных типов молодежи Республики Саха (Якутия) в целом, а также в этнических группах русских, саха (якутов) и группы КМНС (коренных малочисленных народностей Севера). В разработанную авторами специальную анкету включены вопросы об отношении к обществу, о важности осознания своей национальной принадлежности, знания языка и культуры своего народа. По ответам на эти вопросы были выделены следующие типы:
(1) коллективист с развитым этническим самосознанием, (2) коллективист с менее развитым этническим самосознанием, (3) коллективист с неразвитым этническим самосознанием, (4) индивидуалист с развитым этническим
самосознанием, (5) индивидуалист с менее развитым этническим самосознанием, (6) индивидуалист с неразвитым этническим самосознанием.
Если интересоваться отношением индивида к обществу и отвлечься от этнического самосознания, то естественно объединить типы 1 - 3 в общий тип (коллективист), а типы 4-6 объединить в общий тип (индивидуалист). Для этнических групп русских, сах (якутов) и коренных малочисленных народов Севера (КМНС) доли коллективистов к моменту окончания анкетирования (2010 год) были соответственно равны 0.71, 0.75, 0.87. Статистических данных для определения переходных вероятностей нет. Поэтому рассмотрим различные гипотетические варианты течения процесса.
Предположим, что за интервал времени между моментами £ и С + 1 коллективисты и индивидуалисты с одинаковой вероятностью продолжают
1
придерживаться своих взглядов: V ~ Ч № < р, ц < 1) Тогда 2
и доли коллективистов и индивидуалистов постепенно становятся равными. Обозначим Т отношение вероятностей для коллективистов и для
индивидуалистов за интервал времени между моментами
1 - р
и
изменить
г =
свои взгляды на противоположные:
коллективистов постепенно становятся равной индивидуалистов постепенно становятся равной
= 0.4 т____а = 0.2,
Ч. В этом случае доля 1 + г, а доля
1 - Ъ =■
.
Пусть,
Г =3
Ь = - = 0.75 4
например, Р ~ 0.8, q — и.4 Тогда
Будем считать такую долю коллективистов нормальной. Так как — (0-2)*
быстро убывает, то признанная нормальной доля коллективистов в этнической группе русских при а ~ 0.71 быстро достигается:
р\ = 0.742 р12 = 0.7484 р1 = 0.7497 р£ = 0.7499 р1 = 0.75
В этнической группе саха (якутов) при а = 0.75 и таких переходных вероятностях доля коллективистов стационарна: р1=0.75 для всех ^ — 0 . В этнической групп КМНС при й = 0.87 и таких переходных вероятностях доля коллективистов постепенно размываетсядо признанной нормальной:
= 0.774 р\ = 0.7548 р1 = 0.751 = 0.7502, р1 = 0.75
Г и Г + 1
Формально интервалы времени между моментами 1 и 1 А могут быть любыми. Выбор этих интервалов определяется только адекватностью @ переходов от состояния в момент £ к состоянию в
описания матрицей
момент
С + 1
4. Марковская модель с 4 состояниями. Рассмотрим марковскую последовательность ? случайных переменных . с множеством
значений С = (1, 2,3, 4) (номера состояний), начальными вероятностями
Тувинский государственный университет
аг = = 'Т и переходными вероятностями ^ч = ^ + Ч = / КМ = Ч
составляющими 1х4-матрицу (строку) А = и 4Х4 . матрицу ^ ~~ (ч^} Здесь, как и прежде, аг - вероятность появления в начальный момент состояния с номером г . а - вероятность появления в момент £ + 1 состояния с номером ) при условии, что в момент ^ появилось состояние с номером ¿(1^1, ]< 4) Можно выписать общие формулы для*?' и ?■ — -1ч', но они громоздкие. Проще рассматривать нужные частные случаи.
Если исследовать отношение индивида к обществу и его этническое самосознание, не выделяя отношение к языку и культуре, то естественно объединить типы 2 - 3 и точно так же объединить типы 5 - 6 . В итоге получаются 4 типа: - коллективист с развитым этническим самосознанием, -коллективист с недостаточно развитым этническим самосознанием, -
индивидуалист с развитым этническим самосознанием. - индивидуалист с недостаточно развитым этническим самосознанием. Занумеруем эти классы в алфавитном порядке: ~~ 1' АВ — 2, ВА — 3, ВВ — 4. Будем считать коллективизм и индивидуализм признаками ( или на первом месте), а этническое самосознание измерять уровнями ( или на втором месте).
Для этнических групп русских, саха (якутов) и КМНС доли рассматриваемых классов к моменту окончания анкетирования (2010 год) представлены в таблице 1.
Таблица 1
Распределение социокультурных типов среди учащейся молодежи в этнических группах Республики Саха (Якутия)
Национальность Социокультурные типы
- коллективист с развитым этническим самосознанием - коллективист с недостаточно развитым этническим самосознанием - индивидуалист с развитым этническим самосознанием - индивидуалист с недостаточно развитым этническим самосознанием
Русские 0.42 0.29 0.16 0.13
Саха (якуты) 0.63 0.12 0.19 0.06
КМНС 0.67 0.2 0.08 0.05
Статистических данных для определения переходных вероятностей нет. Поэтому рассмотрим различные гипотетические варианты течения процесса.
Сделаем упрощающие предположения. Будем считать, что вероятность перехода от признака или уровня к признаку или уровню равна V- а вероятность обратного перехода равна Я . Одновременную смену признака и уровня будем считать невозможной. При этих условиях матрица переходных вероятностей имеет вид 0" = Р + '/):
Рассмотрим произвольную строку безиндексными обозначениями А ~ {К' Все ненулевые элементы матрицы Q
и строки 1
0 < р. Я < -
начальных /?, У, & больше 0 и меньше 1.
вероятностей с 1 -а - /? - у}
Будем предполагать, что 2 . Естественно выделяется частный
случай Ч ~ V , когда переходы отк и обратно одинаково вероятны. В этом случае матрица? симметрична и вычисления проще.
Стационарным является распределение
легко проверяется. По стандартным правилам <2 <1 = (1-г)2 (1- 2г).
Равенство ВО = В вычисляется определитель матрицы Q : = (1 — г)2 (1 — 27')- Степени
множителей 1 —т,1 — 2Г входят в выражения для разности распределения случайной величины ? М и стационарного. Если эти множители достаточно малы, то стационарное распределение быстро достигается. В случае ? ~ Р , когда переходы от к и обратно одинаково вероятны, верны равенства
р=я=1 В = \ЬХ-,-]
г г 2 и стационарное распределение равномерно: ¡.4 ' 4' 4-' 4.). Доли
всех рассматриваемых классов становятся одинаковыми.
Рассмотрим сначала распределение для этнической группы русских. Пусть, например, V = 0-1- 9=0-4. Тогда 1 -г = 0.5, 1-2 г = о
Ь = {0.64, 0.16, 0.16, 0.04} Будем считать такие доли
рассматриваемых классов нормальными. Так как множители 1 ~ 1 — 2Г
малы, то признанные нормальными доли классов достигаются быстро и при сохранении переходных вероятностей в дальнейшем не меняются. Они выписаны в таблице 2.
При тех же переходных вероятностях динамику распределений в этнических группах саха (якутов) и КМНС описывают таблицы 3 и 4.
Таблица 2
Гипотетическая динамика изменения социокультурных типов у русских до признанного нормальным при заданных переходных вероятностях Р = 0.1, Ц — 0.4
Время 1- Социокультурные типы
- коллективист с развитым этническим самосознанием - коллективист с недостаточно развитым этническим самосознанием - индивидуалист с развитым этническим самосознанием - индивидуалист с недостаточно развитым этническим самосознанием
1 0.52 0.24 0.17 0.07
2 0.58 0.2 0.17 0.06
Е + 3 0.61 0.18 0.16 0.05
Е + 4 0.62 0.17 0.16 0.04
Е + 5 0.63 0.16 0.16 0.04
Е + 6 0.64 0.16 0.16 0.04
При тех же переходных вероятностях динамику распределений в этнических группах саха (якутов) и КМНС описывают таблицы 3 и 4.
Таблица 3
Гипотетическая динамика изменения социокультурных типов у саха до признанного нормальным при заданных переходных вероятностях Р = 0.1, = 0.4
Время 1- Социокультурные типы
- коллективист с развитым этническим самосознанием - коллективист с недостаточно развитым этническим самосознанием - индивидуалист с развитым этническим самосознанием - индивидуалист с недостаточно развитым этническим самосознанием
t+l 0.63 0.15 0.18 0.04
Е + 2 0.63 0.15 0.17 0.04
Г + З 0.64 0.16 0.17 0.04
4 0.64 0.16 0.16 0.04
Таблица 4
Гипотетическая динамика изменения социокультурных типов у коренных малочисленных народов Севера до признанного нормальным при заданных переходных вероятностях
Время Социокультурные типы
- коллективист с развитым этническим самосознанием - коллективист с недостаточно развитым этническим самосознанием - индивидуалист с развитым этническим самосознанием - индивидуалист с недостаточно развитым этническим самосознанием
1 0.65 0.19 0.13 0.04
1 + 2 0.64 0.17 0.14 0.04
Е + 3 0.64 0.17 0.15 0.04
Е + 4 0.64 0.16 0.16 0.04
Вычисления всюду проводились с округлением до сотых.
Замечание. Из таблиц 2-4 видно, что у саха и КМНС при заданных одинаковых переходных вероятностях время достижения стационарного распределения социокультурных типов, условно признанных нормальными, совпадает, а у русских - больше на два периода. Это косвенно может говорить о близости личных, социальных и культурных характеристик у саха и КМНС.
5. Марковская модель с 8 состояниями. Рассмотрим марковскую последовательность f случайных переменных fM. t^O с множеством значении (номера состояний), начальными вероятностями
! L jj и переходными вероятностями !ч
составляющими 1х8-матрицу (строку) ^ ~ и 8х8-матрицу ? —
Здесь, как и прежде, ai - вероятность появления в начальный момент состояния с номером ' , а 'ч - вероятность появления в момент t + 1 состояния с номером J при условии, что в момент ^ появилось состояние с номером : ' Д — - ; — э . Будем рассматривать подходящий для выбранного примера частныи случаи такои модели с достаточно простыми выкладками.
В статье [10, с. 96] приведено распределение группы молодежи Якутии по личным, социальным и культурным признакам, в каждом из которых выделены два уровня. По этой классификации группа была разделена на 8 классов .-i.-i.-i .4.4^ А-- А , .4.-1 А-- г А.-- г г .4 г г г . Классы занумерованы в алфавитном порядке:
AAA — 1, ААВ — 2, ABA — 3, BAA — 4, ABB — 5, BAB — 6, ВВА — 7, BBB — 8 Отступая от принятой в [9, с. 82] терминологии, можно так описать эти классы.
По личным качествам выделялись индивиды с высоким уровнем активности (буква на первом месте в обозначении класса). По социальным качествам выделяются индивиды (коллективисты) с высоким уровнем социального взаимодействия (буква на втором месте в обозначении класса). По культурным качествам выделяются индивиды с высоким уровнем этнического самосознания (важно осознавать свою национальную принадлежность, знать язык и культуру своего народа) (буква на третьем месте в обозначении класса). Альтернативы отмечаются буквой на соответствующем месте.
В качестве начального выберем распределение для молодежи 20 - 24 лет:
Достаточных статистических данных для определения переходных вероятностей нет. Поэтому рассмотрим различные гипотетические варианты течения процесса.
Сделаем упрощающие предположения. Будем считать, что вероятность перехода от признака или уровня к признаку или уровню равна V , а вероятность обратного перехода равна Ч . Одновременную смену признака и уровня будем считать невозможной. При этих условиях матрица переходных вероятностей имеет вид (.2 р + q = г, р + 2 q = s) ■
Все ненулевые элементы матрицы 0 больше 0 и меньше 1. Будем
предполагать, что ^ 3 . Естественно выделяется частный случай
Ч = V , когда переходы от к и обратно одинаково вероятны. В этом случае матрица 0 симметрична и вычисления проще. Стационарным является распределение
Ъ = ~
где Р + Ч По стандартным
с =
V + Я Равенство ~ В легко проверяется, правилам вычисляется определитель матрицы
(1= (1-и)3(1-2и)3(1-Зи)
где
и
.
Степени
множителеи 1-й, 1
2и,
зи
входят в выражения для разности
распределения случайной величины СМ и стационарного. Если эти множители достаточно малы, то стационарное распределение быстро достигается. В случае
Ч ~ Р , когда переходы от к и обратно одинаково вероятны, верны равенства
1
Ъ = с = -
I и стационарное распределение равномерно:
П 1 1 1 1 1 1 П
В =
■1 1 1 1 1 1 1 1|
У 8 8 8 8 8 . Доли всех рассматриваемых классов
одинаковыми.
Пусть V = Я =1/5.
Тогда и = '
1 - и =
3/5'
2 и = 1/5,
и стационарное распределение & равномерно. Умеренно малый размер множителей 1 — и, 1 — 2и, 1 — Зи определяет умеренную скорость
приближении распределения к равномерному. Эти распределения в моменты 1 = ^ 1 - 1 ^ с округлением результатов до тысячных
выписаны в таблице 5.
Таблица 5
Гипотетическая динамика изменения социокультурных типов у молодежи возрастной группы 20-24 лет до приближения распределения к равномерному при заданных переходных вероятностях
1
Время t Социокультурные типы
AAA ААВ ABA BAA ABB BAB ВВА ввв
3 0.116 0.132 0.136 0.101 0.15 0.114 0.117 0.134
6 0.123 0.127 0.127 0.12 0.13 0.123 0.123 0.127
12 0.125 0.125 0.125 0.124 0.126 0.125 0.125 0.125
15 0.125 0.125 0.125 0.125 0.125 0.125 0.125 0.125
Приведем также пример неравномерного стационарного распределения. 1 1 Ч =
Пусть Р 12'
6 .
13 11
и = —,1 - и = —,1 - 2 и = —,1 - 3 и
4 4 2 4
Тогда ч ч ¿4 и
В = {0.30,0.15, 0.15,0.15, 0.07,0.07, 0.07,0.04} при округлении до сотых. Будем считать такие доли рассматриваемых классов нормальными. Умеренно
малый размер множигелеи
1 - и, 1 - 21/, 1 - 311
определяет умеренную
скорость стабилизации распределения. Эти распределения в моменты - = - ! с округлением результатов до сотых выписаны в таблице 6.
Таблица 6
Гипотетическая динамика изменения социокультурных типов у молодежи возрастной группы 20 -24 лет до стацонарного распределения при заданных
1 1
переходных вероятностях
V =
12'
Я =
Время t Социокультурные типы
AAA ААВ ABA BAA ABB BAB ВВА ВВВ
5 0.23 0.14 0.16 0.12 0.11 0.08 0.08 0.05
10 0.28 0.15 0.15 0.14 0.08 0.08 0.08 0.04
15 0.29 0.15 0.15 0.15 0.08 0.07 0.07 0.04
20 0.3 0.15 0.15 0.15 0/07 0.07 0.07 0.04
Можно провести аналогичные вычисления и в общей марковской модели с 8 состояниями, не выводя общих формул и не делая упрощающих предположений о переходных вероятностях. Существующие математические программы позволяют производить достаточно точные вычисления и с матрицами больших порядков.
Заключение. Если имеющихся статистических данных достаточно для определения переходных вероятностей, то можно выбрать стохастическую марковскую модель, с нужной адекватностью описывающую динамику исследуемой общественной группы по индивидуальным, социальным и культурным признакам. Марковская модель дает возможность прогнозировать
Тувинский государственный университет
эту динамику на близкую, среднюю и далекую перспективу. Использование аппарата неоднородных марковских цепей и последовательное изменение переходных вероятностей позволяет учитывать изменение условий. Прогнозы по удачно выбранным марковским моделям хорошо согласуются с реальностью.
Приведенные в статье примеры хотя и используют реальные статистические данные, но имеют иллюстративный характер. Отсутствие достаточных статистических данных для определения переходных вероятностей позволяет рассматривать только их различные гипотетические варианты. Если определить желательную по социальным соображениям стабильную структуру рассматриваемой группы населения, то можно выбрать специальные переходные вероятности так, чтобы они обеспечивали постепенный переход к стационарному распределению по выделенным классам за приемлемое время. Переходные вероятности могут давать представление о направлениях усилий общества для стабилизации, рассматриваемой структуры. В частности, речь может идти о воспитании личности, этнического сознания, общей и национальной культуры. Библиографический список
1. Кемени Дж., Снелл Дж. Конечные цепи Маркова. Москва, Наука, 1970. C. 271.
2. Пытьев Ю.П. Возможность. Элементы теории и применения. Москва, УРСС, 2000. С.
190.
3. Жданок А.И. Основы теории конечных цепей Маркова. Теория и практика. Кызыл, ГУП «Тываполиграф», 2008. С. 168.
4. Савельев Л., Балакин С. Конечные марковские цепи и серии (теория и приложения). Saarbrucken, LAMBERT Academic Publishing, 2012. С. 150.
5. Савельев Л.Я. Элементарная теория вероятностей, 1-2. Новосибирск, НГУ, 2005. С. 158,
190.
6. Толстова Ю.Н. Математико-статистические модели в социологии. Москва, ГУ ВШЭ, 2008. С. 243.
7. Толстова Ю.Н. Измерение в социологии. Москва, 2007. С. 287.
8. Феллер В. Введение в теорию вероятностей и ее приложения, 1-2, Москва, Мир, 1984. С. 528, 752.
9. Загоруйко Н.Г., Савельев Л.Я. Относительная мощность измерительных шкал. Вычислительные системы, вып. 101. Новосибирск, 1984. С. 11 4-129.
10. Абрамова М.А., Костюк В.Г., Гончарова Г.С. Жизненные планы и адаптация различных социокультурных типов молодежи. Журнал социологии и социальной антропологии. Т. 16, № 1, 2013. С. 78-99.
Bibliograflcheskij spisok
1. Kemeni Dzh., Snell Dzh. Konechnye tsepi Markova. Moskva, Nauka, 1970. C. 271.
2. Pytev Yu.P. Vozmozhnost. Elementy teorii i primeneniya. Moskva, URSS, 2000. S. 190.
3. Zhdanok A.I. Osnovy teorii konechnykh tsepej Markova. Teoriya i praktika. Kyzyl, GUP "Tyvapoligraf", 2008. S. 168.
4. Savelev L., Balakin S. Konechnye markovskie tsepi i serii (teoriya i prilozheniya). Saarbrucken, LAMBERT academic publishing, 2012. S. 150.
5. Savelev L.Ya. Elementarnaya teoriya veroyatnostej, 1-2. Novosibirsk, NGU, 2005. S. 158, 190.
6. Tolstova Yu.N. Matematiko-statisticheskie modeli v sotsiologii. Moskva, GU VSHE, 2008. S.
243.
7. Tolstova Yu.N. Izmerenie v sotsiologii. Moskva, 2007. S. 287.
8. feller v. vvedenie v teoriyu veroyatnostej i ee prilozheniya, 1-2, moskva, mir, 1984. s. 528, 752.
9. Zagorujko N.G., Savelev L.Ya. Otnositelnaya moschnost izmeritelnykh shkal. Vychislitelnye sistemy, vyp. 101. Novosibirsk, 1984. S. 114-129.
10. Abramova M.A., Kostyuk V.G., Goncharova G.S. Zhiznennye plany i adaptatsiya razlichnykh sotsiokulturnykh tipov molodezhi. Zhurnal sotsiologii i sotsialnoj antropologii. T. 16, № 1, 2013. S. 78-99.
Савельев Лев Яковлевич - кандидат физико-математических наук, профессор Новосибирского государственного университета, профессор Тувинского государственного университета, НИЛ «Теории вероятностей и ее приложений», E-mail: [email protected]
Гончарова Галина Савитовна - научный сотрудник Института философии и права СО РАН, E-mail: [email protected]
Savelyev Lev - PhD Candidate of Physics and Mathimatics, professor of Novosibirsk State University, Professor of Tuvan State University, Laboratory "Theory of Probability and its Applications", E-mail: [email protected]
Goncharova Galina - Researcher, Institute of Philosophy and Law, Russian Academy of Sciences, E-mail: [email protected].
УДК 519.217.2 +314.3
МАРКОВСКАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДИНАМИКИ КОЛИЧЕСТВА ДЕТЕЙ В СЕМЬЯХ ТУВЫ
Хурума А.К.
Тувинский государственный университет, НИЛ «Теории вероятностей и ее приложений», Кызыл
THE MARKOV MATHEMATICAL MODEL THE DYNAMICS OF THE NUMBER OF CHILDREN IN FAMILIES
OF TUVA
Khuruma A.K.
Tuvan State University, Laboratory "Theory of Probability and its Applications", Kyzyl
В работе строится и исследуется Марковская математическая модель динамики количества детей в семьях Республики Тыва. Дается полное количественное описание параметров и характеристик построенной цепи Маркова, основанное на данных Тывастата РТ. Доказывается, что построенная цепь Маркова является неразложимой и эргодической, найдено ее стационарное распределение и дана оценка скорости сходимости к нему ежегодных абсолютных распределений вероятностей. Произведено прогнозирование динамики изменения количества детей в семьях на ряд лет и сделано сравнение с уже имеющимися статистическими данными.
Ключевые слова: Конечные цепи Маркова, Марковская модель, эргодичность, стационарное распределение, структура семьи, прогнозирование
This article presents the constructed mathematical model of Markov dynamics of the number of children in Tuvan families. Completed quantitative description of the parameters and characteristics of the constructed Markov chain was given, which based on databases of Tuvan Statistics Center of Republic of Tuva. It is proved that the constructed Markov chain is indecomposable and ergodic. The stationary distribution was found and evaluated the rate of convergence to it annual absolute probability distributions. The forecasting dynamics of