曲线图论文剂量效应曲线的累积正态分布拟合

1、曲线图论文|剂量效应曲线的累积正态分布拟合    本文由写手联盟-论文代写代发表 收集整理发布，本站提供相关职称论文发表服务。相关专业提示：正态分布曲线 正态分布曲线面积表 正态分布钟形曲线 excel 正态分布曲线 正态分布曲线图 标准正态分布曲线 正态分布曲线的特点 正态分布累积 正态分布 累积函数 曲线拟合/ D$ |7 I& d9 i. Y$ e# j7 9 Y7 c1 n  ?% A* ?! f8 E摘要：1 R/ q+ 3 p6 C) X6 F  n3 q目的：研究S型剂量效应曲线的累积正态分布

2、函数拟合。方法：构建累积正态分布函数的适当表达形式，对作者实验室所获得的实验数据进行拟合，同时与常用的四参数方程拟合进行比较。结果：两种拟合方法具有高度的一致性。结论：正态分布模型是一种传统的数学模型，在此处的应用为实验数据拟合提供了另一个可能的选择，并为药效学S型剂量效应模型的研究分析提供了另一种思路。特别地，在常用于药效学S型剂量效应曲线拟合的数学模型中，其Hill系数被认为没有任何生理学含义，而累积正态分布函数中的相应参数却有明确的生理学含义，从数学上体现了药理与药效的内在联系。' P8 W9 W: V# W( 6 p- m1 L,   g1 A关键词：S型剂

