稳定性研究分析目的稳定性研究用于评价试剂、校准品或质控品在规定储存、运输或使用条件下检测结果随时间变化的程度。本文演示如何使用ivdtools包建立稳定性分析计划完成节点偏差、条件间及 相对首节点的回归趋势、达到允许限值的时间预测、阿伦纽斯Arrhenius加速稳定性分析及平均动力学温度meankinetic temperatureMKT计算。::: celllibrary(ivdtools) library(readr):::本文数据均为确定性的教学示例仅用于验证分析流程。示例中的 5%允许漂移、10% 加速降解限值、目标温度和活化能不是通用接受标准正式研究必须使用方案、风险评估或适用标准预先规定的参数。选择浓度值或信号值作为响应值应根据方案确定。函数概述函数 主要用途 关键输入或输出stability_plan()规划回归系列的最少时间点 允许漂移、变异、重复数、预期漂移和效能stability_bias()计算各条件相对首节点的绝对/相对偏差 节点汇总、条件内偏差、条件间偏差和限值状态stability_regression()单条件趋势或试验/参考条件配对偏差回归 斜率、单侧置信限、节点结果和诊断信息stability_time()预测点估计或单侧置信限达到允许偏差的时间 预测时间、是否外推及搜索范围arrhenius()根据多个升高温度的降解速率外推目标温度稳定期 动力学级数、速率、活化能和预测稳定期mkt()计算等间隔或时间加权的平均动力学温度 各探头及合并 MKTstability_regression()支持两种模式mode self对单一条件的测量结果随时间回归mode compare在匹配时间点计算试验条件相对参考条件的偏差再对偏差随时间回归。分析前应明确条件映射、时间单位、偏差尺度、预期变化方向和允许限。只有在方案提供接受标准时才能据此作出是否可接受的结论。示例一建立分析计划假设研究方案规定允许相对漂移为 5%预期漂移为 2.5%重复性变异为1.5%目标效能为 90%。比较每个时间点测量 1、2 或 3次时所需的最少时间点::: cellplan - ivdtools::stability_plan( allowable_drift 5, variability 1.5, replicates c(1, 2, 3), expected_drift 2.5, power 0.9, mode simple, bias_type relative ) print(plan):::示例二偏差分析读取和核验原始数据::: cellbias_data - read_csv( ./data/stability-bias-example.csv, show_col_types FALSE ) bias_data:::::: celldim(bias_data) str(bias_data):::节点偏差计算以2-8C为参考条件按各条件自身首节点归一化评价相对偏差是否在示例性的 ±5%限值内::: cellbias_result - ivdtools::stability_bias( data bias_data, condition Condition, time Day, value Result, reference_condition 2-8C, limit 5, bias_type relative, time_unit d ) print(bias_result):::条件内偏差以各条件自己的首节点为基线条件间偏差则在相同时间点比较试验条件和参考条件。::: cellplot(bias_result):::示例三不同保存条件比较的回归分析读取和核验原始数据::: cellregression_data - read_csv( ./data/stability-regression-example.csv, show_col_types FALSE ) regression_data:::试验条件与参考条件比较在每个匹配时间点计算25C相对2-8C的相对偏差并对偏差随时间进行回归。示例预先指定结果随时间降低因此使用direction decrease::: cellregression_result - ivdtools::stability_regression( data regression_data, time Day, value Result, condition Condition, mode compare, test_condition 25C, reference_condition 2-8C, bias_type relative, direction decrease, conf.level 0.95, limit 5, time_unit d ) print(regression_result):::::: cellplot(regression_result):::模型使用单侧置信限评价最不利方向。本示例识别出一个 Cook距离较大的节点但函数不会自动删除应回查原始记录、实验过程和方案规定并在有充分理由时进行包含/不包含该节点的敏感性分析。示例四与首节点比较的回归分析单个条件相对首节点分析除示例三的试验条件/参考条件配对回归外mode self可分析单个条件自身随时间的变化。 以下仍使用回归示例数据但仅选择25C条件各节点的observed_bias相对于该条件首节点 均值计算::: cellself_regression_result - ivdtools::stability_regression( data regression_data, time Day, value Result, condition Condition, mode self, test_condition 25C, bias_type relative, direction decrease, conf.level 0.95, limit 5, time_unit d ) print(self_regression_result):::::: cellplot(self_regression_result):::首节点观测偏差固定为0模型估计值则由全部时间点共同拟合不要求回归截距恰好通过首节点。单侧置信限用于评价预先指定的最不利方向示例性的 ±5%限值仍不能替代研究方案规定。示例五预测达到限值的时间stability_time()使用示例三的回归对象分别计算点估计和单侧置信限首次达到 -5%偏差的时间::: celltime_result - ivdtools::stability_time( regression_result, limit 5, max_time 90 ) print(time_result):::::: cellplot(time_result):::本例最长实际观察时间为 42天而预测达到限值的时间超过了观察范围因此输出标记为EXTRAPOLATED。外推结果依赖线性趋势在未观察区间继续成立的假设不能等同于实测稳定性证据应显著报告外推距离和不确定性并优先用更长的实际观察数据确认。示例六阿伦纽斯Arrhenius加速稳定性分析读取和核验原始数据::: cellarrhenius_data - read_csv( ./data/stability-arrhenius-example.csv, show_col_types FALSE ) arrhenius_data:::动力学和阿伦纽斯拟合示例使用 30、37 和 45 °C 的加速数据外推至 5 °C允许相对变化10%并预先指定 降解方向为降低::: cellarrhenius_result - ivdtools::arrhenius( data arrhenius_data, temperature Temperature, time Day, value Result, target_temp 5, order auto, direction decrease, limit 10, bias_type relative, temp_unit C, time_unit d, conf.level 0.95 ) print(arrhenius_result):::::: cellplot(arrhenius_result):::order auto会在零级和一级动力学候选中按实现规则选择模型。正式分析应在实验设计和机理知识支持下确定候选级数并核对每个温度的降解方向、速率显著性、降解程度和模型残差。从升高温度外推至目标储存温度可能跨越很大范围预测稳定期及其区间必须与模型假设、实际时间数据和验证研究一起解释。示例七MKT 计算读取和核验原始温度数据::: cellmkt_data - read_csv( ./data/stability-mkt-example.csv, show_col_types FALSE ) mkt_data$Timestamp - as.POSIXct( mkt_data$Timestamp, format %Y-%m-%d %H:%M:%S, tz UTC ) mkt_data:::时间加权 MKT使用两个温度探头和实际时间间隔温度单位为摄氏度活化能采用示例值83.144 kJ/mol::: cellmkt_result - ivdtools::mkt( data mkt_data, temp_cols c(Probe_A, Probe_B), time Timestamp, temp_unit C, ea 83.144 ) print(mkt_result):::MKT不是温度的简单算术平均值高温暴露会因阿伦纽斯关系获得更大权重。正式计算应确认探头校准状态、时间顺序、缺失区间、温度单位、曝光持续时间以及活化能假设。MKT只能概括 温度暴露不能替代产品特异的稳定性研究。多样本数据拆分当不同浓度水平、样本类型或批次需要分别评价时应保留样本标识先拆分数据再对每个样本使用相同且预先规定的参数分析。读取并核验多样本原始数据::: cellmultiple_samples_data - read_csv( ./data/stability-multiple-samples-example.csv, show_col_types FALSE ) multiple_samples_data:::使用split()建立按样本命名的数据列表并先查看其中一个样本::: cellsample_data_list - split( multiple_samples_data, multiple_samples_data$Sample ) names(sample_data_list) sample_a_data - sample_data_list[[Sample_A]] sample_a_data:::逐个样本执行分析。结果判读要点方案优先分析前明确条件、节点、重复数、时间单位、方向、偏差尺度和接受限。原始数据只读报告缺失、重复、非数值记录和每一步排除行不自动删除异常记录。配对正确条件比较必须在匹配时间点进行并明确试验条件和参考条件。模型诊断报告斜率、置信限、拟合方法、残差、影响点及是否发生 WLS回退。外推醒目超过观察范围的稳定时间预测应明确标注不作为实测证据。温度假设明确阿伦纽斯分析和 MKT必须记录摄氏/开尔文单位、时间单位和活化能。统计与接受性分开统计显著不等于产品可接受只能依据预设允许限作出判断。独立复核临床、监管或放行决策应由合适的专业人员独立复核。