IVRT 体外释放等效判定与统计分析

要解决什么问题

局部起效化学仿制药的 IVRT（体外释放） 等效研究，监管上通常要同时看曲线形态与释放速率：既要算 f2 相似因子，也要按 NMPA 指导原则做 T/R 斜率配对比值的非参数排序（第一阶段 36 个比值取第 8、第 29 个，界值 75%～133.33%），有时还要用 ANOVA 看制剂、吸收池、时间等因素。实验常见 6×6 TS/RS 配对（或简化 n×n），输入可能是原始浓度，也可能是已算好的释放速率或累计释放曲线。手工在 Excel 里拼统计、画图、写结论，容易在单位换算、排序分位和报告格式上出错。

技能如何工作

本技能对应脚本 ivrt_equivalence.py，按数据形态选择 --format：

• raw（推荐）：从各时间点浓度（μg/mL）经 Q_n = C_n × V、Q_n/A 等步骤得到单位面积释放与 %释放，再进入等效统计；可配合 --V、--A 及 --ts-cells / --rs-cells、含量参数。
• release-rate：已有各池释放速率（如 Higuchi Q/A = k_H·√t + b 的斜率）时，直接做 T/R 比值与排序判定。
• profile：已有平均释放曲线时，侧重 f2 等曲线相似性评价。

核心输出包括：f2 ≥ 50 判定、90% CI（排序法） 及阶段信息（36 或 324 个比值）、混合效应 ANOVA 表，以及 综合等效判决（主要标准：f2 与排序 CI 同时满足）。图表含释放速率对比、累计释放曲线、T/R 比值直方图、各池拟合线、6×6 配对热力图，另有 TS vs RS 等效判定仿真 GIF。最终生成 HTML 报告；SKILL.md 要求报告开头的「使用的统计方法与公式说明」及表图下方的公式标注卡片需在脚本跑完后人工注入，脚本不会自动生成该章节。

适用边界与注意点

• 法规锚点面向 NMPA IVRT/IVPT 指导原则及常见通则思路；不做实验设计咨询、样品制备或临床等效，仅 体外释放等效。
• f2 适用需足够时间点数且 CV 满足文档条件（早期时间点 CV 可放宽至 20%）。
• 释放速率单位约定为 ug/cm²/√h，T/R 比值以 百分单位 表达；浓度测量值保留 4 位小数。
• 与 ivrt-logan918-calc 分工：前者偏 Logan 918-12 透过量→累计释放与细胞动画；本技能偏 统计比较与等效判决；raw 路径可把计算管线内联，无需分步调用。