Медицинские радиологические последствия среди участниц ликвидации аварии на Чернобыльской АЭС. Часть 1. Прогнозные оценки

В настоящее время, когда прошло более 25 лет после аварии на Чернобыльской АЭС, по-прежнему одним из наиболее актуальных остаётся вопрос о её медицинских радиологических последствиях. Более 25 лет функционирует Национальный радиационно-эпидемиологический регистр (НРЭР), созданный по постановлению правительства РФ на базе Федерального государственного бюджетного учреждения «Медицинский радиологический научный центр» Минздравсоцразвития России.

Представленная работа посвящена оценке медицинских радиологических последствий среди участниц ликвидации аварии на Чернобыльской АЭС. В настоящее время НРЭР располагает медико-дозиметрической информацией на 6577 женщин-ликвидаторов [1, 2]. Впервые делается попытка крупномасштабного радиационно-эпидемиологического анализа информации НРЭР об участницах ликвидации последствий аварии на Чернобыльской АЭС. Планируемая работа включает две части. На первом этапе будет сделан прогноз радиологических последствий на основе современных моделей радиационного риска НКДАР ООН [4] и МКРЗ [3]. На втором этапе планируется оценить радиологические последствия по фактическим данным НРЭР и дать сравнение этой оценки с выводами прогноза.

Методы расчёта

В разделе строятся и сравниваются распределения пожизненных атрибутивных рисков у женщин-ликвидаторов по заболеваемости и смертности от солидных раков, а также отдельно

Иванов В . К . – Председатель РНКРЗ, зам. директора по научн. работе, член-корр. РАМН; Меняйло А . Н . * – аспирант; Кащеев В . В . – ст. научн. сотр., к.б.н.; Чекин С . Ю . – ст. научн. сотр.; Максютов М . А . – зав. лаб., к.т.н.; Пряхин Е . А . – программист; Туманов К . А . – научн. сотр. ФГБУ МРНЦ Минздравсоцразвития России.

по локализациям: лёгкие, желудок и молочная железа. Для расчетов использовалась база данных НРЭР по ликвидаторам Чернобыльской АЭС. В ней из 6577 женщин-ликвидаторов, у 1842 имеется информация о полученных дозах. Здесь использовались только данные о женщинах с известными дозами. Чтобы подготовить эти данные к анализу, а затем обработать их для получения искомого результата, было разработано несколько программных средств.

В работе рассматриваются две когорты женщин: когорта А – все 1842 женщины-ликвидаторы, для которых имеются сведения о полученных дозах, и когорта Б – только те 863 женщины из когорты А, которые приехали в зону аварии в 1986 г. Средний возраст женщин на момент прибытия в зону аварии в когорте А составляет 37,2 года, а в когорте Б – 37,5 лет.

