药代动力学非房室分析（NCA）：AUC、半衰期、清除率在线计算

非房室分析（NCA）不假定任何房室模型，直接从浓度-时间点算出暴露量和 处置参数。它是生物等效性和毒代研究的默认方法，也是绝大多数人用 Excel 手搓梯形面积的那件事。

公式

AUC(0→t_last)  = Σ (Ci + Ci+1)/2 × Δt                     线性梯形
AUC(t_last→∞)  = C_last / λz                              外推尾部
AUMC(0→t_last) = Σ (ti·Ci + ti+1·Ci+1)/2 × Δt
AUMC(t_last→∞) = C_last·t_last/λz + C_last/λz²

λz   = 末端若干点对 ln C 做线性回归，取斜率的相反数（至少 3 个点）
t½   = ln2 / λz
MRT  = AUMC(0→∞) / AUC(0→∞)
CL   = 剂量 / AUC(0→∞)            血管外给药时得到的是 CL/F
Vz   = 剂量 / (λz × AUC(0→∞))     血管外给药时是 Vz/F
Vss  = CL × MRT                    只对静脉给药有意义

四条面积公式出自 Gabrielsson & Weiner, Non-compartmental analysis, Methods Mol Biol 2012;929:377-89 的式 (1)–(4)。

为什么默认不是纯线性梯形

线性梯形法用直线连接相邻两点，而真实的消除相是往下凹的指数曲线—— 直线永远在曲线上方，所以下降段一律高估；上升段同理一律低估。 这不是什么微小误差，取决于你采样有多稀。

默认的 linear-up/log-down 只在上升段用直线，下降段改用对数梯形：

下降段（Ci+1 < Ci）：ΔAUC = (Ci − Ci+1) × Δt ÷ ln(Ci / Ci+1)

这个式子对单指数段是精确的，所以消除相估得准得多。

可以验算

用单指数静注数据 C(t) = 100 × e^(−0.2t)（mg/L），剂量 1000 mg， 采样 0 / 1 / 2 / 4 / 8 / 12 / 24 h。这组数据每一项都有解析真值， 不是从哪本书上抄的数：

点开工具就是这组数据，可以直接核。右边那一列就是不换方法的代价： 清除率会低估 5.6%，而清除率是决定给药剂量的那个参数。

残余的不到 0.1% 来自浓度只填了两位小数，不是方法误差。

适用边界

• λz 拟合点数会改变所有下游参数。本工具默认取末端 4 点，可以改。 自动模式按校正 R² 最大挑点数——这是通行做法，但它选出来的点数 未必是你看图会选的那几个，拟合用了哪些点工具会列出来，自己看一眼。
• 外推占比 > 20% 时这个 AUC(0→∞) 不可靠，工具会告警。 生物等效性指南（EMA / FDA）通行要求 AUC(0→t) 至少覆盖 AUC(0→∞) 的 80%。 超了说明采样时间不够长，不是算法问题，补采样才能解决。
• CL 和 Vz 在血管外给药时带着生物利用度 F，得到的是 CL/F、Vz/F。 把 CL/F 当 CL 用会把清除率系统性高估 1/F 倍。工具按你选的途径标注。
• Vss = CL × MRT 只对静脉给药成立。口服时 MRT 里混进了吸收时间 （MAT），算出来的不是真正的稳态分布容积。
• 不做多剂量累积、不做稳态参数（AUCτ、Cav、波动度），不做统计学比较。 单次给药的单个个体，一次一条曲线。
• 浓度低于定量下限（BLQ）的点怎么处理会明显影响结果，而规则因研究而异 （首个 BLQ 计 0、全部剔除、按 LLOQ/2 代入……）。本工具按你填的数算， 不替你做 BLQ 规则，填进来的 0 会让该段退回线性梯形。