中藏药长期稳定性全自动分析技能
自动化分析中藏药长期稳定性（或加速稳定性）的含量测定数据，生成符合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：统计计算
对每个条件下的每批数据执行：
1. 描述性统计 — 均值、标准差、RSD、变化率
2. 线性回归分析 — 斜率、截距、R²、标准误
3. 95%置信区间 — 回归线CI + 预测区间PI
4. 批间一致性 — 斜率极差、截距差异评估
5. 合并回归分析 — 三批联合回归 + 有效期推断
步骤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个子章节）：
1. 实验设计合规性审查
2. 含量数据状态评估
3. 趋势统计学评价
4. 批间一致性评价
5. 审评风险识别
6. 模拟CDE综合意见
步骤5：HTML报告输出
自包含HTML文件，内嵌Base64编码图表，可直接用浏览器打开。
注意事项
-
依赖安装：首次使用需安装Python包：
bash pip3 install numpy matplotlib openpyxl --only-binary=:all: -
中文字体：图表标签使用英文以防CJK字体缺失，报告正文使用中文HTML显示
-
规格限：当前默认含量规格下限为14.0%，可根据实际产品规格在脚本中修改
spec_lower变量 -
数据局限性：分析结果仅基于含量数据。若需完整稳定性评价，建议补充有关物质、溶出度、水分等多指标数据
-
外推限制：有效期外推超过实际考察周期的1.5倍时（ICH Q1E §4.2.4），应谨慎解释
输出文件
分析完成后，在指定输出目录生成 稳定性分析报告_<原始文件名>.html。
