Метод Спирмена—Кёрбера

Метод Спирмена—Кёрбера — это статистический метод, используемый в биологии, медицине и фармакологии для расчёта средних эффективных доз (например, ED₅₀ — доза, вызывающая эффект у 50 % особей) и средних летальных доз (LD₅₀ — доза, вызывающая гибель 50 % подопытных животных). Метод основан на интерполяции кумулятивных частот эффекта в зависимости от дозы и позволяет оценить медианную дозу по ограниченному числу экспериментальных точек, не требуя полного построения кривой «доза—эффект».

История

Метод был разработан независимо двумя учёными в начале XX века. Английский психолог и статистик Чарльз Спирмен (Charles Spearman) в 1908 году предложил способ оценки пороговых доз на основе пробит-анализа, однако его работа не получила широкого распространения. Немецкий фармаколог Фридрих Кёрбер (Friedrich Körber) в 1910 году модифицировал подход, адаптировав его для токсикологических экспериментов. Впоследствии метод был объединён под названием «метод Спирмена—Кёрбера» и стал одним из стандартных в биологических исследованиях.

В СССР и России метод активно применялся с 1930-х годов, особенно в токсикологии и фармакологии, где требовалось быстрое определение LD₅₀ для новых химических соединений. С развитием компьютерных технологий в конце XX века метод был частично вытеснен более точными методами, такими как пробит-анализ и логистическая регрессия, но остаётся востребованным в ситуациях, когда число экспериментальных групп ограничено (например, при работе с редкими или дорогими веществами).

Принцип метода

Метод Спирмена—Кёрбера применяется к данным, полученным в экспериментах с фиксированными дозами (обычно 4–6 доз), каждая из которых вводится группе подопытных животных. Для каждой дозы регистрируется доля особей, у которых проявился эффект (например, гибель или терапевтический ответ). Метод предполагает, что зависимость между логарифмом дозы и вероятностью эффекта имеет S-образную форму (сигмоидальную кривую), а медианная доза соответствует точке, где эффект наблюдается у 50 % особей.

Алгоритм расчёта

1. Логарифмирование доз: Все дозы переводятся в логарифмическую шкалу (обычно по основанию 10 или e).
2. Вычисление разности логарифмов: Определяется шаг между соседними дозами (\(d\)) — разность логарифмов соседних доз.
3. Расчёт кумулятивных частот: Для каждой дозы вычисляется доля особей с эффектом (\(p_i\)).
4. Интерполяция: Медианная доза (\(ED_{50}\) или \(LD_{50}\)) рассчитывается по формуле:

\[ \log_{10}(ED_{50}) = \log_{10}(D_{\text{max}}) - \frac{d}{2} \cdot \left( \sum_{i=1}^{k} p_i \right) \]

где:

• \(D_{\text{max}}\) — максимальная доза в эксперименте,
• \(d\) — шаг между логарифмами доз,
• \(p_i\) — доля особей с эффектом в каждой дозовой группе,
• \(k\) — количество дозовых групп.

Формула может варьироваться в зависимости от того, используется ли логарифмирование по основанию e или 10, а также от того, включается ли в расчёт нулевая доза (контрольная группа).

Условия применимости

• Дозы должны быть равноудалены в логарифмической шкале (геометрическая прогрессия).
• Число доз — не менее 4, оптимально 5–6.
• В каждой дозовой группе — не менее 5–10 особей.
• Эффект должен быть бинарным (например, «жив»/«мёртв») или квантильным.
• Кривая «доза—эффект» должна быть монотонной (с ростом дозы доля эффекта увеличивается).

Преимущества и недостатки

Преимущества

• Простота расчёта: Не требует сложного программного обеспечения; может быть выполнен вручную или в электронной таблице.
• Экономичность: Позволяет получить оценку LD₅₀ при относительно небольшом числе животных (что соответствует принципам гуманного обращения с животными).
• Устойчивость к выбросам: В отличие от некоторых итеративных методов, нечувствителен к единичным аномальным значениям.

Недостатки

• Ограниченная точность: Даёт лишь приближённую оценку, особенно при малом числе доз или неравномерном распределении эффекта.
• Неприменимость для несимметричных кривых: Метод предполагает симметрию кривой «доза—эффект» относительно медианной точки, что не всегда выполняется.
• Отсутствие доверительных интервалов: Классическая версия метода не позволяет рассчитать стандартную ошибку или доверительный интервал для ED₅₀ (хотя существуют модификации, например, метод Спирмена—Кёрбера с бутстреп-оценкой).

Применение

Метод Спирмена—Кёрбера широко используется в следующих областях:

• Токсикология: Определение LD₅₀ для оценки острой токсичности химических веществ, пестицидов, лекарственных препаратов.
• Фармакология: Расчёт ED₅₀ для оценки эффективности анальгетиков, анестетиков, противосудорожных средств.
• Микробиология: Определение летальных доз для микроорганизмов (например, при оценке эффективности антисептиков).
• Экология: Оценка токсичности загрязнителей для водных организмов (например, рыбы или дафний).

В России метод включён в ряд методических указаний по токсикологическим исследованиям, в частности, в «Руководство по экспериментальному (доклиническому) изучению новых фармакологических веществ» (2005) и ГОСТ 32644-2014 «Методы испытаний по воздействию химической продукции на организм человека. Оценка острой токсичности».

Модификации

Существует несколько модификаций метода, направленных на повышение точности:

• Метод Спирмена—Кёрбера с поправкой на контроль: Если в контрольной группе (нулевая доза) наблюдается спонтанный эффект, в расчёт вводится поправочный коэффициент (формула Эбботта).
• Метод Спирмена—Кёрбера для пробит-анализа: Комбинация с пробит-преобразованием для линеаризации кривой.
• Бутстреп-версия: Позволяет оценить доверительные интервалы путём многократной ресэмплинга данных.

Интересные факты

• Метод Спирмена—Кёрбера иногда называют «методом средних доз» или «методом интерполяции по кумулятивным частотам».
• В англоязычной литературе метод известен как Spearman–Kärber method (иногда с вариантом написания «Karber»).
• Несмотря на появление более точных методов, метод Спирмена—Кёрбера остаётся стандартом в ряде регламентирующих документов, например, в рекомендациях OECD (Организация экономического сотрудничества и развития) по тестированию химических веществ (OECD TG 423).

Источники

• Spearman, C. (1908). «The method of 'right and wrong cases' (constant stimuli) without Gauss's formulae». British Journal of Psychology, 2(3), 227–242.
• Körber, F. (1910). «Über die Bestimmung der mittleren letalen Dosis». Archiv für experimentelle Pathologie und Pharmakologie, 63(1), 1–14.
• Finney, D. J. (1971). Probit Analysis (3rd ed.). Cambridge University Press.
• Руководство по экспериментальному (доклиническому) изучению новых фармакологических веществ / Под ред. Р. У. Хабриева. — М.: Медицина, 2005. — 832 с.
• ГОСТ 32644-2014 «Методы испытаний по воздействию химической продукции на организм человека. Оценка острой токсичности». — М.: Стандартинформ, 2014.