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 /( μ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%，改善了傳統的利用最小二乘法手工求算參數的方法，可替代專業性的藥學軟件實現數據求算，既簡單快捷又經濟有效，是一種適合藥學初學者使用的簡便方法。

【Excel 處理藥代動力學數據的討論論文】相關文章：

如何建立和運用Excel數據清單01-22

EXCEL表格中數據透視表怎么用05-24

《EXCEL中數據的排序、篩選與分類匯總》教學設計06-05

excel函數07-30

Oracle數據庫無響應故障處理方式12-05

excel函數[優選]07-31

《初識Excel》說課稿11-14

Excel使用心得05-26

excel表格培訓心得06-22

常用的Excel函數公式11-23