РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПРОГНОЗИРОВАНИЯ ЗАБОЛЕВАЕМОСТИ КЛЕЩЕВЫМ ЭНЦЕФАЛИТОМ > Полезные советы
Тысяча полезных мелочей    

РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПРОГНОЗИРОВАНИЯ ЗАБОЛЕВАЕМОСТИ КЛЕЩЕВЫМ ЭНЦЕФАЛИТОМ

РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПРОГНОЗИРОВАНИЯ ЗАБОЛЕВАЕМОСТИ КЛЕЩЕВЫМ ЭНЦЕФАЛИТОМ

Цокова Т.Н. Козлов Л.Б. Разработана математическая модель прогнозирования инфекционной заболеваемости на модели природно-очаговой инфекции, возбудителем которой является вирус клещевого энцефалита. Математическая модель представлена в виде аддитивного временного ряда, включающая тренд, случайные компоненты и сезонные составляющие, имеющие разную периодичность: менее года, 3 года и многолетнюю. Статья в формате PDF 130 KB

Математическое моделирование эпидемических процессов (ЭП) позволяет осуществлять качественный долгосрочный прогноз инфекционной заболеваемости, выявлять факторы, влияющие на динамику ЭП. В 70-х годах ХХ столетия Л.Н. Большев с соавторами [2, 3] разработали математическую модель прогноза заболеваемости клещевым энцефалитом (КЭ). Для построения модели авторы провели предварительный анализ влияния частоты присасывания клещей и иммунной прослойки населения всех возрастов на частоту заболеваний КЭ. Е.И.Болотин с соавторами [1] разработали методику факторного временного прогнозирования эпидемических проявлений очагов КЭ с использованием выявленных критических уровней заболеваемости.

Нами для разработки математического моделирования ЭП использована многолетняя динамика заболеваемости КЭ на юге Тюменской области и составлены временные ряды в климатических подзонах. Математическая модель позволила без предварительного анализа факторов, влияющих на ЭП, осуществлять прогноз заболеваемости КЭ.

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

Данный метод использован нами для изучения динамики ЭП, вызванного вирусом КЭ. Сложность возникла при идентификации временного ряда. При монотонном тренде происходит устойчивое нарастание или убывание значений временного ряда. Анализировать такой ряд обычно нетрудно. В качестве экстраполирующей функции ряда чаще всего выбирают известные математические зависимости - линейные, параболические, экспоненциальные и др., которые часто не отражают динамику ЭП. Прогнозируемые значения такого ряда будут являться оценочными и не обеспечивают точности расчётов.

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

Для прогнозирования заболеваемости КЭ нами использована аддитивная модель и временной ряд был представлен функциональной зависимостью:

y(t) = f(t) + s(t) + ε(t),                       (1)

где y(t) - значение временного ряда в момент времени t; f(t) - основная составляющая ряда (тренд); s(t) - сезонная составляющая (отражает повторяемость показателей заболеваемости в течение сезона года); ε(t) - случайная компонента (отражает неучтённые и случайные факторы, влияющие на ЭП). Временные ряды могут содержать одновременно все перечисленные компоненты или их различные комбинации. Такой ряд называют временным с сезонной составляющей.

Для получения тренда использовано уравнение полинома Чебышева высокого порядка, что позволило снизить погрешности в прогнозе заболевания. Аппроксимирующая функция f(t) теоретически может быть выражена многочлeном любой степени m, например:

                     (2)

При каждом повышении порядка полинома требуется определение не только нового параметра, αm+1, но ввиду изменения системы "нормальных уравнений", проводят пересчет всех остальных параметров: от αo до αm.

Рассмотрев общий случай использования метода наименьших квадратов, П. Л. Чебышев разработал метод вычисления уравнения регрессии, позволяющий определять добавляемый параметр без пересчета найденных ранее параметров, ограничиваясь лишь вычислением нового параметра. Добавляемый члeн имеет вид: , где  определяется по общей формуле:

                          (3)

 В этом случае уравнение регрессии принимает вид [4]:   

                     (4)

После нахождения уравнения тренда вычисляют остаточную дисперсию по формуле:

,                   (5)

где n-m-1 - степень свободы, m - степень полинома, n - объём выборки.

Для получения Sm применяют формулу:

,                     (6)

