Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

3.1 n 元语法

先从计算 P(wh)P(w|h) 这一任务开始,也就是给定某段历史 hh 时词 ww 的概率。假设历史 hh 是“The water of Walden Pond is so beautifully”,我们希望知道下一个词是 blue 的概率:

P( blue  The water of Walden Pond is so beautifully )(3.1)P (\text { blue } | \text { The water of Walden Pond is so beautifully })\tag{3.1}

一种估计方法是直接使用相对频次:取一个非常大的语料库,统计 The water of Walden Pond is so beautifully 出现了多少次,再统计它后面跟着 blue 的次数。这相当于回答:“在看到历史 hh 的所有情况中,有多少次后面跟着词 ww?”

P(blue|The water of Walden Pond is so beautifully)=C(The water of Walden Pond is so beautifully blue)C(The water of Walden Pond is so beautifully)(3.2)\begin{array}{c} P (\text {blue|The water of Walden Pond is so beautifully}) = \\ \frac {C (\text {The water of Walden Pond is so beautifully blue})}{C (\text {The water of Walden Pond is so beautifully})} \end{array}\tag{3.2}

如果语料库足够大,就可以计算这两个计数,并依据式 3.2 估计概率。然而,即使整个网络也不足以让我们准确估计完整句子的计数。原因在于语言具有创造性:新句子不断产生,我们无法期望对整个句子这样庞大的对象取得准确计数。因此,要估计给定历史 hh 时词 ww 的概率,或者整个词序列 WW 的概率,还需要更巧妙的方法。

首先约定一些记法。本章仍会沿用“词”这一说法,尽管实践中的语言模型通常在第 2 章介绍的 BPE 词元等词元上计算。对于随机变量 XiX_i 取值为“the”的概率 P(Xi="the")P(X_i="the"),简写为 P(the)P(the)。由 nn 个词组成的序列写作 w1wnw_1\ldots w_nw1:nw_{1:n}。因此,w1:n1w_{1:n-1} 表示字符串 w1,w2,,wn1w_1,w_2,\ldots,w_{n-1};我们还会使用等价记法 w<nw_{<n},读作“从 w1w_1 到并包括 wn1w_{n-1} 的所有 ww 元素”。序列中各词分别取特定值的联合概率 P(X1=w1,X2=w2,,Xn=wn)P(X_1=w_1,X_2=w_2,\ldots,X_n=w_n),简写为 P(w1,w2,,wn)P(w_1,w_2,\ldots,w_n)

怎样计算 P(w1,w2,,wn)P(w_1,w_2,\ldots,w_n) 这样的完整序列概率?一种做法是使用概率的链式法则(chain rule)分解它:

P(X1...Xn)=P(X1)P(X2X1)P(X3X1:2)P(XnX1:n1)=k=1nP(XkX1:k1)(3.3)\begin{array}{l l} P (X _ {1}... X _ {n}) & = P (X _ {1}) P (X _ {2} | X _ {1}) P (X _ {3} | X _ {1: 2}) \dots P (X _ {n} | X _ {1: n - 1}) \\ & = \prod_ {k = 1} ^ {n} P (X _ {k} | X _ {1: k - 1}) \end{array}\tag{3.3}

把链式法则应用于词,得到:

P(w1:n)=P(w1)P(w2w1)P(w3w1:2)P(wnw1:n1)=k=1nP(wkw1:k1)(3.4)\begin{array}{l l} P (w _ {1: n}) & = P (w _ {1}) P (w _ {2} | w _ {1}) P (w _ {3} | w _ {1: 2}) \dots P (w _ {n} | w _ {1: n - 1}) \\ & = \prod_ {k = 1} ^ {n} P (w _ {k} | w _ {1: k - 1}) \end{array}\tag{3.4}

链式法则揭示了计算序列联合概率与计算给定先前词语时某个词的条件概率之间的联系。式 3.4 表明,可以把许多条件概率相乘,估计整个词序列的联合概率。然而,使用链式法则似乎并没有真正解决问题!我们仍不知道怎样计算给定一长串先前词语时某个词的精确概率 P(wnw1:n1)P(w_n|w_{1:n-1})。如前所述,不能简单统计语料库中每个长字符串后面每个词的出现次数,因为语言具有创造性,某个特定上下文可能从未出现过。

3.1.1 Markov 假设

n 元模型背后的直觉是:不再依据一个词的完整历史计算其概率,而只用最后几个词来近似这段历史。

