用电子表格Excel计算药物的有效期

用电子表格Excel计算药物的有效期

文档由作者亲自提供，药学研究人员必备！

摘要：[目的]计算出药物有效期。[方法]运用电子表格Excel软件处理药物稳定性试验数据。[结果]计算结果与文献一致。[结论]与图解法相比，该方法快速，简便，并可避免由于观察读取图上刻度所致的误差，能精确地算出有效期。

关键词：有效期；稳定性试验；Excel软件

Calculating Term Of Validity Of Pharmaceuticals By Microsoft Excel

Liu Xing，Li Jin-Wei，Liu jin-yi

(Beijing Tri-Prime Genetic Engineering Co., Ltd. Beijing 102600,China) ABSTRACT: [AIM] To Calculate the Term Of Validity Of Pharmaceuticals . [METHODS] The Data Of Experiment Were Treat With Microsoft Excel . [RESULTS] The Calculation by Microsoft Excel were identical with those of literature.[CONCLUSION] Compare With Graphic Method, It is faster and more convenient. It prevents errors due to reading on graphs, and can calculates Term Of Validity precisely.

KEY WORDS: Validity Of Pharmaceuticals ;Stability Test; Microsoft Excel

新药原料药或制剂，在投产前必须进行稳定性试验，其中有效期的确定很重要。试验方法是采用长期考察，在接近药品实际贮存条件下，取成品包装的3批样品，在温度(25±2)℃，相对湿度(60±5)%条件下，放置不同的时间取样，分析检验稳定性考察各项目。其中含量(相当于标示量的百分数)测定数据，可作为定量指标，规定方法［1］是以含量对时间(月)计算得回归方程式，求出各时间(月)点的含量计算值，再计算在P=0.95的可信区间±Z，将各点连接可得分布于回归直线两侧的曲线。取质量标准规定的含量低限(一般为90%)，在纵坐标90%处作水平线，此线与置信区间下限线-Z相交点，再从横坐标上找与该交点相对应的时间(月)点，即为产品的有效期。

基于以上所述，本文介绍一种应用Excel软件计算有效期的方法，该方法简便，准确。

1 运算原理及公式

按照试验测定的数据，以标示量（%）对时间进行直线回归，得回归方程，

求出计算标示量的表达式，然后再求出标示量95%单侧可信限的下界限的表达式。利用Excel单变量求解功能，即对于y1=f(x1)，给定y1值，反过来求解x1求出有效期x1。

计算标示量的表达式[1]：

=b*x+a

b=(n*Σxy-(Σx)*(Σy))/ (n*Σx2-(Σx)2)

a=(Σy-(Σx)*b)/ n

Z=t n-2*s*SQRT(1/n+(x1-Σx/n)2/（Σ(x-Σx/n)2）)

s= SQRT(Q/(n-2))

Q=Lyy-b*Lxy=(Σy2-(Σy)2/n)-(n*Σxy-(Σx)*(Σy))/ (n *Σx2-(Σx)2)*(Σxy-(Σx)*(Σy)/ n)

95%单侧可信限的下界限的曲线方程表达式：

y1= b*x1+a-Z

y1=((n*Σxy-(Σx)*(Σy))/ (n*Σx2-(Σx)2))*x1+(Σy-(Σx)*( (n*Σxy-(Σx)*(Σy))/ (n*Σx2-(Σx)2)))/n-(t n-2*( SQRT(((Σy2-(Σy)2/n)-(n*Σxy-(Σx)*(Σy))/ (n*Σx2-(Σx)2)*(Σxy-(Σx)*(Σy)/n))/(n-2)))*SQRT(1/n+(x1-Σx/n)2/（Σ(x-Σx/n)2）)

上述x1为自变量为C2，y1为因变量，x和y分别为时间和标示量，n为自变量的数。

2 输入公式

打开Excel，在单元格A1中输入“时间（月）”，在单元格B1中输入“标示量（%）”，在单元格C1中输入“有效期（月）”。以A2:A100和B2:B100分别为时间和标示量数据输入区。在单元格D2中输入95%单侧可信限的下界限的曲线方程表达式：

“=((COUNT(A2:A100)*SUMPRODUCT(A2:A100,B2:B100)-SUM(A2:A100) *SUM(B2:B100))/(COUNT(A2:A100)*SUMSQ(A2:A100)-SUM(A2:A100)^2))*C2 +((SUM(B2:B100)-((COUNT(A2:A100)*SUMPRODUCT(A2:A100,B2:B100)-SUM (A2:A100)*SUM(B2:B100))/(COUNT(A2:A100)*SUMSQ(A2:A100)-SUM(A2:A10 0)^2))*SUM(A2:A100))/COUNT(A2:A100))-(TINV(0.1,(COUNT(A2:A100)-2)))*((

SUMSQ(B2:B100)-SUM(B2:B100)^2/COUNT(A2:A100)-(COUNT(A2:A100)*SU MPRODUCT(A2:A100,B2:B100)-SUM(A2:A100)*SUM(B2:B100))/(COUNT(A2:A 100)*SUMSQ(A2:A100)-SUM(A2:A100)^2)*(SUMPRODUCT(A2:A100,B2:B100) -SUM(A2:A100)*SUM(B2:B100)/COUNT(A2:A100)))/((COUNT(A2:A100))-2))^0. 5*((1/(COUNT(A2:A100))+(C2-(SUM(A2:A100)/COUNT(A2:A100)))^2/(DEVSQ( A2:A100))))^0.5”

其数据输入区分别为A2:A100和B2:B100，n为COUNT(A2:A100)，t n-2为TINV(0.1,(COUNT(A2:A100)-2)，Σx和Σy分别为SUM(A2:A100)和SUM(B2:B100)，Σx2和Σy2分别为SUMSQ(A2:A100)和SUMSQ(B2:B100)，Σxy为SUMPRODUCT(A2:A100,B2:B100)，Σ(x-Σx/n)2为DEVSQ(A2:A100)。

2 求解步骤

先从工具菜单选用单变量求解项，弹出单变量求解对话框，如图1所示，在目标单元格编辑框内用鼠标选取单元格D2，该单元格为求解公式所在位置，在目标值编辑框内键入需要得到的目标值90，在可变单元格编辑框内用鼠标选取单元格C2。

图1

单击“确定”按钮，弹出单变量求解状态框，如图2所示，显示当前解和目标值。

图2

继续单击“确定”按钮，在单元格C2内得到求解结果，如图3所示。