Под пожизненным атрибутивным риском заболеваемости/смертности будем понимать отношение пожизненного избыточного абсолютного риска LAR заболеваемости/смертности к общему суммарному наблюдаемому показателю заболеваемости/смертности, т.е.: _______________ LAR inc ( l,S,D g a max

AR inc ( l,s,D g )

LARn(l,s,Dg)+ ^ S c(gf,a,l,s)■ A(a,l,s)

a = gf

AR mort ( l,s,D g ) =

_______________LARmor (l,S,Dg)______________

LARmor(l,s,Dg)+ am^ Smor(gf,a,s)■ ^mort(a,l,s) a=gf

Здесь D g – доза, полученная в возрасте g , g f – возраст первого облучения, a – возраст дожития, l – локализация опухоли, s – пол, S^inc – функция здорового дожития от возраста g f до возраста a , S " ^orr - функция дожития от возраста g f до возраста a , ^ ^c и A ^ort - фоновые показатели заболеваемости и смертности соответственно.

В данной работе указанный пожизненный атрибутивный риск рассчитывался для каждой женщины-ликвидатора Чернобыльской АЭС, для которых имеются сведения о полученных дозах. Таким образом, для каждой женщины-ликвидатора из когорт А и Б имеется набор данных: общая, полученная за время пребывания в зоне аварии, доза; пожизненный атрибутивный риск заболеваемости и пожизненный атрибутивный риск смертности. Используя эти данные в совокупности можно построить эмпирические плотности распределения дозы и пожизненных атрибутивных рисков для когорт А и Б, а также оценить для этих распределений основные статистические характеристики: среднее, медиана, мода, среднеквадратичное отклонение, интерквартильный размах.

Пусть X – означает исследуемую величину (доза или пожизненный атрибутивный риск), тогда xi – исследуемая величина, соответствующая i-й женщине из когорты А или Б. Пусть так- же всего женщин в исследуемой когорте будет N. Тогда среднее значение этой величины X вычисляется как среднее арифметическое:

N

X = -1 Ух. .

Nh '

Медиана Me(X) – это такое значение исследуемой величины X , когда 50 % женщин из рассматриваемой когорты имеют значение величины x_i < Me(X) .

Для того, чтобы построить эмпирическую плотность распределения величины X, весь набор значений x, в когортах А и Б разбивается на одинаковые интервалы шириной Ах. Затем вы- nj

.

числяется эмпирическая плотность вероятности для каждого интервала, т.е.: f =----—

j N • Ax

Здесь n j – количество женщин, у которых значение величины x i попадает в j -й интервал разбиения. Тогда можно определить наиболее вероятный интервал, т.е. интервал с максимальной f j . Модой будем считать середину этого интервала. Эти три параметра: среднее значение, медиана и мода являются характеристиками среднего для распределения величины X . Рассмотрим также для этого распределения и характеристики разброса.

Среднеквадратичное отклонение S показывает среднее отклонение от величины X и определяется как:

s =

i = 1

N -1

Интерквартильный размах E определяется как разность:

E = x о ,75 - x о , 25 , где x₀ ,75 - это такое значение исследуемой величины X , когда 75 % женщин из рассматриваемой когорты имеют значение величины x i < x 0,75 . И аналогично, x 0,25 – это такое значение исследуемой величины X , когда 25 % женщин из рассматриваемой когорты имеют значение величины x i < x 0,25 .

Дозовое распределение

При использовании базы данных НРЭР и разработанных программных средств, были получены эмпирические плотности дозового распределения среди женщин-ликвидаторов в когортах А и Б. Графики этих плотностей изображены на рисунках 1а и 1б. Здесь используется логарифмический масштаб по оси доз.

Для того, чтобы построить эмпирическую плотность дозового распределения, весь набор доз в когортах А и Б разбивался на одинаковые интервалы шириной в 1 мГр.

В когорте А наиболее вероятной оказалась доза в интервале от 1 до 2 мГр (средняя 1,5 мГр). Число женщин в этом интервале доз составляет 403, т.е. почти 21,9 % от общего числа женщин в когорте А. У 9 женщин из когорты А доза меньше наиболее вероятной (0-1 мГр), у остальных доза выше (вплоть до 830 мГр), поэтому средняя доза в когорте составляет 45,7 мГр. В то же время только 28,3 % женщин из когорты А имеют дозу выше средней, а 50 % – выше 7,4 мГр. Среднеквадратичное отклонение в когорте А составляет 78,4 мГр (в 1,7 раза выше среднего), а интерквартильный размах – 48,7 мГр.

• а)    Когорта А

  M=45.7 мГр; Me=7.4 мГр; Mo=1.5 мГр; S=78.4 мГр

  
  
  

  M=41.7 мГр; Me=7.5 мГр; Mo=1.5 мГр; S=65.1 мГр

  
  
  

  M=46.4 мГр; Me=7.0 мГр; Mo=1.5 мГр; S=79.9 мГр

  
  
  

  M=47.6 мГр; Me=8.0 мГр; Mo=1.5 мГр; S=85.7 мГр

  

  Рис . 1 а . Дозовое распределение среди женщин-ликвидаторов последствий аварии на Чернобыльской АЭС (когорта А).

  
  

В когорте Б, также как и в когорте А, наиболее вероятная доза принадлежит интервалу от 1 до 2 мГр (средняя 1,5 мГр). Женщин в этом интервале доз 111, т.е. 12,9 % от их общего числа. Видно, что в когорте Б процент людей с более высокими дозами больше, чем в когорте А, поэтому в когорте Б по сравнению с когортой А увеличивается средняя полученная доза и составляет 64,6 мГр. У 3 женщин в когорте Б доза ниже наиболее вероятной (0-1 мГр). Только 31,3 % женщин из когорты Б имеют дозу выше средней, а 50 % – выше 18 мГр. Среднеквадратичное отклонение в когорте А составляет 91,9 мГр (в 1,4 раза выше среднего), а интерквартильный размах – 96,9 мГр.

• б)    Когорта Б

  M=64.6 мГр; Me=18.0 мГр; Mo=1.5 мГр; S=91.9 мГр

  
  
  

  M=58.1 мГр; Me=22 мГр; Mo=1.5 мГр; S=72.7 мГр

  
  
  

  M=64.4 мГр; Me=17.4 мГр; Mo=1.5 мГр; S=95.0 мГр

  

  Рис . 1 б . Дозовое распределение среди женщин-ликвидаторов последствий аварии на Чернобыльской АЭС (когорта Б).

  
  

  M=69.3 мГр; Me=16.0 мГр; Mo=1.5 мГр; S=100.9 мГр

  
  
  

Для более детального понимания дозового распределения среди женщин-ликвидаторов на рисунках 1а и 1б также приведены для этих когорт графики эмпирической плотности распределения доз для отдельных возрастных групп. В когорте А к возрастной группе 18-29 лет принадлежат 476 женщин-ликвидаторов (25,8 % от общей когорты А), к возрастной группе 30-44 года – 834 женщины-ликвидатора (45,2 % от общей когорты), к возрастной группе 45-60 лет – 535 женщин-ликвидаторов (29 % от общей когорты). В когорте Б к возрастной группе 18-29 лет принадлежат 221 женщина-ликвидатор (25,8 % от общей когорты А), к возрастной группе 30-44 года – 378 женщин-ликвидаторов (45,2 % от общей когорты), к возрастной группе 45-60 лет – 263 женщины-ликвидатора (29 % от общей когорты). Данные в разных возрастных подкогортах могут незначительно пересекаться, так как один и тот же человек мог быть в Чернобыльской зоне в разных возрастах.

Характеристики дозового распределения для обеих когорт представлены в таблице 1.

Таблица 1

Характеристики дозового распределения у женщин - ликвидаторов в когортах А и Б для различных возрастных групп

