充分利用EXCEL数据处理功能巧算Weibull分布位置参数

交流沙龙

　464　2004,Vol.28,No.10　

药学进展

　2004年第28卷　第10期　第　464页　

ProgressinPharmaceuticalSciences

标准方面是根据中国药典2000年版二部附录对注射剂的要求,对该药品的性状、鉴别、溶液的澄清度与颜色、酸度、有关物质、水分、热原、无菌、含量等项指标进行了研究,采用HPLC法测定本品主药含量,并进行了专属性、精密度、回收率等方面的研究。结果显示,本法重现性好、专属性强,3批样品含量均在9010%～11010%规定限度内。而且该注射剂稳

定性好,室温贮存期长,有利于保管及临床使用。

[参考文献]

[1]许铁南.新喹诺酮类抗菌剂AM1155[J].国外医药抗生

素分册,1998,19(4):30223051

[2]中华人民共和国药典委员会.中华人民共和国药典(二

部附录)[M].北京:化学工业出版社,20001

交流沙龙

L数据处理功能巧算

Weibull分布位置参数

邱家学3

(中国药科大学商学院,江苏南京210009)

[摘　要]　目的:利用Excel的数据处理功能进行Weibull参数的确定。方法:利用Excel的模拟运算表、规划求解和

相关回归分析功能提出了两种计算方法进行药物溶出度Weibull分布参数确定中的难点———位置参数α的确定,从而实现对所有分布参数的确定。结果:本文所提两种方法的计算结果与文献结果一致。结论:本文提出的计算方法操作简便、易于实现。

[关键词]　Excel软件;药物溶出度;Weibull分布

[中图分类号]　TP391.13　　[文献标识码]　B　　[文章编号]　1001-5094(2004)01-0464-03

CalculatingtheParametersofWeibull’sDistributionbyExcelDataProcessing

QIUJia2xue3

(ChinaPharmaceuticalUniversity,Nanjing210009,China)

[Abstract]　Objective:TocalculatetheparametersofWeibull’sdistrbutionbyExcel.Methods:α,thepo2sitionparameterintheprocessofdeterminingtheparametersofWeibull’sdistribution,isconfirmedbythesimulativeoperationtable,programmingsolution,correlationandregressionofExcel,andthenalltheotherdistributionparametersarecalculated.Results:ThecalculationbyExcelisidenticalwiththatofliterature.Conclusion:ThemethodisconvenientandsuitableofdataprocessingindeterminationofWeibull’sdistribu2tionfordrugdissolutionrate.

[Keywords]　Excelsoftware,Drugdissolutionrate,Weibull’sdistribution

[接受日期]　2004205231

3

　通讯作者:　邱家学,副教授;研究方向:　医药产业经济;

Tel:025283271450;　E2mail:qjx@http://www.77cn.com.cn

　2004年第28卷　第10期　第　465页　

ProgressinPharmaceuticalSciences

药学进展

A518单元格。

交流沙龙

　2004,Vol.28,No.10　　465

对于药物释放的各种曲线,可以用有关方程进行拟合,求出符合该曲线的方程,从而实现对诸如

td、t50等参数的确定。常用的模型有零级释放、一级

114　在B19单元格内输入公式“=correl(D4:D9,

F4:F9)”,得α=0101时的相关系数01999718963。1.5　选择A19～B518区域,打开“数据”菜单,选择

释放、Higuchi方程和Weibull分布函数等[1]。许多文献对Weibull分布中参数的计算进行了讨论[2,3],尤其是文献[4]提出了借助Excel强大的计算功能进行Weibull分布参数的求解。实际上,Weibull分布

参数计算中最为关键之处在于位置参数α的确定,一旦位置参数α确定,其他参数(如形状参数m、尺度参数β及td、t50等)就可轻易计算出来,而文献[4]对位置参数α的确定并未作出明确的交代。本文利用Excel的模拟运算表、功能提出了两种确定位置参数α1　“模拟运算表”命令,在弹出的对话框中“输入引用,

列的单元格”处输入“＄A＄1”并”确定”,从而得出位置参数α取0101～4199时的各相关系数值。1.6　整体,,选择A20～B518区域,的“选择性粘确定”。对D20～E518的数,使最大的相关系数及其对应的位

“浮出”,从而完成位置参数α的确定。α

1.7　以求出的位置参数α之值回代到A1单元格,在D4～D9和F4～F9中的数据自动进行修改。以D4～D9和F4～F9分别作为回归分析的因变量和自变量,利用Excel提供的Intercept、Slope函数求出回归方程的截距、斜率。在A16、B16单元格内输入“=intercept(D4:D9,F4:F9)”、“=slope(D4:D9,F4:F9)”即可实现。1.8　对求出的截距、斜率进行变量变换,得到尺度B17单元格内输入参数β和形状参数m。在A17、“=EXP(2A16)”、“=B16”即可,从而完成了Weibull分布的参数确定。

119　求出Weibull分布参数后,在A18、B18单元格内输入公式“=a1-(-A173ln(1-01632))^B17”、

以文献[4]提的某药物在慢速搅拌

(30r/min)下于不同的溶出介质中的累计溶出百分比-时间数据为例进行计算。

具体计算过程如下。

111　在Excel的A3～F3单元格内分别输入“取样时