где - значение ординаты, рассчитанное по уравнению полинома m - степени. Переход к многочлeнам более высокого порядка производят до тех пор, пока остаточная дисперсия продолжает уменьшаться. Если остаточная дисперсия при выравнивании по многочлeну m+1 порядка по сравнению с остаточной дисперсией, полученной для уравнения порядка m, уменьшается незначительно, переход к уравнениям более высокого порядка следует прекратить и аппроксимацию считать достаточной.

Ошибка расчёта (среднеквадратическое отклонение полученной функции от экспериментальных точек) должно быть одного порядка с погрешностью введённых табличных данных, так как среднеквадратическое отклонение зависит от у, n и вида выбранной функции y*. Вычисляют среднеквадратическое отклонение  полученной теоретической кривой y* от экспериментальной у:

                        (7)

и сравнивают с погрешностью эксперимента ε. При этом возможно:

1) δn> ε- аппроксимация слишком грубая, степень полинома необходимо увеличить;

2) δn < ε - аппроксимация физически недостоверна, истинная функция "сплющена", старшие степени полинома лишние и, следовательно, надо уменьшить степень полинома;

3)  δn ≈ ε- степень полинома оптимальна.

Погрешность эксперимента рассчитывали по формуле:

,                 (8)

где yi - наблюдаемые значения,  - среднее значение введённых параметров, n - объём выборки.

Аппроксимирующая функция для строго периодической сезонной составляющей s(t) в уравнении (1) находили, используя тригонометрическую функцию, приподнятую над осью t и сдвинутую вдоль неё, в виде:

s(t) = ao + a1  + b1                (9)

В нашем случае, будучи периодической и заданной в виде таблицы, сезонная составляющая не являлась синусоидой, но была близка к ней. Изменив масштаб по оси t в  раз ( - число полупериодов сезонной составляющей) получили:

s(t) = ao + a1 cos(z) + b1sin(z)                                     (10)

Сжатие графика по оси t не изменит табличные значения s(t) и коэффициентов ao, a1, b1. Коэффициенты ao, a1, b1 определялись из уравнения (10) методом наименьших квадратов. Была получена система из трёх уравнений, после решения которой, получены формулы для расчёта коэффициентов:


b1=

a1=

ao=

zi=                     (11)

Функция s(t) определена на отрезке [1, N] не симметрично точке t=0. Поэтому, после вычисления коэффициентов по уравнениям (11), следует сдвинуть функцию s(t), заменив аргумент  в формуле (9) на выражение:

                                                (12)

Аппроксимация тригонометрическим многочлeном менее точная, чем например степенная иди другая, но к ней прибегают если функция строго периодическая. Для более точного расчёта ao, a1, b1 необходимо точно определить период сезонной составляющей вводимого временного ряда.

Условие δn ≈ ε - степень полинома оптимальна, выполнялось для уравнения тренда с учётом сезонной составляющей.

Весь алгоритм расчётов по формулам 2 - 12 можно автоматизировать. Нами была составлена компьютерная программа, которая позволяет провести полиномиальную аппроксимацию высокого порядка (8 порядка и выше) значений функций заданных таблично. Используя полученное уравнение тренда, проводили прогнозирование заболеваемости. Программы составлены для ПК, работающие в среде Microsoft Windows на программном языке Visual Basik.

Для составления данной программы использована база данных «Tumklech» [8]. Фактическая заболеваемость КЭ была переведена в показатели заболеваемости на 100 тыс. населения. Анализ тренда проводили в два этапа: определяли наличие тренда и выделяли тренд. Для определения наличия тренда использован критерий Стьюдента, позволяющий выявить различие выборочных средних двух половинок временного ряда. Если различие значимо, то гипотеза о наличии тренда не отвергалась. Для выделения тренда использовали модель простого статистического ряда и аддитивную модель временного ряда (временной ряд с сезонной составляющей).

В динамике заболеваемости КЭ хаpaктерна трехлетняя цикличность и в течение сезона года (с 1.04 по 30.09) также наблюдались колебания в показателях заболеваемости. Усреднив значения ряда в соответствии с трёхлетней цикличностью, получим сезонные колебания заболеваемости для каждого года трехлетнего цикла [7]. Однако для долгосрочного прогноза, предложенный способ прогноза будет недостаточно точным.