Характеристики распределения	Когорта А				Когорта Б
	все возрасты	18-29 лет	30-44 года	45-60 лет	все возрасты	18-29 лет	30-44 года	45-60 лет
Среднее (мГр)	45,7	41,7	46,4	47,6	64,6	58,1	64,4	69,3
Медиана (мГр)	7,4	7,5	7,0	8,0	18,0	22	17,4	16,0
Мода (мГр) Среднеквадратичное	1,5	1,5	1,5	1,5	1,5	1,5	1,5	1,5
отклонение (мГр) Интерквартильный	78,4	65,1	79,9	85,7	91,9	72,7	95,0	100,9
размах (мГр) Минимальное	48,7	48,0	48,9	49,0	96,9	85,9	94,2	97,9
значение (мГр) Максимальное	0,1	0,1	0,1	1,0	0,1	0,1	0,1	1,0
значение (мГр)	830	250	800	830	830	250	800	830

Распределение пожизненных атрибутивных рисков

Основной задачей данного исследования является прогнозирование отдалённых последствий облучения женщин-ликвидаторов аварии на ЧАЭС. Для этого существуют международные модели МКРЗ и НКДАР ООН. При использовании этих моделей были рассчитаны величины пожизненных атрибутивных рисков заболеваемости и смертности от всех солидных раков и опухолей некоторых других локализаций для когорт А и Б женщин-ликвидаторов. Затем были построены эмпирические плотности распределения этих пожизненных атрибутивных рисков. На основе полученных распределений можно судить о величинах отдалённых последствий облучения при аварии на ЧАЭС.

Для того, чтобы построить эмпирическую плотность распределения пожизненных атрибутивных рисков, весь набор рисков разбивался на одинаковые интервалы. Для когорты А этот интервал составляет 0,005 %, для когорты Б – 0,007 %. Полученные распределения пожизненных атрибутивных рисков приведены на рисунках 2-5. Здесь представлены атрибутивные риски как смертности, так и заболеваемости для когорт А и Б по всем солидным типам рака и отдельно по трём локализациям опухоли: лёгкие, желудок и молочная железа.

Рассмотрим заболеваемость и смертность от рака отдельно. Сначала сравним между собой оценки пожизненных атрибутивных рисков, рассчитанные по модели МКРЗ и НКДАР ООН для когорт А и Б. Затем сравним между собой оценки, полученные для когорт А и Б, чтобы определить величину различий пожизненных атрибутивных рисков раковых заболеваний в этих когортах.

Заболеваемость а) Когорта А

МКРЗ: M=1.1495 %; Me=0.1872 %; Mo=0.0275 %; S=2.0172 %

НКДАР: M=1.1948 %; Me=0.1864 %; Mo=0.0225 %; S=2.1723 %

МКРЗ: M=1.6472 %; Me=0.3209 %; Mo=0.0375 %; S=2.5278 %

НКДАР: M=1.8994 %; Me=0.3752 %; Mo=0.0425 %; S=2.9213 %

МКРЗ: M=1.1999 %; Me=0.1941 %; Mo=0.0275 %; S=2.0398 %

НКДАР: M=1.2111 %; Me=0.1908 %; Mo=0.0225 %; S=2.0810 %

Рис . 2 а . Распределение пожизненных атрибутивных рисков заболеваемости от солидных типов рака среди женщин-ликвидаторов последствий аварии на Чернобыльской АЭС для различных возрастных групп (когорта А).

МКРЗ: M=0.6451 %; Me=0.1063 %; Mo=0.0175 %; S=1.1787 %

НКДАР: M=0.5605 %; Me=0.0905 %; Mo=0.0125 %; S=1.0425 %

б) Когорта Б

МКРЗ: M=1.5887 %; Me=0.4494 %; Mo=0.0175 %; S=2.3277 %

НКДАР: M=1.6457 %; Me=0.4567 %; Mo=0.0245 %; S=2.5038 %

МКРЗ: M=2.2919 %; Me=0.9442 %; Mo=0.0385 %; S=2.8135 %

НКДАР: M=2.6432 %; Me=1.1073 %; Mo=0.0525 %; S=3.2515 %

МКРЗ: M=1.6421 %; Me=0.4233 %; Mo=0.0245 %; S=2.3872 %

НКДАР: M=1.6543 %; Me=0.4431 %; Mo=0.0245 %; S=2.4329 %

Рис . 2 б . Распределение пожизненных атрибутивных рисков заболеваемости от солидных типов рака среди женщин-ликвидаторов последствий аварии на Чернобыльской АЭС для различных возрастных групп (когорта Б).

МКРЗ: M=0.9573 %; Me=0.2225 %; Mo=0.0175 %; S=1.4699 %

НКДАР: M=0.8352 %; Me=0.1866 %; Mo=0.0175 %; S=1.3118 %

Смертность а) Когорта А

МКРЗ: M=1.5397 %; Me=0.2519 %; Mo=0.0325 %; S=2.6799 %

НКДАР: M=1.2158 %; Me=0.1915 %; Mo=0.0275 %; S=2.1716 %

МКРЗ: M=1.6091 %; Me=0.2647 %; Mo=0.0325 %; S=2.7052 %

МКРЗ: M=2.2112 %; Me=0.4392 %; Mo=0.0525 %; S=3.3657 %

НКДАР: M=1.8576 %; Me=0.3683 %; Mo=0.0425 %; S=2.8545 %

