作者：潘玉善,李文清,王伟,马素英,胡功政
【摘要】 目的 计算 药物溶出度weibull分布参数。方法 采用sas软件非线性回归nlin程序包中的levenbergmarquardt法进行迭代计算。结果 本方法计算的结果与 文献 的记录相吻合。结论 本方法简便、快速、拟合结果准确，适用于药学工作者对药物溶出度weibull分布的计算。
【关键词】 药物溶出度;sas系统;weibull分布
abstract:objective to calculate the weibull′s distribution parameters for drug dissolution.methods the parameters of drug dissolution rate were treated with nonlinear regression by sas software and iterated by levenbergmarquardt algorithm (lma).results the calculation result by sas was identical with that in literature.conclusion the method was simple, convenient, calculated rapidly and the fitting result was accurate. it was suitable for data processing in calculation of weibull′s distribution for drug dissolution.
key words:drug dissolution; sas; weibull’s distribution
药物体外溶出度试验是一种控制药物制剂质量的体外检测方法，在对溶出试验数据进行处理时，其评价方法常用的有：对数曲线法、指数模型法、机率单位法、weibull法、gompertz法[1]，这些方法中又以weibull法最常用。WwW.11665.coM许多文献对weibull分布中参数的计算进行了探讨，目前对weibull分布参数的求解多采用excel软件计算[2～4]，虽然比较常用，但这种方法计算过程复杂而且计算得到的数据也不够直观，需要进一步的处理。此外，比较实用的计算机处理软件有spss、matlab[5，6]，本文提出采用sas系统(statistical analysis system)计算药物溶出度weibull分布参数时，方法简单，过程简练，计算速度快，为药物溶出度weibull分布参数的求解提供了一种准确、快速的方法。
1 数据资料
本文溶出数据来自文献[2]，它提供的某药物在慢速搅拌(30 r·min-1)的条件下，于不同的溶出介质中的累积溶出度-时间数据，见表1。
2 方法与结果
2.1 数据输入与编辑程序
启动sas(8.1)系统，打开工作界面，在显示的主窗口内，包含若干个子窗口，其中有三个最重要的窗口：editor程序编辑窗口，log运行记录窗口，output输出记录窗口。数据输入过程是在editor程序编辑窗口中进行。程序编辑窗口如图1。表1 不同时间下药物累积溶出百分比
程序说明：data语句产生临时数据集weibull(可任意命名，但不能超过八个字节)，表明数据步的开始;input语句指明读取变量t、f，@@表示读入一条观测值后不换行，连续读入数据;cards语句表明以下为数据行，数据行以下的“;”表示数据行结束;proc语句表明调用nlin(非线性回归)过程，选用的迭代计算法是marquardt法;model语句指定了所要拟合的非线性回归模型的weibull表达式，这里的a、b、m分别是文献[7]中的α、1·β-1及m，其中“**”表示乘方;parms语句设定weibull 参数的初始值，a=3，b=0.05，m=1;bounds语句表示weibull参数a的条件是大于或等于0;run语句表示程序结束。编辑完毕后，点击submit图标，提交编辑窗口中的程序或用菜单命令run→submit。
2.2 计算结果
点击最大化按钮，全屏显示output输出记录窗口，如图2，在此图中显示和记录了sas运行后的计算结果。
从output结果输出窗口中可以看出，经过5次迭代达到收敛标准，可以得到最优化解。在0.01 mol·l-1 hcl介质的参数为a=0.7895，b=0.0274(β=36.4564)，m=1.2424;通过相似的操作可以求得药物在0.10 mol·l-1 hcl的条件下，溶出度分布参数的最优化解：a=0.9028，b=0.0418(β=23.9234)，m=0.9109。
2.3 数据分析
sas软件计算结果除三个参数的求算外还有非线性最小二乘法的统计 总结 ，用于对模型的显著性检验和判断。在0.01 mol·l-1 hcl介质中，模型建立的决定系数r squared=1-residual ss/corrected ss=0.9999，说明建立的模型拟合优度令人满意。为求得特征参数td、t50，将上述结果代入weibull公式计算可得，笔者又编辑一个sas程序，如图3，能快速求得特征参数td=18.8750 min，将f值换成0.5即得t50=14.2582 min。在0.10 mol·l-1 hcl介质中，用同样的方法求得模型建立的决定系数r2=0.9949，td=33.5268 min，t50=22.7276 min，所得结果与文献[2，3，6]的结果非常接近。
3 讨论
药物溶出度是在规定溶剂中与一定条件下，从片剂或胶囊剂等固体制剂中溶出的速度和程度。对某些药品尤其是一些水溶性差的药物、缓释制剂，体外溶出度在某种程度上与体内生物利用度呈现一定的相关性，因此可以用体外溶出度试验估计其在体内药物动力学和药物生物利用度的特性[8]。weibull分布函数是处理药物制剂体外溶出度数据最常用的一种方法。为获得准确、可靠的药物制剂溶出参数，除试验条件正确、可靠以外，正确理解weibull分布参数的意义，快速、准确 计算 weibull分布参数和药物溶出度的特征参数也是十分重要的。weibull分布的函数方程为f(t)=1-e-(t-α)**m/β，包含3个未知参数：α：位置参数，m：形状参数，β：尺度参数[7]。 通过建立的weibull模型可以求得特征参数td(药物溶出63.2%所需的时间)，有时也 会计 算t50(药物溶出50%所需的时间)等。在一些 文献 [2～4，7]中提到位置参数a的选取是一个复杂的问题，给参数的计算带来一定的困难。本文介绍的方法通过编写程序，调用非线性回归过程marquardt法进行迭代运算，使用bounds限定语句表示weibull分布中位置参数a的条件是大于或等于0，可以一次性求得weibull分布的3个未知参数和模型的方差分析。此外，可以利用 weibull分布参数的计算结果，编写sas程序快速计算特征参数td、t50的值。本文编写的程序能以模板形式保存，在进行类似的运算时只需调出程序，更改数据步里的相应试验数据即可求解weibull分布参数和药物溶出度的特征参数。
综上所述，本文介绍的sas系统求解药物溶出度weibull分布参数和td、t50等值的方法，编程语言近似于 自然 语言，简单易学，使用者可以将自己的试验数据替换程序中的数据即可进行计算。本方法是建立在数据统计分析的基础上进行数据处理，过程简单、快速、精确度高，是求解药物溶出度weibull分布参数的简便、实用方法，适用于药学工作者对药物溶出度weibull分布的快速计算。
【 参考 文献】
[1] 夏锦辉,刘昌孝.固体药物制剂的体外溶出度的统计学评价分析[j].
