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

　　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].中国药学杂志,2002,35(2):130-131.

　　[2] 张莉,夏运岳.用电子表格Excel计算药物溶出度Weibull分布参数[J].药学进展,2002,26(1):48-50.

　　[3] 邱家学.充分利用Excel数据处理功能巧算Weibull分布位置参数[J].药学进展,2004,28(10):464-466.

　　[4] 吴娟,黄蓓琳,王洪泉.Excel在溶出度试验数据处理中的应用[J].药学服务与研究,2002,12(2):311-313.

　　[5] 黄献,刘裕恒,莫志江.用SPSS拟合药物溶出度Weibull参数[J].中国药房,2006,17(14):1079-1081.

　　[6] 曹俊涵,郭晓波.用MATLAB曲线拟合工具箱计算药物溶出度Weibull分布参数[J].药学进展,2006,30(12):556-558.

　　[7] 陈幼亭.威布尔分布函数处理溶出数据应注意的问题[J].中国医院药学杂志,1998,18(9):419-420.

　　[8] 杨荣平,杨明,刘小彬.中药固体制剂及其溶出度的研究概况[J].世界科学技术:中医药现代化,2005,7(2):45-49,139.

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说医药类用SAS软件计算药物溶出度Weibull分布参数(2)在线全文阅读。