中藏药长期稳定性全自动分析技能

自动化分析中藏药长期稳定性（或加速稳定性）的含量测定数据，生成符合ICH Q1A/Q1E和《中药稳定性研究技术指导原则》的完整HTML分析报告。

使用时机

当用户提出以下需求时触发：

分析长期稳定性/加速稳定性的含量测定数据
对稳定性数据进行回归分析和趋势评价
推断有效期（shelf life estimation）
评价批间一致性（batch consistency）
评估温度敏感性（temperature sensitivity）
需要从CDE审评专家视角审查稳定性数据

数据格式要求

Excel文件应包含以下结构（参考 /storage/Users/currentUser/Desktop/长期稳定性案例1.xlsx）：

| 条件 | 批号 | 0 | 3个月 | 6个月 | 9个月 | ... | |------|------|---|-------|-------|-------|-----| | 25℃，60%RH | 24032601 | 15.2 | 15.6 | 15.1 | ... | ... | | 20℃，55%RH | 24032601 | 15.2 | 15.4 | ... | ... | ... |

必须包含"条件"列（标注温湿度条件名称）和"批号"列
时间点列用数字或"X个月"格式标记
支持任意数量的条件、批次和时间点（推荐3批×7时间点×2条件）

使用方法

# 基本用法
python3 scripts/analyze_stability.py <Excel文件路径>

# 指定输出目录
python3 scripts/analyze_stability.py <Excel文件路径> <输出目录>

分析流程

步骤1：数据解析

自动读取Excel文件，识别条件、批次、时间点结构，验证数据完整性。

步骤2：统计计算

对每个条件下的每批数据执行：

描述性统计 — 均值、标准差、RSD、变化率
线性回归分析 — 斜率、截距、R²、标准误
95%置信区间 — 回归线CI + 预测区间PI
批间一致性 — 斜率极差、截距差异评估
合并回归分析 — 三批联合回归 + 有效期推断

步骤3：图表生成（11张）

| 图表 | 内容 | 法规依据 | |------|------|---------| | 含量趋势图 | 各批含量随时间变化的折线图 | ICH Q1A | | 线性回归拟合 | 各批散点+回归线，标注R² | ICH Q1E | | ICH Q1E核心图 | 回归+95%CI+预测区间+规格限 | ICH Q1E §4.1 | | 箱线图 | 各时间点批间分布与离群值 | FDA Guidance | | 残差分析 | 残差vs拟合值 + 残差vs时间 | ICH Q1E回归诊断 | | 合并回归+有效期 | 三批联合回归+CI+规格限 | ICH Q1E §4.2 | | 变化率对比 | 24个月变化率的批间/条件间对比 | 中药指导原则 |

步骤4：CDE审评分析生成

自动生成审评视角分析章节（6个子章节）：

实验设计合规性审查
含量数据状态评估
趋势统计学评价
批间一致性评价
审评风险识别
模拟CDE综合意见

步骤5：HTML报告输出

自包含HTML文件，内嵌Base64编码图表，可直接用浏览器打开。

注意事项

依赖安装：首次使用需安装Python包：

pip3 install numpy matplotlib openpyxl --only-binary=:all:

中文字体：图表标签使用英文以防CJK字体缺失，报告正文使用中文HTML显示
规格限：当前默认含量规格下限为14.0%，可根据实际产品规格在脚本中修改spec_lower变量
数据局限性：分析结果仅基于含量数据。若需完整稳定性评价，建议补充有关物质、溶出度、水分等多指标数据
外推限制：有效期外推超过实际考察周期的1.5倍时（ICH Q1E §4.2.4），应谨慎解释

输出文件

分析完成后，在指定输出目录生成 稳定性分析报告_<原始文件名>.html。