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

我的该文章已经在发表在：药学进展，2003；27（3）：176用电子表格 Excel 计算药物的有效期刘星，李巾伟，刘金毅(北京三元基因工程有限公司，北京，102600)摘要：目的 计算出药物有效期。 方法 运用电子表格 Excel 软件处理药物稳定性试验数据。结果计算结果与文献一致。结论 与图解法相比，该方法快速，简便，并可避免由于观察读取图上刻度所致的误差，能精确地算出有效期。关键词：有效期；稳定性试验；Excel 软件Calculating Term Of Validity Of Pharmaceuticals By Microsoft ExcelLiu 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 批样品，在温度(252)，相对湿度(605)% 条件下，放置不同的时间取样，分析检验稳定性考察各项目。其中含量(相当于标示量的百分数)测定数据，可作为定量指标，规定方法 1 是以含量对时间(月)计算得回归方程式，求出各时间(月)点的含量计算值 ，再计算 在 P=0.95 的可信区间 Z，将各点连接可得分布于回归直线两侧的曲线。取质量标准规定的含量低限(一般为 90%)，在纵坐标 90%处作水平线，此线与置信区间下限线 -Z 有相交点，再从横坐标上找与该交点相对应的时间(月)点，即为产品的有效期。基于以上所述，本文介绍一种应用 Excel2000 软件计算有效期的方法，该方法简便，准确。1 运算原理及公式按照试验测定的数据，以标示量（%）对时间进行直线回归，得回归方程，求出计算标示量的表达式，然后再求出标示量 95%单侧可信限下界限的表达式。利用 Excel 单变量求解功能，即对于 y=f(x)，给定 y 值，反过来求解 x 求出有效期 x。计算标示量的表达式 1：y=b*x+ab=(n*xy-(x)*(y)/ (n*x 2-(x) 2)a=(y-(x)*b)/nn 为自变量的数。95%单侧可信限的下界限的曲线方程表达式：= b*x+a-ZZ=tn-2*s*SQRT(1/n+(x-x/n) 2/（(x-x/n) 2）)s= SQRT(Q/(n-2)Q=Lyy-b*Lxy=(y 2-(y) 2/n)-(n*xy-(x)*(y)/ (n*x 2-(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)-SUM(B2:B100)2/COUNT(A2:A100)-(COUNT(A2:A100)*SUMPRODUCT(A2:A100,B2:B100)-SUM(A2:A100)*SUM(B2:B100)/(COUNT(A2:A100)*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”上述公式的表达式即为 95%单侧可信限的下界限的曲线方程表达式：y1=(n*xy-(x)*(y)/ (n*x 2-(x) 2)*x1+(y-(x)*( (n*xy-(x)*(y)/ (n*x 2-(x) 2)/n-( tn-2*( SQRT(y 2-(y) 2/n)-(n*xy-(x)*(y)/ (n*x 2-(x) 2)*(xy-(x)*(y)/n)/(n-2)*SQRT(1/n+(x1-x/n) 2/(x 2+n*(x/n) 2-2*(x/n)*(x)上述 x1 为自变量为 C2，y1 为因变量，x 和 y 分别为 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)，x 2和 y 2分别为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 所示。图 33 应用实例某药物在温度（252），相对湿度（605）%条件下进行长期试验，得各时间的标示量如下表 1。表 1 样品各时间的标示量（%）时间（月） 0 3 6 9 12 18标示量（%） 99.3 97.6 97.3 98.4 96.0 94.04 讨论和注意事项该实例引自魏树礼教授主编的生物药剂学与药物动力学 2，书中采用作图法求出的药物有效期结果为 25.5 月，采用 Excel 计算的药物有效期结果为25.57 月。采用 Excel 软件计算，所有运算过程均在计算机内部进行，省去了作图，减少了人为因素的干扰，从而减少了运算误差，同时也节省了大量的时间。在 Excel 软件计算程序中，设计 A2:A100 和 B2:B100 分别为时间和其相对应的标示量的输入区，这样为使用者