Научная статья на тему 'К МЕТОДИКЕ ВЕРОЯТНОСТНОЙ ОЦЕНКИ СРЕДНЕЭФФЕКТИВНОГО ВРЕМЕНИ ГИБЕЛИ ЖИВОТНЫХ ОТ СРЕДНЕСМЕРТЕЛЬНОЙ ДОЗЫ'

К МЕТОДИКЕ ВЕРОЯТНОСТНОЙ ОЦЕНКИ СРЕДНЕЭФФЕКТИВНОГО ВРЕМЕНИ ГИБЕЛИ ЖИВОТНЫХ ОТ СРЕДНЕСМЕРТЕЛЬНОЙ ДОЗЫ Текст научной статьи по специальности «Математика»

CC BY
26
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Гигиена и санитария
Scopus
ВАК
CAS
RSCI
PubMed
Область наук
i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Т.В. Пастушенко, Ю.А. Пилипенко, А.А. Жуков, Ф.А. Перейма

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «К МЕТОДИКЕ ВЕРОЯТНОСТНОЙ ОЦЕНКИ СРЕДНЕЭФФЕКТИВНОГО ВРЕМЕНИ ГИБЕЛИ ЖИВОТНЫХ ОТ СРЕДНЕСМЕРТЕЛЬНОЙ ДОЗЫ»

УДК 615.9.015.38.076.9

Т. В. Пастушенко, Ю. А. Пилипенко, А. А. Жуков, Ф. А. Перейма

К МЕТОДИКЕ ВЕРОЯТНОСТНОЙ ОЦЕНКИ СРЕДНЕЭФФЕКТИВНОГО ВРЕМЕНИ ГИБЕЛИ ЖИВОТНЫХ ОТ СРЕДНЕСМЕРТЕЛЬНОЙ ДОЗЫ

Тернопольскин медицинский институт

В токсиколого-гигненичсских исследованиях для оценки кумулятивных свойств химических веществ по результатам острых опытов успешно используется среднее время гибели животных (ЕТ50), значимость которого впервые отметил Б. М. Штабский. В последнее время некоторые авторы (Г. Н. Красовский и соавт.) показали возможность прямого использования ЕТ50 в остром опыте как показателя кумулятивности при прогнозировании хронической токсичности вредных веществ, что имеет большое практическое значение для их ускоренного гигиенического нормирования. При такой высокой информативности данного показателя нет пока единого методического подхода к его определению. Для расчета ЕТ50 в основном используется классический метод пробит-анализа, который позволяет определить его лишь для одной из введенных доз вещества, близкой к LD50 (М. Л. Беленький) или несколько превышающей ее (Б. М. Штабский, 1973). Существенным недостатком расчета ЕТ50 с помощью пробит-анализа, что справедливо отмстили Г. Н. Красовский и соавт. (1982), является отсутствие учета вариаций временного параметра, обусловленных индивидуальной чувствительностью животных к изучаемому веществу. Поэтому авторами был предложен более совершенствованный графический метод определения ЕТ50, основанный на учете гибели животных в остром опыте по суткам или часам наблюдения с дальнейшим определением среднего времени гибели животных от каждой из введенных доз вещества. Для расчета ЕТ50 строится график в логарифмическом масштабе. Из точки на прямой, соответствующей LD50, опускается перпендикуляр до пересечения с осыо абсцисс и определяется ЕТ50 — среднее время гибели животных от среднесмертельной дозы вещества.

Графические методы являются приближенными и обладают меньшей точностью по сравнению с расчетными. Учитывая это, в настоящем сообщении на основании методики графического определения ЕТ50 (Г. Н. Красовский и соавт., 1982) предлагается расчетный метод вероятностной оценки среднеэффективного времени гибели животных от среднесмертельной дозы вещества.

Для расчета среднеэффективного времени гибели животных (ETfSu, ) от среднесмертельной дозы вещества (LD50) в условиях острого опыта используется линейное уравнение регрессии:

