13. 二项式模型

在机器学习领域,应用最广的两个模型,一个是线性回归模型 另一个就是逻辑回归模型。 线性回归模型就是采用标准连接函数的高斯模型, 高斯模型是处理连续值数据的基本模型。 而逻辑回归模型是处理二分类数据的基本模型, 逻辑回归模型就是标准连接函数的二项式回归模型。 二项式回归模型对应的是指数族中的二项式分布, 二项式分布是统计学中最常见的概率分布之一,应用十分广泛。 本章我们讨论 GLM 框架下的二项式回归模型。

13.1. 伯努利分布

如果一个随机变量只有两种可能状态,就可以认为这个随机变量服从伯努利分布(Bernoulli distribution)。 比如,在广告场景中,用户点击广告的行为可以分成点击和不点击两个状态; 投掷一枚硬币,只能是正面向上或者反面向上(排除硬币站立的情况) 。服从伯努利分布的随机变量通常称为伯努利变量, 伯努利变量只有两个不同的状态,因此它是离散随机变量,伯努利分布属于离散概率分布。

通常会用数字 01 分别表示伯努利变量的两种状态, 假设状态为 1 的概率是 π ,那么状态为 0 的概率就是 1π 。伯努利概率分布的概率分布函数通常可以写成

(13.1.1)f(y;π)=πy(1π)1y

Y=1 的概率和 Y=0 的概率分别为

(13.1.2)P(Y=1)=f(y=1;π)=π1(1π)11=πP(Y=0)=f(y=0;π)=π0(1π)10=1π

伯努利分布的期望和方差分别是

(13.1.3)E[Y]=μ=πV(Y)=π(1π)=μ(1μ)

可以看出 π 其实就是分布的期望参数 μ ,并且伯努利分布的方差是期望的一个二次函数。 伯努利分布仅有一个期望参数, 因此它是一个单参数的概率分布 。

伯努利分布是离散变量的概率分布, 离散分布的概率分布函数称为 概率质量函数 ,概率质量函数的值直接就是概率值。 这一点和连续值分布是不同的,连续值分布的概率分布函数叫做 概率密度函数 ,概率密度函数的值并不是概率值,需要积分才能得到概率值。

13.2. 逻辑回归模型

13.2.1. 模型定义

假设响应变量 Y 是一个伯努利变量, 它的概率分布函数如 公式(13.1.1) 所示, 现在把它转化成指数族的形式

(13.2.1)f(y;π)=exp{yln(π1π)+ln(1π)}

由于参数 π 就是分布的期望 μ ,因此我们直接用符号 μ 替换 π

(13.2.2)f(y;μ)=exp{yln(μ1μ)+ln(1μ)}

和指数族的自然形式对比下,可以直接给出各项的内容。

(13.2.3)自然参数  θ=ln(μ1μ)累积函数  b(θ)=ln(1μ)分散函数  a(ϕ)=ϕ=1

它的期望可以通过累积函数的一阶导数求得

(13.2.4)E[Y]=b(θ)=μ

方差函数通过累积函数的二阶导数得到

(13.2.5)ν(μ)=b(θ)=μ(1μ)

分散函数和方差函数的乘积就是分布的方差

(13.2.6)V(Y)=a(ϕ)ν(μ)=μ(1μ)

可以看到伯努利分布的方差不是常量,而是关于期望参数 μ 的二次函数, 显然伯努利分布的方差会受到期望 μ 的影响。

根据标准连接函数的定义,标准连接函数是使得线性预测器 η 等于自然参数 θ 的连接函数,所以对于伯努利分布,其标准连接函数就是

(13.2.7)η=θ=g(μ)=ln(μ1μ)

在统计学中,这个函数称为 logit (/ˈloʊdʒɪt/ LOH-jit) 函数, 逻辑回归模型的标准连接函数就是 logit 函数,它的一阶导数为

(13.2.8)g(μ)=1μ(1μ)

响应函数是连接函数的反函数, logit 函数的反函数为

(13.2.9)logit(μ)1=eη1+eη=logistic(η)