НКДАР: M=1.2396 %; Me=0.2004 %; Mo=0.0275 %; S=2.1159 %

МКРЗ: M=0.8578 %; Me=0.1423 %; Mo=0.0225 %; S=1.5552 %

НКДАР: M=0.6255 %; Me=0.1021 %; Mo=0.0125 %; S=1.1503 %

Рис . 3 а . Распределение пожизненных атрибутивных рисков смертности от солидных типов рака среди женщин-ликвидаторов последствий аварии на Чернобыльской АЭС для различных возрастных групп (когорта А).

б) Когорта Б

МКРЗ: M=3.0742 %; Me=1.2890 %; Mo=0.0525 %; S=3.7399 %

МКРЗ: M=2.1251 %; Me=0.6149 %; Mo=0.0385 %; S=3.0869 %

НКДАР: M=1.6776 %; Me=0.4844 %; Mo=0.0175 %; S=2.5035 %

НКДАР: M=2.5849 %; Me=1.0787 %; Mo=0.0525 %; S=3.1768 %

МКРЗ: M=2.1984 %; Me=0.5775 %; Mo=0.0385 %; S=3.1589 %

НКДАР: M=1.6949 %; Me=0.4458 %; Mo=0.0245 %; S=2.4746 %

Рис . 3 б . Распределение пожизненных атрибутивных рисков смертности от солидных типов рака среди женщин-ликвидаторов последствий аварии на Чернобыльской АЭС для различных возрастных групп (когорта Б).

МКРЗ: M=1.2719 %; Me=0.2995 %; Mo=0.0245 %; S=1.9374 %

НКДАР: M=0.9297 %; Me=0.2159 %; Mo=0.0175 %; S=1.4396 %

Заболеваемость

Лёгкие:

МКРЗ: M=3.5734 %; Me=0.6651 %; Mo=0.0975 %; S=5.5906 %

НКДАР: M=3.9667 %; Me=0.7461 %; Mo=0.1075 %; S=6.1573 %

Желудок:

МКРЗ: M=1.3259 %; Me=0.2211 %; Mo=0.0325 %; S=2.3043 %

НКДАР: M=1.1609 %; Me=0.1975 %; Mo=0.0275 %; S=1.9254 %

• а)    Когорта А

• Молочная железа:

МКРЗ: M=1.9767 %; Me=0.1517 %; Mo=0.0125 %; S=1.8243 %

НКДАР: M=2.4260 %; Me=0.4247 %; Mo=0.0725 %; S=3.9800 %

Средний в интервале атрибутивный риск заболеваемости раком молочной железы (%)

Лёгкие:

б) Когорта Б

Желудок:

МКРЗ: M=4.9692 %; Me=1.6096 %; Mo=0.0945 %; S=6.4569 %

НКДАР: M=5.5131 %; Me=1.8230 %; Mo=0.1085 %; S=7.1122 %

МКРЗ: M=1.8333 %; Me=0.5311 %; Mo=0.0315 %; S=2.6576 %

НКДАР: M=1.6251 %; Me=0.4694 %; Mo=0.0245 %; S=2.2463 %

Молочная железа:

МКРЗ: M=1.3451 %; Me=0.3645 %; Mo=0.0105 %; S=2.1046 %

НКДАР: M=3.3560 %; Me=1.0393 %; Mo=0.0435 %; S=4.5886 %

Лёгкие:

Смертность а) Когорта А

Желудок:

МКРЗ: M=4.0149 %; Me=0.7596 %; Mo=0.1075 %; S=6.2232 %

НКДАР: M=2.2335 %; Me=0.3956 %; Mo=0.0625 %; S=3.6661 %

МКРЗ: M=1.4045 %; Me=0.2349 %; Mo=0.0325 %; S=2.4361 %

НКДАР: M=0.8897 %; Me=0.1412 %; Mo=0.0175 %; S=1.5899 %

Молочная железа:

МКРЗ: M=1.5564 %; Me=0.2464 %; Mo=0.0175 %; S=2.8587 %

НКДАР: M=1.9702 %; Me=0.3160 %; Mo=0.0475 %; S=3.5195 %

• б)    Когорта Б

  Лёгкие:

  
  

  МКРЗ: M=5.5750 %; Me=1.8364 %; Mo=0.1085 %; S=7.1747 %

  НКДАР: M=3.1465 %; Me=0.9423 %; Mo=0.0595 %; S=4.3051 %

  
  
  
0.001      0.01       0.1        1         10

Желудок:

МКРЗ: M=1.9415 %; Me=0.5641 %; Mo=0.0315 %; S=2.8085 %

НКДАР: M=1.2303 %; Me=0.3539 %; Mo=0.0175 %; S=1.8360 %

Все возрасты

m

Молочная железа:

МКРЗ: M=2.1398 %; Me=0. 5915 %; Mo=0.0175 %; S=3.2899 %

НКДАР: M=2.7052 %; Me=0. 7677 %; Mo=0.0455 %; S=4.0435 %

Заболеваемость