где а и Ь — коэффициенты регрессии; у — доза (Д) или логарифм дозы (^Д); х— среднее время гибели животных (Т) или логарифм среднего времени гибели животных (^Т).

Коэффициенты регрессии а и Ь определяются методом наименьших квадратов по формулам:

Ь =

n 2 ух — 2 у 2 х я 2 хг ~ (2 *)а

2у — Ы.х

а = ■

(2) (3)

1/ = а + Ьх,

где п — число испытанных эффективных (вызывающих гибель животных) доз.

Зная а и Ь, легко найти х по уравнению (1), решая его относительно х при любом заданном значении у:

Подставляя в уравнение (4) значение у, равное ЬО50 или ^ ЬО50, находим среднеэффектив-ное время гибели животных от среднесмертельной дозы вещества — ЕТ •

Для определения доверительных границ среднеэффективного времени гибели животных от среднесмертельной дозы последовательно подставляем в уравнение (4) значения у, соответствующие нижней и верхней доверительной границе среднесмертельной дозы, таким образом определяя ЕТ 5ого для нижней 11 верхней доверительной границы ЬОго. ЕТ ооа) расчитанное для нижней доверительной границы Ь-Эзо, является верхней доверительной границей среднеэффективного времени гибели животных от среднесмертельной дозы, а ЕТ аоп;» рассчитанное для верхней доверительной границы ЬО50 — нижней доверительной границей среднеэффективного времени гибели животных от среднесмертельной дозы.

Техника расчета рассматривается на

приведенном ниже примере.

Пример 1. Для более наглядного сопоставления полученных результатов определения среднего времени гибели животных от среднесмертельной дозы графическим и предлагаемым расчетным методами мы для примера расчета воспользовались данными острого опыта на белых крысах с 2, 3, 6-трихлортолуолом, которые Г. Н. Красовский и соавт. (1982) использовали для де-

монстрацни техники определения ЕТ50 графическим методом.

При определении ETfoci; за основу берется таблица учета времени и частоты гибели животных в остром опыте по суткам или часам наблюдения с расчетом времени гибели животных от каждой из введенных доз вещества, аналогичная предложенной ранее для этой цели Г. Н. Красов-скнм и соавт. (1982), ввиду чего мы считали нецелесообразным привести ее повторно в нашем примере.

Для вероятностной оценки ЕТ Ко<ц» расчетным путем на основании таблицы данных острого опыта на белых крысах с 2, 3, б-трихлортолуолом (Г. Н. Красовский и соавт., 1982) составляется рабочая таблица для определения коэффициентов регрессии (табл. 1). Дозы (Д), вызывающие гибель животных, и среднее время гибели животных (Т) от каждой из введенных доз вещества, так же как и в работе Г. II. Красовского и соавт. (1982) при расчете ЕТ50 графическим методом, логарифмировались.

LDS0 = 4800 (3902 - 5904) мг/кг», 2 х2 = 116,6572.

После заполнения рабочей табл. 1. рассчитываем по формуле (2) коэффициент регрессии b и по формуле (3) — коэффициент регрессии а.

_ 6-39.36708—21.95-10.8008 _

Ь~ 6 19,839222—116,6572 ——0.367979,

21,95 +0,367979-10,8008 а =-1-1—-= 4,320746.

Согласно формулам (1) и (4) цифровое выражение зависимости доза — время для любого среднего значения Д имеет следующий вид: у = 4,320716 — 0,367979 хил« 4,320746— у * — 0,367979

При Д = LD50 = 4800=у получаем:

4.320746— lg 4800 4.320746 — 3,6812 16 л:-" 0.367979 — 0,367979 —'.738.

Отсюда ЕТ Яо<1» равно 54,7 ч.

Далее находим доверительные границы ET^'od), поочередно подставляя в равенство (4) значения у, соответствующие нижней и верхней доверительной границе LD50, при этом обозначив х'— верхнюю, а х"— нижнюю доверительные границы среднеэффективного времени гибели животных.

