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

我的该文章已经在发表在：药学进展，2003；27（3）：176．

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

刘星，李巾伟，刘金毅

(北京三元基因工程有限公司，北京，102600)

摘要：[目的]计算出药物有效期。[方法]运用电子表格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有相交点，再从横坐标上找与该交点相对应的时间(月)点，即为产品的有效期。

基于以上所述，本文介绍一种应用Excel2000软件计算有效期的方法，该方

法简便，准确。

1 运算原理及公式

按照试验测定的数据，以标示量（%）对时间进行直线回归，得回归方程，求出计算标示量的表达式，然后再求出标示量95%单侧可信限下界限的表达式。利用Excel单变量求解功能，即对于y=f(x)，给定y值，反过来求解x求出有效期x。

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

y=b*x+a

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

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

n为自变量的数。

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

= b*x+a-Z

Z=t

*s*SQRT(1/n+(x-Σx/n)2/（Σ(x-Σx/n)2）)

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)

2 输入公式

打开Excel2000，在单元格A1中输入“时间（月）”，在单元格B1中输入“标示量（%）”，在单元格C1中输入“有效期（月）”。以A2:A100和 B2:B100分别为时间和标示量数据输入区。在单元格D2中输入以下公式：

“=((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:A100)^2))*SUM(A2:A100 ))/COUNT(A2:A100))-(TINV(0.1,(COUNT(A2:A100)-2)))*((SUMSQ(B2:B100)-SU M(B2:B100)^2/COUNT(A2:A100)-(COUNT(A2:A100)*SUMPRODUCT(A2:A100,B2:B10 0)-SUM(A2:A100)*SUM(B2:B100))/(COUNT(A2:A100)*SUMSQ(A2:A100)-SUM(A2:A 100)^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”

上述公式的表达式即为95%单侧可信限的下界限的曲线方程表达式：

y1=((n*Σxy-(Σx)*(Σy))/ (n*Σx2-(Σx)2))*x1+(Σy-(Σx)*( (n*Σ

*( SQRT(((Σy2-(Σy)2/n)-(n*Σxy-(Σx)*(Σy))/ (n*Σx2-(Σx)2)))/n-( t

n-2

xy-(Σx)*(Σy))/ (n*Σx2-(Σx)2)*(Σxy-(Σx)*(Σy)/n))/(n-2)))*SQRT(1/n+(x1-Σx/n)2/(Σx2+n*(Σx/n)2-2*(Σx/n)*(Σx)))上述x1为自变量为C2，y1为因变量，x和y分别为A2:A100和 B2:B100，

为TINV(0.1,(COUNT(A2:A100)-2)，Σx和Σy分别为n为COUNT(A2:A100)，t

n-2

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所示，显示当前解和目标值。