【自然语言处理】最大熵马尔可夫模型
有任何的书写错误、排版错误、概念错误等,希望大家包含指正。
由于这部分的参考资料比较少,网上大部分资料重复且不完整,对于一些关键计算没有推导,所以这里我主要讨论几篇论文和讲义。但是这些论文和讲义之间也有些许差别,讨论的过程中我会加入自己的理解,难免存在错误,欢迎大家讨论。
最大熵马尔可夫模型
最大熵马尔可夫模型(maximum-entropy Markov model,MEMM)又称为条件马尔可夫模型(conditional Markov model,CMM)。单纯顾名思义的话,可能会认为最大熵马尔可夫模型是最大熵模型与马尔可夫模型的融合,但其实,它结合了最大熵模型和隐马尔可夫模型(HMM)的共同特点,被广泛应用于序列标注问题。
模型介绍与对比
我们的讨论是以观测与状态是一对一关系为前提,而不考虑多个观测对应一个状态的情况,比如:Yellow River
被认为是一个专有名词,而不是一个形容词和一个名词。
先对比 MEMM 与 HMM。二者的概率图如下。
图 1 HMM (上)和满足一阶马尔可夫性质的 MEMM (下)
图 111 中的中的 MEMM 概率图为满足一阶马尔可夫性质的 MEMM,二阶 MEMM 的概率图只需要在此基础上加上从 st−2s_{t-2}st−2 到 sts_{t}st 的箭头。
对比 HMM 和 MEMM 的概率图,HMM 认为观测作为输出值与当前时刻的状态有关,而 MEMM 认为当前时刻的状态由之前时刻的状态和全部观测(或者多个时刻的观测值)共同确定,也就是说,在 MEMM 中观测不再独立,甚至作为输入值用于确定每个时刻的状态。这种性质决定了 HMM 是生成式模型,而 MEMM 是判别式模型。生成式模型对概率图表示的联合概率进行建模,再利用联合概率计算后验概率进行解码等操作;判别模型直接对概率图表示的条件概率建模,避免了对联合概率建模,使得解码等操作更加直接。
根据 HMM 的概率图,我们知道 HMM 直接对联合概率 P(S,O)P(S, O)P(S,O) 建模,即
P(S,O)=∏t=1TP(st∣st−1)P(ot∣st)P(S, O) = \prod_{t=1}^T P(s_t\mid s_{t-1}) P(o_t\mid s_t) P(S,O)=t=1∏TP(st∣st−1)P(ot∣st)
为了公式简洁、形式统一,不再特殊考虑 t=1t=1t=1 时的初始概率分布。
根据 MEMM 的概率图,我们对条件概率建模,有
P(S∣O)=∏t=1TP(st∣st−1,O)(1)P(S\mid O) = \prod_{t=1}^T P(s_t\mid s_{t-1}, O) \tag{1} P(S∣O)=t=1∏TP(st∣st−1,O)(1)
也有部分参考资料会进一步简化模型,将 HMM 概率图中状态到观测的箭头反向后的概率图作为 MEMM 的概率图,这相当于他们认为每个时刻的状态 sts_tst 不是与整个观测序列 OOO 有关,而是与当前时刻的观测值 oto_tot 有关。这会影响模型效果,但不会影响我们的推导和理解。
另外,对于每个时刻的状态,MEMM 考虑的是观测序列和之前时刻的状态,而 HMM 仅考虑之前时刻的状态。在序列标注问题中,将文本上下文包含的信息(特征)作为解码输出状态的判断依据显然更加合理,所以 MEMM 考虑整个观测序列是必然选择。如果将整个观测序列视为输入,那么输入值很可能无法完整枚举,因此需要大量的特征来刻画观测序列。如在文本中识别一个未见的公司名字时,除了传统的单词识别方法以外,还需要用到很多特征信息,如大写字母、结尾词、词性、格式、在文本中的位置等信息。也就是说,我们无法直接使用大量不同的观测序列,故需要提取其中的上下文信息(特征)来代表输入的观测序列。HMM 不适合处理用很多特征描述观测序列的情况,因为 HMM 在建模时不会将整个观测序列视为输入值,也就无法提取充足的特征。为此,MEMM 直接采用条件概率模型,从而使观测序列可以用特征表示,借助最大熵框架进行特征选取。
注意,对于某个时刻状态的解码一般不以整个观测序列作为输入,而是选取该时刻周围若干时刻的观测作为输入。之所以上文以“整个观测序列”进行叙述,其一是为了强调 MEMM 相比于 HMM 更加注重整体、上下文,其二是因为如此叙述会使得语言更加简洁。
多篇论文中常提到的一点是,HMM 需要根据状态输出观测,但是直观上,一个状态输出(影响)多个观测是很难实现的,比如,根据某个位置单词的词性判断周围的几个单词具体是什么,显然不合理,这也是 HMM 的观测独立假设的由来之一;反观 MEMM,将观测序列作为输入,可以很好地利用多个相互作用的特征和观测与观测之间的长期依赖关系,比如在定语从句 This zoo is not the only place where the animal can be seen.
中,where
需要根据先行词 zoo
确定,这种依赖关系跨越了半个句子,是 HMM 无法处理的,而 MEMM 直接输入整个句子作为观测,这使 MEMM 可以处理更复杂的情况。
对比 MEMM 与最大熵模型。最大熵模型没有考虑时序关系,以观测序列为条件,每个词是单独进行判别的,状态之间的关系无法得到充分利用,但是最大熵模型的优点在于可以使用任意的复杂特征。具有时序关系的 HMM 可以建立状态之间的马尔可夫性,这是最大熵模型所没有的。一个很自然的想法就是将两者的优势结合起来,这就得到了最大熵马尔可夫模型。
最大熵模型是一种特殊的逻辑回归模型(softmax),所以在有些文章中会看到。
学习问题
文献 [1 McCallum et al. ,2000] 与文献 [3 Michael Collins 讲义] 有些许不同,这里分开讨论。
基本假设,可能的状态数为 NNN,所有可能的状态的集合为 Q={q1,q2,…,qN}Q = \{q_1,q_2,\dots, q_N\}Q={q1,q2,…,qN}。
先讨论文献 [3]。假设训练集中的样本个数为 DDD,样本的序列长度均为 TTT,第 ddd 个样本的状态序列 Sd=(s1(d),s2(d),…,sT(d))S_d=(s^{(d)}_1,s_2^{(d)},\dots, s^{(d)}_T)Sd=(s1(d),s2(d),…,sT(d)),观测序列为 Od=(o1(d),o2(d),…,oT(d))O_d = (o^{(d)}_1,o^{(d)}_2,\dots, o^{(d)}_T)Od=(o1(d),o2(d),…,oT(d))。对于第 ddd 个样本的每个时刻 ttt,可以定义一个三元组
ht(d)=<Sd,Od,t>h_t^{(d)} =<\mathcal S_d,\mathcal O_d,t> ht(d)=<Sd,Od,t>
其中 Sd\mathcal S_dSd 为状态序列 SdS_dSd 的子集,在文献 [3] 中,由于讨论的是二阶马尔可夫性的 MEMM,所以 Sd={st−2(d),st−1(d)}\mathcal S_d = \{s^{(d)}_{t-2},s^{(d)}_{t-1}\}Sd={st−2(d),st−1(d)};Od\mathcal O_dOd 为观测序列 OdO_dOd 的子集,在文献 [3] 中 Od=Od\mathcal O_d = O_dOd=Od;ttt 表示序列的位置。可见,hth_tht 捕获 sts_tst 的条件信息,可以将 hth_tht 与 HMM 中的 oto_tot 类比理解,只不过在 MEMM 中不再简单地以某个时刻的观测值作为条件,而是以一组信息作为条件。另外,因为本文为了统一,仅讨论一阶 MEMM,故 Sd={st−1(d)}\mathcal S_d = \{s_{t-1}^{(d)}\}Sd={st−1(d)};至于 Od\mathcal O_dOd ,在理论推导中是取 OdO_dOd 的真子集还是与之相同都无所谓。
这样一来,特征函数 fff 会将 hhh 和 sss 一并考虑,即
f(h,s)={1,if(h,s)andmeetcertainconditions0,otherwisef (h,s) = \left\{ \begin{array}{ll} 1, & {\textbf {if}} \space\space (h,s) \space\space \textbf{and} \space\space \textbf{meet}\space\space \textbf{certain} \space\space \textbf{conditions}\\ 0, & \textbf{otherwise}\\ \end{array} \right. f(h,s)={1,0,if (h,s) and meet certain conditionsotherwise
式 (1)(1)(1) 可以被写为
P(S∣O)=∏t=1TP(st∣ht)(2)P(S\mid O) = \prod_{t=1}^T P(s_t\mid h_t) \tag{2} P(S∣O)=t=1∏TP(st∣ht)(2)
根据最大熵模型,其中 P(st∣ht)P(s_t\mid h_t)P(st∣ht) 可以表示为
P(st∣ht;w)=exp(w⋅f(ht,st))∑i=1Nexp(w⋅f(ht,s=qi))(3)P(s_t\mid h_t;w) = \frac{\exp\left(w·f(h_t,s_t)\right)}{\sum\limits_{i=1}^N\exp\left( w·f(h_t, s=q_i) \right)} \tag{3} P(st∣ht;w)=i=1∑Nexp(w⋅f(ht,s=qi))exp(w⋅f(ht,st))(3)
其中 www 和 fff 均为向量,w⋅fw·fw⋅f 表示内积运算。
文献 [3] 中直接给出了采用正则化的对数似然函数来学习模型参数 www,
L(w)=log∏d=1D∏t=1TP(st(d)∣ht(d);w)−λ2∥w∥22=∑d=1D∑t=1TlogP(st(d)∣ht(d);w)−λ2∑i=1nwi2\begin{align} L(w) &= \log \prod_{d = 1}^D\prod_{t=1}^T P(s_t^{(d)}\mid h_t^{(d)};w) - \frac{\lambda}{2} \Vert w \Vert_2^2 \notag\\ &= \sum_{d=1}^D\sum_{t=1}^T \log P(s_t^{(d)}\mid h_t^{(d)};w) - \frac{\lambda}{2} \sum_{i=1}^n w_i^2 \tag{4} \end{align} L(w)=logd=1∏Dt=1∏TP(st(d)∣ht(d);w)−2λ∥w∥22=d=1∑Dt=1∑TlogP(st(d)∣ht(d);w)−2λi=1∑nwi2(4)
将全部样本每个时刻的条件概率之积作为似然函数,同时引入正则化项。加入正则化项主要是考虑到最大似然估计可能会遇到特征函数数量特别大的情况,从而导致个别特征函数只有一个样本符合的极端情形。可以想象到,这些特征函数对应的权重会趋于无穷,所以样本对应的条件概率为 111,这使得模型对于这些样本不具有很好的泛化能力。举个简单的例子,训练集中满足特征函数 f1f_1f1 的只有样本 Yellow River
,训练后的特征函数 f1f_1f1 对应的权重 w1w_1w1 趋于无穷,P(st=River∣st−1=Yellow)=1P(s_t = {\rm River} \mid s_{t-1} = {\rm Yellow})=1P(st=River∣st−1=Yellow)=1,也就是说,训练后的模型只要遇到 Yellow
就一定会将下一个时刻的状态解码为 River
,这显然是不合理的,加上正则化项以缓解这类问题。
我认为文献 [3] 中的学习算法,即式 (4)(4)(4) 有些许不合理,因为无法类似于最大熵模型一样证明对偶函数与对数似然函数等价。
再来讨论我认为比较合理的文献 [1] 提出的学习方案。
该文献提出的方案为,在一阶 MEMM 的前提下,将 <st−1(d),st(d),Od><s_{t-1}^{(d)},s_t^{(d)},O_d><st−1(d),st(d),Od> 视为一个单元,那么对于大小为 DDD,各序列长度均为 TTT 的训练集而言,可以拆分出大约 D×TD\times TD×T 个单元。这些单元相当于特征为 ht(d)=<st−1(d),Od,t>h_t^{(d)}=<s_{t-1}^{(d)},O_d,t>ht(d)=<st−1(d),Od,t>,标签为 st(d)s_t^{(d)}st(d) 的新样本,但是为了区别于“样本”,仍然采用“单元”来叙述。为了表示更加简洁,我们为这些单元统一下标,即第 iii 个单元为 <hi,si><h_i, s_i><hi,si>。
我们总共要训练出 NNN 个最大熵模型,NNN 为上面介绍的可能的状态数。每个单元的三元组 hih_ihi 中的前一个时刻状态总共有 NNN 种取值,将相同取值的单元划分为同组,总共可以将单元分为 NNN 组,每组单元训练一个最大熵模型,即每组单元对应一个权重向量。这些最大熵模型的训练方法类似,所以这里只讨论其中一组的训练方法。
假设这组单元对应的前一个时刻状态取值均为 s′s's′,单元个数为 Ks′K_{s'}Ks′,对这组单元重新顺序编号为 1∼Ks′1\sim K_{s'}1∼Ks′,那么这些单元对应的 <h,s><h,s><h,s> 的下标也对应重新编号为 1∼Ks′1\sim K_{s'}1∼Ks′。定义特征函数的经验期望:
EP~(f)=1Ks′∑k=1Ks′f(hk,sk)(5)E_{\tilde P} (f) = \frac{1}{K_{s'}} \sum_{k=1}^{K_{s'}} f (h_k, s_k) \tag{5} EP~(f)=Ks′1k=1∑Ks′f(hk,sk)(5)
使用估计的条件概率计算的特征函数的真实期望:
EP(f)=1Ks′∑k=1Ks′∑i=1NP(s=qi∣s′,hk)f(hk,s=qi)=1Ks′∑k=1Ks′∑i=1NPs′(s=qi∣hk)f(hk,s=qi)(6)\begin{align} E_P(f) &= \frac{1}{K_{s'}} \sum_{k=1}^{K_{s'}} \sum_{i = 1}^N P(s = q_i\mid s', h_k) f(h_k, s=q_i) \\ &= \frac{1}{K_{s'}} \sum_{k=1}^{K_{s'}} \sum_{i = 1}^N P_{s'}(s = q_i\mid h_k) f(h_k, s=q_i) \\ \end{align} \tag{6} EP(f)=Ks′1k=1∑Ks′i=1∑NP(s=qi∣s′,hk)f(hk,s=qi)=Ks′1k=1∑Ks′i=1∑NPs′(s=qi∣hk)f(hk,s=qi)(6)
最大熵模型对应优化问题的约束条件(之一)为
EP~(fi)=EP(f),i=1,2,…,nE_{\tilde P}(f_i) = E_P(f),\space\space\space\space i = 1,2,\dots, n EP~(fi)=EP(f), i=1,2,…,n
其中 nnn 为特征函数个数。最大熵模型的学习转化为约束最优化问题,优化问题对应的目标函数为
maxPs′H(Ps′)=−∑k=1Ks′∑i=1NPs′(hk)Ps′(s=qi∣hk)logPs′(s=qi∣hk)\max_{P_{s'}} \space \space\space\space H(P_{s'}) = -\sum_{k=1}^{K_{s'}} \sum_{i=1}^N P_{s'}(h_k)P_{s'}(s=q_i\mid h_k) \log P_{s'}(s=q_i\mid h_k) Ps′max H(Ps′)=−k=1∑Ks′i=1∑NPs′(hk)Ps′(s=qi∣hk)logPs′(s=qi∣hk)
由于枚举的单元都是被观察到的单元,所以 Ps′(hk)=1P_{s'}(h_k)=1Ps′(hk)=1。故,完整的约束优化问题为
maxPs′H(Ps′)=−∑k=1Ks′∑i=1NPs′(s=qi∣hk)logPs′(s=qi∣hk)s.t.Ks′EP(fi)=Ks′EP~(fi),i=1,2,…,n∑k=1Ks′∑i=1NPs′(s=qi∣hk)=Ks′\begin{array}{l} \max\limits_{P_{s'}} & H(P_{s'}) = -\sum\limits_{k=1}^{K_{s'}} \sum\limits_{i=1}^NP_{s'}(s=q_i\mid h_k) \log P_{s'}(s=q_i\mid h_k) \\ s.t. & K_{s'}E_P(f_i) = K_{s'}E_{\tilde P} (f_i),\space\space\space\space i=1,2,\dots, n \\ & \sum\limits_{k=1}^{K_{s'}} \sum\limits_{i=1}^N P_{s'}(s=q_i\mid h_k) = K_{s'} \end{array} Ps′maxs.t.H(Ps′)=−k=1∑Ks′i=1∑NPs′(s=qi∣hk)logPs′(s=qi∣hk)Ks′EP(fi)=Ks′EP~(fi), i=1,2,…,nk=1∑Ks′i=1∑NPs′(s=qi∣hk)=Ks′
转换为最小化问题后对应的拉格朗日函数为
L(Ps′,w)=−H(Ps′)+w0(Ks′−∑k=1Ks′∑i=1NPs′(s=qi∣hk))+Ks′∑i=1nwi(EP~(fi)−EP(fi))=∑k=1Ks′∑i=1NPs′(s=qi∣hk)logPs′(s=qi∣hk)+w0(Ks′−∑k=1Ks′∑i=1NPs′(s=qi∣hk))+∑i=1nwi∑k=1Ks′(fi(hk,sk)−∑i=1NPs′(s=qi∣hk)fi(hk,s=qi))\begin{align} L(P_{s'},w) &= -H(P_{s'})+ w_0 \Big(K_{s'} - \sum_{k=1}^{K_{s'}}\sum_{i=1}^N P_{s'} (s=q_i\mid h_k)\Big) + K_{s'}\sum_{i=1}^n w_i\big( E_{\tilde P}(f_i) - E_{P}(f_i) \big) \notag\\ &= \sum\limits_{k=1}^{K_{s'}} \sum\limits_{i=1}^N P_{s'}(s=q_i\mid h_k) \log P_{s'}(s=q_i\mid h_k) + w_0 \Big(K_{s'} - \sum_{k=1}^{K_{s'}}\sum_{i=1}^N P_{s'} (s=q_i\mid h_k)\Big) + \sum_{i=1}^n w_i\sum_{k=1}^{K_{s'}}\Big( f_i(h_k, s_k) - \sum_{i=1}^N P_{s'}(s=q_i\mid h_k)f_i(h_k, s=q_i) \Big) \notag \\ \end{align} L(Ps′,w)=−H(Ps′)+w0(Ks′−k=1∑Ks′i=1∑NPs′(s=qi∣hk))+Ks′i=1∑nwi(EP~(fi)−EP(fi))=k=1∑Ks′i=1∑NPs′(s=qi∣hk)logPs′(s=qi∣hk)+w0(Ks′−k=1∑Ks′i=1∑NPs′(s=qi∣hk))+i=1∑nwik=1∑Ks′(fi(hk,sk)−i=1∑NPs′(s=qi∣hk)fi(hk,s=qi))
拉格朗日函数对 Ps′(s∣h)P_{s'}(s\mid h)Ps′(s∣h) 求偏导,其中 sss 和 hhh 均为确切取值,假设这 Ks′K_{s'}Ks′ 个单元中为 <h,s><h,s><h,s> 的单元个数为 KKK,则
∂L∂Ps′(s∣h)=K(logPs′(s∣h)+1)−Kw0−K∑i=1Nwifi(h,s)\frac{\partial L}{\partial P_{s'}(s\mid h)} = K\big(\log P_{s'}(s\mid h)+1\big) - Kw_0 - K\sum_{i=1}^N w_if_i(h, s) ∂Ps′(s∣h)∂L=K(logPs′(s∣h)+1)−Kw0−Ki=1∑Nwifi(h,s)
令偏导为零,利用 ∑i=1NPs′(s=qi∣h)\sum\limits_{i=1}^N P_{s'}(s=q_i\mid h)i=1∑NPs′(s=qi∣h) 可得最大熵模型
Ps′(s∣h)=exp(w⋅f(h,s))Z(h,s′)(7)P_{s'}(s\mid h) = \frac{\exp\big( w·f(h, s) \big)}{Z(h,s')} \tag{7} Ps′(s∣h)=Z(h,s′)exp(w⋅f(h,s))(7)
其中
Z(h,s′)=∑i=1Nexp(w⋅f(h,s=qi))Z(h, s') = \sum\limits_{i=1}^N\exp\big( w·f(h, s=q_i) \big) Z(h,s′)=i=1∑Nexp(w⋅f(h,s=qi))
其中 w⋅fw·fw⋅f 为向量内积。
将式 (7)(7)(7) 代入拉格朗提函数得到对偶函数
ψ(w)=∑k=1Ks′∑i=1NPs′(s=qi∣hk)logPs′(s=qi∣hk)+∑i=1nwi∑k=1Ks′(fi(hk,sk)−∑i=1NPs′(s=qi∣hk)fi(hk,s=qi))=∑k=1Ks′∑i=1NPs′(s=qi∣hk)logPs′(s=qi∣hk)+∑k=1Ks′∑i=1nwifi(hk,sk)−∑k=1Ks′∑i=1NPs′(s=qi∣hk)∑i=1nwifi(hk,s=qi)=∑k=1Ks′∑i=1nwifi(hk,sk)+∑k=1Ks′∑i=1NPs′(s=qi∣hk)logPs′(s=qi∣hk)−∑k=1Ks′∑i=1NPs′(s=qi∣hk)∑i=1nwifi(hk,s=qi)=∑k=1Ks′∑i=1nwifi(hk,sk)+∑k=1Ks′∑i=1NPs′(s=qi∣hk)(logPs′(s=qi∣hk)−∑i=1nwifi(hk,s=qi))=∑k=1Ks′∑i=1nwifi(hk,sk)−∑k=1Ks′∑i=1NPs′(s=qi∣hk)logZ(hk,s′)=∑k=1Ks′∑i=1nwifi(hk,sk)−∑k=1Ks′logZ(hk,s′)∑i=1NPs′(s=qi∣hk)=∑k=1Ks′∑i=1nwifi(hk,sk)−∑k=1Ks′logZ(hk,s′)\begin{align} \psi(w) &= \sum\limits_{k=1}^{K_{s'}} \sum\limits_{i=1}^N P_{s'}(s=q_i\mid h_k) \log P_{s'}(s=q_i\mid h_k) + \sum_{i=1}^n w_i\sum_{k=1}^{K_{s'}}\Big( f_i(h_k, s_k) - \sum_{i=1}^N P_{s'}(s=q_i\mid h_k)f_i(h_k, s=q_i) \Big) \notag \\ &= \sum_{k=1}^{K_{s'}} \sum_{i=1}^N P_{s'} (s=q_i \mid h_k) \log P_{s'} (s = q_i\mid h_k) + \sum_{k=1}^{K_{s'}}\sum_{i=1}^n w_i f_i(h_k,s_k) - \sum_{k=1}^{K_{s'}}\sum_{i=1}^NP_{s'}(s = q_i\mid h_k) \sum_{i=1}^n w_i f_i(h_k, s=q_i) \notag \\ &= \sum_{k=1}^{K_{s'}}\sum_{i=1}^n w_i f_i(h_k,s_k) + \sum_{k=1}^{K_{s'}} \sum_{i=1}^N P_{s'} (s=q_i \mid h_k) \log P_{s'} (s = q_i\mid h_k) - \sum_{k=1}^{K_{s'}}\sum_{i=1}^NP_{s'}(s = q_i\mid h_k) \sum_{i=1}^n w_i f_i(h_k, s=q_i) \notag\\ &= \sum_{k=1}^{K_{s'}}\sum_{i=1}^n w_i f_i(h_k,s_k) + \sum_{k=1}^{K_{s'}} \sum_{i=1}^N P_{s'} (s=q_i \mid h_k) \Big(\log P_{s'} (s = q_i\mid h_k) - \sum_{i=1}^n w_i f_i(h_k, s=q_i) \Big) \notag\\ &= \sum_{k=1}^{K_{s'}}\sum_{i=1}^n w_i f_i(h_k,s_k) - \sum_{k=1}^{K_{s'}} \sum_{i=1}^N P_{s'} (s=q_i \mid h_k) \log Z(h_k, s')\notag\\ &= \sum_{k=1}^{K_{s'}}\sum_{i=1}^n w_i f_i(h_k,s_k) - \sum_{k=1}^{K_{s'}} \log Z(h_k, s')\sum_{i=1}^N P_{s'}(s=q_i\mid h_k) \notag\\ &= \sum_{k=1}^{K_{s'}}\sum_{i=1}^n w_i f_i(h_k,s_k) - \sum_{k=1}^{K_{s'}} \log Z(h_k, s') \notag\\ \end{align} ψ(w)=k=1∑Ks′i=1∑NPs′(s=qi∣hk)logPs′(s=qi∣hk)+i=1∑nwik=1∑Ks′(fi(hk,sk)−i=1∑NPs′(s=qi∣hk)fi(hk,s=qi))=k=1∑Ks′i=1∑NPs′(s=qi∣hk)logPs′(s=qi∣hk)+k=1∑Ks′i=1∑nwifi(hk,sk)−k=1∑Ks′i=1∑NPs′(s=qi∣hk)i=1∑nwifi(hk,s=qi)=k=1∑Ks′i=1∑nwifi(hk,sk)+k=1∑Ks′i=1∑NPs′(s=qi∣hk)logPs′(s=qi∣hk)−k=1∑Ks′i=1∑NPs′(s=qi∣hk)i=1∑nwifi(hk,s=qi)=k=1∑Ks′i=1∑nwifi(hk,sk)+k=1∑Ks′i=1∑NPs′(s=qi∣hk)(logPs′(s=qi∣hk)−i=1∑nwifi(hk,s=qi))=k=1∑Ks′i=1∑nwifi(hk,sk)−k=1∑Ks′i=1∑NPs′(s=qi∣hk)logZ(hk,s′)=k=1∑Ks′i=1∑nwifi(hk,sk)−k=1∑Ks′logZ(hk,s′)i=1∑NPs′(s=qi∣hk)=k=1∑Ks′i=1∑nwifi(hk,sk)−k=1∑Ks′logZ(hk,s′)
这里的最大熵模型同样满足对偶函数与对数似然函数等价。对数似然函数为
L(w)=log∏k=1Ks′Ps′(sk∣hk)=∑k=1Ks′logPs′(sk∣hk)=∑k=1Ks′(∑i=1nwifi(hk,sk)−logZ(hk,s′))=∑k=1Ks′∑i=1nwifi(hk,sk)−∑k=1Ks′logZ(hk,s′)\begin{align} L(w) &= \log \prod_{k=1}^{K_{s'}} P_{s'} (s_k \mid h_k)\notag \\ &= \sum_{k=1}^{K_{s'}} \log P_{s'} (s_k\mid h_k) \notag\\ &= \sum_{k=1}^{K_{s'}} \left( \sum_{i=1}^n w_if_i(h_k,s_k) - \log Z(h_k, s') \right) \notag\\ &= \sum_{k=1}^{K_{s'}}\sum_{i=1}^n w_i f_i(h_k,s_k) - \sum_{k=1}^{K_{s'}} \log Z(h_k, s')\notag \\ \end{align} L(w)=logk=1∏Ks′Ps′(sk∣hk)=k=1∑Ks′logPs′(sk∣hk)=k=1∑Ks′(i=1∑nwifi(hk,sk)−logZ(hk,s′))=k=1∑Ks′i=1∑nwifi(hk,sk)−k=1∑Ks′logZ(hk,s′)
可见,对数似然函数与对偶函数相同。
极大化 ψ(w)\psi(w)ψ(w) 可以采用最大熵模型中提到的最优化算法,比如通用迭代尺度法(GIS)、改进的迭代尺度法(IIS)等。文献 [3] 中介绍到,对每组单元对应的 ψ(w)\psi(w)ψ(w) 分别采用 GIS 进行迭代,最终得到 NNN 组权重向量,对应 NNN 个最大熵模型,这 NNN 个最大熵模型完整地描述了概率分布 P(s∣s′,h)P(s\mid s', h)P(s∣s′,h)。
可以想象,如果满足二阶马尔可夫性,那么就需要训练 N×NN\times NN×N 个最大熵模型来描述概率分布 P(s∣s′,s′′,h)P(s\mid s', s'', h)P(s∣s′,s′′,h),具体过程与上面类似。
对比两个文献提出的学习方法,文献 [1] 的观点更能让我信服,而且可以通过比较合理的推导进行证明。
这部分推导属于我自己研究的,正确性无法保证,仅供参考。
解码问题
Daniel Jurafsky 在《Speech and Language Processing》中提到,MEMM 的解码也可以采用类似于 HMM 中近似算法的贪心思想。使用这种方式解码,需要对输入观测序列从头到尾运行分类器,每次分类都作出当前最好的决定。算法流程大致描述为,从头到尾扫描整个观测序列,对于每个时刻,根据前一个时刻我们的解码得到的状态选择当前时刻对应的最大熵模型学习到的条件概率分布来确定当前时刻的解码状态,依次类推。特别地,第一个时刻的解码状态可以仅依赖于观测。
至于 Viterbi 算法也是与 HMM 一致。定义
δt(i)=maxs1,s2,…,st−1P(O,st=qi,st−1,…,s1∣w),1≤i≤N\delta_t(i) = \max _{s_1,s_2,\dots, s_{t-1}} P(O,s_{t}=q_i,s_{t-1},\dots, s_1\mid w),\space\space\space\space 1\le i\le N δt(i)=s1,s2,…,st−1maxP(O,st=qi,st−1,…,s1∣w), 1≤i≤N
其中 PPP 是 Ps′P_{s'}Ps′ 的简写。对式 (1)(1)(1) 变形可得
δt(i)=maxs1,s2,…,st−1P(st=qi∣st−1,…,s1,O,w)P(st−1,…,s1,O∣w)=max1≤j≤NP(st=qi∣ht,w)δt−1(j)=max1≤j≤N[δt−1(j)]P(st=qi∣ht,w)\begin{align} \delta_t(i) &= \max_{s_1, s_2,\dots, s_{t-1}} P(s_t=q_i\mid s_{t-1},\dots, s_1, O, w) P(s_{t-1},\dots, s_1, O \mid w) \notag\\ &= \max_{1\le j \le N} P(s_t=q_i \mid h_t, w) \delta_{t-1}(j) \notag\\ &= \max_{1\le j \le N} [\delta_{t-1}(j)]P(s_t=q_i \mid h_t, w) \notag\\ \end{align} δt(i)=s1,s2,…,st−1maxP(st=qi∣st−1,…,s1,O,w)P(st−1,…,s1,O∣w)=1≤j≤NmaxP(st=qi∣ht,w)δt−1(j)=1≤j≤Nmax[δt−1(j)]P(st=qi∣ht,w)
如果考虑一阶马尔可夫性,那么其中 ht=<st−1,O,t>h_t=< s_{t-1},O, t>ht=<st−1,O,t>。初始化
δ1(i)=P(s1=qi∣h1,w)\delta_1(i) = P(s_1 = q_i\mid h_1, w) δ1(i)=P(s1=qi∣h1,w)
其中 h1=<O,1>h_1 = <O, 1>h1=<O,1>。
定义在时刻 ttt 状态为 qiq_iqi 的所有单个路径 (s1,s2,…,st−1,st=qi)(s_1,s_2,\dots ,s_{t-1},s_t=q_i)(s1,s2,…,st−1,st=qi) 中概率最大的路径的第 t−1t - 1t−1 个结点为
ψt(i)=argmax1≤j≤N[δt−1(j)],,1≤i≤N\psi_t(i) = {\rm arg} \max_{1\le j \le N} [\delta_{t-1}(j)],\space\space\space\space,1\le i \le N ψt(i)=arg1≤j≤Nmax[δt−1(j)], ,1≤i≤N
ψt(i)\psi_t(i)ψt(i) 保存的是时刻 ttt 状态 st=qis_t=q_ist=qi 由时刻 t−1t-1t−1 的哪种状态转移而来。初始化
ψ1(i)=0\psi_1(i) = 0 ψ1(i)=0
终止
P∗=max1≤i≤NδT(i)sT∗=argmax1≤i≤N[δT(i)]P^* = \max_{1\le i \le N}\delta_T(i) \\ s_T^* = {\rm arg}\max \limits_{1\le i\le N} [\delta_T(i)] P∗=1≤i≤NmaxδT(i)sT∗=arg1≤i≤Nmax[δT(i)]
路径回溯,只需要根据 st∗=ψt+1(st+1∗)s_t^* = \psi_{t+1}(s_{t+1}^*)st∗=ψt+1(st+1∗) 逆序递推即可得到最优路径 S∗=(s1∗,s2∗,…,sT∗)S^*=(s_1^*,s_2^*,\dots, s_T^*)S∗=(s1∗,s2∗,…,sT∗)。
关于“双向 MEMM”、“MEMM 的缺点”等扩展问题在这里不展开,可以见其他文章。
REF
[1] McCallum A, Freitag D, Pereira F C N. Maximum entropy Markov models for information extraction and segmentation[C]//Icml. 2000, 17(2000): 591-598.
[2] Ratnaparkhi A. A maximum entropy model for part-of-speech tagging[C]//Conference on empirical methods in natural language processing. 1996.
[3] fall2014-loglineartaggers.pdf (columbia.edu)
[4] Collins M. Log-linear models[J]. Self-published Tutorial, 2005.
[5] Collins M. Log-linear models, memms, and crfs[J]. Columbia University lecture, 2015.
[6] Lafferty J, McCallum A, Pereira F C N. Conditional random fields: Probabilistic models for segmenting and labeling sequence data[J]. 2001.
[7] Wang H, Fei H, Yu Q, et al. A motifs-based Maximum Entropy Markov Model for realtime reliability prediction in System of Systems[J]. Journal of Systems and Software, 2019, 151: 180-193.
[8]《Speech and Language Processing》Daniel Jurafsky 等著
[9]《百面机器学习》诸葛越等著
[10] 【机器学习】最大熵模型【上】最大熵模型概述与约束最优化问题 - CSDN
[11] 【机器学习】最大熵模型【下】最大熵模型学习的最优化算法 - CSDN
[12] 【自然语言处理】隐马尔可夫模型【Ⅴ】解码问题 - CSDN
[13] (【自然语言处理】隐马尔可夫模型【Ⅱ】隐马尔科夫模型概述 - CSDN
【自然语言处理】最大熵马尔可夫模型相关推荐
- 隐马尔可夫模型,最大熵马尔可夫模型和条件随机场的区别与联系
隐马尔可夫模型(Hidden Markov Model,HMM),最大熵马尔可夫模型(Maximum Entropy Markov Model,MEMM)以及条件随机场(Conditional ...
- 最大熵模型(ME)和最大熵马尔可夫模型(MEMM)
文章目录 最大熵模型(Maximum Entropy Model, MEM) 最大熵思想 最大熵模型表示 最大熵模型学习 最大熵模型求解过程举例 最大熵模型与极大似然估计 最大熵模型与逻辑回归 模型学 ...
- 【机器学习】最大熵马尔科夫模型
点击上方蓝色字体,关注AI小白入门哟 跟着博主的脚步,每天进步一点点 本文介绍了最大熵马尔可夫模型,在隐马尔可夫模型(隐状态序列)的基础上应用最大熵模型思想,将一个概率生成模型转化为概率判别模型,同样 ...
- 【自然语言处理】隐马尔可夫模型【Ⅰ】马尔可夫模型
有任何的书写错误.排版错误.概念错误等,希望大家包含指正. 由于字数限制,分成六篇博客. [自然语言处理]隐马尔可夫模型[Ⅰ]马尔可夫模型 [自然语言处理]隐马尔科夫模型[Ⅱ]隐马尔科夫模型概述 [自 ...
- 马尔可夫模型与条件随机场模型
https://blog.csdn.net/s151506879/article/details/51912595 条件随机场模型是由Lafferty在2001年提出的一种典型的判别式模型.它在观测序 ...
- 第五篇:HMM 隐马尔可夫模型
目录 概览 训练 预测(维特比算法) 概览 HMM本身也是用于序列标注,为什么需要HMM? POS 标记,我们按照一句话为一个类别,也就是将该句子中的每个词的对应的类别,连起来作为一个类别,比如NN_ ...
- 条件随机场(CRF)和隐马尔科夫模型(HMM)最大区别在哪里?CRF的全局最优体现在哪里?
作者:尔总的马甲 链接:https://www.zhihu.com/question/53458773/answer/554436625 来源:知乎 著作权归作者所有.商业转载请联系作者获得授权,非商 ...
- 标记偏置 隐马尔科夫 最大熵马尔科夫 HMM MEMM
转载自:http://blog.csdn.net/xum2008/article/details/38147425 隐马尔科夫模型(HMM): 图1. 隐马尔科夫模型 隐马尔科夫模型的缺点: 1.HM ...
- 马尔可夫模型与隐马尔可夫模型
http://www.cnblogs.com/baiboy/p/hmm2.html 转自:白宁超 摘要:最早接触马尔可夫模型的定义源于吴军先生<数学之美>一本.直到做自然语言处理时,才真正 ...
最新文章
- javascript中this那些事
- verilog中wire和reg的区别,什么时候用wire?什么时候用reg?
- 通过点击切换文本框内容的脚本示例
- JS判断是否是科学计算结果并处理
- HTML5期末大作业:腾讯游戏网站设计——腾讯游戏官网(13页) HTML+CSS+JavaScript web网页设计与开发
- 如何获得SQL Server 2012 MCSA认证
- PCI Express 16X和PCIE 2.0有什么区别?哪个更好?
- Snipaste截屏软件使用
- 中国超市智能储物柜市场需求动态与盈利前景预测报告2022-2027年
- 两种python调包进行拉丁超立方采样(lhs)
- Unity导入模型贴贴图一面有贴图另一面透明的解决方法
- 修改登录页面Login
- 光盘刻录系列之二刻录光盘的程序步骤
- 字体图标之Symbol用法
- 在INTEL KABYALAKE平台上运行COREBOOT 记录
- 基于Python Django Mysql 开发的宠物用品商城
- 北理工计算机贾云,徐畅_北京理工大学计算机学院
- happyhbase
- pytorch实现图像的腐蚀和膨胀
- 应届生,找java后端开发要什么水平才算合格?