Известно, что в природе существуют циклические колебания погоды, связанные с цикличностью излучения Солнца. Большой цикл Солнца длится 11 лет. Он связан с цикличностью образования солнечных пятен. 1979-1982 годы максимальной активности; 1986-1987 годы минимальной активности. В настоящее время Солнечная активность минимальная [6].

Методом дисперсионного анализа нами также установлено наличие многолетнего цикла в динамике заболеваемости КЭ. В 1999г. наблюдалась максимальная заболеваемость КЭ во всех климатических подзонах юга Тюменской области.

Применение многолетней цикличности к нашим расчётам, означает усложнение модели временного ряда ещё одной составляющей, зависящей от многолетней цикличности. В этом случае прогноз модели будет точнее. Условно назовём эту модель - аддитивная модель с двумя составляющими. Уравнение этой модели:

y(t) =f(t) + e1(t) + e2(t),                                                 (13)

где ε1(t) - сезонная составляющая с трёхлетним периодом, ε2(t) -сезонная составляющая с многолетним периодом. Данная модель может быть использована для долгосрочного прогнозирования заболеваемости КЭ.

СПИСОК ЛИТЕРАТУРЫ:

  1. Болотин Е.И., Цициашвили Г.Ш., Голычева И.В. // Паразитология. - 2002. - Т.36.- Вып. 2. - С.89.
  2. Большев Л.Н., Гольдфарб Л.Г. // Вопросы эпидемиологии и профилактики клещевого энцефалита. - М., 1970. - С. 154.
  3. Большев Л.Н., Гольдфарб Л.Г., Круопис Ю.И. // Вопросы эпидемиологии и профилактики клещевого энцефалита. - М., 1970. - С. 171.
  4. Венецкий И.Г., Кильдишев Г.С. Теория вероятности и математическая статистика. - М., «Статистика». - 1975. - 264 с.
  5. Ивашев-Мусатов О.С. Теория вероятности и математическая статистика. - М., «Наука». - 1979. - 254 с.
  6. Ишков В.Н. «Вселенная и мы» http://www.астронет.ru/ / Солнце в текущем 23 цикле солнечной активности.
  7. Козлов Л.Б., Кашуба Э.А., Цокова Т.Н. и др. Способ прогноза заболеваемости клещевыми инфекциями // Патент RU 2294697 С2. - Бюл. № 7 от 10.03.2007г.
  8. Козлов Л.Б., Цокова Т.Н., Огурцов А.А. и др. / Заболеваемость клещевым энцефалитом в Тюменской области «Tumklech» // База данных. - Свидетельство №2007620363 от 18.10.2007. - Правообладатель: ФГУЗ «Центр гигиены и эпидемиологии в Тюменской области».


ГЛУЩЕНКО ЛЮДМИЛА ФЁДОРОВНА

ГЛУЩЕНКО ЛЮДМИЛА ФЁДОРОВНА Статья в формате PDF 175 KB...

17 06 2026 1:14:35

Договор купли продажи

Договор купли продажи Статья в формате PDF 103 KB...

16 06 2026 10:17:55

АВЕРЬЯНОВ ПЕТР ФЕДОРОВИЧ

АВЕРЬЯНОВ ПЕТР ФЕДОРОВИЧ Статья в формате PDF 82 KB...

15 06 2026 12:23:11

МИНИМИЗАЦИЯ АППАРАТУРЫ ДЛЯ ТЕРМИЧЕСКОГО АНАЛИЗА

МИНИМИЗАЦИЯ АППАРАТУРЫ ДЛЯ ТЕРМИЧЕСКОГО АНАЛИЗА Статья в формате PDF 171 KB...

13 06 2026 16:25:17

ФАКТОР НЕКРОЗА ОПУХОЛЕЙ И ГИБЕЛЬ КЛЕТОК ПЕЧЕНИ

ФАКТОР НЕКРОЗА ОПУХОЛЕЙ И ГИБЕЛЬ КЛЕТОК ПЕЧЕНИ Статья в формате PDF 97 KB...

12 06 2026 17:40:44

ПСИХОЛОГИЧЕСКАЯ ВЕКТОРНАЯ МОДЕЛЬ ЗДОРОВЬЕСОЗИДАЮЩЕГО ПОТЕНЦИАЛА ЛИЧНОСТИ

