CC BY-ND
24
ИСПОЛЬЗОВАНИЕ МЕТОДА ПОСТРОЕНИЯ ТАБЛИЦ ЖИЗНИ ДЛЯ ПРОГНОЗИРОВАНИЯ РИСКА РАЗВИТИЯ ПРОФЕССИОНАЛЬНОЙ ПАТОЛОГИИ
Л.Н. Будкарь, И.В. Бугаева, Л.Г. Терешина. А.Г.Столяр, Т.Ю.Обухова, Е.А. Карпова,
USE OF LIFE TABLES FOR PREDICTION OF OCCUPATIONAL PATHOLOGY HAZARD RATE
L.N. Budkar, I.V. Bugaeva, L.G. Tereshina, A.G.Stolar, T.U. Obuchova, E.A. Carpova ФГУН ЕМНЦ ПОЗРПП Роспотребнадзора, г. Екатеринбург
Документальное фиксирование всех этапов развития заболевания в зависимости от пылевого стажа работника в клинике профессиональных болезней позволяет использовать метод построения таблиц жизни. Прогнозируется постановка диагноза силикоз у 50% наблюдаемых работников через 29 лет пылевого стажа. Риск развития профессиональной патологии легких на протяжении 40 летнего контакта с неблагоприятным производственным фактором не превышает 20%.
Documenting all stages of disease progression related to the dust exposure records of the employees observed at the clinic of occupational diseases enabled its specialists to draw the patients' life tables. Prediction of silicosis to be diagnosed after a 29-year exposure to dust in 50% in the employees observed was made. The hazard rate of occupational lung pathology development after a 40-year exposure to dust does not exceed 20%.
Особенность наблюдения за работниками, подвергающихся воздействию неблагоприятных производственных факторов, заключается в том, что на протяжении ряда лет (иногда более 20) пациент систематически обследуется, данные состояния здоровья фиксируются в историях болезни (ИБ). Кроме того, в ИБ собираются сведения (санитарно-гигиеническая характеристика, копия трудовой книжки), в которых представлены характеристики условий труда работника. В результате для каждого работника документируются сроки и динамика развития заболевания — начала контакта с вредным фактором производства, появления первых признаков патологии, формирования диагностических изменений, время установления окончательного диагноза, декомпенсации со-
стояния и в некоторых случаях сроков инвали-дизации. Таким образом, в клинике профессиональных болезней имеется уникальная возможность документально отследить все этапы развития патологии для каждого пациента, что служит основой для применения метода построения таблиц жизни (ТЖ).
Во всех случаях наблюдений, когда исследователя интересует время и/или вероятность наступления определенного, критического события (такого как постановка окончательного диагноза, декомпенсации и т. д.), можно использовать метод построения таблиц жизни (анализ выживаемости, дожития). Анализ выживаемости [1, 2] позволяет учитывать в исследовании данные, касающиеся тех объектов, которые не участвуют в конечном событии, так
зе
декабрь №12 (201)
как они тоже могут внести свой вклад в расчет вероятностей наступления конечного события. Случаи могут входить в исследование в различное время. Техники анализа выживаемости уникальны в том, что они используют информацию о тех случаях, когда терминальное событие не наступило в течение времени наблюдения. Например, за время наблюдения не у всех пациентов наступила ремиссия. К этим наблюдениям применяется термин цензуриро-ванные случаи.
Таблицы жизни используются для анализа выживаемости при возможности классифицирования (разделения) наблюдения на значимые равные интервалы времени, такие как секунды, месяцы, годы и т. д. Можно применять ТЖ для расчета возможности наступления конечного (критического) события во время любого временного интервала на протяжении исследования. Техника ТЖ определяет вероятность конечного события для каждого интервала путем деления числа случаев, испытавших это событие на протяжении интервала, на число случаев, не испытавших конечное событие в это же время, после того как знаменатель согласован с цензурированными наблюдениями.
Кумулятивная выживаемость (КВ) — то есть вероятностная выживаемость к концу любого определяемого (известного) временного интервала.
Расчет
кумулятивной выживаемости по ТЖ
Для расчета КВ в ТЖ приводятся следующие данные.
1. Время начала (старта) интервала — старт каждого интервала в годах, месяцах и т. д. Временем начала наблюдения служит время начала контакта с неблагоприятным производственным фактором.
2. Число вошедших в этот интервал — это число, не имеющих окончательного диагноза профессионального заболевания на начало данного временного отрезка.
3. Число изъятых в течение интервала — это число цензурированных случаев.
4). Число подвергнутых риску — это число случаев, вошедших в данный интервал, минус половина от числа изъятых на протяжении интервала. Используется для подсчета вклада цензурированных данных. Этот показатель применяется для вычисления пропорции выживания и пропорции окончания как знаменатель.
5. Число терминальных событий — это число случаев, которые испытали конечное собы-
тие. В данном исследовании — это число случаев подозрения на профпатологию, или число установленных окончательных диагнозов, или число умерших пациентов.
6. Пропорция окончания — это пропорция случаев, испытавших конечное событие на протяжении данного интервала. Это оценка вероятности того, что случай испытает терминальное событие в течение интервала.
7. Пропорция выживания (пропорциональная выживаемость) — это пропорция случаев, переживших данный интервал.
8. Кумулятивная пропорция выживаемости (кумулятивная выживаемость (КВ)) — это вероятность выживания к концу любого определенного интервала. КВ равна произведению пропорции выживания (пропорциональная выживаемость) текущего временного интервала на пропорцию выживания предыдущего интервала.
9. Плотность вероятности — это оценка вероятности потери на протяжении отдельного (определенного) интервала.
10. Степень риска — это пропорция тех, кто дожил до определенного интервала ко всем, кто испытает терминальное событие в этом интервале. Это — коэффициент (процент, доля, степень), и может иметь значение более 1.
11. Стандартная ошибка КВ — это оценка вариабельности кумулятивной пропорции выживаемости.
12. Стандартная ошибка плотности вероятности — это оценка вариабельности плотности вероятности.
13. Стандартная ошибка степени риска — это оценка вариабельности степени риска.
14. Срединное время выживания — это время, в течение которого половина случаев достигнет терминального события. То есть, это время, в течение которого кумулятивная пропорция выживания составляет 0,50.
Функции выживания и риска. Кумулятивная пропорция выживаемости является оценкой (расчетом) функции выживания 8 (1). Она определяется как пропорция случаев переживших дольше определенного интервала 1. Высокая функция выживания отражает высокую вероятность не иметь признаков профпатологии, или окончательного диагноза, или высокую вероятность не умереть. Функция риска — это оценка (пропорция, степень, коэффициент) условной вероятности потери, разделенной временным интервалом. Символ И (1) применяется для представления функции риска. Он индикатор вероятности, что случай подвергнет-
ДШШ №12 (201) ЗНиСО
37
ся терминальному событию в определенный временной период. Высокая функция риска свидетельствует о высокой вероятности диагностики первых признаков профессиональной патологии, вероятности установления окончательного диагноза профессионального заболевания. Она может быть получена из функции выживания 8 (1) и может принимать любые значения от 0 до бесконечности.
Приведем пример использования построения таблиц жизни для прогнозирования риска развития силикоза.
Целью данной работы была оценка динамики развития риска профессиональной патологии на протяжении пылевого стажа у лиц, подвергавшихся воздействию силикозоопасной пыли.
Проанализированы истории болезни 223 работников, имеющих контакт с силикозоопасной пылью. Средний возраст обследуемых составил 56,80 ± 0,52 лет (от 37 до 79 лет). Стаж работы в контакте с пылью кремния диоксида был 26,45 ± 0,44 года (от 8 до 42 лет).
В группе 148 человек страдали силикозом (66,4 %), у 38 пациентов (17,0 %) профессиональный диагноз только подозревался и 37 человек (16,6 %) не имели признаков пылевой патологии.
Рассчитывались следующие характеристики: вероятность не иметь профессионального заболевания и риск развития пылевой патологии легких (пакет прикладных программ SPSS, версия 11,0).
Результаты исследований
Анализ вероятности отсутствия пылевой патологии у лиц, контактирующих с пылью диоксида кремния, показал, что она прогрессивно уменьшается на протяжении контакта с вредным фактором. Так, через 15 лет работы в условиях воздействия силикозоопасной пыли кумулятивная вероятность не иметь диагноза силикоз составила 88,56 ± 0,42 %, через 25 лет — 45,73 ± 1,23 %, через 35 лет - 6,79 ± 2,84 % и к максимально наблюдаемому сроку контакта (40 лет) - 6,79 ± 2,84% (рис. 1).
Срединная длительность стажа работы в условиях воздействия пыли кремния (то есть стаж работы в упомянутых условиях труда, при котором прогнозируется развитие профессионального заболевания у половины наблюдаемых), составила 29,01 года.
Риск развития силикоза неуклонно возрастал и достиг максимального значения 17,78 ±
Рис. 1. Динамика функции вероятности не иметь диагноза силикоз у наблюдаемых пациентов
38
декабрь №12 (201)
4,60 % к 35 годам контакта с неблагоприятным производственным фактором (рис. 2).
К 40 годам стажа работы в условиях воздействия кремнийсодержащей пыли риск развития силикоза равнялся нулю, что, возможно, связано с тем, что все у кого могла развиться профессиональная патологии из наблюдаемых работников все уже получили данное заболевание.
В сводной табл. 1 приведена динамика всех анализируемых характеристик у наблюдаемых
лиц в зависимости от длительности стажа работы в условиях воздействия пыли кремния.
Таким образом, длительность стажа работы в условиях воздействия силикозоопасной пыли, при котором прогнозируется развитие силикоза у половины контактирующих, составляет 29 лет. На протяжении первых 15 лет контакта с пылью вероятность не заболеть сохраняется на уровне 85—100 % (риск развития профпатологии при этом не превышает двух процентов).
Функция риска
0,15"
п а.
¥ о
О. 0,05"
0 10
Стаж до постановки диагноза силикоз, лет
Рис. 2. Динамика риска развития силикоза у наблюдаемых пациентов
Таблица 1. Основные характеристики развития силикоза в зависимости от длительности стажа работы в условиях воздействия кремнийсодержащей пыли
Пылевой стаж, лет Кумулятивная вероятность отсутствия пылевой патологии, % Риск развития пылевой патологии, %
5 100,00 ± 0,00 0,00 ± 0,00
10 96,84 ± 1,18 0,64 ± 0,24
15 88,56 ± 2,15 1,79 ± 0,42
20 67,28 ± 3,27 5,46 ± 0,82
25 45,73 ± 3,67 7,63 ± 1,23
30 22,87 ± 3,43 13,33 ± 2,26
35 8,79 ± 2,84 17,78 ± 4,60
40 8,79 ± 2,84 00,00 ± 0,00
декабрь №12 (201) зНифо
39
Наибольшее число случаев появления профессиональной патологии органов дыхания от воздействия пыли кремния приходится на период 20—25-летнего стажа в условиях воздействия упомянутого фактора (80 случаев постановки диагноза — 35,9 % ). Максимальный риск развития профессиональных изменений в состоянии здоровья у лиц, контактирующих с неблагоприятным производственным фактором, даже на протяжении 35-летнего контакта является незначительным и составляет менее 20,0 %.
Выводы
1. Метод построения таблиц жизни может быть использован для анализа риска развития профессиональной патологии, а также прогнозирования сроков развития и течения профессиональных заболеваний.
2. Максимальная длительность пылевого стажа, при которой вероятность не иметь профессиональной патологии сохраняется более 95,0 %, составляет 10 лет. Срединное время дли-
тельности вредного стажа, при котором прогнозируется 50,0 % вероятность развития силикоза легких, составило 29 лет.
3. Наибольшее число случаев появления профессиональной патологии органов дыхания от воздействия пыли кремния приходится на период 25-летнего вредного стажа. Риск развития профессиональных изменений в состоянии здоровья у лиц, контактирующих с неблагоприятным производственным фактором, даже на протяжении 40-летнего контакта наблюдается менее 20,0 %.
СПИСОК ЛИТЕРАТУРЫ
1. Buhl A., Zofel P. SPPS: Methoden fur die Markt-und Meinungsforschung. München, 2000. P. 580
2. Бююль А. SPPS: искусство обработки информации. Анализ статистических данных и восстановление скрытых закономерностей. Пер. с нем. [Текст] / А. Бююль, П. Це-фель //СПб.: ООО «ДиасофтЮП»б 2002.608 с.
