S 曲线
XiaO / 2023-10-29
S 曲线
S 曲线因其函数图像形状像字母 S 得名。S 型函数是有界、可微的实函数,在实数范围内均有取值,且导数恒为非负,有且只有一个拐点。
所有连续非负的凸形函数的积分都是 S 型函数,因此许多常见概率分布的累积分布函数会是 S 型函数,如逻辑函数,古德曼函数,误差函数,广义逻辑函数。
逻辑 S 曲线的的参数意义
斜率因子和不对称系数对 S 曲线的整体形状影响最大:
- 斜率因子控制曲线的陡度。斜率因子值越大,下渐近线和上渐近线之间的过渡越陡峭。斜率因子值越小,曲线越平缓。
- 不对称系数 (s) 控制曲线的倾斜度或偏移度。
- 如果对称因子 s=1,则曲线是完全对称的,即左右两侧的变化速度相同。
- 如果 s<1,称为左偏(left skewed),偏离左边。
- 如果 s>1,称为右偏(right skewed),偏离右边。
其他参数对曲线的影响如下:
- 下渐近线控制最小 y 值。
- 上渐近线控制最大 y 值。
- 拐点使整条曲线在 x 轴上向左或向右移动。
因此渐近线和拐点并不改变 S 曲线的形状。
S 曲线的不同表达式
1. Five-parameter logistic curve
$$ y\ =b\ +\ \frac{t-b}{\left(1+\left(2^{\frac{1}{s}}-1\right)\cdot\left(\frac{ES_{50}}{x}\right)^{h}\right)^{s}} $$
$$ y’\ =\ -{ES_{50}}^{h}hs\left(2^{\frac{1}{s}}-1\right)\left(b-t\right)x^{-\left(h+1\right)}\left(\frac{{ES_{50}}^{h}\left(2^{\frac{1}{s}}-1\right)}{x^{h}}+1\right)^{-\left(s+1\right)} $$
(1). 参数 t: 系统所能达到的最大效应值,与 y 的单位相同。
(2). 参数 b: 系统所产生的最低效应值,与 y 的单位相同。
(3). $ES_{50}$
: 刺激在达到其引发的效应从最小到最大的一半时的值,与 x 单位相同。
- 当我们谈论
$ES_{50}$
时,我们通常说的是刺激在达到其引发的效应从最小到最大的一半时的值。然而,这个"一半"并不总是指最大效应的一半,即 50%。它依赖于系统具体的上限与下限,也就是说,$ES_{50}$
所代表的"一半效应"对应于实际最大效应与最低效应的中间点的效应值。譬如,一个具体系统所能产生的最大效应为 40%,最低效应为 2%,此时的$ES_{50}$
对应的效应值其实是 21%。
(4). h: Hill’s Slope 曲线的陡峭程度,无单位。
- 在五参数 S 曲线中,Hill 斜率(h)是用来描述
$ES_{50}$
点(或曲线中点)附近响应变化趋势的一个参数,也就是说,它描述了在$ES_{50}$
点附近,效应随刺激变化的速度。Hill 斜率并不是$ES_{50}$
点的实际斜率,而是描述了在这个点附近,曲线形状或效应变化的一个度量。 - 进一步来说,曲线的 Hill 斜率反应了该系统对于外界刺激的敏感程度。如果 Hill 斜率较大(绝对值大),说明在
$ES_{50}$
点附近反应变化速度快,曲线在这一区间陡峭;相反,如果 Hill 斜率小(绝对值小),说明在$ES_{50}$
点附近反应变化速度慢,曲线在这一区间较平缓。
(5). s: symmetry parameter,对称因子,无单位,反映了 S 曲线的对称性。
2. Four-parameter logistic curve from USP
即上述对称因子等于 1 时的情况:
$$ y\ =b\ +\ \frac{t-b}{1+\left(\frac{ES_{50}}{x}\right)^{h}} $$
$$ y’=-\frac{{ES_{50}}^{h}h\left(b-t\right)x^{-\left(h+1\right)}}{\left(\frac{{ES_{50}}^{h}}{x^{h}}+1\right)^{2}} $$
$$ y’’\ =\ \frac{\left({ES_{50}}^{h}ht-{ES_{50}}^{h}bh\right)\left(hx^{\left(h-1\right)}\left(x^{\left(2h+1\right)}+2{ES_{50}}^{h}x^{\left(h+1\right)}+{ES_{50}}^{\left(2h\right)}x\right)-x^{h}\left(\left(2h+1\right)x^{2h}+2{ES_{50}}^{h}\left(h+1\right)x^{h}+{ES_{50}}^{2h}\right)\right)}{\left(x^{\left(2h+1\right)}+2{ES_{50}}^{h}x^{\left(h+1\right)}+{ES_{50}}^{\left(2h\right)}x\right)^{2}} $$
数学意义上的理解:
- 累积变化量对时间的一阶导数,就是累积变化量随时间的变化率,其物理意义就是速度; 累积变化量对时间的二阶导数,就是累积变化量随时间变化率随时间的变化率,也就是速度随时间的变化率,其物理意义就是加速度。加速度是由作用在物体上的外力和物体的质量决定的。
- S 曲线的拐点处,曲线斜率达到最大值,一阶导数(速度)有极值。
- 如果满足二阶可导条件,则该拐点二阶导数(加速度)为零。
3. Four-parameter logistic curve from 欧洲药典 (PhEur)
欧洲监管机构给出的公式明确将对数的底数 L 设置为 e (即它们总是使用自然对数)。该公式为:
$$
y=\delta+\frac{\alpha-\delta}{1+e^{-\beta(x-\gamma)}}
$$
这里有符号的明显变化。最显著的变化是使用 γ 表示 $ln(ES_{50})$
,并定义 $x = ln(z)$
。庆幸的是,对于其他参数来说都是直接替换。假设 β 为正:
上渐近线: $\alpha$
下渐近线: $\delta$
斜率参数: $-\beta$
$ln(ES_{50})$
: $\gamma$