例如,二元模型只使用给定前一个词时的条件概率 P(wnwn1)P(w_n|w_{n-1}),近似给定所有先前词语时某个词的概率 P(wnw1:n1)P(w_n|w_{1:n-1})。换句话说,我们不计算:

P( blue  The water of Walden Pond is so beautifully )(3.5)P (\text { blue } | \text { The water of Walden Pond is so beautifully })\tag{3.5}

而使用下面的概率进行近似:

P( blue  beautifully )(3.6)P (\text { blue } | \text { beautifully })\tag{3.6}

因此,使用二元模型预测下一词的条件概率时,会作出如下近似:

P(wnw1:n1)P(wnwn1)(3.7)P (w _ {n} | w _ {1: n - 1}) \approx P (w _ {n} | w _ {n - 1})\tag{3.7}

一个词的概率只取决于前一个词,这项假设称为 Markov 假设(Markov assumption)。Markov 模型(Markov model)是一类概率模型,它们假设无需向过去回看太远,就能预测未来某个单位的概率。二元模型回看一个词;我们可以把它推广为回看两个词的三元模型,并进一步推广为回看 n1n-1 个词的 n 元模型。

下面给出使用 n 元语法近似序列中下一词条件概率的一般公式。这里用 NN 表示 n 元语法的大小,因此 N=2N=2 表示二元语法,N=3N=3 表示三元语法。于是,给定完整上下文时某个词的概率近似为:

P(wnw1:n1)P(wnwnN+1:n1)(3.8)P (w _ {n} | w _ {1: n - 1}) \approx P (w _ {n} | w _ {n - N + 1: n - 1})\tag{3.8}

利用单个词概率的二元假设,把式 3.7 代入式 3.4,即可计算完整词序列的概率:

P(w1:n)k=1nP(wkwk1)(3.9)P (w _ {1: n}) \approx \prod_ {k = 1} ^ {n} P (w _ {k} | w _ {k - 1})\tag{3.9}

3.1.2 如何估计概率

怎样估计二元或 n 元概率?一种直观方法称为最大似然估计(maximum likelihood estimation,MLE)。为了得到 n 元模型参数的 MLE 估计,我们从语料库取得计数,再把计数归一化到 0 与 1 之间。对概率模型而言,归一化就是除以某个总计数,使所得概率位于 0 与 1 之间,且总和为 1。

例如,要计算给定前一个词 wn1w_{n-1} 时词 wnw_n 的某个二元概率,可以统计二元语法 C(wn1wn)C(w_{n-1}w_n) 的次数,再除以共享同一个首词 wn1w_{n-1} 的所有二元语法之总数:

P(wnwn1)=C(wn1wn)wC(wn1w)(3.10)P (w _ {n} | w _ {n - 1}) = \frac {C (w _ {n - 1} w _ {n})}{\sum_ {w} C (w _ {n - 1} w)}\tag{3.10}

由于以给定词 wn1w_{n-1} 开头的所有二元语法计数之和,必然等于该词的一元计数 C(wn1)C(w_{n-1}),所以可以把公式简化为:

P(wnwn1)=C(wn1wn)C(wn1)(3.11)P (w _ {n} | w _ {n - 1}) = \frac {C (w _ {n - 1} w _ {n})}{C (w _ {n - 1})}\tag{3.11}

下面用一个只含三个句子的微型语料库完成示例。首先,需要在每个句子开头加入特殊符号 s\langle s\rangle,为第一个词提供二元上下文;还要在句末加入特殊符号 /s\langle/s\rangle[1]

sI am Sam/ssSam I am/ssI do not like green eggs and ham/s\begin{array}{l} \langle s \rangle \text {I am Sam} \langle / s \rangle \\ \langle s \rangle \text {Sam I am} \langle / s \rangle \\ \langle s \rangle \text {I do not like green eggs and ham} \langle / s \rangle \end{array}

下面计算这个语料库中的若干二元概率:

P(I<s>)=23=0.67P(Sam<s>)=13=0.33P(amI)=23=0.67P(</s>Sam)=12=0.5P(Samam)=12=0.5P(doI)=13=0.33\begin{array}{l l l} P (\mathsf {I} | <\mathsf {s}>) = \frac {2}{3} = 0.67 & P (\mathsf {Sam} | <\mathsf {s}>) = \frac {1}{3} = 0.33 & P (\mathsf {am} | \mathsf {I}) = \frac {2}{3} = 0.67 \\ P (</\mathsf {s}> | \mathsf {Sam}) = \frac {1}{2} = 0.5 & P (\mathsf {Sam} | \mathsf {am}) = \frac {1}{2} = 0.5 & P (\mathsf {do} | \mathsf {I}) = \frac {1}{3} = 0.33 \end{array}