Когорта А . Для начала рассмотрим в общем все солидные типы рака в когорте А. Как видно из рисунка 2а, пожизненный атрибутивный риск заболеваемости солидными типами рака, рассчитанный по модели МКРЗ, практически не отличается от аналогичного риска, рассчитанного по модели НКДАР. Средний пожизненный атрибутивный риск заболеваемости для всех возрастов по модели МКРЗ равен 1,1495 %, при этом у 26,5 % женщин риск выше среднего. Для модели НКДАР аналогичный средний пожизненный атрибутивный риск равен 1,1948 % и у 25,8 % этот риск выше среднего. Медиана распределения пожизненного атрибутивного риска по модели МКРЗ равна 0,1872 %, а по модели НКДАР – 0,1864 %. Таким образом, по среднему значению риска модель НКДАР даёт оценку выше, чем модель МКРЗ всего в 1,04 раза, а по медиане распределений оценки практически не отличаются. Наиболее вероятным для обеих моделей оказался пожизненный атрибутивный риск заболеваемости солидными типами рака в интервале от 0,02 % до 0,025 %. Число женщин в этом интервале рисков для модели МКРЗ составляет 86 человек (4,6 % от всей когорты), а для модели НКДАР – 79 человек (4,2 % от всей когорты). Число интервалов пожизненного атрибутивного риска, находящихся левее наиболее вероятного, равно 5 для модели МКРЗ и 4 для модели НКДАР. Количество женщин, у которых пожизненные атрибутивные риски соответствуют этим интервалам, равно 223 (12,1 % от всей когорты) по модели МКРЗ и 175 (9,5 % от всей когорты) по модели НКДАР. Для модели МКРЗ среднеквадратичное отклонение составляет 2,0172 % (в 1,75 раза выше среднего), а интерквартильный размах – 1,2731 %. Соответствующие цифры для модели НКДАР равны 2,1723 % (в 1,81 раза выше среднего) и 1,1895 %. Как и в случае со средними величинами, величины, характеризующие разброс пожизненного атрибутивного риска, также незначительно различаются между собой для обеих моделей (по среднеквадратичному отклонению модель НКДАР даёт оценку в 1,08 раза выше, чем модель МКРЗ, а по интерквартильному размаху – в 1,07 раза ниже, чем МКРЗ).

Из вышесказанного можно сделать вывод, что в когорте А пожизненный атрибутивный риск заболеваемости солидными типами рака незначительно различается при расчёте по двум разным моделям МКРЗ и НКДАР в величинах среднего и величинах разброса атрибутивного риска. Для более детального понимания распределения пожизненного атрибутивного риска заболеваемости солидными типами рака на рисунке 2а приведены также графики эмпирической плотности распределения для различных возрастных групп в исследуемой когорте. Характеристики этих возрастных групп такие же, какие использовались при анализе дозового распределения.

В таблице 2 приведены основные характеристики распределения пожизненных атрибутивных рисков заболеваемости солидными раками для этих возрастных групп.

Для более детального понимания распределения рисков среди женщин-ликвидаторов кратко рассмотрим также распределение пожизненных атрибутивных рисков заболеваемости в когорте А для нескольких отдельных локализаций опухоли: лёгкие, желудок, молочная железа. Графики распределений атрибутивных рисков заболеваемости раком этих локализаций приведены на рисунке 4а, а в таблице 3 приведены основные параметры этих распределений для всех возрастов в исследуемой когорте.

Таблица 2

Характеристики распределения пожизненных атрибутивных рисков заболеваемости и смертности солидными типами рака у женщин - ликвидаторов последствий аварии на Чернобыльской АЭС для различных возрастных групп в когортах А и Б

Заболеваемость
Характеристики распределения	Модель МКРЗ								Модель НКДАР ООН
	Все возрасты		18-29 лет		30-44 года		45-60 лет		Все возрасты		18-29 лет		30-44 года		45-60 лет
	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б
Среднее	1,1495	1,5887	1,6472	2,2919	1,1999	1,6421	0,6451	0,9573	1,1948	1,6457	1,8994	2,6432	1,2111	1,6543	0,5605	0,8352
Медиана	0,1872	0,4494	0,3209	0,9442	0,1941	0,4233	0,1063	0,2225	0,1864	0,4567	0,3752	1,1073	0,1908	0,4431	0,0905	0,1866
Мода	0,0275	0,0175	0,0375	0,0385	0,0275	0,0245	0,0175	0,0175	0,0225	0,0245	0,0425	0,0525	0,0225	0,0245	0,0125	0,0175
Среднеквадратичное отклонение	2,0172	2,3277	2,5278	2,8135	2,0398	2,3872	1,1787	1,4699	2,1723	2,5038	2,9213	3,2515	2,0810	2,4329	1,0425	1,3118
Интерквартильный размах	1,2731	2,1252	2,0543	3,2973	1,3368	2,2739	0,7437	1,3531	1,1895	2,1607	2,3843	3,7404	1,3791	2,2310	0,6480	1,1686
Минимальное значение	0,00223	0,0032	0,0039	0,0043	0,0032	0,0032	0,0073	0,0094	0,0014	0,0034	0,0045	0,0051	0,0034	0,0034	0,0056	0,0075
Максимальное значение	16,6522	16,6522	10,3128	10,3128	16,6522	16,6522	12,8774	12,8774	17,5526	17,5526	12,1739	12,1739	17,5526	17,5526	11,7911	11,7911
Смертность
Характеристики распределения	Модель МКРЗ								Модель НКДАР ООН
	Все возрасты		18-29 лет		30-44 года		45-60 лет		Все возрасты		18-29 лет		30-44 года		45-60 лет
	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б