3、量效应曲线，累积正态分布，四参数方程，一致性检验。- ?8 M5 v+ o( e% B9 K; a% p/ |+ l& S# T  T1 E前言：. K" m; w' e" T4 B" T( |7 F& n  K2 , y! w6 F& B: F剂量效应曲线在药效学研究中占有重要的地位。药物的作用强度一般随剂量增加而增大，但符合简单直线关系的较少，一般均为各种曲线关系。如以剂量的对数对药效作图，多呈S型曲线。在拟合过程中，因计算量方面的考虑，历史上总是想办法转化为直线进行处理，如将S型曲线的1

4、5%85%近似当作直线处理，或用各种数学方法将实验数据的曲线关系转化为直线关系，再利用对最小二乘法解偏微方程所获的直线拟合公式，得出拟合方程参数。随着计算机的普及和性能提高，计算量已不再作为一个主要的考虑因素，现在可根据最小二乘法的原理，用迭代法直接对残差平方和求极小值，避免复杂的数学变换所带来的误差。用于拟合S型曲线的函数方程有多种，但未见以累积正态分布函数进行拟合，现以细胞因子药效测定数据对此进行探讨。( B( T3 2 i+ q1 q$ A' ?" A& X4 a! a% & c1实验部份：( - z1 g& y. F$ Q/ u" ,

5、 q$ 8 u! f. O1.1材料和方法：; x# e0 u( a# w5 s3 U+ T0 x. % L6 g4 u' 1.1.1材料：G-CSF、IL-6活性测定实验数据，由本室提供。  U$ Z3 j" K( ?* g! G3 ( z: _5 b4 L4 Y1.1.2软件：MSEXCELXP。( c4 V* h- Y+ l$ X; o: m% C: A, I3 ! t9 ?6 3 - J+ y1.1.3方法：- z' v, g. l: W" + x( G& c, o+ o/ I7 rMTT法：细胞因子活性测定法参考文献1、

6、2、3、4。5 y7 U( j+ z; N$ U" K3 w& '' n( S# * 4 y9 i最小二乘法：以实验测定结果中的样品浓度（取对数）或稀释度单位为x，实测光密度值为y，组成一个实验点。将每点的x值代入相应的数学拟合模型，在假定的方程参数初值下得到y计算值，与各点的实测y值相减得残差。求得所有点的残差平方和，变动方程参数值使残差平方和最小。用MSEXCEL软件，以规划求解加载宏自动调整参数值，获得最终的拟合参数结果5。  N3 i% Q% y- R+ T+ V/ N: U" B' t6 J2 R) k 

7、; X! p一致性检验：拟合S型剂量效应曲线的目的在于获得相应拟合参数，用于评价被测药物的药效学特征。因此，比较本文中两种数学模型的拟合效果，可以用比较其获得的具有相同生理意义的拟合参数，评价其拟合效果的一致性。8 S8 K( v4 ; D: Z8 f4 L. M+ 7 f1.2建立数学模型：- Q/ C  v" * q* l3 m7 m: n6 u& Y8 h# W( s0 k" E1.2.1在药效实验中，药物剂量与药效的关系是一个重要的研究内容，如图1所示。其中，以药物剂量的对数对药效作图常呈S型，称为S型剂量效应曲线。+ ?% N

8、5 J# K: L% - x) _3 s1 G: o4 Z5 & n2 B8 U9 0 1.2.2在细胞因子类生物药品如G-CSF、IFN、IL-6、IL-11、GM-CSF、EPO等的生物学活性测定中，常采用MTT法。如以细胞因子浓度的对数（或者倍比稀释度）对光密度读数作图，均呈典型的S型曲线。内含微处理器或用外接微机处理数据的现代酶标仪，可以直接给出按四参数方程（式1）拟合得到的拟合参数，也可以按最小二乘法对实验数据进行拟合以获得拟合参数，拟合曲线如图2。四参数方程是一个广泛采用的S型曲线拟合模型，具有数据处理功能的酶标仪多数具备，本文中以其作为对照模型。, |) q, 3 / u

9、1 e) Y" e. H5 A6 _! H( b' 1 y" q1.2.3正态分布曲线及累积正态分布曲线如图3，其中（式2）为大家所熟知的正态分布方程，对（式2）进行积分得（式3），为累积正态分布方程。注意图2中的四参数方程曲线非常类似于如图3中的累积正态分布曲线。# N: R# L; A/ k6 7 o$ a% F3 U4 P' _2 E6 t. e2 D+ X  ! t1.2.4在细胞因子类生物药品活性的实际测定中，常以倍比稀释度（浓度呈对数单位递减）做横座标，以光密度读数做纵座标，绘制S型剂量效应曲线如图4，函数关系上与图2等价。,

10、 l$ " a+ ?" S; a6 Z. d& l5 B5 U; R在式（1）中，参数A、D分别代表最大（当浓度极大时）最小（当浓度极小时）吸收（药效指标），也称为曲线的上下平台；C为中点浓度，当样品浓度等于C时的光密度等于(A+D)/2，正处于曲线的拐点处，此点为中点稀释度，浓度为ED50（EC50、TC50、TD50、LC50、LD50）；B表示曲线的陡峭程度，称斜率因子，一般认为无明确的生理学含义。7 H. A! g& X3 y3 q. j5 d( f/ q+ $ Q3 Z; - j/ A& E2 W5 M, u5 l1.2.5在（式3）表示的

11、累积正态分布函数中，累积概率F(X)的值域为0100%，除正态分布的两个特征参数均值u、方差σ外，引入两个参数A、D，将值域转化为最大、最小吸收（药效），构建用于拟合S型剂量效应曲线的累积正态分布函数：: o. U" y  c5 R5 Q$ o) C3 g9 p% K! R$ a; E在笔者采用的MS-EXCEL软件中，函数NORMDIST(x,μ,σ,TRUE)可获得相应的累积正态分布值，因此函数（式4）可表达为y=NORMDIST(LN(x),μ,σ,TRUE)*(A-D)

12、+D。并且以为横座标，以y为纵座标，可得适用于拟合处理S型剂量效应曲线实验数据的累积正态分布函数模型。其中LN(x)x为药品浓度，对x进行对数化处理可用自然对数或以其它底数的对数；为药效指标。+ S( Q1 5 l6 X2 X. ': j# S- C# '9 W% I0 G1.2.6在细胞因子类生物药品的活性测定中，常以倍比稀释度做横座标，则细胞因子浓度随稀释度增加而降低，代表细胞因子药效作用的光密度读数也随稀释度增加而降低，此时用于拟合的函数（式4）应进行相应的变换同样，函数（式5）可表达为y=(1?NORMDIST(LN(x),μ,σ,TR

13、UE)*(A-D)+D，其中x为相应实验点的药品稀释倍数，获得如图4拟合细胞因子活性测定的累积正态分布曲线。在图4中，四参数方程曲线（式1）参数取值为A=0.8，D=0.1，C=0.01，B=-1.4，此时的中点稀释倍数为6.6439(稀释倍数为2)；以累积正态分布曲线（式5）对之拟合得参数值为A=0.7938，D=0.1066，μ=6.6420，σ=1.6846时，两曲线基本重合（12等距点残差平方和为9.4E-05），并且具有相同生理学含义的参数值也基本相等。5 - u! F4 ) B  n% _; T& Y1 H2 |

14、0; e2 '0 f# E' E" O7 '1.3拟合处理：3 d7 s, q( h, e0 a! : _3 I, '0 A* W" |+ w以作者实验室中的部分细胞因子类生物药品的活性测定实验数据，分别用四参数方程（式1）与累积正态分布（式5）两种拟合函数模型，采用MS-EXCEL软件按最小二乘法原理以前述方法进行拟合，获得各样品测试的拟合参数结果。- c8 M$   m$ E  " a& % / a6 C/ q6 z& B, M: m; P' B1.4拟合

15、效果比较：. C  h0 R" w* d" p* V! D: ?8 m: d2 J* : Y# 6 H& l残差平方和的最小值是评价拟合好坏的直接指标。从实际应用考虑，如上述两种拟合方法无差别，则获得的具有相同生理含义的相应拟合参数应相等，将多个样品测试结果的相应拟合参数成对作图应为一直线，以直线方程拟合，则有截距为0、斜率为1、相关系数为1的理论期望。因此，以四参数方程拟合的相应参数值为横座标，以累积正态分布拟合的相应参数值为纵座标，对相应的拟合参数结果作拟合效果一致性分析。( p6 F% u; Y! r$ H3 / * r 

16、0;J1 t+ d3 q" # N! 1.5结果：! Y! J, & F, ?/ P$ h8 x  e8 - W3 _  E" P, ( F7 s1.5.1对139例G-CSF活性测定实验数据进行分析，两拟合方法获得的具有相同生理含义的相应拟合参数及残差平方和的最小值极为接近，因篇幅原因不能给出全部的数据。考虑到中点稀释度是评价药效的指标，是最主要的参数，139例数据中的第77例样品该指标相对误差最大，按倍比稀释的12孔两平行列原始测定数据为：0.404、0.4、0.354、0.261、0.167、0.113、0.07、0.0

17、53、0.034、0.014、0.007、0和0.413、0.403、0.363、0.26、0.173、0.111、0.062、0.043、0.029、0.018、0.009、0；按（式1）拟合所得的最大吸收A0.446954、最小吸收D0.011734、斜率因子B-1.15479、中点浓度C0.045794换算为中点稀释度4.448687、残差平方和0.00112，按（式5）拟合所得的最大吸收A0.437945、最小吸收D0.015487、方差S2.003633、均值U或中点稀释度4.495227、残差平方和0.001474。表1给出全部139例所得的各项可比较指标，按上述一致性检验方法的分

18、析结果，均与理论期望非常接近。+ V# , |' F7 N! C. D+ m- N7 n5 V5 Z: H8 - / g6 S$ w1.5.2对47例IL-6活性测定实验数据进行分析，结果见表2，两拟合方法同样有很好的一致性。$ K- K( z& U& R6 : E1 ?0 / m1 ?' l7 4 a6 J2讨论：0 a3 . g7 I, r8 x- |8 a: B3 . M% A8 U( T. a2.1用于描述药物浓度与药效学关系的数学模型，一般有线性、对数线性、双曲线、S曲线等6、7、8。最常见的S型剂量效应曲线模型，其数学表达式为（式6），变换后可得（式

19、7）。( D; A% E& '$ c' h/ y( G0 T* t: q! g: t3 5 ?5 R可以看出，当E=y,c=x,E=D,E=A?D,EC=C,r=?B0max50时，（式7）则变为（式1），即两个不同表达的数学模型是完全等价的。换言之，用于细胞因子活性测定实验数据拟合，常见于现代酶标仪内置处理程序的四参数方程，等同用于药效学实验数据拟合的S型剂量效应曲线方程。特别是当1,0max0E=E=（以药效百分数表达）时，则方程更为简捷。( r. F+ G2 E8 t) a! c  d2 ' 8 C( 8 q9 u! j6 V 

20、; B2.2在药理学教科书中，还简单提及了质反应量效曲线如图5与统计分布的关系，文献7中也明确提及了质反应的“正态累积曲线”，但无进一步的介绍。本文中的数据及分析显示（式3）与（式1）的曲线极为相似，从而采用累积正态分布拟合实验数据；作者认为量反应与质反应的S型量效曲线均可能采用此法拟合，虽没有从数学上证明（式3）可完全替代（式1），也因数据来源的局限没有讨论适用的范围，但为拟合方程的选择提出了另一种可能，是否可用（式3）拟合其它情况下的S型量效曲线，应视具体的实验数据确定。# E% , B$ s9 I/ M: Q- I% M6 E- M1 _5 w2.3（式6）中的r称为Hill系

21、数，以及等价函数（式1）中的参数B，一般认为没有任何生理学含义，只表示曲线的陡峭程度。如采用本文介绍的累积正态分布函数（式3），对S型剂量效应曲线进行拟合，除相应参数与（式1）、（式6）中相应参数同样具有生理学含义外，特别是参数σ，代表正态分布方差，也就是测试药物在实验个体上药效分布的数学描述，具有明确的生理学含义。对照两种拟合方法，有理由相信，上述表示曲线陡峭程度的参数其生理学含义，同样是代表测试药物在实验个体上的药效分布，不同处在于它们并不能直接代表分布方差。/ $ : v: P5 e0 A6 c( X0 g7 5 8 . r  d2.4药效学实验数据采用（式3）拟合，则比用（式6）或等价的（式1）拟合可提供更多的有价值信息。如在细胞因子类生物药品的活性测定中，有可能通过有目的调整细胞培养条件，调整检测细胞对细胞因子敏感性的分布方差，从而调整拟合曲线的陡峭程度，以获得更好的