用SAS软件计算药物溶出度Weibull分布参数 【摘要】目的计算药物溶出度Weibull分布参数。方法采用SAS软件非线性回归NLIN 程序包中的Levenberg Marquardt法进行迭代计算。结果本方法计算的结果与文献的记录相吻合。结论本方法简便、快速、拟合结果准确,适用于药学工作者对药物溶出度Weibull 分布的计算。 【关键词】药物溶出度;SAS系统;Weibull分布 Calculation of Weibull’s distribution parameters for d rug dissolution with SAS PAN Yu shan1,LI Wen qing2,WANG Wei1, MA Su ying3,HU Gong zheng1 (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 Marquardt 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 da ta 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 计算结果 点击最大化按钮,全屏显示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]的结果非常接近。 本文来源：https://www.wddqw.com/doc/42d0686bd8ef5ef7ba0d4a7302768e9950e76e4a.html