Среднее	1,5397	2,1251	2,2112	3,0742	1,6091	2,1984	0,8578	1,2719	1,2158	1,6776	1,8575	2,5849	1,2396	1,6949	0,6255	0,9297
Медиана	0,2519	0,6149	0,4392	1,2890	0,2647	0,5775	0,1423	0,2995	0,1915	0,4844	0,3683	1,0787	0,2004	0,4458	0,1021	0,2159
Мода	0,0325	0,0385	0,0525	0,0525	0,0325	0,0385	0,0225	0,0245	0,0275	0,0175	0,0425	0,0525	0,0275	0,0245	0,0125	0,0175
Среднеквадратичное отклонение	2,6799	3,0869	3,3657	3,7399	2,7052	3,1589	1,5552	1,9374	2,1716	2,5035	2,8545	3,1768	2,1159	2,4746	1,1503	1,4396
Интерквартильный размах	1,7068	2,8727	2,7909	4,4522	1,8163	3,0680	0,9984	1,8070	1,2888	2,2360	2,3233	3,7025	1,3952	2,2897	0,7302	1,3080
Минимальное значение	0,0029	0,0044	0,0054	0,0059	0,0044	0,0044	0,0096	0,0124	0,0018	0,0034	0,0044	0,0050	0,0034	0,0034	0,0068	0,0088
Максимальное значение	21,4558	21,4558	13,6117	13,6117	21,4558	21,4558	16,6505	16,6505	17,5434	17,5434	11,8397	11,8397	17,5434	17,5434	12,7324	12,7324

Радиация и риск . 2011. Том 20. № 3                                                          Научные статьи

Таблица 3

Заболеваемость
Характеристики распределения	Лёгкие				Желудок				Молочная железа
	Модель МКРЗ		Модель НКДАР		Модель МКРЗ		Модель НКДАР		Модель МКРЗ		Модель НКДАР
	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б
Среднее	3,5734	4,9692	3,9667	5,5131	1,3259	1,8333	1,1609	1,6251	0,9767	1,3451	2,4260	3,3560
Медиана	0,6651	1,6096	0,7461	1,8230	0,2211	0,5311	0,1975	0,4694	0,1517	0,3645	0,4247	1,0393
Мода	0,0975	0,0945	0,1075	0,1085	0,0325	0,0315	0,0275	0,0245	0,0125	0,0105	0,0725	0,0435
Среднеквадратичное отклонение	5,5906	6,4569	6,1573	7,1122	2,3043	2,6576	1,9254	2,2463	1,8243	2,1046	3,9800	4,5886
Интерквартильный размах	4,3619	7,5692	4,8745	8,4829	1,5053	2,4514	1,3348	2,3606	0,9559	1,7292	2,9225	4,7743
Минимальное значение	0,0095	0,0096	0,0102	0,0103	0,0029	0,0037	0,0029	0,0029	0,0013	0,0027	0,0071	0,0071
Максимальное значение	42,7853	42,7853	46,0192	46,0192	18,6451	18,6451	17,9800	17,9800	14,6312	14,6312	32,4002	32,4002
Смертность
	Лёгкие				Желудок				Молочная железа
Характеристики	Модель		Модель		Модель		Модель		Модель		Модель
распределения	МКРЗ		НКДАР		МКРЗ		НКДАР		МКРЗ		НКДАР
	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б	Ког. А	Ког. Б
Среднее	4,0149	5,5750	2,2335	3,1465	1,4045	1,9415	0,8897	1,2303	1,5564	2,1398	1,9702	2,7052
Медиана	0,7596	1,8364	0,3956	0,9423	0,2349	0,5641	0,1412	0,3539	0,2464	0,5915	0,3160	0,7677
Мода	0,1075	0,1085	0,0625	0,0595	0,0325	0,0315	0,0175	0,0175	0,0175	0,0175	0,0475	0,0455
Среднеквадратичное отклонение	6,2232	7,1747	3,6661	4,3051	2,4361	2,8085	1,5899	1,8360	2,8587	3,2899	3,5195	4,0435
Интерквартильный размах	4,9544	8,5526	2,7208	4,6270	1,5984	2,6017	0,9643	1,6234	1,5575	2,7958	2,0204	3,5775
Минимальное значение	0,0109	0,0109	0,0040	0,0042	0,0030	0,0039	0,0014	0,0024	0,0021	0,0044	0,0024	0,0059
Максимальное значение	46,0945	46,0945	33,6207	33,6207	19,5897	19,5897	13,1397	13,1397	21,9355	21,9355	26,9674	26,9674

Радиация и риск . 2011. Том 20. № 3                                                          Научные статьи

Когорта Б. Теперь рассмотрим в общем все солидные типы раков в когорте Б. Как и в когорте А, пожизненный атрибутивный риск заболеваемости солидными раками, рассчитанный по модели МКРЗ, практически не отличается от аналогичного риска, рассчитанного по модели НКДАР. Средний пожизненный атрибутивный риск заболеваемости для всех возрастов по модели МКРЗ равен 1,5887 %, при этом у 30,5 % женщин – риск выше среднего. Для модели НКДАР аналогичный средний пожизненный атрибутивный риск равен 1,6457 % и также у 30,5 % этот риск выше среднего. Медиана распределения пожизненного атрибутивного риска по модели МКРЗ равна 0,4494 %, а по модели НКДАР – 0,4567 %. Таким образом, точно так же, как и в когорте А, по среднему значению риска модель НКДАР даёт оценку выше, чем модель МКРЗ всего в 1,04 раза, а по медиане распределений оценки практически не отличаются. Наиболее вероятным для модели МКРЗ оказался пожизненный атрибутивный риск заболеваемости солидными типами рака в интервале от 0,014 % до 0,021 %, а для модели НКДАР – от 0,021 % до 0,028 %. Количество женщин в этом интервале рисков для модели МКРЗ составляет 35 человек (4,1 % от всей когорты Б), а для модели НКДАР – 38 человек (4,4 % от всей когорты Б).

