用SAS软件计算药物溶出度Weibull分布参数

【摘要】 目的 计算药物溶出度Weibull分布参数。方法 采用SAS软件非线性回归NLIN程序包中的Levenberg睲arquardt法进行迭代计算。结果 本方法计算的结果与文献的记录相吻合。结论 本方法简便、快速、拟合结果准确，适用于药学工作者对药物溶出度Weibull分布的计算。

【关键词】 药物溶出度;SAS系统;Weibull分布

Calculation of Weibull’s distribution parameters for drug dissolution with SAS

　　PAN Yu瞫han1,LI Wen瞦ing2,WANG Wei1, MA Su瞴ing3，HU Gong瞶heng1

　　(Henan Agricultural University 1.College of Animal Husbandry and Veterinary Sciences; 2.College of Life Sciences, Zhengzhou，Henan 450002, China; 3.College of Pharmacy, Henan Xinxiang Medical College, Xinxiang，Henan 453003, China)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 Levenberg睲arquardt 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法最常用。许多文献对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 计算结果