ПСИХОЛОГИЧЕСКАЯ ВЕКТОРНАЯ МОДЕЛЬ ЗДОРОВЬЕСОЗИДАЮЩЕГО ПОТЕНЦИАЛА ЛИЧНОСТИ В данной статье осуществлены анализ и обощение зарубежных психологических концепций, объясняющих активность человека в отношении своего здоровья, и на этой основе дано авторское определение понятия «здоровьесозидающий потенциал личности». Особое внимание авторы уделяют рассмотрению структуры здоровьесозидающего потенциала, описанию психологических механизмов его формирования и выявлению закономерностей его развития в различные возрастные периоды. Авторами впервые представлена и научно обоснована векторная модель здоровьесозидающего потенциала личности, показаны её возможности при выявлении психологических детерминант, влияющих на показатели целостного здоровья человека. ...

11 06 2026 1:27:51

РАЗРАБОТКА СОСТАВА, ТЕХНОЛОГИЯ И СТАНДАРТИЗАЦИЯ ОФТАЛЬМОЛОГИЧЕСКИХ ЛЕКАРСТВЕННЫХ ФОРМ С ОРТОФЕНОМ

РАЗРАБОТКА СОСТАВА, ТЕХНОЛОГИЯ И СТАНДАРТИЗАЦИЯ ОФТАЛЬМОЛОГИЧЕСКИХ ЛЕКАРСТВЕННЫХ ФОРМ С ОРТОФЕНОМ Предложены офтальмологические лекарственные формы с ортофеном - глазные лекарственные пленки и пролонгированные глазные капли. Разработан их состав, технология длч производства в аптечных условиях, стандартизация. Проведены подробные биофармацевтические исследования in vitro по выбору оптимальных вспомогательных компонентов. Выбран способ количественного анализа ортофена в разработанных лекарственных форм - с помощью спектрофотометрии. ...

06 06 2026 16:33:43

ЭТНОЭПИДЕМИОЛОГИЧЕСКИЕ ПРОБЛЕМЫ ДАЛЬНЕГО ВОСТОКА

ЭТНОЭПИДЕМИОЛОГИЧЕСКИЕ ПРОБЛЕМЫ ДАЛЬНЕГО ВОСТОКА Статья в формате PDF 102 KB...

03 06 2026 18:51:51

РАЗБИЕНИЕ СТРУКТУРИРОВАННОГО 3D ПРОСТРАНСТВА НА МОДУЛЯРНЫЕ ЯЧЕЙКИ И МОДЕЛИРОВАНИЕ НЕВЫРОЖДЕННЫХ МОДУЛЯРНЫХ СТРУКТУР

РАЗБИЕНИЕ СТРУКТУРИРОВАННОГО 3D ПРОСТРАНСТВА НА МОДУЛЯРНЫЕ ЯЧЕЙКИ И МОДЕЛИРОВАНИЕ НЕВЫРОЖДЕННЫХ МОДУЛЯРНЫХ СТРУКТУР Обсуждаются разбиения 3D прострaнcтва на модулярные ячейки с целью последующего конструирования невырожденных модулярных 3D структур кристаллов. ...

02 06 2026 9:59:42

МЕТАБОЛИЧЕСКАЯ РЕКОНСТРУКЦИЯ КОРРЕКЦИИ КАРДИОСКЛЕРОЗА МЕДЛЕННЫМ УПРАВЛЕНИЕМ ДИНАМИЧЕСКИМ ХАОСОМ

МЕТАБОЛИЧЕСКАЯ РЕКОНСТРУКЦИЯ КОРРЕКЦИИ КАРДИОСКЛЕРОЗА МЕДЛЕННЫМ УПРАВЛЕНИЕМ ДИНАМИЧЕСКИМ ХАОСОМ В отличие от известной методики математика Лоскутова и кардиолога Ардашева по лечению тяжелых аритмий хирургической абляцией зон хаотизации сердца, в работе предложен в форме биоинформационной и математической моделей подход для терапевтической, противовоспалительной методики снижения хаотизации. Проведена метаболическая реконструкция патобиохимии кардиосклероза и его коррекции. Кардиосклероз рассматривается, как аутовоспалительный процесс на базе медленного (недели, месяцы) «неправильного» взаимодействия депо углеводов и жиров. Расчеты показывают, что при медленных (годы) сценариях тренировки сердца и защите его от свободных радикалов и воспалений при стрессе цитопротекторами и пептидотерапией, могут возникать мультициклы, обеспечивающие снижение хаоса. Это создает условия прекондиционирования, тесно связанные с условиями для обновления клеток в сердце. ...