logit 函数的反函数就是我们熟知 logistic 函数, logistic 函数中文叫做 逻辑函数 ,它是标准连接的伯努利回归模型的响应函数 ,因此一般把伯努利回归模型叫做 逻辑回归模型(logistic regression model)

(13.2.10)响应函数y^=μ^=r(η)=eη1+eη

备注

很多人把 logistic 函数称为 sigmoid,这是不准确的。 sigmoid 定义是:拥有S形状的一类函数。sigmoid 是一类函数的统称,并不是特指某一个函数, logistic 函数是 sigmoid 中的一例,其它的还有 Arctangent 函数、 Hyperbolic tangent 函数 、Gudermannian 函数等等。

最后整理下逻辑回归模型的关键组件

(13.2.11)标准连接函数:η=g(μ)=logit(μ)=lnμ1μ=ln(μ)ln(1μ)响应函数:μ=r(η)=eη1+eη方差函数:ν(μ)=μ(1μ)分散函数:a(ϕ)=1连接函数导数:g=1μ(1μ)

13.2.2. 参数估计

大部分有关逻辑回归模型的资料中,都是采用完全最大似然法估计模型的参数, 比如梯度下降法、牛顿法等等。 然而逻辑回归模型是可以纳入到 GLM 框架中的, 因此逻辑回归模型也是可以用 IRLS 算法进行的参数估计的。

逻辑模型的对数似然函数为

(13.2.12)(μ;y)=i=1N{yiln(μi1μi)+ln(1μi)}=i=1N{yilnμi+(1yi)ln(1μi)}

IRLS 算法中权重矩阵 W 和工作响应矩阵 Z 的计算公式分别为:

(13.2.13)Wii=1a(ϕ)ν(μ^i)(gi)2=μ^i(1μ^i)
(13.2.14)Zi=(yiμ^i)gi+ηi=(yiμ^i)μ^i(1μ^i)+ηi

偏差统计量为:

(13.2.15)D=2{(y;y)(μ^;y)}=2i=1N{yilnyiμ^i+(1yi)ln(1yi1μ^i)}

13.2.3. odds 与 logit

GLM 中连接函数的作用是把线性预测器 η 和响应变量的期望值 μ 连接在一起,本质上就是把 η 的空间和 μ 的空间进行映射, 并且这种映射关系必须是 双射 ,也就是对于任意一个 η 都可以得到一个唯一的 μ 与之对应,反过来也成立,对于任意一个 μ 都可以得到一个唯一的 η 与之对应, 这就要求连接函数必须是单调可逆的。 在很多有关 logistic 回归模型的资料中都会提到一个概念,odds , 为了令读者对二项式回归模型理解的更透传, 这里我们介绍下 odds 与标准连接函数 logit 的关系。

学过基础数学技能的人都知道,概率(probability)是用来描述事件发生的可能性的。 概率一般是通过频次来计算的,比如投掷一枚骰子 n 次,其中点数 1 的次数是 a ,那么点数 1 概率为

(13.2.16)p(1)=an

点数不是 1 的概率为

(13.2.17)p(¬1)=1an=nan

用概率来描述事件发生可能性是符合人的直觉的, 因此在日常生活中概率的应用是广泛的。 然而在统计学中,除了概率以外,还可以用 几率(odds) 来描述事件发生的可能性。在英语里,odds 的意思就是指几率、可能性。 odds 指的是 事件发生的概率不发生的概率 之比。

(13.2.18)odds=probability of eventprobability of no event=p1p

在上面的例子中,点数为 1odds

(13.2.19)odds(1)=a/n(na)/n=ana=frequency of eventfrequency of no event

可以看到事件总次数 n 是可以被抵消掉的, 因此 odds 也可以看做是频次之比。 由于 odds 是概率或者频次的比值,显然 odds 的取值范围是 [0,)odds 的值越大,事件发生的可能性就越大。 概率的值域范围是 [0,1] ,而 odds 的值域范围是 [0,) ,从概率到 odds 的转变,实现了值域的改变。

如果对 odds 取自然对数,就得到了 logit

(13.2.20)logit(odds)=ln(odds)=lnp1p

odds 的自然对数就称为 logitlogitlog-it 的简写。 odds 取自然对数后,输出值的范围就变成了 (,) ,正好和线性预测器 η 的值域范围变得一致了。 从概率到 logit 值域范围的演变过程为

