本节摘要:剂量-反应关系是毒理学的"量刑标尺"。本节讲清两类曲线(个体剂量-效应、群体剂量-反应)、S 型曲线的生物学来源、LD50/ED50/TD50/治疗指数与 NOAEL/LOAEL 的定义和算法,并给出一段可复算的概率单位回归示例,最后列出四个最常见的误读。
上一节解决了"够不够格立案",本节给罪行定量。你在 1.1 看到的 WHO 分级表,其数值全部产自本节要讲的曲线拟合。
把一批实验动物分组,每组给不同剂量,记录结局,就得到剂量-反应数据。群体反应画出来几乎总是 S 型:低剂量段平缓、中段陡升、高段饱和。这不是数学家的偏好,而是生物系统的共性——个体敏感性近似对数正态分布,累加起来就是对称的 S 型。

曲线上有三个必须会读的部位:中点(LD50 或 ED50)、斜率、平台。中点是"平均量刑";斜率反映群体异质性——陡峭的曲线意味着剂量稍变、反应剧变(如氰化物),平缓的曲线意味着个体差异大;平台则提示机制饱和(受体占满、代谢酶饱和)。
四个核心参数:LD50(半数致死量)、ED50(半数有效量,药理学用)、TD50(半数中毒量)。药物安全性的经典指标治疗指数 TI 等于 TD50 除以 ED50,越大越安全;更稳健的替代是安全范围,用 ED99 与 TD01 之比,因为它照顾到曲线两端。
急性参数(LD50)回答"一次多少致命",但现代毒理学的日常是慢性低剂量问题。这里的主角换成 NOAEL(未观察到有害作用剂量)与 LOAEL(观察到有害作用的最低剂量)。它们不是拟合出来的,而是从试验剂量分组中直接读出:某剂量组与对照组相比、经统计检验无显著差异的最高剂量即 NOAEL。
从 NOAEL 到人群安全剂量,中间还要除以不确定系数(通常种内 10、种间 10,共 100 倍),得到 ADI(每日允许摄入量)。这条换算链是第 8 章的主场,这里先记住:NOAEL 是整条链的证据地基,而它天生依赖试验设计的剂量分组。
用一组经典式数据演示 LD50 的算法。某化学物经口染毒小鼠,每组 10 只:
| 剂量(mg/kg) | 死亡数 | 死亡比例 |
|---|---|---|
| 20 | 0 | 0.00 |
| 40 | 2 | 0.20 |
| 60 | 5 | 0.50 |
| 80 | 8 | 0.80 |
| 100 | 10 | 1.00 |
概率单位(probit)回归把"死亡比例"变换为正态等价的分数,再对对数剂量做线性回归。用统计软件一行代码的事,手工版过程如下(展示思路,便于核对软件输出):
# 概率单位回归求 LD50(R 语言) dose <- c(20, 40, 60, 80, 100) dead <- c(0, 2, 5, 8, 10) n <- rep(10, 5) fit <- glm(cbind(dead, n - dead) ~ log(dose), family = binomial(link = "probit")) summary(fit) # 看斜率与截距的显著性 ld50 <- exp(-coef(fit)[1] / coef(fit)[2]) ld50 # 输出约为 60 mg/kg # 95% 置信区间可用 MASS 包的 dose.p 函数: # library(MASS); dose.p(fit, p = 0.5)
输出 LD50 约 60 mg/kg,对照 1.1 的 WHO 表落在"中等毒性"档(50–500)。解读要点:0% 与 100% 的组对拟合贡献的信息有限(probit 变换把它们推向无穷),所以试验设计要让中间组落在 20%–80% 区间;置信区间比点估计更重要——本例区间较宽,因为每组只有 10 只动物,这正是 1.3 节统计学的伏笔。
⚠️ 常见坑:比较两种物质的毒性时只报 LD50 点值、不报置信区间与物种——两份 LD50 相差不到置信区间宽度的比较没有意义。
定罪需要定量,定量依赖证据的可靠性。下一节处理最后一块立案拼图:哪些变量会让同样的剂量产生悬殊的后果,以及统计学如何防止冤案。
问:为什么 LD50 的经典地位在下降? 因为监管实际需要的是毒性分级与治疗参考,而不是精确到个位数的半数致死量;分级法与限值试验用少得多的动物就能给出分级结论(这是第 7 章 3R 的入口)。LD50 仍在比较学与教学中使用,但不再是大宗必需品。
问:S 型曲线的斜率有什么实用价值? 斜率告诉你"剂量稍变、反应剧变"的风险程度:陡曲线(氰化物类)意味着抢救时间窗以分钟计,平缓曲线意味着个体差异大、需要更宽的安全边际。斜率还影响 BMD 置信区间的宽度——第 8 章的推导会用到。
问:NOAEL 和 ADI 是一回事吗? 不是。NOAEL 是动物试验读出的剂量,ADI 是 NOAEL(或 BMDL)除以不确定系数后的人群管理参考值。前者是证据,后者是决策——这个区分贯穿第 8 章。
从 U 型曲线到毒物兴奋效应:本节的 S 型框架默认"剂量越高效应越强",但带反馈调节的体系(内分泌、应激反应)可能出现低剂量激活、高剂量抑制的倒 U 型。判读这类数据的第一原则是终点要客观、剂量点要足够密——非单调结论比单调结论需要更强的数据支撑,这也是第 4 章内分泌干扰争议的方法论根源。
把本章参数连成一句话:LD50 给急性分级,NOAEL 给慢性起点,BMDL 给更稳的起点,治疗指数给药物窗口——它们共享同一条剂量-反应逻辑,只是终点与时间尺度不同。抓住这条线,第 8 章的推导只剩除法。