《药物稳定性试验统计分析方法》由会员分享，可在线阅读，更多相关《药物稳定性试验统计分析方法（3页珍藏版）》请在金锄头文库上搜索。
1、. .药物稳定性试验统计分析方法在确定有效期的统计分析过程中，一般选择可以定量的指标进展处理，通常根据药物含量变化计算，按照长期试验测定数值，以标示量对时间进展直线回归，获得回归方程，求出各时间点标示量的计算值y,然后计算标示量y95单侧可信限的置信区间为yz，其中：12-21式中，tN-2概率0.05，自由度N-2的t单侧分布值见表12-4，N为数组；X0给定自变量；自变量X的平均值；12-22式中，；Lyyy的离差平方和，；Lxyxy的离差乘积之和；b直线斜率。将有关点连接可得出分布于回归线两侧的曲线。取质量标准中规定的含量低限根据各品种实际规定限度确定与置信区间下界限相交点对应的时间，即
2、为药物的有效期。根据情况也可拟合为二次或三次方程或对数函数方程。此种方式确定的药物有效期，在药物标签及说明书中均指明什么温度下保存，不得使用“室温之类的名词。例：某药物在温度252，相对温度6010%的条件下进展长期实验，得各时间的标示量如表12-4。表12-4 供试品各时间的标示量时间/月03691218标示量/%99.397.697.398.496.094.0以时间为自变量x，标示量%y为因变量进展回归，得回归方程y= 99.180.26x，r=0.8970，查T单侧分布表，当自由度为4，P=0.05得tN-2=2.132当X0=0时，即0月= 2.1320.9297= 1.356按回归方
3、程计算0月时的y值得99.18%，那么y值置信区间yz，即：99.18 + 1.356=100.5499.181.356=97.82其他各时间3、6、9、12、18月的y及置信区间按同法计算，结果见表12-5。表12-5 稳定性数据表时间/月实测标示量y/%计算标示量y/%下界值yz上界值y+z099.399.1897.82100.54397.698.4097.3499.45697.397.6296.7798.47998.496.8496.0297.661296.096.0695.0897.041894.094.5092.9296.092492.9490.6195.273091.3888.27
4、94.493689.8285.9193.72用时间与y、y、yz、yz作图，得图12-6，从标示量90%处划一条直线与置信区间下界限相交，自交点作垂线于时间轴相交处，即为有效期，本例有效期为25.5个月。图12-6图12-6 药品产品有效期估算图六、经典恒温法前述实验方法主要用于新药申请，但在实际研究工作中，也可考虑采用经典恒温法，特别对水溶液的药物制剂，预测结果有一定的参考价值。经典恒温法的理论依据是前述Arrhenius的指数定律K=Ae-E/RT，其对数形式为12-23以logK对1/T作图得一直线，此图称Arrhenius图，直线斜率为-E/2.303R，由此可计算出活化能E。假设将直
5、线外推至室温，就可求出室温时的速度常数K25。由K25可求出分解10%所需的时间即T0.9或室温贮藏假设干时间以后剩余的药物的浓度。实验设计时，除了首先确定含量测定方法外，还要进展预试，以便对该药的稳定性有一个根本的了解，然后设计实验温度与取样时间。方案好后，将样品放入各种不同温度的恒温水浴中，定时取样测定其浓度或含量，求出各温度下不同时间药物的浓度变化。以药物浓度或浓度的其他函数对时间作图，以判断反响级数。假设以logc对T作图得一直线，那么为一级反响。再由直线斜率求出各温度的速度常数，然后按前述方法求出活化能和T0.9。要想得到预期的结果，除了精心设计实验外，很重要的问题是对实验数据进展正
6、确的处理。化学动力学参数如反响级数、K、E、T1/2的计算，有图解法和统计学方法，后一种方法比拟准确、合理，故近来在稳定性的研究中广泛应用。下面介绍线性回归法。例如某药物制剂，在40、50、60、70四个温度下进展加速实验，测得各个时间的浓度，确定为一级反响，用线性回归法求出各温度的速度常数，结果见表12-6。表12-6 动力学数据表T/1/T103K105/h-1lgK403.1922.664.575503.0947.944.100603.00122.383.650702.91356.503.248将上述数据logK对1/T进展一元线性回归，得回归方程：logK=T+10.64E=4765.982.3038.319=91309.77J/mol=91.31kJ/mol除经典恒温法外，还有线性变温法，Q10法，活化能估算法等，在研究工作中，有时可以应用。. v .