(13.2.21)probability: [0,1]odds: [0,)logit: (,)

逻辑回归模型的期望值 μ 就表示一个概率值, 其取值范围是 [0,1] 。线性预测器 η 的取值范围是 (,)logit 作为逻辑(二项式)回归模型的标准连接函数, 其作用就是实现 μη 的映射。

13.3. 二项式分布

在英语语境中,会把随机变量的单次采样称为一次实验, 连续多次独立实验的结果形成的序列,称为一次 trial 。 如果把伯努利变量进行多次独立取样, 就得到一个伯努利状态序列, 如果把这个序列中状态为 1 的次数看做一个随机变量, 这个变量就是一个二项式(binomial)变量, 二项式变量的概率分布称为二项式分布(binomial distribution)。 二项式分布表示进行 n 次伯努利实验,状态 1 的次数的概率分布。 假设响应变量 Y 服从二项式分布, 则概率分布函数为

(13.3.1)f(y;n,π)=(ny)πy(1π)ny

其中符号 π 表示单次伯努利实验状态为 1 的概率,也就是伯努利分布的期望参数。 符号 n 是进行的实验次数,通常是已知的常量。 可以看出二项式分布的概率分布函数就是在伯努利概率分布函数的基础上加了组合数 (ny) ,之所以是组合数,而不是排列数,是因为二项式随机变量 Y 表示的次数,与顺序无关,因此是组合数。

二项式分布的期望和方差分别是

(13.3.2)E[Y]=μ=nπV(Y)=nπ(1π)

通常实验次数 n 是已知的常量,并不是未知参数, 二项式分布的未知参数只有 π ,这和伯努利分布是同一个参数。 事实上,伯努利分布是二项式分布的一个特例, 二项式分布中,当 n=1 时,就是退化成了伯努利分布。 因此,有些资料中,把伯努利分布也称作二项式分布,这是合理的。

13.4. 二项式回归模型

13.4.1. 模型定义

假设响应变量 Y 服从二项式分布,其概率分布函数为

(13.4.1)f(y;n,π)=(ny)πy(1π)ny

其中 π 是单次实验成功的概率,n 是实验的总次数, y 是成功的次数,这个分布函数中 π 是唯一的参数。 转化成自然指数族的形式为

(13.4.2)f(y;n,π)=exp{yln(π)+nln(1π)yln(1π)+ln(ny)}=exp{yln(π1π)+nln(1π)+ln(ny)}

和指数族的自然形式对比下,可以直接给出各项的内容。

(13.4.3)自然参数  θ=ln(π1π)累积函数  b(θ)=nln(1π)分散函数  a(ϕ)=ϕ=1

累积函数的一阶导数和二阶导数分别为

(13.4.4)b(θ)=bππθ=n1ππ(1π)=nπb(θ)=2bπ2(πθ)2+bπ2pθ2=n(1π)2(1π)2π2+n1π(1π)p(12π)=nπ2+nπ(12π)=nπ(1π)

通过累积函数的导数可以分别得到分布的期望和方差,

(13.4.5)E[Y]=μ=b(θ)=nπVar(Y)=a(ϕ)b(θ)=a(ϕ)ν(μ)=nπ(1π)=μ(1μn)

同样,二项式分布的方差是关于期望的一个函数,方差会受到期望的影响。

二项式分布的自然参数和伯努利分布的自然参数是完全一样, 因此二项式模型的标准连接函数也是 logit 函数。

(13.4.6)η=θ=g(μ)=ln(π1π)=ln(μnμ)

连接函数的导数为

(13.4.7)g(μ)=nμ(nμ)

二项式模型的响应函数同样也是 logistic 函数。

(13.4.8)y^=μ^=r(η)=neη1+eη=nπ^=nlogistic(η)

对比下伯努利回归模型与二项式回归模型, 可以看出,无论是连接函数还是响应函数,仅仅只是差了一个常量 n 而已,而 n 是一个已知的常量,并且当 n=1 时,二者就完全一样了。 连接函数都是映射的线性预测器 ηπ 的关系,这和 n 无关。 因此可以认为 伯努利模型和二项式模型是同一个模型,二项式回归模型也是逻辑回归模型