, 4,320746—lg 3902 4,320746 — 3.5913

Ig* ~ 0,367979 — 0,367979 ~ 1,982;

„ 4,320746—lg 5904 4,320746 — 3.7712

lg* _ 0,367979 — 0,367979 —'.493.

1 Доверительные границы ЬЭи для 2,3,6-трихлортолуо-ла, рассчитаны нами, так как в работе Г. Н. Красовского и соавт. (1982) они не приводятся.

Таблица 1

Расчет коэффициентов регрессии а и Ь для вероятностного определения ЕТ;?д,, .2,3,6-трихлортолуола

D, иг/кг Т, сут lgD IgT XII *•

У X

6000 28 3,7782 1,4472 5,467811 2,094387

5500 36 3,7-104 1,5563 5,821184 2,422069

5000 48 3.6990 1.6812 6,218758 2,826433

4500 84 3,6532 1,9243 7,029852 3.702930

4000 108 3,6021 2,0334 7,324510 4,134715

3000 144 3,4771 2,1584 7.504972 4,658690

21,95 10.8008 39,36708 19.83922

Отсюда нижняя доверительная граница ETsom равняется 31,1 ч, а верхняя — 95,9 ч.

В итоге ЕТ|о(|) =54,7 (31,1—95,9) ч.

Таким образом, следует отмстить, что установленное нами ЕТ 50(i) для 2, 3, 6-трихлортолуола предлагаемым методом (54,7 ч) незначительно отличается от ЕТ50 рассчитанного Г. Н. Красовским и соавт. (1982) графическим методом (56 ч).

Нами выполнено более 50 контрольных расчетов ЕТ sod) в 4 координатах, т. е. в 4 вариантах: 1-й — когда Д логарифмировалась, а Т выражалось в натуральных величинах (lgD—Т), 2-й — когда Д выражалась в натуральных величинах, а Т логарифмировалось (Д — lg Т), 3-й— когда Д и Т логарифмировались (lg D — IgT), 4-й — когда Д и Т выражались в натуральных величинах (Д—Т). Для каждого из 4 вариантов расчета ЕТ so(i) доверительные границы определялись в 2 вариантах, первый из которых приведен выше. Второй же заключался в том, что сначала, последовательно подставляя в уравнение (4) значения у, равные LD|6 и LDs4, находили среднее эффективное время гибели животных для LDi6 и LD84 — ЕТ sod) и ЕТ Еоо» соответственно. Далее рассчитывали ст, т, ml (t — критерий Стьюдента). При этом 2o = ETso(i) —ЕТ sJd), а средняя ошибка (т) среднеэффективного времени гибели животных (ЕТ 50(i>) равна: 2а

m = ~yW'

где N — общее число животных в группах, которые использовались для расчета ETsSo) (доверительные границы ЕТ по(Г) находят при Р = 0,05 для числа степеней свободы f=N — 1).

