充分利用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@

　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单元格 …… 此处隐藏：3572字，全部文档内容请下载后查看。喜欢就下载吧 ……