虽然严格来说二项式分布的期望是 μ=nπ ,但是由于 n 是已知常量,仅仅起到一个倍数的作用, 所以在强调 期望参数 时,可以只考虑 π , 而忽略 n 。这一点请铭记,否则在看某些资料时会迷惑!

最后,我们汇总下二项式回归(逻辑回归)模型的一些关键组件。

(13.4.9)标准连接函数:η=g(μ)=logit(μ)=lnμnμ=ln(μ)ln(nμ)响应函数:μ=r(η)=neη1+eη=nπ方差函数:ν(μ)=μ(1μn)分散函数:a(ϕ)=1连接函数导数:g=nμ(nμ)

13.4.2. 参数估计

二项式回归模型的对数似然函数为

(13.4.10)(μ^;y)=i=1N{yiln(π^i1π^i)+niln(1π^i)+ln(niyi)}=i=1N{yiln(μ^inμ^i)+niln(niμ^i)niln(ni)+ln(niyi)}

IRLS 算法中 WZ 的计算公式分别为

(13.4.11)Wii=1a(ϕ)ν(μ^i)(gi)2=niμ^i(niμ^i)
(13.4.12)Zi=(yiμ^i)gi+ηi=ni(yiμ^i)μ^i(niμ^i)+ηi

偏差统计量为

(13.4.13)D=2{(y;y)f(μ^;y)m}=2i=1N{yiln(yiniyi)+niln(niyi)yiln(μ^iniμ^i)niln(niμ^i)}=2i=1N{yilnyiyiln(niyi)+niln(niyi)yilnμ^i+yiln(niμ^i)niln(niμ^i)}=2i=1N{[yilnyiyilnμ^i][yiln(niyi)yiln(niμ^i)]+[niln(niyi)niln(niμ^i)]}=2i=1N{yilnyiμ^i+(niyi)ln(niyiniμ^i)}

我们用符号 oik 表示从样本中 观测(observed) 到的 k 个状态的次数, 比如 oi0 表示失败的次数 niyioi1 表示成功的次数 yi 。用符号 eik 表示模型 拟合(fitted) 的结果, 比如 ei0 表示模型预测的失败的次数 niμ^iei1 表示模型预测的成功的次数 μ^i 。 则偏差统计量可以简写为

(13.4.14)D=2i=1Nkoiklnoikeik

二项式模型的偏差统计量是不包含冗余参数的,比如分散参数 ϕ ,所以可以直接用它的渐近分布进行假设检验。

(13.4.15)Dχ2(Np)

注意对于 ni 比较小的数据,这个近似的效果会比较差。

二项式模型的皮尔逊卡方统计量为

(13.4.16)χ2=i=1N(yiμ^i)2μ^i(1μ^in)=i=1N(yiniπ^i)2niπ^i(1π^i)

二项式模型的卡方统计量 χ2 和偏差统计量 D 是近似相等,可以用泰勒级数进行证明,这里省略证明。 当 ni 比较小时,卡方统计量会比偏差更准确一些, 但是但 ni 非常小时,无论偏差统计量还是卡方统计量都不在准确。

13.5. 其它连接函数

我们已经很清楚连接函数的作用了, 它的的作用就是把线性预测器 η 和模型的期望参数 μ 进行可逆的映射 。二项式模型的期望参数 μ=π (注意,这里我们忽略 n ) 是一个概率值,它的合理取值范围是 [0,1] 。因此,能用来做二项式模型响应函数(连接函数的反函数)的函数, 它的输出域就必须是 [0,1] 。 而在统计学中,有一类函数是符合这个特点的,那就是概率分布的 累积分布函数 。 概率分布的累积分布函数,满足单调性,并且输出域是 [0,1] ,因此累积分布函数的反函数是可以作为二项式模型的连接函数的。

假设函数 f(s) 是一个概率密度(质量)函数,其累积分布函数 F(s) 就是对 f(s) 的积分。

(13.5.1)F(s)=tf(s)ds