Рассчитанные величины ЕТ 5o(i> в 4 координатах сопоставляли между собой и ЕТ50, рассчитанными графическим методом Г. Н. Красовского и соавт. (1982). Также сопоставлялись рассчитанные предлагаемым нами методом в 2 вариантах доверительные границы ЕТ Ь'ош . Результаты такого сопоставления для некоторых химических веществ представлены в табл. 2.

3*

— 67 —

Таблица 2

Сопоставление расчетных величины ЕТ®® (1) н их доверительных границ с ЕТ50, определенными графическим методом

Статистический показатель ЕТко. рассчитанные графи-ческим методом Г. Н. Красовского и соавт. (1982) „50 bTjodj. рассчитанные предлагаемым методом в координатах

IRD-T D-luT IgD-lgT D-T

М

Доверительные границы:

1-й вариант

2-й вариант

М

Доверительные границы:

1-й вариант

2-й вариант

М

Доверительные границы:

1-й вариант

2-й вариант

М

Доверительные границы:

1-й вариант

2-й вариант

Метил-N- (2-бензимидозолил) карбамат

36,0

Бензотриазол 37,0

Э тилкаРбитол 54,0

35,7

29,2—42,2 30,5—40,9

37,8

34,2—41,7 33,8-41,9

54,1

46,9-61,8 47,6—60,6

Этиловый эфир, тризтиленгликоля

70,0

70,7

61,2—81,0 61,4—80,0

35,9

29,5—41,9 31,0—40,8

37,7

34,2—41,6 33,6—41,7

53,3

46,6—61,2 49,8—56,8

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

69,7

60,4-80,5 60,2—79,1

34,7

28,8—41,6 29,5—39,8

37,1

33,7—41,3 32,9—41,4

52,7

46,0—60,9 46,0—59,4

68,7

59,6-80,0 58,8—78,6

36,9

29,9—42,5 31,8—42,0

38,3

34,7—42,0 34,4—42,2

54,8

47,4—62,1 48,4—61,2

71,8

62,0—81,5 62,7—80,8

Как видно из табл. 2, установленные ЕТ so(i> предлагаемым методом показатели ЕТ fod) в различных координатах практически не различаются. Нет существенной разницы между этими величинами и ЕТ50, определенными графическим методом Г. Н. Красовского и соавт. (1982). При расчете доверительных границ ЕТ fo<n в 2 вариантах полученные величины также почти одинаковы.

На основании анализа полученных результатов расчета ЕТ в различных координатах предлагаемым методом мы рекомендуем для этой цели пользоваться 4-м вариантом, когда D и Т выражаются в натуральных величинах, что при несущественной разнице с полученными с помощью

I—i гтч 50

других вариантов расчета показателями hi sou) намного упрощает технику расчетов и сокращает время определения. По той же причине доверительные границы ЕТ г'o(i) мы предлагаем рассчитывать по доверительным границам LD50, т. е. пользоваться 1-м вариантом.

В заключение следует отметить, что при хорошем соответствии расчетных оценок ЕТ sou) и приближенных графических ЕТ50, принципиально следует предпочитать расчет, который прост, требует минимума времени и лишен субъективизма. Предлагаемый метод расчета среднеэффективного времени гибели животных от среднесмертельной дозы вещества позволяет устанавливать

доверительные границы ETsoti) • Полученные таким образом вероятностные величины ЕТ &S<i» дадут возможность более точно прогнозировать хроническую токсичность химических веществ путем расчета максимально недействующих доз (МНД) в хронических опытах по предложенному для этой цели уравнению Г. Н. Красовского и соавт. (1982), рассчитывая при этом не только средние МНД по средним LD50 и ЕТ sod) ,по и доверительные интервалы МНД по доверительным интервалам LD50 и ЕТ ¿от .

Предложенный нами простой методический прием вероятностной оценки среднеэффективного времени гибели животных от среднесмертельной дозы (ЕТ to<i> ) может найти широкое применение в токсиколого-гигиенических исследованиях для прогнозирования кумулятивности и хронической токсичности вредных веществ с целью их гигиенического регламентирования в воде.

Литература. Беленький М. JI. Элементы количественной оценки фармакологического эффекта. Л., 1963, с. 117—131.

Красовский Г. Н„ Егорова Н. А., Жолдакова 3. И. и др. —

Гиг. и сан., 1982, № 7, с. 12—14. Красовский Г. Н., Жолдакова 3. И., Егорова Н. А. — В кн.: Проблема лороговости в токсикологии, М., 1979, с. 27—50.

Штабский Б. М. — Гиг. и сан., 1973, № 8, с. 24—27. Штабский Б. М, — Гиг. труда, 1974, № 1, с. 23-27.

Поступила 2I.03.S3

i Надоели баннеры? Вы всегда можете отключить рекламу.