对数正态模型及饮酒房室模型在血药 酒精浓度曲线模拟中的应用.docVIP

对数正态模型及饮酒房室模型在血药 酒精浓度曲线模拟中的应用 【摘要】 将对数正态模型和饮酒房室模型用于血药酒精浓度随时间变化曲线的模拟，对数正态模型和房室模型均能较好的模拟酒精血药浓度变化过程。配对 t 检验显示，对数正态模型与房室模型的模拟结果相一致;分别 计算 两种模拟结果的峰浓度和达峰时间，结果基本一致。 【关键词】 对数正态模型; 饮酒房室模型; 配对t检验 　　饮酒房室模型实指一房室模型［1］，主要用来 研究 血管外给药包括口服、肌内注射等多种给药方式的药动学过程。对数正态模型［2～3］常被用于急性心肌梗塞后心肌酶曲线的模拟，计算心肌酶的释放量，并以此估计梗塞面积。本研究将对数正态模型 应用 于血药酒精浓度的曲线模拟，并与一房室模型模拟结果进行比较。 　　 　　1数学模型提出 　　1.1饮酒房室模型［1］饮酒血液中酒精浓度C与时间T的关系如下:C= k0 Fx0V（ka -k） （e-kT-e-k aT）ka 为饮酒的一级吸收速率常数;k为酒精的消除速率常数;F为饮酒的生物利用度;x0 为初始的饮酒量;V为人体血液体积。 对于一个人体ka 、k、F、V均为常数，初始的饮酒量x0 也为常数，故模型可以变为:C=α（e-kT-e-k aT）其中α、k、ka 为待定参数。 　　1.2对数正态模型［2］ 将酒精浓度C与时间T表示为如下关系:C= b T ·EXP （ ln T-c）2-2d2 其中b、c、d为待定参数。 　　2曲线模拟 　　利用以上两个模型采用最小二乘法原理［4］对饮酒驾车的血药酒精浓度C （mg/ml ）随时间T（ h ）变化的数据进行曲线模拟，见表1。 　　表1 血药酒精浓度C随时间T变化数据（略） 　　2.1房的可决系数R2 分别为0.9829、 0.9884， 可决系数均大于0.95，拟合效果较好。 　　模拟图与可决系数显示，两个模型均能较好的模拟酒精血药浓度的变化过程。 　　3配对t检验 　　3.1利用模拟方程分别 计算 酒精浓度的估计值 　　表2 酒精浓度估计值（略） 　　3.2t检验考察两个模型模拟的结果有无显著性差异，将两个模型的估计值C1 、C2 进行配对t检验［6］。Di =C1 -C2 ;μd =∑n i=1Di /n H0 :μd =0;H1 :μd ≠0在假设H0 :μd =0成立下的统计量为: T= |μd -0| Sd / nt（n-1）Sd =∑n i=1D2 i - ∑n i=1Di 2 /n /（n（n-1））根据以上公式结合表2数据，计算得到假设成立下的统计量:T= |μd -0| Sd / n =0.3987由 SPSS 计算该统计量下的P=0.9959，按α=0.05的水平，Pα接受假设H0 ，即认为两个模型的模拟结果相一致，无显著性差异。 　　3.3峰浓度和达峰时间的计算峰浓度指饮酒后血液中酒精的最大浓度Cmax;达峰时间指酒精浓度达到最大时所经历的时间Tmax。利用拟合曲线方程结合软件分别计算出两个模型的峰浓度和达峰时间为:房室模型:Cmax=81.34Tmax=1.30对数正态模型:Cmax=83.71Tmax=1.26比较两个模型模拟结果的峰浓度和达峰时间基本相一致。 　　4结论 　　对数正态模型是统计学中 应用 的数学模型，本 研究 将其应用于酒精血药浓度变化的模拟，获得了较好的效果，与房室模型的结果相一致，比较两个模型的峰浓度和达峰时间也基本一致。因此，其可以应用于一般血管外给药的药动学的研究。 【 参考 文献 】 　　1 毛宗秀.高等数学.第3版. 北京:人民卫生出版社，1999，325～327. 　　2 张恒久，于忠卓，孙庆丰，等.推断急性心肌梗塞大小的一种数学 方法 .中华心血管杂志，1987，15（2）:128. 　　3 Sobel BE，Roberts，Larson KB.Estimation of infract size from ser-um MB creatine phosphokinase activity;applica-tions and limita-tions. Am J Cardiol， 1976，37:474. 　　4 姜启源，邢文训，谢金星，等. 大学数学实验. 北京:清华大学出版社，2005，2:318～319. 　　5 刘定远.医药数理统计方法.第3版.北京:人民卫生出版社，1999，200～201. 　　6 刘筱娴.医学统计学.北京: 科学 出版社，2000，10:43～50. 1