close

Вход

Забыли?

вход по аккаунту

?

Патент BY3546

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
(19)
BY (11) 3546
(13)
C1
(51)
(12)
6
A 61B 5/08,
G 06F 17/00,
G 06F 159:00
ГОСУДАРСТВЕННЫЙ ПАТЕНТНЫЙ
КОМИТЕТ РЕСПУБЛИКИ БЕЛАРУСЬ
54)
CПОСОБ ПРОГНОЗИРОВАНИЯ ОБОСТРЕНИЙ БРОНХИАЛЬНОЙ
АСТМЫ В ЗАВИСИМОСТИ ОТ ПОГОДНО-КЛИМАТИЧЕСКИХ
УСЛОВИЙ
(21) Номер заявки: 960487
(22) 1996.10.01
(46) 2000.09.30
(71) Заявители: Витебский филиал НИИ радиационной медицины Минздрава РБ, Минский
Республиканский эндокринологический диспансер
НИИ радиационной медицины Минздрава РБ
(BY)
(72) Авторы: Доценко Э.А., Прищепа И.М., Доценко
М.Л., Василевская Т.А. (BY)
(73) Патентообладатели: Витебский
филиал
НИИ
радиационной медицины Минздрава РБ,
Минский
Республиканский
эндокринологический диспансер НИИ радиационной
медицины Минздрава РБ (BY)
(57)
Способ прогнозирования обострений бронхиальной астмы в зависимости от погодно-климатических условий, включающий ретроспективный анализ силы влияния комплексных погодных показателей температуры, атмосферного давления и влажности на больных бронхиальной астмой, математического моделирования
и заключения о характере дня, отличающийся тем, что ретроспективный анализ проводят в течение не менее 24 месяцев, вычисляют дискриминантные функции и координаты групповых центроидов для благоприятной и неблагоприятной для астматиков погоды, используют показатели конкретного дня, рассчитывают
новые значения дискриминантных функций, соотносят их с полученными координатами групповых центроидов для астмоопасных и астмобезопасных дней, при этом выделяют интервал “зимние месяцы” и при
значении дискриминантной функции изучаемого дня (-0,414) и менее определяют день как астмоопасный,
при значении более (-0,414) - астмобезопасный, выделяют интервал “весенние месяцы” и при значении дискриминантной функции изучаемого дня (+1,057) и более определяют день как астмоопасный, при значении
менее (+1,057) - астмобезопасный, выделяют интервал “летние месяцы” и при значении дискриминантной
функции изучаемого дня (-0,284) и менее определяют день как астмоопасный, при значении более (0,284) - астмобезопасный, выделяют интервал “осенние месяцы” и при значении дискриминантной функции
изучаемого дня (-0,627) и менее определяют день как астмоопасный, при значении более (-0,627) - астмобезопасный.
(56)
1. Федосеев Г.Б., Убайдуллаев А.М. Ранняя диагностика и первичная профилактика бронхиальной астмы.
- Ташкент: Медицина, 1989. - С. 76-87.
Изобретение относится к медицине, в частности к пульмонологии и аллергологии, и может быть использовано для прогноза обострений бронхиальной астмы.
В качестве прототипа выбрана математическая модель обострений БА в зависимости от погодных условий, предложенная Г.Б. Федосеевым и А.М. Убайдуллаевым [1]. Способ основан на построении модели
множественной регрессии, где в качестве входных факторов использованы погодные показатели, а выходного - число обращений больных БА за экстренной медицинской помощью. Общими с прототипом признаками
являются:
1. Объект исследования - бронхиальная астма;.
2. Формирование когорты наблюдения на основе данных обращаемости за экстренной медицинской помощью;
3. Использование методов многомерной статистики;
4. Общность показателей физического состояния атмосферы.
BY 3546 C1
Описанный способ имеет следующие недостатки:
1. Небольшой срок наблюдения (до 1 года), что не позволяет обнаружить глобальные закономерности;
2. Отсутствует оценка воспроизводимости моделей;
3. Построение моделей для оценки течения, а не прогноза БА, так как не введено понятие “пик” обострений БА и невозможно оценить - много обострений БА или мало;
4. Проведение анализа по годичному интервалу, в то время как известно, что в отдельные сезоны (зима,
весна, лето, осень) влияние погодных условий на течение болезни различно;
5. Для построения моделей берутся погодные показатели в тот же день, что и обострения БА, что неприменимо в практике, поскольку не позволяет проводить профилактические лечебные мероприятия.
Задачей изобретения является повышение эффективности прогноза течения БА в зависимости от погодных условий.
Предлагаемый способ основан на нескольких этапах.
1 этап. Формирование компьютерной базы данных, которая включает:
1) число обращений больных БА за экстренной медицинской помощью в течение не менее чем 2-х лет;
2) погодно-климатические показатели, включающие ежедневно оцениваемые: средняя температура воздуха (X1); максимальная температура воздуха (Х6); минимальная температура воздуха (Х7); количество
осадков в 2100 (Х9); количество осадков в 900 (X10); количество осадков за сутки (X11); скорость ветра в 2100
(X12); скорость ветра в 900 (Х13); средняя относительная влажность (Х18); минимальная относительная
влажность (Х22); атмосферное давление (Х24); направление ветра (Х28); изменения абсолютных показателей: градиент средней температуры (Х3); модуль градиента средней температуры (Х4); разница между максимальной и минимальной температурой воздуха в пределах одного дня (Х8); градиент изменения скорости
ветра за ночь (X14); модуль изменения скорости ветра за ночь (Х15); градиент изменения скорости ветра за
день (X16); модуль изменения скорости ветра за день (X17); градиент изменения средней влажности (Х20);
модуль градиента средней влажности (Х21); градиент изменения атмосферного давления за сутки (Х25); модуль градиента изменения атмосферного давления (Х26), комплексный коэффициент жесткости погоды
(Х23) [1], в котором учитываются среднесуточная температура воздуха, среднесуточная скорость ветра, относительная влажность и суточная изменчивость температуры воздуха, индекс патогенности температуры
воздуха (Х2); индекс патогенности влажности (X19); индекс патогенности изменения атмосферного давления
(Х27); индекс патогенности скорости ветра (Х30); индекс патогенности изменения температуры воздуха (Х5);
индекс патогенности метеорологических условий (Х31) вычисляется путем суммации частных индексов патогенности. Учитывались также погодные показатели за 1 (Х41-Х431), 2 (Х51-Х531), 3 (X61-X631), 4 (X71-X731),
5 (X1-X831), 6 (X91-X931) дней до вызова скорой помощи.
Формирование базы проводят в интервале не менее чем 24 месяца, для того, чтобы включить в исследование
года с различными погодно-климатическими режимами, например 1991 и 1993 гг.
2 этап. Проводят однофакторный дисперсионный анализ силы влияния погодных показателей на частоту
обострений БА и отбирают те из них, влияние которых статистически значимо. Рассматривают погодные
факторы в дни, на 3-е суток предшествующих обострению БА. При анализе указанной базы были отобраны
показатели Х61, Х63, Х64, Х66, Х68 (температурные показатели), Х618, Х620 (показатели влажности),
Х624, Х625, Х626 (показатели давления) со сдвигом на 3 дня.
3 этап. Строят классификационный ряд обострений БА, для чего рассчитывают среднегодичное число
обострений БА за 2 года, с которым арифметически суммируют 1,5S (S - среднеквадратичное отклонение).
Дни, в которых количество вызовов было меньше (Хср+1,5S), считают благоприятными (“1”), больше - неблагоприятными для астматиков (“2”) (Приложение 1).
4 этап. С помощью прикладной статистической программы (например, “Statgrafics 3.1”) рассчитывают
коэффициенты дискриминантных функций (Приложение 2), координаты групповых центроидов и срединной
точки, делящей расстояние между ними строго пополам (Приложение 3), и эффективность отнесения дней
по классам (Приложение 4).
5 этап. В уравнения дискриминантных функций подставляют значения погодных факторов и рассчитывают ее числовое значение. В зависимости от приближения полученного результата к тому или иному центроиду (что оценивают путем его сравнения со срединной точкой, делящей расстояние между центроидами
строго пополам) выносят суждение о принадлежности дня к астмобезопасному (1) и астмоопасному (2).
п. 1. Для повышения точности из указанной базы выделяют осенние периоды, что позволяет существенно
повысить эффективность прогнозирования.
1 этап. Из общей базы данных выделен осенний интервал (сентябрь, октябрь, ноябрь 1991 и 1993 гг.).
2 этап. Проведен однофакторный дисперсионный анализ силы влияния погодных факторов на число обострений БА. Статистически значимы оказались показатели Х61, Х63, Х64, Х66, Х68 (температурные показатели), Х618, Х620 (показатели влажности), Х624, Х625, Х626 (показатели давления) со сдвигом на 3 дня.
3 этап. Строим классификационный ряд: среднее число вызовов составило 18,5, S – 7,4, X+1,5S – 30. При
превышении числа обращений за экстренной медицинской помощью данного значения дню присваивается
номер “2”, при меньшем - “1”.
2
BY 3546 C1
4 этап. Рассчитывают коэффициенты дискриминантных функций для осеннего периода (Приложение 2),
координаты групповых центроидов и срединной точки (Приложение 3) и эффективность отнесения к классам (Приложение 4).
5 этап. В уравнение дискриминантной функции для осеннего периода подставляют количественные значения погодных факторов и получают ее значение. При приближении его к координатам центроида “1” день
отмечают как благоприятный, “2” - неблагоприятный.
Пример 1.
День 6.09.91 г. В осенний сезон дискриминантная функция имеет следующий вид.
Ф=0,280*Х61-0,028*Х63-0,063*Х64-0,365*Х66-0,007*Х68-0,025*Х618-0,018*Х6200,051*Х624+0,071Х625-0,071*Х626+42,416
Величины показателей с 3-х дневным лагом (то есть 3.09.91 г.) были следующие: Х61 - 17,2, Х63 - 2,8,
Х64 - 2,8, Х66 - 23, X68 - 11,1, X618 - 68,5, X620 - 2, X624 - 749, X625 - (-6,5), X626 - 6,5. Вычисленная дискриминантная функция составила (-2,366), что меньше, чем (-0,627) и наиболее близко к центроиду (-1,426);
следовательно, предсказанный класс для 6.09.91 г. равен 2, это астмоопасный день. Истинное число вызовов
составило 31, что является превышением Xcp+1,5S и относится к астмоопасному дню.
Пример 2.
День 16.09.91 г. В осенний сезон дискриминантная функция имеет следующий вид.
Ф=0,280*Х61-0,028*Х63-0,063*Х64-0,365*Х66-0,007*Х68-0,025*Х618-0,018*Х6200,051*Х624+0,071Х625-0,071*Х626+42,416
Величины первичных показателей с 3-дневным лагом (то есть 13.09.91 г.) были следующие: Х61 - 8,3,
Х63 - -2,6, Х64 - 2,6, Х66 - 12,4, X68 - 7,2, X618 - 55, X620 - (-24), X624 - 751, X625 - 0, X626 - 0. Вычисленная дискриминантная функция составила (0,829), что больше (-0,627) и наиболее близко к центроиду
(+0,171); следовательно, предсказанный класс для 16.09.91 г. равен 1, это астмобезопасный день. Истинное
число вызовов составило 24, что не является превышением Xcp+1,5S и относится к астмобезопасному дню.
п.2. Для повышения точности описанный анализ проводят в течение зимних периодов по указанной базе
данных, что позволяет существенно повысить эффективность прогнозирования.
1 этап. Из общей базы данных выделен зимний интервал (декабрь, февраль, январь 1991 и 1993 гг.).
2 этап. Проведен однофакторный дисперсионный анализ силы влияния погодных факторов на число обострений БА. Статистически значимы оказались показатели Х61, Х63, Х64, Х66, X68 (температурные показатели), X618, X620 (показатели влажности), X624, X625, X626 (показатели давления) со сдвигом на 3 дня.
3 этап. Строим классификационный ряд: среднее число вызовов составило 23,7, S-4,39, X+1,5S-30. При
превышении числа обращений за экстренной медицинской помощью данного значения дню присваивается
номер “2”, при меньшем - “1”.
4 этап. Рассчитывают коэффициенты дискриминантных функций для зимнего периода (Приложение 2),
координаты групповых центроидов и срединной точки (Приложение 3) и эффективность отнесения к классам (Приложение 4).
5 этап. В уравнение дискриминантной функции для зимнего периода подставляют количественные значения погодных факторов и получают ее значение. При приближении его к координатам центроида “1” день
отмечают как благоприятный, “2” - неблагоприятный.
Пример 3.
День 5.02.91 г. В зимний сезон дискриминантная функция имеет следующий вид.
Ф=0,334*Х61-0,038*X63-0,083*X64-0,181*X66-0,049*Х68+0,028*X618+
0,031*X620+0,076*X624+0,044Х625+0,163*X626-58,966
Величины показателей с 3-дневным лагом (то есть 2.02.91 г.) были следующие: Х61 - (-13,3), Х63 - 2,4, Х64 2,4, Х66 - (-9,2), X68 - (7,9), Х618 - 85,5, X620 - 2,5, X624 - 762, X625 - 1, X626 - 1. Вычисленная дискриминантная функция составила (-1,83), что меньше (-0,414) и наиболее близко к центроиду (-0,924); следовательно,
предсказанный класс для 5.02.91 г. равен 2, это астмоопасный день. Истинное число вызовов составило 31,
что является превышением Xcp+1,5S и относится к астмоопасному дню.
Пример 4.
День 8.01.91 г. В зимний сезон дискриминантная функция имеет следующий вид.
Ф=0,334*Х61-0,038*X63-0,083*X64-0,181*X66-0,049*Х68+0,028*X618+
0,031*X620+0,076*X624+0,044Х625+0,163*X626-58,966
Величины первичных показателей с 3-дневным лагом (то есть 5.01.91 г.) были следующие: Х61 - 0,6, Х63
- (-1,2), Х64 - 1,2, Х66 - 1,6, X68 - 2,6, Х618 - 91,5, X620 - 4,5, X624 - 739, X625 - (-1), X626 - 1. Вычисленная
дискриминантная функция составила (-0,252), что больше (-0,414) и наиболее близко к центроиду (0,095);
следовательно, предсказанный класс для 8.01.91 г. равен 1, это астмобезопасный день. Истинное число вызовов составило 20, что не является превышением Xcp+1,5S и относится к астмобезопасному дню.
п.3. Для повышения точности описанный анализ проводят в течение весенних периодов по указанной базе данных, что позволяет существенно повысить эффективность прогнозирования.
3
BY 3546 C1
1 этап. Из общей базы данных выделен весенний интервал (март, апрель, май 1991 и 1993 гг.).
2 этап. Проведен однофакторный дисперсионный анализ силы влияния погодных факторов на число обострений БА. Статистически значимы оказались показатели Х61, Х63, Х64, Х66, X68 (температурные показатели), Х618, X620 (показатели влажности), X624, X625, X626 (показатели давления) со сдвигом на 3 дня.
3 этап. Строим классификационный ряд: среднее число вызовов составило 18,3, S – 5,2, X+1,5S - 26. При
превышении числа обращений за экстренной медицинской помощью данного значения дню присваивается
номер “2”, при меньшем - “1”.
4 этап. Рассчитывают коэффициенты дискриминантных функций для весеннего периода (Приложение 2),
координаты групповых центроидов и срединной точки (Приложение 3) и эффективность отнесения к классам (Приложение 4).
5 этап. В уравнение дискриминантной функции для весеннего периода подставляют количественные значения погодных факторов и получают ее значение. При приближении его к координатам центроида “1” день
отмечают как благоприятный, “2” - неблагоприятный.
Пример 5.
День 15.03.91 г. Весной дискриминантная функция имеет следующий вид.
Ф=0,375*Х61-0,011*Х63-0,438*Х64-0,435*Х66+0,191*Х68-0,011*Х618+ 0,009*Х620-0,025*Х6240,148*Х625-0,019*Х626+20,590
Величины показателей с 3-дневным лагом (то есть 12.03.91 г.) были следующие: Х61 - 0,1, Х63 - 1,3, Х64
- 1,3, Х66 - 1,5, X68 - 16,5, Х618 - 85, X620 - 38, X624 - 746, X625 - (-4), X626 - 4. Вычисленная дискриминантная функция составила (3,8158), что больше (+1,057) и наиболее близко к центроиду (2,255); следовательно, предсказанный класс для 15.03.91 г. равен 2, это астмоопасный день. Истинное число вызовов составило 44, что является превышением Xcp+1,5S и относится к астмоопасному дню.
Пример 6.
День 8.03.91 г. Весной дискриминантная функция имеет следующий вид.
Ф=0,375*Х61-0,011*Х63-0,438*Х64-0,435*Х66+0,191*Х68-0,011*Х618+ 0,009*Х620-0,025*Х6240,148*Х625-0,019*Х626+20,590
Величины первичных показателей с 3-дневным лагом (то есть 5.03.91 г.) были следующие: Х61 - (-0,2),
Х63 - 2, Х64 - 2, Х66 - 1,6, X68 - 2,6, X618 - 85, X620 - 24, X624 - 750, X625 - (-6), X626 – 6. Вычисленная
дискриминантная функция составила (0,7226), что меньше (+1,057) и наиболее близко к центроиду (-0,110);
следовательно, предсказанный класс для 8.03.91 г. равен 1, это астмобезопасный день. Истинное число вызовов составило 21, что не является превышением Хср+1,5S и относится к астмобезопасному дню.
п.4. Для повышения точности, описанный анализ проводят в течение летнего периода по указанной базе
данных, что позволяет существенно повысить эффективность прогнозирования.
1 этап. Из общей базы данных выделен летний интервал (июнь, июль, август 1991 и 1993 гг.).
2 этап. Проведен однофакторный дисперсионный анализ силы влияния погодных факторов на число обострений БА. Статистически значимы оказались показатели Х61, Х63, Х64, Х66, X68 (температурные показатели), X618, X620 (показатели влажности), X624, X625, X626 (показатели давления) со сдвигом на 3 дня.
3 этап. Строим классификационный ряд: среднее число вызовов составило 18,5, S – 5,5, X+1,5S - 27. При
превышении числа обращений за экстренной медицинской помощью данного значения дню присваивается
номер “2”, при меньшем - “1”.
4 этап. Рассчитывают коэффициенты дискриминантных функций для летнего периода (Приложение 2),
координаты групповых центроидов и срединной точки (Приложение 3) и эффективность отнесения к классам (Приложение 4).
5 этап. В уравнение дискриминантной функции для летнего периода подставляют количественные значения погодных факторов и получают ее значение. При приближении его к координатам центроида “1” день
отмечают как благоприятный, “2” - неблагоприятный.
Пример 7.
День 15.07.91 г. В летний сезон дискриминантная функция имеет следующий вид.
Ф=0,582*Х61+0,253*Х63-0,168*Х64-0,480*Х66+0,103*Х68+0,028*Х618-0,003*Х6200,070*Х624+0,068Х625-0,031*Х626-54,574
Величины показателей с 3-дневным лагом (то есть 12.07.91 г.) были следующие: Х61 - 15,8, Х63 - (-0,7),
Х64 - 0,7, Х66 - 21, X68 - 8,4, X618 - 78, X620 - 9, X624 - 746, X625 - 0, X626 - 0. Вычисленная дискриминантная функция составила (-0,5109), что меньше (-0,284) и наиболее близко к центроиду (-0,634); следовательно, предсказанный класс для 15.07.91 г. равен 2, это астмоопасный день. Истинное число вызовов составило 32, что является превышением Хср+1,5S и относится к астмоопасному дню.
4
BY 3546 C1
Пример 8.
Рассмотрим 22.06.91г. В летний сезон дискриминантная функция имеет следующий вид.
Ф=0,582*Х61+0,253*Х63-0,168*Х64-0,480*Х66+0,103*Х68+0,028*Х618-0,003*Х6200,070*Х624+0,068Х625-0,031*Х626-54,574
Величины первичных показателей с 3-дневным лагом (то есть 19.06.91 г.) были следующие: Х61 - 20,3,
Х63 - 0,3, Х64 - 0,3, X66-25, X68-7,1, X618 - 68, X620 - (-16), X624 - 747, X625 - 1, X626 - 1. Вычисленная
дискриминантная функция составила (0,2764), что больше (-0,284) и наиболее близко к центроиду (0,066);
следовательно, предсказанный класс для 22.06.91 г. равен 1, это астмобезопасный день. Истинное число вызовов составило 25, что не является превышением Xcp+1,5S и относится к астмобезопасному дню.
Приложение 1
Средние и стандартные отклонения (S) числа обострений БА в зависимости от базы данных
Сезон
Зима
Весна
Лето
Осень
Хср
23,7
18,3
18,5
18,5
S
Границы классов
1
0-30
0-26
0-27
0-30
4.4
5.2
5.5
7.4
2
>30
>26
>27
>30
Приложение 2
Коэффициенты дискриминантных функций (по первичным погодным показателям)
Показатель
Х61
Х63
Х64
Х66
Х68
Х618
Х620
Х624
Х625
Х626
Константа
Зима
0,334
-0,038
-0,083
-0,181
-0,049
0,028
0,031
0,076
0,044
0,163
-58,966
Весна
0,375
-0,011
-0,438
-0,435
0,191
-0,011
0,009
-0,025
-0,148
-0,019
20,590
Лето
0,582
0,253
-0,168
-0,480
0,103
0,028
-0,003
0,070
0,068
-0,031
-54,574
Осень
0,280
-0,028
-0,063
-0,365
-0,007
-0,025
-0,018
-0,051
0,071
-0,071
42,416
Приложение 3
Координаты групповых центроидов (входные параметры - первичные погодные показатели)
Сезон
Зима
Весна
Лето
Осень
Координата центроида
1
2
0,095
-0,924
-0,110
2,255
0,066
-0,634
0,171
-1,426
5
Координата срединной
точки
-0,414
+1,057
-0,284
-0,627
BY 3546 C1
Приложение 4
Эффективность отнесения к классам обострения с учетом общего числа наблюдений (анализ проведен на основе первичных погодных показателей)
Истинный класс
Истинный класс
1
Зима
2
1
Весна
2
1
Лето
2
1
Осень
2
Предсказанный класс
1
2
68,2
22,3
(73)
(24)
2,7
6,8
(3)
(7)
84,2
11,1
(91)
(12)
1,9
2,8
(2)
(3)
62,4
28,2
(73)
(33)
4,3
5,1
(5)
(6)
69,6
19,6
(78)
(22)
3,8
7,1
(4)
(8)
Государственный патентный комитет Республики Беларусь.
220072, г. Минск, проспект Ф. Скорины, 66.
6
Число (%) общих совпадений
75,0
87,1
67,5
76,7
Документ
Категория
Без категории
Просмотров
0
Размер файла
156 Кб
Теги
by3546, патент
1/--страниц
Пожаловаться на содержимое документа