Число интервалов пожизненного атрибутивного риска, находящихся левее наиболее вероятного, равно 2 для модели МКРЗ и 3 для модели НКДАР. Количество женщин, у которых пожизненные атрибутивные риски соответствуют этим интервалам, равно 12 (1,4 % от всей когорты Б) по модели МКРЗ и 50 (5,8 % от всей когорты Б) по модели НКДАР.

В когорте Б для модели МКРЗ среднеквадратичное отклонение составляет 2,3277 % (в 1,47 раза выше среднего), а интерквартильный размах – 2,1252 %. Соответствующие цифры для модели НКДАР равны 2,5038 % (в 1,52 раза выше среднего) и 2,1607 %. Как и в когорте А, величины, характеризующие разброс пожизненного атрибутивного риска, также не сильно различаются между собой для обеих моделей.

Поэтому можно сделать вывод, что в когорте Б, как и в когорте А, пожизненный атрибутивный риск заболеваемости солидными раками незначительно различается при расчёте по двум разным моделям МКРЗ и НКДАР в величинах среднего и величинах разброса этого риска.

Из рисунков 1 и 2 видно, что единственный интервал дозового распределения в когортах А и Б, находящийся левее наиболее вероятного, привёл к появлению нескольких интервалов левее наиболее вероятного в распределении пожизненных атрибутивных рисков заболеваемости солидными раками в этих когортах.

Как и в случае с когортой А, рассмотрим заболеваемость раком нескольких отдельных локализаций в когорте Б. Графики распределений пожизненного атрибутивного риска заболеваемости раком этих локализаций приведены на рисунке 4б, а в таблице 3 также приведены основные параметры этих распределений для всех возрастов в исследуемой когорте.

Из этого следует вывод, что различие в пожизненных атрибутивных рисках заболеваемости раком рассмотренных локализаций, полученные по моделям МКРЗ и НКДАР ООН, примерно одинаковы для обеих когорт.

Смертность

Когорта А . Аналогично тому, как это было при рассмотрении заболеваемости, сначала сравним основные характеристики распределений пожизненных атрибутивных рисков смертности от всех солидных типов рака по моделям МКРЗ и НКДАР в когорте А.

Средний пожизненный атрибутивный риск смертности для всех возрастов по модели МКРЗ составляет 1,5397 % и у 26,5 % женщин этот риск выше среднего. Аналогичный риск по модели НКДАР составляет 1,2158 %, и у 26,2 % женщин этот риск выше среднего. Отсюда видно, что модель МКРЗ по этому параметру превосходит модель НКДАР в 1,27 раза. Медианы распределений пожизненных атрибутивных рисков по моделям МКРЗ и НКДАР соответственно равны 0,2519 % и 0,1915 %. Здесь модель МКРЗ превосходит модель НКДАР в 1,32 раза. Наиболее вероятным оказывается атрибутивным риск смертности, лежащий в интервале 0,03-0,035 % для модели МКРЗ и 0,025-0,03 % для модели НКДАР. Число женщин, находящихся в наиболее вероятном интервале рисков, равно 61 (3,3 % от всей когорты) для модели МКРЗ и 80 (4,3 % от всей когорты) для модели НКДАР. Количество людей с пожизненным атрибутивным риском меньшим, чем наиболее вероятный, равно 178 (9,6 % от всей когорты) для модели МКРЗ и 224 (12,2 % от всей когорты) для модели НКДАР.

Для модели МКРЗ среднеквадратичное отклонение составляет 2,6799 % (в 1,74 раза выше среднего), а интерквартильный размах – 1,7068 %. Соответствующие цифры для модели НКДАР равны 2,1716 % (в 1,79 раза выше среднего) и 1,2888 %. Таким образом, и по этим параметрам модель МКРЗ превосходит модель НКДАР. По отношению к модели МКРЗ это различие в величинах разброса пожизненного атрибутивного риска смертности составляет 19 % для среднеквадратичного отклонения и 24,5 % для интерквартильного размаха.

Из вышесказанного следует, что все рассматриваемые характеристики распределения пожизненного атрибутивного риска смертности от солидных раков по модели МКРЗ несколько выше аналогичных характеристик для модели НКДАР.

Как и в случае с пожизненным атрибутивным риском заболеваемости солидными раками в таблице 2 приведены основные характеристики распределения пожизненного атрибутивного риска смертности для различных возрастных групп.

Наибольший средний атрибутивный риск смертности в когорте А соответствует раку лёгких (2,2-4,0 %).

