3.3. Оценка избыточного относительного риска смертности (err)
При оценке радиационных рисков лиц, подвергшихся облучению, НКДАР ООН рекомендует использовать модели UNSCEAR-94.
Пусть ожидаемое количество смертей есть m0, а наблюдаемое количество смертей – m. Тогда наблюдаемое количество смертей можно выразить через ожидаемое в аддитивной форме:
m = m0 + EAR, (1)
где EAR – превышение количества наблюдаемых случаев смертей над ожидаемым, вызванное облучением. Взаимосвязь m и m0 может быть выражена и через относительное превышение числа смертей ERR или в мультипликативной форме:
m = m0×(1 + ERR). (2)
В общем случае избыточные риски EAR и ERR зависят от дозы облучения, возраста при облучении, текущего возраста и пола.
Из (1) и (2) следует, что:
EAR = m - m0, (3)
ERR = (m - m0)/m0. (4)
Из (3) и (4) следует очевидная связь между величинами, характеризующими превышение естественного уровня смертности:
EAR = m0×ERR. (5)
Согласно общепринятой терминологии EAR – избыточный абсолютный риск (Excess Absolute Risk), ERR – избыточный относительный риск (Excess Relative Risk).
В работе оценивался избыточный относительный риск (excess relative risk, ERR), равный оценке параметра в функции (6) для дозы:
f(D)=0 (1+ D). (6)
Подбор параметра осуществлялся с помощью максимизации функции правдоподобия.
Оценивалась значимость зависимости риска от дозы. Статистическим тестом для этих целей служит тест отношения правдоподобий.
Для оценки заболеваемости использовалась мультипликативная модель (ERR). Модель для оценки ERR имеет следующий вид:
, (7)
где
– фоновая заболеваемость (заболеваемость
людей, если бы они не подверглись
облучению), зависящая от календарного
года c,
региона r
и достигнутого возраста a;
– функция, характеризующая относительные
изменения в заболеваемости, связанные
с дозой внешнего облучения d.
Аналогично исследованию стандартизованных показателей смертности, анализ зависимости смертности от дозы внешнего облучения проводился для трех когорт:
Ликвидаторы из всех регионов (полная когорта). Использовались данные только о тех ликвидаторах, на которых имеется информация о дозе внешнего облучения. Для класса болезней системы кровообращения (МКБ10: I00-I99), средняя доза внешнего облучения для данной когорты (85608 человека) равняется 0,113 Гр, максимальная доза – 1,41 Гр и минимальная – 0,0001 Гр. Общее количество случаев смертности в период с 1986 по 2012 гг. составило 10432. Для класса цереброваскулярных болезней (МКБ10: I60-I69), средняя доза внешнего облучения для данной когорты (85608 человека) равняется 0,113 Гр, максимальная доза – 1,41 Гр и минимальная – 0,0001 Гр. Общее количество случаев смертности в период с 1986 по 2012 гг. составило 1859.
Ликвидаторы из всех регионов, въехавшие в зону аварии в первые 365 дней (до 25.04.1987). Для класса болезней системы кровообращения (МКБ10: I00-I99), средняя доза внешнего облучения для данной когорты (42563 человек) равняется 0,164 Гр, максимальная доза – 1,4 Гр и минимальная – 0,0001 Гр. Общее количество случаев смертности в период с 1986 по 2012 гг. составило 5398. Для класса цереброваскулярных болезней (МКБ10: I60-I69), средняя доза внешнего облучения для данной когорты (42563 человека) равняется 0,164 Гр, максимальная доза – 1,4 Гр и минимальная – 0,0001 Гр. Общее количество случаев смертности в период с 1986 по 2012 гг. составило 949.
Ликвидаторы из всех регионов, въехавшие в зону аварии в первые 480 дней (до 31.01.1988). Для класса болезней системы кровообращения (МКБ10: I00-I99), средняя доза внешнего облучения для данной когорты (69242 человек) равняется 0,132 Гр, максимальная доза – 1,4 Гр и минимальная – 0,0001 Гр. Общее количество случаев смертности в период с 1986 по 2012 гг. составило 8697. Для класса цереброваскулярных болезней (МКБ10: I60-I69), средняя доза внешнего облучения для данной когорты (69242 человека) равняется 0,132 Гр, максимальная доза – 1,4 Гр и минимальная – 0,0001 Гр. Общее количество случаев смертности в период с 1986 по 2012 гг. составило 1565.
В таблицах 12 - 14 приводятся оценки величины ERR с 95% доверительным интервалом для различных групп локализаций смертности по рассматриваемым когортам, в таблицах также представлено количество случаев (исключая людей, имеющих онкологическое заболевание на всем периоде исследования) и величина p-value характеризующая статистическую значимость исследуемого параметра (ERR).
Таблица 12 - Значения ERR смертности для полной когорты
Класс заболеваний |
Случаев |
ERR, 1/Гр |
95% ДИ |
p-value |
Болезни системы кровообращения |
9453 |
0.4884 |
(0.2180; 0.7588 |
0.001 |
Цереброваскулярные болезни |
1664 |
0.6196 |
(-0.04278; 1.282) |
0.06675 |
Таблица 13 - Значения ERR смертности для первой референтной когорты
Класс заболеваний |
Случаев |
ERR, 1/Гр |
95% ДИ |
p-value |
Болезни системы кровообращения |
4832 |
0.6836 |
(0.2671; 1.100) |
0.001295 |
Цереброваскулярные болезни |
832 |
1.653 |
(0.37; 2.94) |
0.01172 |
Таблица 14 - Значения ERR смертности для второй референтной когорты
Класс заболеваний |
Случаев |
ERR, 1/Гр |
95% ДИ |
p-value |
Болезни системы кровообращения |
7864 |
0.3686 |
(0.06761; 0.6697) |
0.01639 |
Цереброваскулярные болезни |
1394 |
0.5584 |
(-0.1891; 1.306) |
0.1432 |