其中 f(s)0,f(s)ds=1 ,累计分布函数表示的是随机变量 s 大于等于 t 的概率, P(st) , 显然,F(s) 的输出范围是 [0,1]

13.5.1. 恒等连接函数

首先我们看下均匀分布的累积分布函数, 假设概率分布函数 f(s) 是均匀分布的概率密度函数, 随机变量 s 的范围是 [c1,c2] ,均匀分布 f(s) 的概率密度函数为

(13.5.2)f(s)={1c2c1c1sc20otherwise

均匀分布的累积概率分布函数 F(x)

(13.5.3)F(x)=c1xf(s) ds=xc1c2c1for c1xc2

β1=c1c2c1β2=1c2c1 ,此时 F(x) 等价于

(13.5.4)F(x)=β1+β2x

F(x) 就是响应变量的期望 μ ,线性部分就是 η , 连接函数 g(μ) 就是恒等函数, 恒等连接函数就相当于是均匀分布的累积概率分布函数。

但是恒等连接函数有个限制,就是只有 x 在区间 [c1,c2] 才有意义,此时参数 β 也被约束在某个区间内。 然而,GLM 的通用参数估计算法 IRLS 并不能解决带约束的参数估计问题, 因此恒等连接函数在二项式模型中并不常用。

13.5.2. probit 回归

现在我们看下正态分布的累积分布函数。 假设概率分布 f(s) 是标准正态分布 N(0,1) ,其CDF又称为累积正态分布函数(cumulative normal distribution function,CNDF), 累积正态分布函数是对标准正态分布 N(0,1) 概率密度函数的积分。 习惯上用符号 ϕ(x) 表示标准正态分布的概率密度函数, 用符号 Φ 表示累积正态分布函数,注意不要和指数族的分散参数搞混,在这里不是尺度参数。 累积正态分布函数为

(13.5.5)Φ(x)=xϕ(x)dx=1σ2πexp[12(sμσ)2]ds=Φ(xμσ)

它的反函数通常称为 probit 函数,用符号 Φ1(x) 表示。

(13.5.6)probit=Φ1(x)

采用 probit 函数作为连接函数二项式回归模型称为 probit 回归模型, probitlogit 非常的类似,二者的(反)函数图像都是 S 型曲线, 并且以点 (0,0.5) 为对称点呈现对称结构,二者只是在曲率上稍微有些差别。

../../../_images/binomial_02.jpg

图 13.5.1 logitprobit 函数曲线对比。

对二值或分组的二项式数据使用 probit 回归模型通常会产生与逻辑回归相似的输出。 但是 logit 模型可以解释为胜率比(odds ratio),而 probit 没有这样的解释。 但是,如果线性关系中涉及正态性(通常在生物学领域中就是如此), 则 probit 可能是合适的模型。 当研究人员对赔率不感兴趣而对预测或分类感兴趣时,也可以使用它。 然后,如果 probit 模型的偏差显着低于相应的 logit 模型的偏差,则首选前者。 当然,比较二项式家族中的任何连接函数时,偏差小的模型永远是最优的选择。

我们可以非常容易的应用 IRLS 算法对 probit 模型进行参数估计, 只需要替换算法中的连接函数 g 、响应函数 r , 以及连接函数的导数 g

(13.5.7)连接函数 η=g(μ)=Φ1(μ)响应函数 μ=g1(η)=Φ(η)连接函数一阶导 g(μ)=ϕ(η)

probit 连接函数相比 logit 连接函数复杂了很多, 它不如 logit 更加的流行。 但是,当特征数据存在正态性时,probit 可能更合适。

13.5.3. log-log 和 clog-log

logitprobit 的曲线都是以点 (0,0.5) 对称的”S”型曲线, 所以,二分类 logitprobit 模型假定响应数据中为 01 的比例是相同的。 然而,当响应数据中 01 的比例相差巨大时, clog-loglog-log 可能会提供更好的模型效果, 因为它们具有非对称性,这些不对称的连接函数有时会更适合特殊的数据情况。

clog-loglog-log 是不对称的”S”型, 对于 clog-log 模型,S曲线的上部比 logitprobit 更大或更长; 而 log-log 模型恰好相反,S曲线的底部向左拉长或倾斜。

../../../_images/binomial_04.jpg

图 13.5.2 clog-loglog-log 函数

上图展示了 clog-loglog-log 连接函数的非对称性, 非对称结构使得我们可以更关注其中的一类。 假设在二分类场景中, y 表示样本是正类的概率, y¯=1y 表示样本是负类的概率。 如果我们使用对称的连接函数,则可以拟合一个同时适用于两个类别的模型, 两个类别的模型仅仅是系数的正负号不同。 此时, probit(y)probit(y¯) 是互补的, logit(y)logit(y¯) 是互补的。 但这在非对称的链接模型中,将不再适用, 在非对称连接函数 log-logclog-log 的情况下, loglog(y) 不再和 loglog(y¯) 互补, cloglog(y) 不再和 cloglog(y¯) 互补, 而是 loglog(y)cloglog(y¯) 互补。 下面用公式展示出各自的互补关系。

(13.5.8)probit(y¯)=Φ1y¯=Φ1(1y)=Φ1(y)=probit(y)logit(y¯)=ln(y¯1y¯)=ln(1y1(1y)))=ln(1yy))=ln(y1y))=logit(y)loglog(y¯)=ln{ln(y¯)}=ln{ln(1y)}=[ln{ln(1y)}]=cloglog(y)cloglog(y¯)=ln{ln(1y¯)}=ln[ln{1(1y)}]=[ln{ln(y)}]=loglog(y)