对于 MLE n 元参数估计的一般情形(N=2,3,N=2,3,\ldots):

P(wnwnN+1:n1)=C(wnN+1:n1wn)C(wnN+1:n1)(3.12)P (w _ {n} | w _ {n - N + 1: n - 1}) = \frac {C (w _ {n - N + 1 : n - 1} w _ {n})}{C (w _ {n - N + 1 : n - 1})}\tag{3.12}

式 3.12 与式 3.11 一样,用某个特定序列的观测频次除以其前缀的观测频次来估计 n 元概率。这个比值称为相对频次(relative frequency)。前文已经指出,这种使用相对频次估计概率的方法是最大似然估计的一例。在 MLE 中,所得参数集合会使给定模型 MM 时训练集 TT 的似然最大,即 P(TM)P(T|M) 最大。例如,假设 Chinese 在一个含 100 万词的语料库中出现 400 次。从另一段同样含约 100 万词的文本中随机选择一个词,它是 Chinese 的概率是多少?MLE 概率估计为 400/1,000,000=0.0004400/1,000,000=0.0004。0.0004 并不是 Chinese 在所有情境中出现概率的最佳估计;在其他语料库或上下文中,它可能极少出现。但正是这个概率最有可能使 Chinese 在一个百万词语料库中出现 400 次。第 3.6 节将介绍如何稍微修改 MLE 估计,得到更好的概率估计。

下面考察一个真实但很小的语料库示例。它取自现已停止运行的 Berkeley Restaurant Project;这是上个世纪的一个对话系统,可以回答有关加利福尼亚州伯克利餐馆数据库的问题(Jurafsky et al., 1994)。下面是几条用户查询样例;文本已经过规范化,转换为小写并去除标点(网站上提供了 9,332 个句子的样本):

can you tell me about any good cantonese restaurants close by

tell me about chez panisse

i’m looking for a good place to eat breakfast

when is caffe venezia open during the day

图 3.1 展示根据经过文本规范化的 Berkeley Restaurant Project 句子构建的二元语法中,部分词的二元计数。请注意,大多数值为零。事实上,这里特意选择了彼此能够形成搭配的样例词;如果从八个随机词中选择,矩阵还会更加稀疏。

iwanttoeatchinesefoodlunchspend
i5827090002
want2060816651
to204686206211
eat0020162420
chinese100008210
food1501501400
lunch20000100
spend10100000

图 3.1 Berkeley Restaurant Project 的 9,332 句语料库中,词表 V=1446V=1446 的八个词之二元计数。零计数以灰色表示。每个单元格给出列标签词出现在行标签词之后的次数。因此,i 行 want 列的单元格表示 want 在语料库中跟随 i 出现了 827 次。

图 3.2 展示归一化后的二元概率,即把图 3.1 中的每个单元格除以对应行的一元计数;所使用的一元计数如下:

iwanttoeatchinesefoodlunchspend
253392724177461581093341278
iwanttoeatchinesefoodlunchspend
i0.0020.3300.00360000.00079
want0.002200.660.00110.00650.00650.00540.0011
to0.0008300.00170.280.0008300.00250.087
eat000.002700.0210.00270.0560
chinese0.006300000.520.00630
food0.01400.01400.000920.003700
lunch0.005900000.002900
spend0.003600.00360000

图 3.2 Berkeley Restaurant Project 的 9,332 句语料库中八个词的二元概率。零概率以灰色表示。

下面是另外几个有用的概率:

P(i<s>)=0.25P(english|want)=0.0011P(food|english)=0.5P(</s>food)=0.68\begin{array}{l l} P (\mathrm{i} | <\mathrm{s}>) = 0.25 & P (\text {english|want}) = 0.0011 \\ P (\text {food|english}) = 0.5 & P (</\mathrm{s}> | \text {food}) = 0.68 \end{array}

现在,只需把相应二元概率相乘,就能计算 I want English food 或 I want Chinese food 等句子的概率:

P(<s>i want english food</s>)=P(i<s>)P(wanti)P(englishwant)P(foodenglish)P(</s>food)=0.25×0.33×0.0011×0.5×0.68=0.000031\begin{array}{r l} P (<\mathrm{s}> \text {i want english food} </\mathrm{s}>) & \\ = P (\mathrm{i} | <\mathrm{s}>) P (\text {want} | \mathrm{i}) P (\text {english} | \text {want}) & \\ P (\text {food} | \text {english}) P (</\mathrm{s}> | \text {food}) & \\ = 0.25 \times 0.33 \times 0.0011 \times 0.5 \times 0.68 & \\ = 0.000031 & \end{array}

练习 3.2 要求计算 i want chinese food 的概率。

这些二元统计捕捉了哪些语言现象?上面的一些二元概率编码了通常被视为纯句法的事实,例如 eat 后面通常是名词或形容词,to 后面通常是动词。另一些概率可能反映个人助理任务,例如句子以 I 开头的概率很高。还有一些甚至可能属于文化事实而非语言事实,例如人们寻找中餐的概率高于寻找英餐的概率。

3.1.3 处理大型 n 元模型的规模问题

实践中的语言模型可能非常大,从而引发一系列实际问题。

对数概率 语言模型的概率总是以对数概率(log probability)形式存储,并在对数空间中计算。这是因为概率按定义小于或等于 1,相乘的概率越多,乘积就越小。把足够多的 n 元概率相乘会造成数值下溢。在对数空间中相加等价于在线性空间中相乘,因此我们通过加法组合对数概率。对数概率相加所得结果不会像概率乘积那么小。所有计算和存储都在对数空间中进行;只有最终需要报告概率时,才对对数概率取指数,转换回普通概率:

p1×p2×p3×p4=exp(logp1+logp2+logp3+logp4)(3.13)p _ {1} \times p _ {2} \times p _ {3} \times p _ {4} = \exp (\log p _ {1} + \log p _ {2} + \log p _ {3} + \log p _ {4})\tag{3.13}

本书在没有指定底数时,用 log 表示自然对数 ln

更长的上下文 为便于教学,前文只介绍了二元模型;当训练数据充足时,通常会使用以前两个词为条件的三元模型,或 4 元、5 元模型。使用这些更大的 n 元语法时,需要在句子左右边界假定额外上下文。例如,为计算句子最开头的三元概率,会为第一个三元语法使用两个伪词,即 P(Iss)P(I|\langle s\rangle\langle s\rangle)

研究者已经创建了一些大型 n 元数据集。例如,从精心整理、包含 10 亿词的美国英语语料库 COCA(Corpus of Contemporary American English)中抽取的最高频 100 万个 n 元语法(Davies, 2020);从 1 万亿词英语网络文本中抽取的 Google Web 5-gram 语料库(Franz and Brants, 2006);以及 Google Books Ngrams 语料库,其中包含来自汉语、英语、法语、德语、希伯来语、意大利语、俄语和西班牙语的 8000 亿个词元(Lin et al., 2012a)。

甚至可以使用距离极远的 n 元上下文。无限元语法(infini-gram,∞-gram)项目(Liu et al., 2024)允许 n 元语法具有任意长度。其核心思想是避免提前计算代价高昂(占用大量空间和时间)的巨型 n 元计数表;系统转而使用一种称为后缀数组(suffix array)的高效表示,在推理时快速计算任意 nn 的 n 元概率。这样,就能在含 5 万亿词元的巨大语料库上计算任意长度的 n 元语法。

构建大型 n 元语言模型时,效率十分重要。标准做法包括:只使用 4~8 位量化概率,而不是 8 字节浮点数;把词语字符串存储在磁盘上,在内存中只用 64 位哈希表示;使用“反向 trie”等特殊数据结构表示 n 元语法。语言模型也经常进行剪枝(pruning),例如只保留计数超过某个阈值的 n 元语法,或者用熵剪除不太重要的 n 元语法(Stolcke, 1998)。KenLM 等高效语言模型工具包(Heafield, 2011; Heafield et al., 2013)使用有序数组,并借助归并排序,以尽可能少的语料库遍历次数高效构建概率表。

Footnotes
  1. 我们需要句末符号,才能使二元语法成为真正的概率分布。如果没有句末符号,不再是所有句子的概率之和为 1,而是每一给定长度下所有句子的概率之和分别为 1。这样的模型会定义无穷多个概率分布,每种句子长度对应一个分布。见练习 3.5。