
第9章 双室模型.ppt
99页第九章第九章双双 室室 模模 型型n单单室室模模型型是是最最简简单单、、最最基基础础、、理理论论上上也也相相当当成成熟熟的的一一种种模模型型但但在在应应用用上上有有其其局局限限性性,,不不少少药药物物不符合单室模型动力学特征因为:不符合单室模型动力学特征因为:n体内各部分的血流都有一定的流速;体内各部分的血流都有一定的流速;n血血浆浆中中的的任任何何物物质质向向体体内内各各部部位位分分布布都都需需要要一一定定的时间n而药物是随血流进入各组织、器官和体液,而药物是随血流进入各组织、器官和体液,n因此,绝对符合单室模型的药物是不存在的因此,绝对符合单室模型的药物是不存在的 回顾n如果药物进入体内以后,只能很快进入机体的某如果药物进入体内以后,只能很快进入机体的某些部位,但很难较快地进入另一些部位,药物要些部位,但很难较快地进入另一些部位,药物要完全向这些部位分布,需要不容忽视的一段时间,完全向这些部位分布,需要不容忽视的一段时间,这时这时从速度论的观点将机体划分为药物分布均匀从速度论的观点将机体划分为药物分布均匀程度不同的两个独立系统程度不同的两个独立系统,即,即““双室模型双室模型””。
n双室模型由于双室模型由于““分布速度分布速度” ” 上的差别,将其分上的差别,将其分为两个隔室,即为两个隔室,即““中央室中央室” ” 和和““周边室周边室” ” Two Compartment Modeln中中央央室室::指指在在双双室室模模型型中中,,一一般般将将血血流流丰丰富富、、分分布布较较快快的的组组织织、、器器官官称称为为中中央央室室如如心心、、肝肝、、脾脾、、肺肺、、肾等n周周边边室室::指指在在双双室室模模型型中中,,一一般般将将血血液液供供应应较较少少、、药药物物分分布布缓缓慢慢的的组组织织、、器器官官称称为为周周边边室室如如肌肌肉肉、、骨骼、脂肪等骨骼、脂肪等Two Compartment Modeln双双室室模模型型在在药药物物动动力力学学分分析析中中应应用用的的很很广广,,它它的的理理论论也也相相当当成成熟熟,,随随着着电电子子计计算算机机以以及及测测定定手手段段的的不不断断精精确确化化,,双双室室模模型型有有了了较较深深入入的的研研究究,,并并能较精确地反映药物的体内过程能较精确地反映药物的体内过程n它它既既阐阐明明了了药药物物进进入入机机体体与与离离开开机机体体的的规规律律,,又又描述了药物在系统内各隔室之间的转运情况。
描述了药物在系统内各隔室之间的转运情况Two Compartment Modeln药药物物经经中中央央室室进进入入各各组组织织、、器器官官,,并并从从中中央央室室消除,消除,n在在中中央央室室与与周周边边室室之之间间药药物物进进行行着着可可逆逆性性的的转转运,运,n因因而而周周边边室室的的作作用用好好似似一一个个与与中中央央室室相相连连接接的的贮库n本章主要讨论药物从中央室消除的双室模型本章主要讨论药物从中央室消除的双室模型Two Compartment Model 本章主要内容n§1 §1 双室模型静脉注射给药双室模型静脉注射给药 n§2 §2 双室模型静脉滴注给药双室模型静脉滴注给药 n§3 §3 双室模型血管外给药双室模型血管外给药 第一节第一节双室模型静脉注射给药双室模型静脉注射给药Intravenous injection of two Compartment Model单单室室模模型型静静脉脉注注射射给给药药药药物物动动力力学学参参数数可采用以下两种方法可采用以下两种方法:(一)血药浓度法(一)血药浓度法(二)尿药浓度法(二)尿药浓度法 1 1 速率法速率法 2 2 亏量法亏量法一、血药浓度一、血药浓度法法n1. 1. 模型的建立模型的建立n2 2.血药浓度与时间函数关系式的确定.血药浓度与时间函数关系式的确定n3 3.求基本药物动力学参数.求基本药物动力学参数n4. 4. 求模型参数求模型参数n5. 5. 求其他药物动力学参数求其他药物动力学参数1 1、模型的建立、模型的建立 一般来说,双室模型药物其体内过程可这样设想:一般来说,双室模型药物其体内过程可这样设想: 药药物物经经静静脉脉注注射射后后,,进进入入中中央央室室,,然然后后再再逐逐渐渐向向周边室转运。
周边室转运 同时,周边室也把其中一部分药物返回中央室同时,周边室也把其中一部分药物返回中央室 另外,仅中央室内有一部分药物同时被消耗掉另外,仅中央室内有一部分药物同时被消耗掉 如图:如图: 药物动态变化包括三个部分:药物动态变化包括三个部分:①① 药物从中央室向周边室转运一部分(出);药物从中央室向周边室转运一部分(出);②② 药物从中央室消除一步分(出);药物从中央室消除一步分(出);③③ 药物从周边室向中央室返回一部分(进入)药物从周边室向中央室返回一部分(进入)对中央室对中央室药物动态变化包括两个部分:药物动态变化包括两个部分:①①药物从中央室向周边室转运一部分(进入);药物从中央室向周边室转运一部分(进入);②② 药物从周边室向中央室返回一部分(出)药物从周边室向中央室返回一部分(出)见模型图)(见模型图)对周边室对周边室静脉注射给药的模型图静脉注射给药的模型图双室双室模型模型X Xc c;;C C;;V Vc cX Xp p;;C Cp p;;V Vp pX X0 0ivk12k21k10中央室中央室周边室周边室 假设药物所有的转运过程均服从一级动力学过程,假设药物所有的转运过程均服从一级动力学过程,即药物的转运速率与该时刻体内药物量成正比,那即药物的转运速率与该时刻体内药物量成正比,那么,模型内各室药物转运可用下列微分方程组描述:么,模型内各室药物转运可用下列微分方程组描述:((1 1))((2 2))n初始条件:当初始条件:当t t==0 0时,时,n中央室药物量中央室药物量XcXc=X=X0 0,,周边室药物量周边室药物量X XP P=0=0,,药物全药物全部在中央室。
部在中央室n将(将(1 1)和()和(2 2)式联立得一阶线性齐次微分方程)式联立得一阶线性齐次微分方程组,在满足初始条件的情况下,采用拉氏变换法组,在满足初始条件的情况下,采用拉氏变换法或微积分代入法求解得中央室和周边室的药物量或微积分代入法求解得中央室和周边室的药物量与时间的关系结果如下与时间的关系结果如下:: 2 2、、血药浓度与时间的函数关系血药浓度与时间的函数关系中央室和周边室药物量如下:中央室和周边室药物量如下:((3 3))((4 4))中央室和周边室药物浓度如下:中央室和周边室药物浓度如下:((5 5))((6 6))3 3、求基本参数(、求基本参数(A A、、B B、、αα、、ββ))设则式中:式中:A A、、B B称为经验常数;称为经验常数;αα为分布相的速度常数或称为快配置系数;为分布相的速度常数或称为快配置系数;ββ称为消除相的混杂参数或称为慢配置系数称为消除相的混杂参数或称为慢配置系数7 7))((8 8))((9 9)) A A、、B B、、αα、、ββ是由模型参数(是由模型参数( k k1212、、k k2121、、k k1010 ))构成构成 的,的, αα、、ββ分别可用下式表示分别可用下式表示: :A A、、B B、、αα、、ββ与与k k1212、、k k2121、、k k1010的关系的关系((1010))((1111))A A、、B B、、αα、、ββ与与k k1212、、k k2121、、k k1010的关的关系系((1212))((1313)) 因因为为分分布布相相的的速速度度比比消消除除相相的的速速度度快快,,即即αα>>>>ββ,,当当时时间间t→∞t→∞((或或充充分分大大))时时,,e e-αt-αt <<<< e e- -βtβt,, e e-αt-αt则则与与e e-βt-βt相相比比是是相相对对小小量量,,可可以以忽忽略略不不计计。
即即,,当当t t充充分分大大时时,, e e-αt-αt→0→0,,而而e e-βt-βt仍仍保保持持一定的数值则(一定的数值则(9 9)式可化简为:)式可化简为:应用残数法求应用残数法求A A、、B B、、αα、、ββ((1414))n以以lgClgC1 1→t→t作图可得到一条直线,由直线的斜率作图可得到一条直线,由直线的斜率((b b1 1),),即可求出即可求出ββ,,由由ββ可求出消除相的生物可求出消除相的生物半衰期半衰期t t1/2(β)1/2(β) 两边取对数,得:两边取对数,得:n以以lgClgC1 1→t→t作图可得到一条直线,直线外推至与纵作图可得到一条直线,直线外推至与纵轴相交,得截距,由截距(轴相交,得截距,由截距(lgBlgB),),即可求出即可求出B Bn其其中中::C C1 1为为中中央央室室实实测测浓浓度度,, BeBe-βt-βt为为外外推推浓浓度度,,(( C C1 1 -Be -Be-βt-βt ))为残数浓度,设残数浓度为为残数浓度,设残数浓度为CrCr移项得移项得设设取对数得取对数得残数残数浓度浓度n以以lgCrlgCr→t→t作图亦为一条直线即残数线作图亦为一条直线即残数线, ,根据残数根据残数线的斜率线的斜率b b2 2和截距和截距a a2 2分别可求出分别可求出αα和和A A以及分布相以及分布相的生物半衰期的生物半衰期t t1/2(α )1/2(α )。
n实际工作中应借助于更精确的计算机程序法,直实际工作中应借助于更精确的计算机程序法,直接对血药浓度-时间数据回归分析求四个基本参接对血药浓度-时间数据回归分析求四个基本参数n在设计实验时应注意,在分布相时间内,若取样在设计实验时应注意,在分布相时间内,若取样太迟太少,可能看不到分布相,而将双室模型当太迟太少,可能看不到分布相,而将双室模型当成单室模型成单室模型n因为分布相速度很快,这一点在实验设计时必须因为分布相速度很快,这一点在实验设计时必须考虑小结小结4 4、求模型参数(、求模型参数(k k1212、、k k2121、、k k1010)) 因为:因为:当时间当时间t t==0 0时时所以:所以:其中:其中:X X0 0 -- --静脉注射给药剂量;静脉注射给药剂量; C C0 0 --t --t等于等于0 0时的浓度;时的浓度; VcVc -- --中央室的表观分布容积中央室的表观分布容积又因为:又因为:所以:所以:中央室中央室求求V V1 1∵∴∵求求k k2121n∴n∵n∴求求k k1212和和k k1010((1 1)) 血药浓度血药浓度- -时间曲线下面积(时间曲线下面积(AUCAUC))((2 2)) 总体清除率总体清除率((TBCLTBCL))((3 3)) 总表观分布容积(总表观分布容积(V Vββ ))((4 4)) 周边室表观分布容积(周边室表观分布容积(V V2 2))5 5、求其他药物动力学参数、求其他药物动力学参数((1 1)血药浓度)血药浓度- -时间曲线下面积时间曲线下面积积分法梯形法求梯形法求AUCAUC梯形法求静脉注射给药梯形法求静脉注射给药AUCAUC示意图示意图((2 2)总体清除率()总体清除率(TBCLTBCL)) n总体清除率等于单位时间内从体内清除药物的表总体清除率等于单位时间内从体内清除药物的表观分布容积数,即表示单位时间内流出的容积。
观分布容积数,即表示单位时间内流出的容积n式式中中V Vββ为为双双室室模模型型总总表表观观分分布布容容积积而而我我们们讨讨论论的的模模型型只只从从中中央央室室消消除除,,所所以以总总体体清清除除率率的的公式可以写成:公式可以写成:又因为又因为又因为又因为所以所以((3 3)总表观分布容积()总表观分布容积(V Vββ )) ((4 4)周边室表观分布容积()周边室表观分布容积(V V2 2)) n∵n∴n性肾排泄动力学中,原形药物的肾排泄速率性肾排泄动力学中,原形药物的肾排泄速率与中央室内的药量之间的关系成正比关系,用公与中央室内的药量之间的关系成正比关系,用公式表示为:式表示为:二、尿药排泄数据的处理二、尿药排泄数据的处理n 为为原原形形药药物物的的瞬瞬时时排排泄泄速速率率;;k ke e为为一一级级肾肾排排泄速度常数;泄速度常数;x x1 1为中央室任一时间的药量为中央室任一时间的药量中央室的药量与时间的关系如下:中央室的药量与时间的关系如下:两边同时求一阶导数得:两边同时求一阶导数得: 此式含义:双室模型静脉注射给药,尿中原形此式含义:双室模型静脉注射给药,尿中原形药物排泄的瞬时速率与时间的函数关系。
药物排泄的瞬时速率与时间的函数关系设则 对对志志愿愿者者而而言言,,尿尿中中原原形形药药物物的的瞬瞬时时速速率率很很难难求求出出,,实实际际工工作作中中常常用用某某段段集集尿尿期期间间的的平平均均速速率率对对该该段段时时间间的的中中点点时时间间作作图图,,其其半半对对数数图图形形为为一一条条二二项项指指数数曲曲线线,,即即用用△△XuXu/△t/△t代代替替dXudXu/ /dt dt ,,求求出出A’A’和和B’B’及及αα和和ββ值尿药速率法求药物动力学参数尿药速率法求药物动力学参数 所以代入一级肾排一级肾排泄速度常泄速度常数数 ∵ 求k21、k12、k10 与单室模型类似,双室模型药物也可采用亏量法,与单室模型类似,双室模型药物也可采用亏量法,根据尿中原形药物的量来求算药物动力学参数根据尿中原形药物的量来求算药物动力学参数 亏量法求药物动力学参数亏量法求药物动力学参数 n当当t→∞t→∞时,则尿中原型药物排泄的总量为:时,则尿中原型药物排泄的总量为:n则亏量:则亏量:含义:双室模型静脉注射给药尿中原型药物含义:双室模型静脉注射给药尿中原型药物的亏量与时间的函数关系的亏量与时间的函数关系设则采用残数法来求采用残数法来求A’’A’’、、B’’B’’、、αα、、ββ四个基本参数,四个基本参数,再求其他药物动力学参数。
再求其他药物动力学参数第二节第二节双室模型静脉滴注给药双室模型静脉滴注给药Intravenous injection of two Compartment Model血药浓度血药浓度法法n1. 1. 模型的建立模型的建立n2 2.血药浓度与时间函数关系式的确定.血药浓度与时间函数关系式的确定n3 3.求基本药物动力学参数.求基本药物动力学参数n4. 4. 求模型参数求模型参数n5. 5. 求其他药物动力学参数求其他药物动力学参数1 1、模型的建立、模型的建立中央室药物动态变化包括四个方面中央室药物动态变化包括四个方面::n药物从体外以恒速滴入,以补充中央室内的药药物从体外以恒速滴入,以补充中央室内的药物量;物量;n药物不断从中央室以药物不断从中央室以k k1212的速度向周边室转运;的速度向周边室转运;n药物不断从周边室以药物不断从周边室以k k2121的速度向中央室转运;的速度向中央室转运;n药物以药物以k k1010的速度从中央室消除的速度从中央室消除对中央室对中央室对周边室对周边室周边室内的药物的动态变化包括二个方面:周边室内的药物的动态变化包括二个方面:n药物以药物以k k1212的速度从中央室进入周边室;的速度从中央室进入周边室;n药物以药物以k k2121的速度从周边室返回中央室。
的速度从周边室返回中央室见下图静脉滴注给药的模型图静脉滴注给药的模型图双室双室模型模型X Xc c;;C C;;V Vc cX Xp p;;C Cp p;;V Vp p药物药物k0k12k21k10中央室中央室周边室周边室由模型图可得微分方程式如下:n中央室:中央室:n周边室:周边室: X Xc c、、X Xp p 分别为中央室和周边室的药量;分别为中央室和周边室的药量; V Vc c、、V Vp p分别为中央室和周边室的表观分布容积;分别为中央室和周边室的表观分布容积; k k12 12 为从周边室向中央室转运的一级速度常数;为从周边室向中央室转运的一级速度常数; k k21 21 为从中央室向周边室转运的一级速度常数;为从中央室向周边室转运的一级速度常数; k k10 10 为从中央室消除的一级速度常数;为从中央室消除的一级速度常数; k k0 0 为为零级静脉滴注速度常数零级静脉滴注速度常数Po.Po.符号的含义如下符号的含义如下::n对上述微分方程组,应用拉氏变化等方法可求得对上述微分方程组,应用拉氏变化等方法可求得中央室的药量及药物浓度与时间的函数关系。
中央室的药量及药物浓度与时间的函数关系2.2.血药浓度与时间的函数关系血药浓度与时间的函数关系n在满足模型的初始条件:当在满足模型的初始条件:当t=0t=0时,时,X Xc c= =X Xp p=0=0,,解上解上述微分方程组,得:述微分方程组,得:中央室药量中央室药量-时间时间中央室药物浓度中央室药物浓度-时间时间C=Xc/Vc12n双双室室模模型型静静脉脉滴滴注注给给药药中中央央室室的的药药量量((或或药物浓度)与时间的函数关系药物浓度)与时间的函数关系nk0为为零零级级静静脉脉滴滴注注速速度度常常数数,,α、、β分分别别为为分分布布相相和和消消除除相相的的混混杂杂参参数数Vc为为中中央央室室的的表表观观分分布布容容积积,,T为为静静脉脉滴注时间滴注时间n中央室是双室模型中的主要隔室,求动力学参中央室是双室模型中的主要隔室,求动力学参数主要以中央室为主数主要以中央室为主公式含义公式含义((1 1)滴注期间的血药浓度)滴注期间的血药浓度- -时间过程时间过程滴注期间:滴注期间:t=Tt=T32345((2 2)静滴停止后的血药浓度)静滴停止后的血药浓度- -时间过程时间过程滴注停止后:滴注停止后:t=T+t’t=T+t’623 3、求基本药物动力学参数、求基本药物动力学参数设n根据上式,按残数法求出α、β、R和S,然后再分别求其它参数如:Vc、k12、k21、k10和Vβ等。
第三节第三节双室模型血管外给药双室模型血管外给药Outside blood vessel administration of two Compartment ModelPo.Po.n1. 模型的建立模型的建立n2 2.血药浓度与时间函数关系式的确定.血药浓度与时间函数关系式的确定n3 3..残残数数法法求求基基本本药药物物动动力力学学参参数数((kaka、、αα和和β)β)n4 4..求其它药物动力学参数求其它药物动力学参数血药浓度法血药浓度法一、模型与函数关系式的建立一、模型与函数关系式的建立Po.Po.n血管外给药包括口服、肌肉注射或皮下注射,透皮血管外给药包括口服、肌肉注射或皮下注射,透皮给药,粘膜给药等给药,粘膜给药等n与血管内给药相比,血管外给药,药物有一个吸收与血管内给药相比,血管外给药,药物有一个吸收过程,药物过程,药物逐渐进入逐渐进入血液循环,而静脉给药,药物血液循环,而静脉给药,药物几乎几乎同时进入同时进入血液循环血液循环n血管外给药,药物的吸收和消除常用一级速率过程血管外给药,药物的吸收和消除常用一级速率过程描述,这种模型称之为一级吸收模型描述,这种模型称之为一级吸收模型。
1 1.模型的建立.模型的建立Po.Po.血管外给药的模型图血管外给药的模型图双室双室模型模型X Xc c;;C C;;V Vc cX Xp p;;C Cp p;;V Vp pX Xa akak12k21k10给药给药部位部位中央室中央室周边室周边室x0 X X0 0 、、X X 分别表示给药剂量和任意时刻给药部位的药量;分别表示给药剂量和任意时刻给药部位的药量; X Xc c、、X Xp p 分别表示中央室和周边室的药量;分别表示中央室和周边室的药量; C C、、C C p p 分别表示中央室和周边室的药物浓度分别表示中央室和周边室的药物浓度; ; V Vc c 、、V Vp p分别表示中央室和周边室的表观分布容积;分别表示中央室和周边室的表观分布容积; k k12 12 表示从周边室向中央室转运的一级速度常数;表示从周边室向中央室转运的一级速度常数; k k21 21 表示从中央室向周边室转运的一级速度常数;表示从中央室向周边室转运的一级速度常数; k k10 10 表示从中央室消除的一级速度常数;表示从中央室消除的一级速度常数; k ka a 表示一级吸收速度常数。
表示一级吸收速度常数Po.Po.符号的含义如下符号的含义如下::Po.Po.根据模型图建立微分方程如下根据模型图建立微分方程如下::对给药部位:对给药部位:对中央室:对中央室:对周边室:对周边室:n三式组成的方程组称为三式组成的方程组称为三维线性常数齐次微分三维线性常数齐次微分方程组方程组2.2.血药浓度与时间的函数关系血药浓度与时间的函数关系n在满足模型的初始条件:当在满足模型的初始条件:当t=0t=0时,时,X=XX=X0 0,,X X1 1=X=X2 2=0=0,,解上述微分方程组,得:解上述微分方程组,得:中央室药量中央室药量-时间时间中央室药物浓度中央室药物浓度-时间时间C=Xc/Vcn双双室室模模型型血血管管外外给给药药中中央央室室的的药药量量((或或药药物浓度)与时间的函数关系物浓度)与时间的函数关系nka为吸收相的以及消除速度常数,为吸收相的以及消除速度常数,nα、、β分别为分布相和消除相的混杂参数分别为分布相和消除相的混杂参数nVc为中央室的表观分布容积,为中央室的表观分布容积,nF为生物利用度为生物利用度n如下图:如下图:公式含义公式含义n从血药浓度-时间曲线可以看出,药物浓度先上升,后下降,最后平稳的减少,可将曲线分为a、b、c三个部分,即三个时相。
吸收相吸收相分布相分布相消除相消除相3 3、残数法估计基本参数、残数法估计基本参数 假设假设 n中央室是双室模型中的主要隔室,求动力学参中央室是双室模型中的主要隔室,求动力学参数主要以中央室为主数主要以中央室为主则中央室药物浓度则中央室药物浓度- -时间的函数关系式变为:时间的函数关系式变为:n因为:因为:kaka>>>>α α >>>> β βn根根据据上上式式,,应应用用两两次次残残数数法法,,即即可可求求出出基基本本参参数数N N、、L L、、M M、、kaka、、αα、、ββ基本公式如下:基本公式如下:求得求得M M和和ββ求得求得L L和和αα求得求得N N和和k ka ak ka a>>>>αα>>>>ββ,,当当t→∞t→∞时时k ka a>>>>αα当当t→∞t→∞时时见图见图3 3、、4 4图图3.3.双室双室模型模型血管血管外给外给药血药血药浓药浓度度- -时间时间图图图图4.4.双室双室模型模型血管血管外给外给药血药血药浓药浓度度- -时间时间半对半对数图数图末端末端直线直线第一条第一条残数线残数线第二条第二条残数线残数线第四节第四节模型的确定模型的确定Compartment Model前 言n经典药物动力学研究的理论基础是隔室经典药物动力学研究的理论基础是隔室模型理论,其计算结果的可靠性,首先模型理论,其计算结果的可靠性,首先依赖于模型的准确性。
在药物动力学研依赖于模型的准确性在药物动力学研究中,实验测定了不同时间的血药浓度,究中,实验测定了不同时间的血药浓度,然后进行数据处理然后进行数据处理n此时,首先遇到的问题是,此时,首先遇到的问题是,该药物属于该药物属于几室模型几室模型? ?是单室模型还是双室模型或其是单室模型还是双室模型或其它多室模型它多室模型?只有确定模型以后,才可?只有确定模型以后,才可能正确求出有关药物动力学参数,对该能正确求出有关药物动力学参数,对该药物的动力学特征作出正确的评价药物的动力学特征作出正确的评价前 言n一般情况下,可先用半对数纸作图,进一般情况下,可先用半对数纸作图,进行初步判断若所得为一直线,则可能行初步判断若所得为一直线,则可能是单室模型;若不是直线,则可能是多是单室模型;若不是直线,则可能是多室模型,要作进一步判断室模型,要作进一步判断前 言模型的判别方法模型的判别方法n一、作图法一、作图法n二、残差平方和(二、残差平方和(SMNSMN))法法n三、拟合度(三、拟合度(r r2 2))法法n四、四、AICAIC法法n五、五、F F检验法检验法一、作图法n例如,某药物静脉注射给药,于不同时间例如,某药物静脉注射给药,于不同时间测得血药浓度数据,见下表测得血药浓度数据,见下表1。
n 表表1 静脉注射给药血药浓度时间数据静脉注射给药血药浓度时间数据 t(h) 0.033 0.25 0.5 1.0 1.5 2.0 C(mg/ml) 7.10 5.80 5.40 4.00 3.40 2.95 t(h) 3.0 4.0 6.0 12.0 C(mg/ml) 2.75 2.20 1.90 1.56n根根据据表表中中的的数数据据,,以以血血药药浓浓度度的的对对数数((lgC)对时间(对时间(t))作图,结果见图作图,结果见图1n根据图根据图1 初步判断:该药不是单室模型,因初步判断:该药不是单室模型,因其半对数图不呈直线,因此可能是双室或其半对数图不呈直线,因此可能是双室或三室模型三室模型 二、残差平方和(SUM)判断n式中,式中,Ci为实测浓度;为实测浓度; 为根据设定的模为根据设定的模型计算的浓度型计算的浓度残差平方和公式为:残差平方和公式为:((1))n将上述数据按双室模型静脉注射给药模型将上述数据按双室模型静脉注射给药模型处理,得药物动力学方程如下:处理,得药物动力学方程如下:((2))根据残差平方和公式求得:根据残差平方和公式求得: SUM=0.2428按双室静注按双室静注处理处理按三室静注按三室静注处理处理n将上述数据按三室模型静脉注射给药模型将上述数据按三室模型静脉注射给药模型处理,得药物动力学方程如下:处理,得药物动力学方程如下:((3))根据残差平方和公式求得:根据残差平方和公式求得: SUM=0.4196n因为:因为:0.2428<<0.4196 对对同同一一组组实实验验数数据据,,求求出出SUM值值,,SUMSUM值值小小,,说说明明所所选选择择的的模模型型能能较较好好地地拟拟合合该该药药物的理论物的理论,判别该药物为双室模型。
判别该药物为双室模型结结论论n根据拟合度公式根据拟合度公式n按按双双室室模模型型求求得得拟拟合合度度值值为为:: = 0.998554;;n按三室模型求得拟合度值为:按三室模型求得拟合度值为: = 0.997501三、拟合度( )法进一步判断拟合度公式为拟合度公式为:n因为因为: 0.998554>>0.997501n所所以以,,根根据据对对同同一一组组实实验验数数据据求求出出的的拟拟合合度度值值大大,,说说明明所所选选择择模模型型能能较较好好地地拟拟合合该该药药物物的的理理论论,,判判别别该该药药物物为为双双室室模模型型这这与残差平方和法判断结果一致与残差平方和法判断结果一致结结论论四、四、AIC判别法判别法 n用残差平方和用残差平方和((SMC))法及拟合度(法及拟合度( ))法进行模型确定时,有时不能很好进行判断;法进行模型确定时,有时不能很好进行判断;近年来发展了一种新的判断线性药物动力学近年来发展了一种新的判断线性药物动力学模型的较好方法,即模型的较好方法,即AIC判别法,已被广泛判别法,已被广泛应用n AIC即即(Akaike,s Information Criterion)英英文缩写。
文缩写Akaike等定义等定义AIC为:为: AIC=NlnRe+2pn式中,式中,N为实验数据的个数;为实验数据的个数;p为所设模型为所设模型参数的个数;参数的个数;Re为权重残差平方和为权重残差平方和np和和Re的计算公式为的计算公式为(n为隔室数为隔室数):: p=2n ;式中,式中,Wi为权重系数,为权重系数,n通常取实测值平方的倒数,如果数据在高通常取实测值平方的倒数,如果数据在高浓度的准确性比低浓度大,则浓度的准确性比低浓度大,则Wi=1,,Re可以用可以用SUM计算;计算;Ci为实测浓度;为根据为实测浓度;为根据最小二乘法计算所得的估算浓度最小二乘法计算所得的估算浓度n根根据据不不同同模模型型的的AIC值值,, 可可以以确确定定最最佳佳的的模模型型 AICAIC值值愈愈小小,,则则认认为为该该模模型型拟拟合合愈愈好好特特别别是是当当两两种种模模型型的的残残差差平平方方和和值值相相近时,近时,AIC值较小的模型被认为较合适值较小的模型被认为较合适结论五、五、F检验法判别检验法判别nF检验检验(F test)法亦可用于模型判别,但需要法亦可用于模型判别,但需要查阅查阅F值表。
值表n式中,式中,Rel和和Re2分别由第分别由第1种和第种和第2种模型得到的种模型得到的加权残差平方和;加权残差平方和;df为自由度,即各自的实验总为自由度,即各自的实验总的数目的数目(采血次数采血次数)减去参数的数目减去参数的数目 例例如如::某某实实验验共共有有7 7个个实实验验点点,,若若按按单单室室参参数数的的数数目目为为2 2,,双双室室为为4 4,,则则dfdf分分别别为为5 5与与3 3F F值值显显著著性性可可与与F F值值表表中中的的自自由由度度为为( (dfdf1 1- -dfdf2 2) )及及dfdf2 2的的F F界值进行比较进行判定界值进行比较进行判定 n例如,某实验单室模型例如,某实验单室模型ReRel l为为0.004020.00402,自,自由度由度dfdf1 1 =5=5;;双室模型双室模型ReRe2 2为为0.0004470.000447,自,自由度由度dfdf2 2=3=3,,则:则:n查查F F表表,,相相应应自自由由度度的的F F界界值值( (概概率率5 5%%) )为为9.959.95,因此,双室模型更好因此,双室模型更好。












