Excel 处理药代动力学数据的讨论.docVIP

Excel 处理药代动力学数据的讨论 　　在临床药物代谢动力学研究中，药动学参数的求解、血药浓度的计算、数据分布的设计、给药方案的设计及图表绘制、组织管理等，都要求进行大量且繁琐的数据处理及绘制各式图形，需要借助各种价格昂贵的专业配套软件才可完成。因此，本文向大家介绍一款常见的适合实际操作的办公数据处理软件，即利用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 /( mu;g /ml) 、lnc、外推浓度Me - beta;t /( mu;g /ml) 、残数浓度cr 1 /( mu;g /ml) 、Le - at /( mu;g /ml) 、cr 2( mu;g /ml) 、lncr 2。在A3 ～ A8 栏中分别输入时间和血药浓度。A17 ～ J17分别输入beta;、M、t1 /2( beta;) 、alpha;、L、Ka、N、t1 /2( alpha;) 、K21。 　　2. 2 计算以文献数据为例，根据公式( 1-5) ，在Excel 工作台上，对lnc-t 作” X、Y 散点图”。选中A3 ～ A14的t /h 和C3 ～C14的lnc，单击菜单栏中的” 插入” 选择” 图表” ，在弹出的对话框中选择”X、Y 散点图” ，作出血药浓度与时间的曲线。根据公式( 1-5) ，结合中线性拟合出的公式，可以看出斜率即为- beta;，参数beta; = 0. 1081，则半衰期t1 /2( beta;) = 0. 693beta; = 6.41h。图中截距就是lnM，在单元格A18中输入公式” = EXP( 1.8338) “ ，然后按” enter” 键，所输出的结果即为参数M 的值，M =6. 26。在已知beta; 和M 的情况下，可以求出外推浓度Mebeta;t /( mu;g /ml) ，在单元格D3中输入公式” = 6. 26 * EXP ( - 0. 1081 *A3) “ ，然后按” enter” 键即可显示结果，将鼠标光标放在单元格的右下角，当变成实心十字形状时选中单元格下拉至D14，利用电子表格所具有的自动填充功能可将余下数值计算出来。为求算alpha;，可先调整c = Ne - Kat + Le - at + e- beta;t，得: c - M- beta;t =Ne - kat + Le - at，因为Kalpha; gt; alpha;，当t 足够大时，e - katrarr;0，故得到: ln( c - M- beta;t )= alpha;2. 303 + ln Lc 为实测浓度，Me - beta;t /( mu;g /ml) 为外推浓度，( c - - beta;t ) 为残数浓度1( cr 1) ，则有:lncr 1 = - alpha;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 散点图” ，连续选择” 下一步” 到” 完成”。由图中的方程式可以清晰的看出斜率即是- alpha;，则参数alpha; =0. 3282，残数半衰期t1 /2( alpha;) = 2. 11h，图中截距就是lnL，在单元格E18中输入公式” = EXP( 0. 8240) “ ，然后按” enter” 键，所输出的结果即为参数L 的值，L = 2. 28。在已知alph