药代动力学16 第十三章 曲线拟合以及常用的药代动力学软件.pdfVIP

第十三章 药动学数据的曲线拟合以及常用软件 计算药动学参数，是药代动力学进一步应用的基础。如何测定有关的动力 学参数呢？常用的方法是：首先在用药后的若干不同时间，采取血样（或尿样）， 测定其血药浓度值或尿中药量（这些数值称为实测值或观察值，用C 表示），这 i 样就有了药物浓度经时曲线数据；然后，依据半对数坐标图，选定一种模型方程 ) C （是时间t 的曲线函数）计算理论估算值(用 表示) ，按照观察值和理论估算值之 i 差的平方和（即残差平方和）或加权残差平方和(均用Re表示)最小的原则，采用 适当的算法，求出有关的动力学参数。这种方法，在数学上称为曲线拟合(fitting a curve) 。由于所采用的线性药代动力学的模型方程是多指数项之和的函数形式， 并且是所含动力学参数的非线性函数，所以这种曲线拟合方法称为非线性最小二 乘法。 一、最小二乘法的一般原理 设y 是变量x 的函数，含有m个待定参数a ，a ，…，a 。记为 1 2 m y ＝f （x ；a ，a ，…，a ） 1 2 m 若对x和y 作n次观察，测得观察值(x ，y ) ，(x ，y ) ，…，(x ，y ) 。根据这样 1 1 2 2 n n 一组二维数据，即平面上的若干点，要求确定这个一元函数y ＝f （x ；a ，a ，…， 1 2 am ）；(i=1，2 ，…，n) ，即一条曲线，使这些点与曲线总体来说尽量接近。并使y ) 的观察值y 与理论估算值y f （x ；a ，a ，…，a ）；(i=1，2 ，…，m) 的误差 i i i 1 2 m n n 平方和，即残差平方和 ) 2 2 取得最小值， R ∑(y =−y ) ∑(y =−f (x ,a ,a ,…,a )) e i i i i 1 2 m i 1 i 1 n n 或者加权残差平方和 ) 2 2 取得 R ∑(y =−y ) /w ∑(y =−f (x ,a ,a ,…,a )) /w e i i i