直观的来讲,就是在对称链接的情况下,有 r(η)+r(η)=1 , 而在非对称链接的情况下 r(η)+r(η)1clog-loglog-log 的关系是对称的, rloglog(η)+rcloglog(η)=1

log-log

(13.5.9)连接函数:η=g(μ)=ln[ln(μ)]响应函数:μ=r(η)=exp[exp(η)]连接函数导数:g=μlnμ

Clog-log

(13.5.10)连接函数:η=g(μ)=ln[ln(1μ)]响应函数:μ=r(η)=1exp[exp(η)]连接函数导数:g=(μ1)ln(1μ)

13.6. 分组数据与比例数据

伯努利模型对应着二分类的场景,也就是一条数据可以有0或1两个类别,此时响应变量的值 yi 就是类别, 模型预测每条数据所属的类别。 而二项式模型中,一条数据就表示一组实验结果,相比于伯努利数据,每条数据中多了一个表示实验次数的 ni, 响应变量的值 yi 不再是0或1的二分类值,而是表示 ni 次实验中成功的次数,其取值范围是 0yini

举个实际的例子说明下,假设有两个赌徒想预测一个篮球运动员的投篮命中率。 一个赌徒的做法是, 收集了这个球员 每次投篮时 的身体状态信息、天气状态、队员状态、对手状态等等信息,以及本次投篮行为的结果, 进球还是没进球。然后训练了一个 伯努利模型(二分类模型) 预测球员 单次投篮进球的概率 。 另一个赌徒的做法是,收集这个球员 每场比赛 的信息,这个球员在每场比赛中投了几次,进了几次,以及其它一些信息。 然后训练了一个 二项式模型 预测球员的 一场比赛中进球率,即在投篮n次的情况下进球几次。

在二项式模型下,一条数据表示一组实验的结果,并且多了一个表示试验次数的 ni(xi,ni,yi) , 此时样本数据称为分组数据(grouped data),一条数据相当于一个组。 然而有时数据中的 ni 是缺失的,并且 yi 不再是成功次数,而是成功的比例 yi=yi/ni ,此时数据样本变成 (xi,yi) , 这时称为比例(proportional)数据 。

有些时候数据中可能缺少实验次数、成功次数这样的数据,而仅有一个比例数据,即成功率, 响应变量 Y 的值是一个 [0,1] 的比例值。 此时,可以使用线性回归模型拟合这个比例值, 但是,线性回归模型的输出值可能会超出 [0,1] 的区间范围。 这时,我们可以先把 y 值转换一下 。

(13.6.1)ynew=logit(y)=ln(y1y)

经过 logit 函数转换后,新的 ynew 值就是实数域范围了,然后再应用线性回归模型。 但是注意 logit 函数不能处理0和1的样本值。