01 06 2026 17:34:34

МАГНИТНЫЕ ПОДРЕШЕТКИ, ИНДУЦИРОВАННЫЕ КАТИОННЫМИ ВАКАНСИЯМИ (НА ПРИМЕРЕ ФЕРРИМАГНИТНОГО ПИРРОТИНА)

МАГНИТНЫЕ ПОДРЕШЕТКИ, ИНДУЦИРОВАННЫЕ КАТИОННЫМИ ВАКАНСИЯМИ (НА ПРИМЕРЕ ФЕРРИМАГНИТНОГО ПИРРОТИНА) На основе анализа s-d обменного взаимодействия в структурах типа NiAs с частично вакантными катионными позициями, моделировались различного рода зависимости результирующей намагниченности от температуры нестехиометрических ферримагнетиков. На основе исследований пирротина методами ЯГР и РФА доказано, что двухподрешеточный ферримагнетик, содержащий в структуре катионные вакансии, должен рассматриваться, при определенном типе распределения вакансий, как ферримагнетик с четырьмя магнитными подрешетками. В данном случае, дополнительные магнитные подрешетки можно рассматривать как подрешетки, индуцированные хаpaктером распределения катионных вакансий в структуре. Квантово-механические расчеты в рамках модели молекулярного поля температурных изменений намагниченности отдельно для каждой из подрешеток, а также анализ результирующей термокривой намагниченности, объясняют ряд экспериментально полученных кривых зависимости намагниченности от температуры нестехиометрического пирротина с различной плотностью вакансий в структуре. ...

31 05 2026 19:16:20

ФЕРМЕНТНЫЕ ПРЕПАРАТЫ ИЗ ВНУТРЕННОСТЕЙ ПРУДОВЫХ РЫБ

ФЕРМЕНТНЫЕ ПРЕПАРАТЫ ИЗ ВНУТРЕННОСТЕЙ ПРУДОВЫХ РЫБ Статья в формате PDF 259 KB...

30 05 2026 18:59:53

ФОРМА И ТОПОГРАФИЯ ПОДЖЕЛУДОЧНОЙ ЖЕЛЕЗЫ У КРЫСЫ

ФОРМА И ТОПОГРАФИЯ ПОДЖЕЛУДОЧНОЙ ЖЕЛЕЗЫ У КРЫСЫ Поджелудочная железа белой крысы имеет три основные части – головка (дуоденальная часть), тело (пилорическая часть) и хвост (желудочно-селезеночная часть). По сравнению с человеком, она отличается большей рыхлостью, изогнутостью, разветвленностью. Встречаются два крайних варианта формы (в виде молотка или трилистника) и топографии поджелудочной железы у белой крысы. ...

29 05 2026 23:25:20

НЕКОТОРЫЕ АСПЕКТЫ НАРКОПРЕСТУПНОСТИ

НЕКОТОРЫЕ АСПЕКТЫ НАРКОПРЕСТУПНОСТИ Статья в формате PDF 251 KB...

22 05 2026 10:49:36

ОЧИСТКА СТОЧНЫХ ВОД НА АВТОМОЙКЕ

ОЧИСТКА СТОЧНЫХ ВОД НА АВТОМОЙКЕ Статья в формате PDF 254 KB...

13 05 2026 0:44:45

СИСТЕМНАЯ МЕДИЦИНА В САНАТОРНО-КУРОРТНОЙ ПРАКТИКЕ

СИСТЕМНАЯ МЕДИЦИНА В САНАТОРНО-КУРОРТНОЙ ПРАКТИКЕ Статья в формате PDF 144 KB...

12 05 2026 21:12:20

МИГРАЦИЯ ТЯЖЕЛЫХ МЕТАЛЛОВ ЧЕРЕЗ ВРЕДНЫЕ ПРИВЫЧКИ

МИГРАЦИЯ ТЯЖЕЛЫХ МЕТАЛЛОВ ЧЕРЕЗ ВРЕДНЫЕ ПРИВЫЧКИ Статья в формате PDF 110 KB...

11 05 2026 22:17:44

Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::