Когорта Б . Теперь рассмотрим в общем все солидные типы рака в когорте Б. Средний пожизненный атрибутивный риск смертности для всех возрастов по модели МКРЗ равен 2,1251 %, при этом у 30,1 % женщин он выше среднего. Для модели НКДАР аналогичный средний пожизненный атрибутивный риск равен 1,6776 % и также у 30,1 % – выше среднего. Медиана распределения пожизненного атрибутивного риска смертности по модели МКРЗ равна 0,6149 %, а по модели НКДАР – 0,4844 %. Таким образом, точно так же, как и в когорте А, по среднему значению риска модель МКРЗ даёт оценку выше, чем модель НКДАР, всего в 1,27 раза. По медиане распределений относительное различие немного меньше, чем в когорте А: оценка медианы по МКРЗ превосходит такую же оценку по НКДАР в 1,27 раза. Наиболее вероятным для модели МКРЗ оказался пожизненный атрибутивный риск смертности от солидных типов рака в интервале от 0,035 % до 0,042 %, а для модели НКДАР – от 0,014 % до 0,021 %. Число женщин в этом интервале рисков для модели МКРЗ составляет 28 человек (3,2 % от всей когорты Б), а для модели НКДАР – 31 (3,6 % от всей когорты Б).

Число интервалов атрибутивного риска, находящихся левее наиболее вероятного, равно 5 для модели МКРЗ и 2 – для модели НКДАР. Количество женщин, у которых пожизненные атрибутивные риски соответствуют этим интервалам, равно 64 (7,4 % от всей когорты Б) по модели МКРЗ и 17 (2,0 % от всей когорты Б) по модели НКДАР.

В когорте Б для модели МКРЗ среднеквадратичное отклонение пожизненного атрибутивного риска смертности от солидных типов рака составляет 3,0869 % (в 1,45 раза выше среднего), а интерквартильный размах – 2,8727 %. Соответствующие цифры для модели НКДАР равны 2,5035 % (в 1,49 раза выше среднего) и 2,2360 %. Как и в когорте А, величины, характеризующие разброс пожизненного атрибутивного риска, также не сильно различаются между собой для обеих моделей.

Поэтому можно сделать вывод, что в когорте Б, так же, как и в когорте А, все характеристики распределения пожизненного атрибутивного риска смертности от солидных раков по модели МКРЗ несколько выше аналогичных характеристик для модели НКДАР.

Видно, что как и в случае заболеваемости солидными типами рака, единственный интервал дозового распределения в когортах А и Б, находящийся левее наиболее вероятного, привёл к появлению нескольких интервалов левее наиболее вероятного в распределениях пожизненного атрибутивного риска смертности от солидных раков в этих когортах.

Как и в случае с когортой А, рассмотрим смертность от рака нескольких отдельных локализаций в когорте Б. Графики распределений пожизненного атрибутивного риска смертности от рака этих локализаций приведены на рисунке 5б. В таблице 3 приведены основные параметры этих распределений для всех возрастов в исследуемой когорте.

Сравним теперь между собой оценки рисков, полученные для когорт А и Б. Из таблицы 3 видно, что относительное различие в характеристиках пожизненного атрибутивного риска смертности от рака всех рассматриваемых локализаций (все солидные типы рака в целом, лёгкие, желудок, молочная железа) практически такое же, как и относительное различие в характеристиках пожизненного атрибутивного риска заболеваемости раком этих локализаций. Т.е., как и в случае с заболеваемостью, среднее значение пожизненного атрибутивного риска смертности от рака всех рассматриваемых локализаций, полученное для когорты Б, в 1,38-1,40 раза больше, чем аналогичное среднее значение, полученное для когорты А. Медиана распределения пожизненного атрибутивного риска смертности, рассчитанная для когорты Б, практически так же, как и в случае заболеваемости, в 2,38-2,51 раза больше, чем эта же медиана, рассчитанная для когорты А. Это справедливо и для среднеквадратичного отклонения пожизненного атрибутивного риска смертности: среднеквадратичное отклонение, рассчитанное для когорты Б, как и в случае заболеваемости больше в 1,15-1,17 раза, чем аналогичная величина, рассчитанная для когорты А.

Выводы

Различие между пожизненными атрибутивными рисками, рассчитанными по модели МКРЗ и НКДАР, значительно зависит от локализации опухоли. В случае заболеваемости оно наиболее выражено для рака молочной железы, а в случае смертности – для рака лёгких. В целом для всех солидных типов рака пожизненные атрибутивные риски заболеваемости, полученные по моделям МКРЗ и НКДАР ООН, незначительно различаются друг от друга, а пожизненный атрибутивный риск смертности, полученный по модели МКРЗ, несколько выше аналогичного риска, рассчитанного по модели НКДАР.

Наибольший средний пожизненный атрибутивный риск получился для смертности от рака лёгких в когорте Б (5,6 %). В целом по всем солидным типам рака наибольший средний пожизненный атрибутивный риск соответствует смертности в когорте Б (1,7-2,1 %). Эти результаты показывают, что в отдалённой перспективе число радиационно-индуцированных опухолей у женщин-ликвидаторов последствий аварии на Чернобыльской атомной станции будет невелико. Однако следует заметить, что средние пожизненные атрибутивные риски, полученные для когорты Б, в 1,38-1,4 раза выше, чем для когорты А, поэтому женщин, которые приехали в зону аварии в 1986 г., можно считать группой риска с повышенной вероятностью получить раковое заболевание. Тем не менее, очень важно проверить на прямых статистических исследованиях правильность этих оценок.