【www.ahwmw.com--组织人事】

　　在临床药物代谢动力学研究中，药动学参数的求解、血药浓度的计算、数据分布的设计、给药方案的设计及图表绘制、组织管理等，都要求进行大量且繁琐的数据处理及绘制各式图形，需要借助各种价格昂贵的专业配套软件才可完成。因此，本文向大家介绍一款常见的适合实际操作的办公数据处理软件，即利用VBA 编译的Excel 软件，进行药代动力学参数的计算，该方法具有简便、快速适合操作等优点。

　　1 材料

　　1. 1 计算机奔腾II 以上，内存256M 以上，硬盘3. 2G 以上。

　　1. 2 软件Windows 7 或Windows XP，Microsoft Excel 2003( MicrosoftCorporation) 或其他版本均可。

　　1. 3 数据及数据处理基本公式选自文献中口服某药物后在不同时间点测得的血药浓度时间数据。

　　2 Excel 的操作过程

　　2. 1 数据输入运行程序，选定工作表在A 2 ～ R2 栏中分别输入取样时间t /h、实测血药浓度c /( μg /ml) 、lnc、外推浓度Me - βt /( μg /ml) 、残数浓度cr 1 /( μg /ml) 、Le - at /( μg /ml) 、cr 2( μg /ml) 、lncr 2。在A3 ～ A8 栏中分别输入时间和血药浓度。A17 ～ J17分别输入β、M、t1 /2( β) 、α、L、Ka、N、t1 /2( α) 、K21。

　　2. 2 计算以文献数据为例，根据公式( 1-5) ，在Excel 工作台上，对lnc-t 作" X、Y 散点图"。选中A3 ～ A14的t /h 和C3 ～C14的lnc，单击菜单栏中的" 插入" 选择" 图表" ，在弹出的对话框中选择"X、Y 散点图" ，作出血药浓度与时间的曲线。根据公式( 1-5) ，结合中线性拟合出的公式，可以看出斜率即为- β，参数β = 0. 1081，则半衰期t1 /2( β) = 0. 693β = 6.41h。图中截距就是lnM，在单元格A18中输入公式" = EXP( 1.8338) " ，然后按" enter" 键，所输出的结果即为参数M 的值，M =6. 26。在已知β 和M 的情况下，可以求出外推浓度Meβt /( μg /ml) ，在单元格D3中输入公式" = 6. 26 * EXP ( - 0. 1081 *A3) " ，然后按" enter" 键即可显示结果，将鼠标光标放在单元格的右下角，当变成实心十字形状时选中单元格下拉至D14，利用电子表格所具有的自动填充功能可将余下数值计算出来。为求算α，可先调整c = Ne - Kat + Le - at + e- βt，得: c - M- βt =Ne - kat + Le - at，因为Kα > α，当t 足够大时，e - kat→0，故得到: ln( c - M- βt )= α2. 303 + ln Lc 为实测浓度，Me - βt /( μg /ml) 为外推浓度，( c - - βt ) 为残数浓度1( cr 1) ，则有:lncr 1 = - αt + ln L若要对lncr 1-t 作线性拟合图示，首先要知道cr 1，在Excel 工作台上，选中单元格E3，输入公式" = B3-D3" ，然后按" enter"键，结果便显示于单元格E3中，鼠标移至单元格E3的右下角，待光标变成实心十字的时候选中单元格下拉至E14，利用自动填充功能计算出余下数值。计算lncr 1，则在单元格F3中输入公式" = LN( E3) " ，按回车键，显示结果后，再用自动填充功能算出余下数值，但由于lncr 1中cr 1的数值不能为负数，所以填充的结果中有部分数值并不存在。在对lncr 1-t 作"X、Y 散点图" 时，与作出图2 的方法相同，选中F7 ～ F11的t /h 和C7 ～ C11的lncr 1，单击菜单栏中的" 插入"选择" 图表" ，在弹出的对话框中选择" X、Y 散点图" ，连续选择" 下一步" 到" 完成"。由图中的方程式可以清晰的看出斜率即是- α，则参数α =0. 3282，残数半衰期t1 /2( α) = 2. 11h，图中截距就是lnL，在单元格E18中输入公式" = EXP( 0. 8240) " ，然后按" enter" 键，所输出的结果即为参数L 的值，L = 2. 28。在已知α 和L 的情况下，可以求出Le - αt /( μg /ml) ，在单元格G3中输入公式" = 3. 23* EXP( - 0. 4* A3) " ，然后按" enter" 键即可显示结果，将鼠标光标放在单元格的右下角，当变成实心十字形状时选中单元格下拉至G14，利用电子表格所具有的自动填充功能可将余下数值计算出来。计算kα及N，t1 /2( α) ，为此先将公式调整为:c - Me - βt - Le - αt = Ne - kat ;因为c - Me - βt = cr 1，上式可进一步整理为:cr 1 - Le - αt = Ne - kat或Le - at - cr 1 = - Ne - kat ;令Le - at - cr 1 = cr 2，即残数浓度2，Le - at 为残数浓度1 上的外推浓度。故取对数得:lncr 2 = - ka t + ln( - N)由以上作图的经验，可以看出，若要求算ka，只需对lncr 2 - t作图，求出其斜率即可。但首先要求出cr 2和lncr 2。在单元格H3中输入公式" = G3-E3" ，按回车键显示结果，用自动填充功能算出余下数值。在单元格I3中输入公式" = LN( H3) " ，按"enter" 键显示出结果，下拉自动填充柄至I14算出其他数值，但由于lncr 2中的cr 2不能为负值，所以求算的结果当中有部分数值并不存在，将显示" #NUM! " 的字样。对lncr 2 - t 作"X、Y 散点图"。选中A3 ～ A8的t /h 和I3 ～I8的lncr 2，单击菜单栏中的" 插入" 选择" 图表" ，在弹出的对话框中选择" X、Y 散点图" ，连续选择" 下一步" 到" 完成"。

　　2. 3 结果在求算出所需全部参数之后，便可以求算出动力学方程，该实践所求得到的动力学方程如下:c = 2. 28e - 0. 33t + 6. 26e - 0. 11t - 18. 87e - 2. 07t按此式计算的理论值与实验值进行比较，说明用上式来表征此实验结果是合适的。

　　3 结论

　　本论文主要是对Excel 求算数据的实践过程进行研究，以药物动力学参数为主线，讨论电子表格在药代动力学参数中的应用，并对用Excel 做出的线性拟合的直线结果进行分析，与最小二乘法结合，优缺点互补，达到最理想最实用的效果。实践结果表明，利用最小二乘法与添加趋势线的拟合直线的方法可以大大简化数据的计算过程，精确率可达到98%，改善了传统的利用最小二乘法手工求算参数的方法，可替代专业性的药学软件实现数据求算，既简单快捷又经济有效，是一种适合药学初学者使用的简便方法。

本文来源：https://www.ahwmw.com/jianghuafayan/72280/