间(min)、累计溶出百分比F(t)、ln(1/(1-f(t)))、ln(ln(1/(1-F(t))))、)”t-α、ln(t-α字样作为数据的α所取的最小值,表头,并在A1单元格内输入0101(

实际上可以取(0,t1)中的任一数值),在A4～A9、B4～B9单元格内分别输入所测的时间、累计溶出百分比(以小数表示)数据。

表1　不同时间下药物累计溶出百分比F(t)(%)介质

5

0101mol/LHCL141990110mol/LHCL13156

10

15

20

25

30

35135521456515676125831992719437146441985115860199

时间t(min)

“=a1-(-A173ln(1-015))^B17”,可以方便地计

算出td、t50。至此,利用模拟运算表和相关回归分析功能求解位置参数α从而实现对Weibull分布有关参数的确定。

2　利用规划求解和相关回归分析功能计算

1.2　在C4～F4单元格内分别输入公式“=ln(1/(1-B4))”、“=ln(C4)”、“=A4-＄A＄1”、“=ln(E4)”,计算出第一个测量时间t1(本例为t1=5)和

基本思路是,为位置参数α任给一个初值(取

)和因变(0,t1)中的数值),可以确定自变量ln(t-α

量ln(ln(1/(1-F(t))))的数据,通过回归分析建立回归方程,将自变量诸值代入回归方程,求出因变量的估计值。利用“规划求解”求出使得残差(因变量的原始值与估计值之差)平方和最小时的位置参数α,此即为所求。在此基础上计算出Weibull分布各参数值。具体操作过程如下。2.1　数据输入及初步数据变换同111和112所述。

)、2.2　以ln(t-αln(ln(1/(1-F(t))))分别作为自

变量和因变量进行回归分析,在A20、A21单元格内

累计溶出百分比F1(本例F1=011499)所对应的相关数值,然后利用填充柄求出其他测量时间和累计溶出百分比所对应的数值。

1.3　在A20单元格内输入α所取的最小值(本例取0101),并以一定的步长值(本例取0101)、终止值(应

取t12步长值,本例取4199)按列方式产生填充序列。本例中即以0101、0102、……、4199分别占据A20～

交流沙龙

　466　2004,Vol.28,No.10　

药学进展

　2004年第28卷　第10期　第　466页　

ProgressinPharmaceuticalSciences

输入“=intercept(D4:D9,F4:F9)”、“=slope(D4:D9,F4:F9)”即可求出截距和斜率。2.3　在G4单元格输入公式“=F43＄A＄21+＄A＄20”求出因变量的一个估计值,并用填充柄填充G5～G9单元格求出因变量的其他估计值。2.4　用D4～D9与G4～G9对应相减的值作为残差

变单元格”,并“确定”,即可在A1单元格中得到所

需的位置参数α。2.6　对尺度参数β、t50的计算与形状参数m和td、117～119基本一致,不再赘述。3　结果与讨论

3.1　对表1的数据计算出的结果列于下表中,从中

放入H4～H9单元格中,并在I4单元格存放利用公

式“=SUMSQ(H4:H9)”计算出的平方和。

(如当前打开的2.5　打开“数据”,选择“规划求解”Excel无此功能,需通过“工具”中的“加载宏”里的

“规划求解”实现)。在“规划求解”弹出的对话框中,

,其与文献[4]的结

果也是相同的。。3.1,、操

选择I4为“目标单元格”,取“最小值”,A1作为方法

模拟运算表+规划求解+相关回归

0.01mol/LHCL

0.10mol/LHCL0.01mol/LHCL0.10mol/LHCL

r

α

0.861.660.858101.66300

β

35.6284218.8854535.6591618.86557

m1.236670.844281.236920.84428

tdt50

0.999960.998440.999960.99844

18.8364634.1189518.8367734.12279

14.2292522.6963114.2297822.69651

3.3　本文充分利用了Excel的模拟运算表、规划求

2002.42424251

[2]史文欣,张琦岩.片剂溶出度测定数据处理的电子计算

解和相关回归分析功能,解决了原先手工或编写计算机程序“艰难地”求解位置参数α的难题,并且操作简便,Excel又“随机”可见,信手拈来。

[参考文献]

[1]毕殿洲.药剂学(第4版)[M].北京:人民卫生出版社,

机程序[J].沈阳药学院学报,1987,4(1):5291

[3]邱家学.灰色系统理论在求解片剂溶出度Weibull分布

参数中的应用[J].中国药科大学学报,1988,19(1):172

211

[4]张莉,夏运岳.电子表格Excel计算药物溶出度Weibull

分布参数[J].药学进展,2002,26(1):482501

注射用舒巴坦钠与头孢他啶配伍的稳定性

尚　慧13,　景　茹2,　田　丹3,　石桂华4

(11山东省医学影像研究所,山东济南250021;21山东省立医院,山东济南250021;31济南市妇幼保健院,

山东济南250021;41山东济南第五人民医院,山东济南250021)

[摘　要]　目的:建立HPLC法同时测定舒巴坦钠和头孢他啶含量,研究注射用舒巴坦钠与头孢他啶配伍的稳定

性。方法:采用HPLC法测定配伍后0～6h内舒巴坦钠和头孢他啶的含量,同时观察不同温度下配伍液的外观、pH值。结果:配伍液中舒巴坦钠和头孢他啶的含量6h内均无明显变化,配伍液外观澄明,pH值基本稳定。结论:注

[接受日期]　2004205224

3

　通讯作者:　尚慧,主管药师;研究方向:　从事医院药物临床工作;

Tel:053127920656;　E2mail:liuyuan@http://www.77cn.com.cn

本文来源：https://www.bwwdw.com/article/e3ji.html