在臨床藥物代謝動力學研究中，藥動學引數的求解、血藥濃度的計算、資料分佈的設計、給藥方案的設計及圖表繪製、組織管理等，都要求進行大量且繁瑣的資料處理及繪製各式圖形，需要藉助各種價格昂貴的專業配套軟體才可完成。因此，本文向大家介紹一款常見的適合實際操作的辦公資料處理軟體，即利用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%，改善了傳統的利用最小二乘法手工求算引數的方法，可替代專業性的藥學軟體實現資料求算，既簡單快捷又經濟有效，是一種適合藥學初學者使用的簡便方法。
