日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

转:数学之美番外篇:平凡而又神奇的贝叶斯方法 收藏

發布時間:2025/4/5 编程问答 45 豆豆
生活随笔 收集整理的這篇文章主要介紹了 转:数学之美番外篇:平凡而又神奇的贝叶斯方法 收藏 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

為什么80%的碼農都做不了架構師?>>> ??

轉自:http://blog.csdn.net/pongba/archive/2008/09/21/2958094.aspx

?

數學之美番外篇:平凡而又神奇的貝葉斯方法

By 劉未鵬(pongba)

C++的羅浮宮( http://blog.csdn.net/pongba )

TopLanguage( http://groups.google.com/group/pongba )

概率論只不過是把常識用數學公式表達了出來。

——拉普拉斯

記得讀本科的時候,最喜歡到城里的計算機書店里面去閑逛,一逛就是好幾個小時;有一次,在書店看到一本書,名叫貝葉斯方法。當時數學系的課程還沒有學到概率統計。我心想,一個方法能夠專門寫出一本書來,肯定很牛逼。后來,我發現當初的那個樸素歸納推理成立了——這果然是個牛逼的方法。

——題記

目錄

0. 前言
1. 歷史
??? 1.1 一個例子:自然語言的二義性
??? 1.2 貝葉斯公式?
2. 拼寫糾正
3. 模型比較與貝葉斯奧卡姆剃刀
??? 3.1 再訪拼寫糾正
??? 3.2 模型比較理論(Model Comparasion)與貝葉斯奧卡姆剃刀(Bayesian Occam’s Razor)
??? 3.3 最小描述長度原則
??? 3.4 最優貝葉斯推理
4. 無處不在的貝葉斯
??? 4.1 中文分詞
??? 4.2 統計機器翻譯????
??? 4.3 貝葉斯圖像識別,Analysis by Synthesis???
??? 4.4 EM 算法與基于模型的聚類??
??? 4.5 最大似然與最小二乘??
5. 樸素貝葉斯方法(又名“愚蠢者的貝葉斯(idiot’s bayes)”)??
??? 5.1 垃圾郵件過濾器??
??? 5.2 為什么樸素貝葉斯方法令人詫異地好——一個理論解釋
6. 層級貝葉斯模型?
??? 6.1 隱馬可夫模型(HMM)
7. 貝葉斯網絡

0. 前言

這是一篇關于貝葉斯方法的科普文,我會盡量少用公式,多用平白的語言敘述,多舉實際例子。更嚴格的公式和計算我會在相應的地方注明參考資料。貝葉斯方法被證明是非常 general 且強大的推理框架,文中你會看到很多有趣的應用。

1. 歷史

托馬斯·貝葉斯(Thomas Bayes)同學的詳細生平在 這里 。以下摘一段 wikipedia 上的簡介:

所謂的貝葉斯方法源于他生前為解決一個“逆概”問題寫的一篇文章,而這篇文章是在他死后才由他的一位朋友發表出來的。在貝葉斯寫這篇文章之前,人們已經能夠計算“正向概率”,如“假設袋子里面有N個白球,M個黑球,你伸手進去摸一把,摸出黑球的概率是多大”。而一個自然而然的問題是反過來:“如果我們事先并不知道袋子里面黑白球的比例,而是閉著眼睛摸出一個(或好幾個)球,觀察這些取出來的球的顏色之后,那么我們可以就此對袋子里面的黑白球的比例作出什么樣的推測”。這個問題,就是所謂的逆概問題。

實際上,貝葉斯當時的論文只是對這個問題的一個直接的求解嘗試,并不清楚他當時是不是已經意識到這里面包含著的深刻的思想。然而后來,貝葉斯方法席卷了概率論,并將應用延伸到各個問題領域,所有需要作出概率預測的地方都可以見到貝葉斯方法的影子,特別地,貝葉斯是機器學習的核心方法之一。這背后的深刻原因在于,現實世界本身就是不確定的,人類的觀察能力是有局限性的(否則有很大一部分科學就沒有必要做了——設想我們能夠直接觀察到電子的運行,還需要對原子模型爭吵不休嗎?),我們日常所觀察到的只是事物表面上的結果,沿用剛才那個袋子里面取球的比方,我們往往只能知道從里面取出來的球是什么顏色,而并不能直接看到袋子里面實際的情況。這個時候,我們就需要提供一個猜測(hypothesis,更為嚴格的說法是“假設”,這里用“猜測”更通俗易懂一點),所謂猜測,當然就是不確定的(很可能有好多種乃至無數種猜測都能滿足目前的觀測),但也絕對不是兩眼一抹黑瞎蒙——具體地說,我們需要做兩件事情:1. 算出各種不同猜測的可能性大小。2. 算出最靠譜的猜測是什么。第一個就是計算特定猜測的后驗概率,對于連續的猜測空間則是計算猜測的概率密度函數。第二個則是所謂的模型比較,模型比較如果不考慮先驗概率的話就是最大似然方法。

1.1 一個例子:自然語言的二義性

下面舉一個自然語言的不確定性的例子。當你看到這句話:

The girl saw the boy with a telescope.

你對這句話的含義有什么猜測?平常人肯定會說:那個女孩拿望遠鏡看見了那個男孩(即你對這個句子背后的實際語法結構的猜測是:The girl saw-with-a-telescope the boy )。然而,仔細一想,你會發現這個句子完全可以解釋成:那個女孩看見了那個拿著望遠鏡的男孩(即:The girl saw the-boy-with-a-telescope )。那為什么平常生活中我們每個人都能夠迅速地對這種二義性進行消解呢?這背后到底隱藏著什么樣的思維法則?我們留到后面解釋。

1.2 貝葉斯公式

貝葉斯公式是怎么來的?

我們還是使用 wikipedia 上的一個例子:

一所學校里面有 60% 的男生,40% 的女生。男生總是穿長褲,女生則一半穿長褲一半穿裙子。有了這些信息之后我們可以容易地計算“隨機選取一個學生,他(她)穿長褲的概率和穿裙子的概率是多大”,這個就是前面說的“正向概率”的計算。然而,假設你走在校園中,迎面走來一個穿長褲的學生(很不幸的是你高度近似,你只看得見他(她)穿的是否長褲,而無法確定他(她)的性別),你能夠推斷出他(她)是男生的概率是多大嗎?

一些認知科學的研究表明(《決策與判斷》以及《 Rationality for Mortals 》第12章:小孩也可以解決貝葉斯問題),我們對形式化的貝葉斯問題不擅長,但對于以頻率形式呈現的等價問題卻很擅長。在這里,我們不妨把問題重新敘述成:你在校園里面 隨機游走 ,遇到了 N 個穿長褲的人(仍然假設你無法直接觀察到他們的性別),問這 N 個人里面有多少個女生多少個男生。

你說,這還不簡單:算出學校里面有多少穿長褲的,然后在這些人里面再算出有多少女生,不就行了?

我們來算一算:假設學校里面人的總數是 U 個。60% 的男生都穿長褲,于是我們得到了 U * P(Boy) * P(Pants|Boy) 個穿長褲的(男生)(其中 P(Boy) 是男生的概率 = 60%,這里可以簡單的理解為男生的比例;P(Pants|Boy) 是條件概率,即在 Boy 這個條件下穿長褲的概率是多大,這里是 100% ,因為所有男生都穿長褲)。40% 的女生里面又有一半(50%)是穿長褲的,于是我們又得到了 U * P(Girl) * P(Pants|Girl) 個穿長褲的(女生)。加起來一共是 U * P(Boy) * P(Pants|Boy) + U * P(Girl) * P(Pants|Girl) 個穿長褲的,其中有 U * P(Girl) * P(Pants|Girl) 個女生。兩者一比就是你要求的答案。

下面我們把這個答案形式化一下:我們要求的是 P(Girl|Pants) (穿長褲的人里面有多少女生),我們計算的結果是 U * P(Girl) * P(Pants|Girl) / [U * P(Boy) * P(Pants|Boy) + U * P(Girl) * P(Pants|Girl)] 。容易發現這里校園內人的總數是無關的,可以消去。于是得到

P(Girl|Pants) = P(Girl) * P(Pants|Girl) / [P(Boy) * P(Pants|Boy) + P(Girl) * P(Pants|Girl)]

注意,如果把上式收縮起來,分母其實就是 P(Pants) ,分子其實就是 P(Pants, Girl) 。而這個比例很自然地就讀作:在穿長褲的人( P(Pants) )里面有多少(穿長褲)的女孩( P(Pants, Girl) )。

上式中的 Pants 和 Boy/Girl 可以指代一切東西,所以其一般形式就是:

P(B|A) = P(A|B) * P(B) / [P(A|B) * P(B) + P(A|~B) * P(~B) ]

收縮起來就是:

P(B|A) = P(AB) / P(A)

其實這個就等于:

P(B|A) * P(A) = P(AB)

難怪拉普拉斯說概率論只是把常識用數學公式表達了出來

然而,后面我們會逐漸發現,看似這么平凡的貝葉斯公式,背后卻隱含著非常深刻的原理。

2. 拼寫糾正

經典著作《人工智能:現代方法》的作者之一 Peter Norvig 曾經寫過一篇介紹如何寫一個拼寫檢查/糾正器的文章(原文在 這里 ,徐宥的翻譯版在 這里 ,這篇文章很深入淺出,強烈建議讀一讀),里面用到的就是貝葉斯方法,這里我們不打算復述他寫的文章,而是簡要地將其核心思想介紹一下。

首先,我們需要詢問的是:“問題是什么?

問題是我們看到用戶輸入了一個不在字典中的單詞,我們需要去猜測:“這個家伙到底真正想輸入的單詞是什么呢?”用剛才我們形式化的語言來敘述就是,我們需要求:

P(我們猜測他想輸入的單詞 | 他實際輸入的單詞)

這個概率。并找出那個使得這個概率最大的猜測單詞。顯然,我們的猜測未必是唯一的,就像前面舉的那個自然語言的歧義性的例子一樣;這里,比如用戶輸入: thew ,那么他到底是想輸入 the ,還是想輸入 thaw ?到底哪個猜測可能性更大呢?幸運的是我們可以用貝葉斯公式來直接出它們各自的概率,我們不妨將我們的多個猜測記為 h1 h2 .. ( h 代表 hypothesis),它們都屬于一個有限且離散的猜測空間 H (單詞總共就那么多而已),將用戶實際輸入的單詞記為 D ( D 代表 Data ,即觀測數據),于是

P(我們的猜測1 | 他實際輸入的單詞)

可以抽象地記為:

P(h1 | D)

類似地,對于我們的猜測2,則是 P(h2 | D)。不妨統一記為:

P(h | D)

運用一次貝葉斯公式,我們得到:

P(h | D) = P(h) * P(D | h) / P(D)

對于不同的具體猜測 h1 h2 h3 .. ,P(D) 都是一樣的,所以在比較 P(h1 | D) 和 P(h2 | D) 的時候我們可以忽略這個常數。即我們只需要知道:

P(h | D) ∝ P(h) * P(D | h) (注:那個符號的意思是“正比例于”,不是無窮大,注意符號右端是有一個小缺口的。)

這個式子的抽象含義是:對于給定觀測數據,一個猜測是好是壞,取決于“這個猜測本身獨立的可能性大小(先驗概率,Prior )”和“這個猜測生成我們觀測到的數據的可能性大小”(似然,Likelihood )的乘積。具體到我們的那個 thew 例子上,含義就是,用戶實際是想輸入 the 的可能性大小取決于 the 本身在詞匯表中被使用的可能性(頻繁程度)大小(先驗概率)和 想打 the 卻打成 thew 的可能性大小(似然)的乘積。

下面的事情就很簡單了,對于我們猜測為可能的每個單詞計算一下 P(h) * P(D | h) 這個值,然后取最大的,得到的就是最靠譜的猜測。

一點注記 :Norvig 的拼寫糾正器里面只提取了編輯距離為 2 以內的所有已知單詞。這是為了避免去遍歷字典中每個單詞計算它們的 P(h) * P(D | h) ,但這種做法為了節省時間帶來了一些誤差。但話說回來難道我們人類真的回去遍歷每個可能的單詞來計算他們的后驗概率嗎?不可能。實際上,根據認知神經科學的觀點,我們首先根據錯誤的單詞做一個 bottom-up 的關聯提取,提取出有可能是實際單詞的那些候選單詞,這個提取過程就是所謂的基于內容的提取,可以根據錯誤單詞的一些模式片段提取出有限的一組候選,非常快地縮小的搜索空間(比如我輸入 explaination ,單詞里面就有充分的信息使得我們的大腦在常數時間內把可能性 narrow down 到 explanation 這個單詞上,至于具體是根據哪些線索——如音節——來提取,又是如何在生物神經網絡中實現這個提取機制的,目前還是一個沒有弄清的領域)。然后,我們對這有限的幾個猜測做一個 top-down 的預測,看看到底哪個對于觀測數據(即錯誤單詞)的預測效力最好,而如何衡量預測效率則就是用貝葉斯公式里面的那個 P(h) * P(D | h) 了——雖然我們很可能使用了 一些啟發法來簡化計算 。后面我們還會提到這樣的 bottom-up 的關聯提取。

3. 模型比較與奧卡姆剃刀

3.1 再訪拼寫糾正

介紹了貝葉斯拼寫糾正之后,接下來的一個自然而然的問題就來了:“為什么? ”為什么要用貝葉斯公式?為什么貝葉斯公式在這里可以用?我們可以很容易地領會為什么貝葉斯公式用在前面介紹的那個男生女生長褲裙子的問題里是正確的。但為什么這里?

為了回答這個問題,一個常見的思路就是想想:非得這樣嗎? 因為如果你想到了另一種做法并且證明了它也是靠譜的,那么將它與現在這個一比較,也許就能得出很有價值的信息。那么對于拼寫糾錯問題你能想到其他方案嗎?

不管怎樣,一個最常見的替代方案就是,選擇離 thew 的 編輯距離 最近的。然而 the 和 thaw 離 thew 的編輯距離都是 1 。這可咋辦捏?你說,不慌,那還是好辦。我們就看到底哪個更可能被錯打為 thew 就是了。我們注意到字母 e 和字母 w 在鍵盤上離得很緊,無名指一抽筋就不小心多打出一個 w 來,the 就變成 thew 了。而另一方面 thaw 被錯打成 thew 的可能性就相對小一點,因為 e 和 a 離得較遠而且使用的指頭相差一個指頭(一個是中指一個是小指,不像 e 和 w 使用的指頭靠在一塊——神經科學的證據表明緊鄰的身體設施之間容易串位)。OK,很好,因為你現在已經是在用最大似然方法了,或者直白一點,你就是在計算那個使得 P(D | h) 最大的 h 。

而貝葉斯方法計算的是什么?是 P(h) * P(D | h) 。多出來了一個 P(h) 。我們剛才說了,這個多出來的 P(h) 是特定猜測的先驗概率。為什么要摻和進一個先驗概率?剛才說的那個最大似然不是挺好么?很雄辯地指出了 the 是更靠譜的猜測。有什么問題呢?既然這樣,我們就從給最大似然找茬開始吧——我們假設兩者的似然程度是一樣或非常相近,這樣不就難以區分哪個猜測更靠譜了嗎?比如用戶輸入tlp ,那到底是 top 還是 tip ?(這個例子不怎么好,因為 top 和 tip 的詞頻可能仍然是接近的,但一時想不到好的英文單詞的例子,我們不妨就假設 top 比 tip 常見許多吧,這個假設并不影響問題的本質。)這個時候,當最大似然不能作出決定性的判斷時,先驗概率就可以插手進來給出指示——“既然你無法決定,那么我告訴你,一般來說 top 出現的程度要高許多,所以更可能他想打的是 top ”)。

以上只是最大似然的一個問題,即并不能提供決策的全部信息。

最大似然還有另一個問題:即便一個猜測與數據非常符合,也并不代表這個猜測就是更好的猜測,因為這個猜測本身的可能性也許就非常低。比如 MacKay 在《Information Theory : Inference and Learning Algorithms》里面就舉了一個很好的例子:-1 3 7 11 你說是等差數列更有可能呢?還是 -X^3 / 11 + 9/11*X^2 + 23/11 每項把前項作為 X 帶入后計算得到的數列?此外曲線擬合也是,平面上 N 個點總是可以用 N-1 階多項式來完全擬合,當 N 個點近似但不精確共線的時候,用 N-1 階多項式來擬合能夠精確通過每一個點,然而用直線來做擬合/線性回歸的時候卻會使得某些點不能位于直線上。你說到底哪個好呢?多項式?還是直線?一般地說肯定是越低階的多項式越靠譜(當然前提是也不能忽視“似然”P(D | h) ,明擺著一個多項式分布您愣是去拿直線擬合也是不靠譜的,這就是為什么要把它們兩者乘起來考慮。),原因之一就是低階多項式更常見,先驗概率( P(h) )較大(原因之二則隱藏在 P(D | h) 里面),這就是為什么我們要用 樣條 來插值,而不是直接搞一個 N-1 階多項式來通過任意 N 個點的原因。

以上分析當中隱含的哲學是,觀測數據總是會有各種各樣的誤差,比如觀測誤差(比如你觀測的時候一個 MM 經過你一不留神,手一抖就是一個誤差出現了),所以如果過分去尋求能夠完美解釋觀測數據的模型,就會落入所謂的數據 過配(overfitting) 的境地,一個過配的模型試圖連誤差(噪音)都去解釋(而實際上噪音又是不需要解釋的),顯然就過猶不及了。所以 P(D | h) 大不代表你的 h (猜測)就是更好的 h。還要看 P(h) 是怎樣的。所謂 奧卡姆剃刀 精神就是說:如果兩個理論具有相似的解釋力度,那么優先選擇那個更簡單的(往往也正是更平凡的,更少繁復的,更常見的)。

過分匹配的另一個原因在于當觀測的結果并不是因為誤差而顯得“不精確”而是因為真實世界中對數據的結果產生貢獻的因素太多太多,跟噪音不同,這些偏差是一些另外的因素集體貢獻的結果,不是你的模型所能解釋的——噪音那是不需要解釋——一個現實的模型往往只提取出幾個與結果相關度很高,很重要的因素(cause)。這個時候觀察數據會傾向于圍繞你的有限模型的預測結果呈 正態分布 ,于是你實際觀察到的結果就是這個正態分布的 隨機取樣 ,這個取樣很可能受到其余因素的影響偏離你的模型所預測的中心,這個時候便不能貪心不足地試圖通過改變模型來“完美”匹配數據,因為那些使結果偏離你的預測的貢獻因素不是你這個有限模型里面含有的因素所能概括的,硬要打腫臉充胖子只能導致不實際的模型,舉個教科書例子:身高和體重的實際關系近似于一個二階多項式的關系,但大家都知道并不是只有身高才會對體重產生影響,物理世界影響體重的因素太多太多了,有人身材高大卻瘦得跟稻草,有人卻是橫長豎不長。但不可否認的是總體上來說,那些特殊情況越是特殊就越是稀少,呈圍繞最普遍情況(胖瘦適中)的正態分布,這個分布就保證了我們的身高——體重相關模型能夠在大多數情況下做出靠譜的預測。但是——剛才說了,特例是存在的,就算不是特例,人有胖瘦,密度也有大小,所以完美符合身高——體重的某個假想的二階多項式關系的人是不存在的,我們又不是歐幾里德幾何世界當中的理想多面體,所以,當我們對人群隨機抽取了 N 個樣本(數據點)試圖對這 N 個數據點擬合出一個多項式的話就得注意,它肯定得是二階多項式,我們要做的只是去根據數據點計算出多項式各項的參數(一個典型的方法就是最小二乘);它肯定不是直線(我們又不是稻草),也不是三階多項式四階多項式.. 如果硬要完美擬合 N 個點,你可能會整出一個 N-1 階多項式來——設想身高和體重的關系是 5 階多項式看看?

3.2 模型比較理論(Model Comparasion)與貝葉斯奧卡姆剃刀(Bayesian Occam’s Razor)

實際上,模型比較就是去比較哪個模型(猜測)更可能隱藏在觀察數據的背后。其基本思想前面已經用拼寫糾正的例子來說明了。我們對用戶實際想輸入的單詞的猜測就是模型,用戶輸錯的單詞就是觀測數據。我們通過:

P(h | D) ∝ P(h) * P(D | h)

來比較哪個模型最為靠譜。前面提到,光靠 P(D | h) (即“似然”)是不夠的,有時候還需要引入 P(h) 這個先驗概率。奧卡姆剃刀就是說 P(h) 較大的模型有較大的優勢,而最大似然則是說最符合觀測數據的(即 P(D | h) 最大的)最有優勢。整個模型比較就是這兩方力量的拉鋸。我們不妨再舉一個簡單的例子來說明這一精神:你隨便找枚硬幣,擲一下,觀察一下結果。好,你觀察到的結果要么是“正”,要么是“反”(不,不是少林足球那枚硬幣:P ),不妨假設你觀察到的是“正”。現在你要去根據這個觀測數據推斷這枚硬幣擲出“正”的概率是多大。根據最大似然估計的精神,我們應該猜測這枚硬幣擲出“正”的概率是 1 ,因為這個才是能最大化 P(D | h) 的那個猜測。然而每個人都會大搖其頭——很顯然,你隨機摸出一枚硬幣這枚硬幣居然沒有反面的概率是“不存在的”,我們對一枚隨機硬幣是否一枚有偏硬幣,偏了多少,是有著一個先驗的認識的,這個認識就是絕大多數硬幣都是基本公平的,偏得越多的硬幣越少見(可以用一個 beta 分布 來表達這一先驗概率)。將這個先驗正態分布 p(θ) (其中 θ 表示硬幣擲出正面的比例,小寫的 p 代表這是 概率密度函數 )結合到我們的問題中,我們便不是去最大化 P(D | h) ,而是去最大化 P(D | θ) * p(θ) ,顯然 θ = 1 是不行的,因為 P(θ=1) 為 0 ,導致整個乘積也為 0 。實際上,只要對這個式子求一個導數就可以得到最值點。

以上說的是當我們知道先驗概率 P(h) 的時候,光用最大似然是不靠譜的,因為最大似然的猜測可能先驗概率非常小。然而,有些時候,我們對于先驗概率一無所知,只能假設每種猜測的先驗概率是均等的,這個時候就只有用最大似然了。實際上,統計學家和貝葉斯學家有一個有趣的爭論,統計學家說:我們讓數據自己說話。言下之意就是要摒棄先驗概率。而貝葉斯支持者則說:數據會有各種各樣的偏差,而一個靠譜的先驗概率則可以對這些隨機噪音做到健壯。事實證明貝葉斯派勝利了,勝利的關鍵在于所謂先驗概率其實也是經驗統計的結果,譬如為什么我們會認為絕大多數硬幣是基本公平的?為什么我們認為大多數人的肥胖適中?為什么我們認為膚色是種族相關的,而體重則與種族無關?先驗概率里面的“先驗”并不是指先于一切經驗,而是僅指先于我們“當前”給出的觀測數據而已,在硬幣的例子中先驗指的只是先于我們知道投擲的結果這個經驗,而并非“先天”。

然而,話說回來,有時候我們必須得承認,就算是基于以往的經驗,我們手頭的“先驗”概率還是均勻分布,這個時候就必須依賴用最大似然,我們用前面留下的一個自然語言二義性問題來說明這一點:

The girl saw the boy with a telescope.

到底是 The girl saw-with-a-telescope the boy 這一語法結構,還是 The girl saw the-boy-with-a-telescope 呢?兩種語法結構的常見程度都差不多(你可能會覺得后一種語法結構的常見程度較低,這是事后偏見,你只需想想 The girl saw the boy with a book 就知道了。當然,實際上從大規模語料統計結果來看后一種語法結構的確稍稍不常見一丁點,但是絕對不足以解釋我們對第一種結構的強烈傾向)。那么到底為什么呢?

我們不妨先來看看 MacKay 在書中舉的一個漂亮的例子:

圖中有多少個箱子?特別地,那棵書后面是一個箱子?還是兩個箱子?還是三個箱子?還是.. 你可能會覺得樹后面肯定是一個箱子,但為什么不是兩個呢?如下圖:

很簡單,你會說:要是真的有兩個箱子那才怪了,怎么就那么巧這兩個箱子剛剛好顏色相同,高度相同呢?

用概率論的語言來說,你剛才的話就翻譯為:猜測 h 不成立,因為 P(D | h) 太小(太巧合)了。我們的直覺是:巧合(小概率)事件不會發生。所以當一個猜測(假設)使得我們的觀測結果成為小概率事件的時候,我們就說“才怪呢,哪能那么巧捏?!”

現在我們可以回到那個自然語言二義性的例子,并給出一個完美的解釋了:如果語法結構是 The girl saw the-boy-with-a-telecope 的話,怎么那個男孩偏偏手里拿的就是望遠鏡——一個可以被用來 saw-with 的東東捏?這也忒小概率了吧。他咋就不會拿本書呢?拿什么都好。怎么偏偏就拿了望遠鏡?所以唯一的解釋是,這個“巧合”背后肯定有它的必然性,這個必然性就是,如果我們將語法結構解釋為 The girl saw-with-a-telescope the boy 的話,就跟數據完美吻合了——既然那個女孩是用某個東西去看這個男孩的,那么這個東西是一個望遠鏡就完全可以解釋了(不再是小概率事件了)。

自然語言二義性很常見,譬如上文中的一句話:

參見《決策與判斷》以及《 Rationality for Mortals 》第12章:小孩也可以解決貝葉斯問題

就有二義性:到底是參見這兩本書的第 12 章,還是僅僅是第二本書的第 12 章呢?如果是這兩本書的第 12 章那就是咄咄怪事了,怎么恰好兩本書都有第 12 章,都是講同一個問題,更詭異的是,標題還相同呢?

注意,以上做的是似然估計(即只看 P(D | h) 的大小),不含先驗概率。通過這兩個例子,尤其是那個樹后面的箱子的例子我們可以看到,似然估計里面也蘊含著奧卡姆剃刀:樹后面的箱子數目越多,這個模型就越復雜。單個箱子的模型是最簡單的。似然估計選擇了更簡單的模型。

這個就是所謂的貝葉斯奧卡姆剃刀(Bayesian Occam’s Razor) ,因為這個剃刀工作在貝葉斯公式的似然(P(D | h) )上,而不是模型本身( P(h) )的先驗概率上,后者是傳統的奧卡姆剃刀。關于貝葉斯奧卡姆剃刀我們再來看一個前面說到的曲線擬合的例子:如果平面上有 N 個點,近似構成一條直線,但絕不精確地位于一條直線上。這時我們既可以用直線來擬合(模型1),也可以用二階多項式(模型2)擬合,也可以用三階多項式(模型3),.. ,特別地,用 N-1 階多項式便能夠保證肯定能完美通過 N 個數據點。那么,這些可能的模型之中到底哪個是最靠譜的呢?前面提到,一個衡量的依據是奧卡姆剃刀:越是高階的多項式越是繁復和不常見。然而,我們其實并不需要依賴于這個先驗的奧卡姆剃刀,因為有人可能會爭辯說:你怎么就能說越高階的多項式越不常見呢?我偏偏覺得所有階多項式都是等可能的。好吧,既然如此那我們不妨就扔掉 P(h) 項,看看 P(D | h) 能告訴我們什么。我們注意到越是高階的多項式,它的軌跡彎曲程度越是大,到了八九階簡直就是直上直下,于是我們不僅要問:一個比如說八階多項式在平面上隨機生成的一堆 N 個點偏偏恰好近似構成一條直線的概率(即 P(D | h) )有多大?太小太小了。反之,如果背后的模型是一條直線,那么根據該模型生成一堆近似構成直線的點的概率就大得多了。這就是貝葉斯奧卡姆剃刀。

這里只是提供一個關于貝葉斯奧卡姆剃刀的科普,強調直觀解釋,更多理論公式請參考 MacKay 的著作 《Information Theory : Inference and Learning Algorithms》第 28 章。

3.3 最小描述長度原則

貝葉斯模型比較理論與信息論有一個有趣的關聯:

P(h | D) ∝ P(h) * P(D | h)

兩邊求對數,將右式的乘積變成相加:

ln P(h | D) ∝ ln P(h) + ln P(D | h)

顯然,最大化 P(h | D) 也就是最大化 ln P(h | D)。而 ln P(h) + ln P(D | h) 則可以解釋為模型(或者稱“假設”、“猜測”)h 的編碼長度加上在該模型下數據 D 的編碼長度。使這個和最小的模型就是最佳模型。

而究竟如何定義一個模型的編碼長度,以及數據在模型下的編碼長度則是一個問題。更多可參考 Mitchell 的 《Machine Learning》的 6.6 節,或 Mackay 的 28.3 節)

3.4 最優貝葉斯推理

所謂的推理,分為兩個過程,第一步是對觀測數據建立一個模型。第二步則是使用這個模型來推測未知現象發生的概率。我們前面都是講的對于觀測數據給出最靠譜的那個模型。然而很多時候,雖然某個模型是所有模型里面最靠譜的,但是別的模型也并不是一點機會都沒有。譬如第一個模型在觀測數據下的概率是 0.5 。第二個模型是 0.4 ,第三個是 0.1 。如果我們只想知道對于觀測數據哪個模型最可能,那么只要取第一個就行了,故事到此結束。然而很多時候我們建立模型是為了推測未知的事情的發生概率,這個時候,三個模型對未知的事情發生的概率都會有自己的預測,僅僅因為某一個模型概率稍大一點就只聽他一個人的就太不民主了。所謂的最優貝葉斯推理就是將三個模型對于未知數據的預測結論加權平均起來(權值就是模型相應的概率)。顯然,這個推理是理論上的制高點,無法再優了,因為它已經把所有可能性都考慮進去了。

只不過實際上我們是基本不會使用這個框架的,因為計算模型可能非常費時間,二來模型空間可能是連續的,即有無窮多個模型(這個時候需要計算模型的概率分布)。結果還是非常費時間。所以這個被看作是一個理論基準。

4. 無處不在的貝葉斯

以下我們再舉一些實際例子來說明貝葉斯方法被運用的普遍性,這里主要集中在機器學習方面,因為我不是學經濟的,否則還可以找到一堆經濟學的例子。

4.1 中文分詞

貝葉斯是機器學習的核心方法之一。比如中文分詞領域就用到了貝葉斯。Google 研究員吳軍在《數學之美》系列中就有一篇是介紹中文分詞的,這里只介紹一下核心的思想,不做贅述,詳細請參考吳軍的文章( 這里 )。

分詞問題的描述為:給定一個句子(字串),如:

南京市長江大橋

如何對這個句子進行分詞(詞串)才是最靠譜的。例如:

1. 南京市/長江大橋

2. 南京/市長/江大橋

這兩個分詞,到底哪個更靠譜呢?

我們用貝葉斯公式來形式化地描述這個問題,令 X 為字串(句子),Y 為詞串(一種特定的分詞假設)。我們就是需要尋找使得 P(Y|X) 最大的 Y ,使用一次貝葉斯可得:

P(Y|X) ∝ P(Y)*P(X|Y)

用自然語言來說就是 這種分詞方式(詞串)的可能性 乘以 這個詞串生成我們的句子的可能性。我們進一步容易看到:可以近似地將 P(X|Y) 看作是恒等于 1 的,因為任意假想的一種分詞方式之下生成我們的句子總是精準地生成的(只需把分詞之間的分界符號扔掉即可)。于是,我們就變成了去最大化 P(Y) ,也就是尋找一種分詞使得這個詞串(句子)的概率最大化。而如何計算一個詞串:

W1, W2, W3, W4 ..

的可能性呢?我們知道,根據聯合概率 的公式展開:P(W1, W2, W3, W4 ..) = P(W1) * P(W2|W1) * P(W3|W2, W1) * P(W4|W1,W2,W3) * .. 于是我們可以通過一系列的條件概率(右式)的乘積來求整個聯合概率。然而不幸的是隨著條件數目的增加(P(Wn|Wn-1,Wn-2,..,W1) 的條件有 n-1 個), 數據稀疏問題 也會越來越嚴重,即便語料庫再大也無法統計出一個靠譜的 P(Wn|Wn-1,Wn-2,..,W1) 來。為了緩解這個問題,計算機科學家們一如既往地使用了“天真”假設:我們假設句子中一個詞的出現概率只依賴于它前面的有限的 k 個詞(k 一般不超過 3,如果只依賴于前面的一個詞,就是2元 語言模型 (2- gram),同理有 3-gram 、 4-gram 等),這個就是所謂的“有限地平線”假設。雖然這個假設很傻很天真,但結果卻表明它的結果往往是很好很強大的,后面要提到的樸素貝葉斯方法使用的假設跟這個精神上是完全一致的,我們會解釋為什么像這樣一個天真的假設能夠得到強大的結果。目前我們只要知道,有了這個假設,剛才那個乘積就可以改寫成: P(W1) * P(W2|W1) * P(W3|W2) * P(W4|W3) .. (假設每個詞只依賴于它前面的一個詞)。而統計 P(W2|W1) 就不再受到數據稀疏問題的困擾了。對于我們上面提到的例子“南京市長江大橋”,如果按照自左到右的貪婪方法分詞的話,結果就成了“南京市長/江大橋”。但如果按照貝葉斯分詞的話(假設使用 3-gram),由于“南京市長”和“江大橋”在語料庫中一起出現的頻率為 0 ,這個整句的概率便會被判定為 0 。 從而使得“南京市/長江大橋”這一分詞方式勝出。

一點注記 :有人可能會疑惑,難道我們人類也是基于這些天真的假設來進行推理的?不是的。事實上,統計機器學習方法所統計的東西往往處于相當表層(shallow)的層面,在這個層面機器學習只能看到一些非常表面的現象,有一點科學研究的理念的人都知道:越是往表層去,世界就越是繁復多變。從機器學習的角度來說,特征(feature)就越多,成百上千維度都是可能的。特征一多,好了, 高維詛咒 就產生了,數據就稀疏得要命,不夠用了。而我們人類的觀察水平顯然比機器學習的觀察水平要更深入一些,為了避免數據稀疏我們不斷地發明各種裝置(最典型就是顯微鏡),來幫助我們直接深入到更深層的事物層面去觀察更本質的聯系,而不是在淺層對表面現象作統計歸納。舉一個簡單的例子,通過對大規模語料庫的統計,機器學習可能會發現這樣一個規律:所有的“他”都是不會穿 bra 的,所有的“她”則都是穿的。然而,作為一個男人,卻完全無需進行任何統計學習,因為深層的規律就決定了我們根本不會去穿 bra 。至于機器學習能不能完成后者(像人類那樣的)這個推理,則是人工智能領域的經典問題。至少在那之前, 聲稱統計學習方法能夠終結科學研究 原文 )的說法 是純粹外行人說的話

4.2 統計機器翻譯

統計機器翻譯因為其簡單,自動(無需手動添加規則),迅速成為了機器翻譯的事實標準。而統計機器翻譯的核心算法也是使用的貝葉斯方法。

問題是什么?統計機器翻譯的問題可以描述為:給定一個句子 e ,它的可能的外文翻譯 f 中哪個是最靠譜的。即我們需要計算:P(f|e) 。一旦出現條件概率貝葉斯總是挺身而出:

P(f|e) ∝ P(f) * P(e|f)

這個式子的右端很容易解釋:那些先驗概率較高,并且更可能生成句子 e 的外文句子 f 將會勝出。我們只需簡單統計(結合上面提到的 N-Gram 語言模型)就可以統計任意一個外文句子 f 的出現概率。然而 P(e|f) 卻不是那么好求的,給定一個候選的外文局子 f ,它生成(或對應)句子 e 的概率是多大呢?我們需要定義什么叫 “對應”,這里需要用到一個分詞對齊的平行語料庫,有興趣的可以參考 《Foundations of Statistical Natural Language Processing》第 13 章,這里摘選其中的一個例子:假設 e 為:John loves Mary 。我們需要考察的首選 f 是:Jean aime Marie (法文)。我們需要求出 P(e|f) 是多大,為此我們考慮 e 和 f 有多少種對齊的可能性,如:

John (Jean) loves (aime) Marie (Mary)

就是其中的一種(最靠譜的)對齊,為什么要對齊,是因為一旦對齊了之后,就可以容易地計算在這個對齊之下的 P(e|f) 是多大,只需計算:

P(John|Jean) * P(loves|aime) * P(Marie|Mary)

即可。

然后我們遍歷所有的對齊方式,并將每種對齊方式之下的翻譯概率 ∑ 求和。便可以獲得整個的 P(e|f) 是多大。

一點注記 :還是那個問題:難道我們人類真的是用這種方式進行翻譯的?highly unlikely 。這種計算復雜性非常高的東西連三位數乘法都搞不定的我們才不會笨到去使用呢。根據認知神經科學的認識,很可能我們是先從句子到語義(一個逐層往上(bottom-up)抽象的 folding 過程),然后從語義根據另一門語言的語法展開為另一門語言(一個逐層往下(top-down)的具體化 unfolding 過程)。如何可計算地實現這個過程,目前仍然是個難題。(我們看到很多地方都有 bottom-up/top-down 這樣一個對稱的過程,實際上有人猜測這正是生物神經網絡原則上的運作方式,對視覺神經系統的研究尤其證明了這一點,Hawkins 在 《On Intelligence》 里面提出了一種 HTM (Hierarchical Temporal Memory)模型正是使用了這個原則。)

4.3 貝葉斯圖像識別,Analysis by Synthesis

貝葉斯方法是一個非常 general 的推理框架。其核心理念可以描述成:Analysis by Synthesis (通過合成來分析)。06 年的認知科學新進展上有一篇 paper 就是講用貝葉斯推理來解釋視覺識別的,一圖勝千言,下圖就是摘自這篇 paper :

首先是視覺系統提取圖形的邊角特征,然后使用這些特征自底向上地激活高層的抽象概念(比如是 E 還是 F 還是等號),然后使用一個自頂向下的驗證來比較到底哪個概念最佳地解釋了觀察到的圖像。

4.4? EM 算法與基于模型的聚類

聚類 是一種 無指導的機器學習 問題,問題描述:給你一堆數據點,讓你將它們最靠譜地分成一堆一堆的。聚類算法很多,不同的算法適應于不同的問題,這里僅介紹一個基于模型的聚類,該聚類算法對數據點的假設是,這些數據點分別是圍繞 K 個核心的 K 個正態分布源所隨機生成的,使用 Han JiaWei 的《Data Ming: Concepts and Techniques》中的圖:

圖中有兩個正態分布核心,生成了大致兩堆點。我們的聚類算法就是需要根據給出來的那些點,算出這兩個正態分布的核心在什么位置,以及分布的參數是多少。這很明顯又是一個貝葉斯問題,但這次不同的是,答案是連續的且有無窮多種可能性,更糟的是,只有當我們知道了哪些點屬于同一個正態分布圈的時候才能夠對這個分布的參數作出靠譜的預測,現在兩堆點混在一塊我們又不知道哪些點屬于第一個正態分布,哪些屬于第二個。反過來,只有當我們對分布的參數作出了靠譜的預測時候,才能知道到底哪些點屬于第一個分布,那些點屬于第二個分布。這就成了一個先有雞還是先有蛋的問題了。為了解決這個循環依賴,總有一方要先打破僵局,說,不管了,我先隨便整一個值出來,看你怎么變,然后我再根據你的變化調整我的變化,然后如此迭代著不斷互相推導,最終收斂到一個解。這就是 EM 算法。

EM 的意思是“Expectation-Maximazation”,在這個聚類問題里面,我們是先隨便猜一下這兩個正態分布的參數:如核心在什么地方,方差是多少。然后計算出每個數據點更可能屬于第一個還是第二個正態分布圈,這個是屬于 Expectation 一步。有了每個數據點的歸屬,我們就可以根據屬于第一個分布的數據點來重新評估第一個分布的參數(從蛋再回到雞),這個是 Maximazation 。如此往復,直到參數基本不再發生變化為止。這個迭代收斂過程中的貝葉斯方法在第二步,根據數據點求分布的參數上面。

4.5 最大似然與最小二乘

學過線性代數的大概都知道經典的最小二乘方法來做線性回歸。問題描述是:給定平面上 N 個點,(這里不妨假設我們想用一條直線來擬合這些點——回歸 可以看作是擬合 的特例,即允許誤差的擬合),找出一條最佳描述了這些點的直線。

一個接踵而來的問題就是,我們如何定義最佳?我們設每個點的坐標為 (Xi, Yi) 。如果直線為 y = f(x) 。那么 (Xi, Yi) 跟直線對這個點的“預測”:(Xi, f(Xi)) 就相差了一個 ΔYi = |Yi – f(Xi)| 。最小二乘就是說尋找直線使得 (ΔY1)^2 + (ΔY2)^2 + .. (即誤差的平方和)最小,至于為什么是誤差的平方和而不是誤差的絕對值和,統計學上也沒有什么好的解釋。然而貝葉斯方法卻能對此提供一個完美的解釋。

我們假設直線對于坐標 Xi 給出的預測 f(Xi) 是最靠譜的預測,所有縱坐標偏離 f(Xi) 的那些數據點都含有噪音,是噪音使得它們偏離了完美的一條直線,一個合理的假設就是偏離路線越遠的概率越小,具體小多少,可以用一個正態分布曲線來模擬,這個分布曲線以直線對 Xi 給出的預測 f(Xi) 為中心,實際縱坐標為 Yi 的點 (Xi, Yi) 發生的概率就正比于 EXP[-(ΔYi)^2]。(EXP(..) 代表以常數 e 為底的多少次方)。

現在我們回到問題的貝葉斯方面,我們要想最大化的后驗概率是:

P(h|D) ∝ P(h) * P(D|h)

又見貝葉斯!這里 h 就是指一條特定的直線,D 就是指這 N 個數據點。我們需要尋找一條直線 h 使得 P(h) * P(D|h) 最大。很顯然,P(h) 這個先驗概率是均勻的,因為哪條直線也不比另一條更優越。所以我們只需要看 P(D|h) 這一項,這一項是指這條直線生成這些數據點的概率,剛才說過了,生成數據點 (Xi, Yi) 的概率為 EXP[-(ΔYi)^2] 乘以一個常數。而 P(D|h) = P(d1|h) * P(d2|h) * .. 即假設各個數據點是獨立生成的,所以可以把每個概率乘起來。于是生成 N 個數據點的概率為 EXP[-(ΔY1)^2] * EXP[-(ΔY2)^2] * EXP[-(ΔY3)^2] * .. = EXP{-[(ΔY1)^2 + (ΔY2)^2 + (ΔY3)^2 + .. ]} 最大化這個概率就是要最小化 (ΔY1)^2 + (ΔY2)^2 + (ΔY3)^2 + .. 。 熟悉這個式子嗎?

5. 樸素貝葉斯方法

樸素貝葉斯方法是一個很特別的方法,所以值得介紹一下。我們用樸素貝葉斯在垃圾郵件過濾中的應用來舉例說明。

5.1 貝葉斯垃圾郵件過濾器

問題是什么?問題是,給定一封郵件,判定它是否屬于垃圾郵件。按照先例,我們還是用 D 來表示這封郵件,注意 D 由 N 個單詞組成。我們用 h+ 來表示垃圾郵件,h- 表示正常郵件。問題可以形式化地描述為求:

P(h+|D) = P(h+) * P(D|h+) / P(D)

P(h-|D) = P(h-) * P(D|h-) / P(D)

其中 P(h+) 和 P(h-) 這兩個先驗概率都是很容易求出來的,只需要計算一個郵件庫里面垃圾郵件和正常郵件的比例就行了。然而 P(D|h+) 卻不容易求,因為 D 里面含有 N 個單詞 d1, d2, d3, .. ,所以P(D|h+) = P(d1,d2,..,dn|h+) 。我們又一次遇到了數據稀疏性,為什么這么說呢?P(d1,d2,..,dn|h+) 就是說在垃圾郵件當中出現跟我們目前這封郵件一模一樣的一封郵件的概率是多大!開玩笑,每封郵件都是不同的,世界上有無窮多封郵件。瞧,這就是數據稀疏性,因為可以肯定地說,你收集的訓練數據庫不管里面含了多少封郵件,也不可能找出一封跟目前這封一模一樣的。結果呢?我們又該如何來計算 P(d1,d2,..,dn|h+) 呢?

我們將 P(d1,d2,..,dn|h+)? 擴展為: P(d1|h+) * P(d2|d1, h+) * P(d3|d2,d1, h+) * .. 。熟悉這個式子嗎?這里我們會使用一個更激進的假設,我們假設 di 與 di-1 是完全條件無關的,于是式子就簡化為 P(d1|h+) * P(d2|h+) * P(d3|h+) * .. 。這個就是所謂的 條件獨立假設 ,也正是樸素貝葉斯方法的樸素之處。而計算 P(d1|h+) * P(d2|h+) * P(d3|h+) * .. 就太簡單了,只要統計 di 這個單詞在垃圾郵件中出現的頻率即可。關于貝葉斯垃圾郵件過濾更多的內容可以參考 這個條目 ,注意其中提到的其他資料。

一點注記 :這里,為什么有這個數據稀疏問題,還是因為統計學習方法工作在淺層面, 世界上的單詞就算不再變多也是非常之多的,單詞之間組成的句子也是變化多端,更不用說一篇文章了,文章數目則是無窮的,所以在這個層面作統計,肯定要被數據稀疏性困擾。我們要注意,雖然句子和文章的數目是無限的,然而就拿郵件來說,如果我們只關心郵件中句子的語義(進而更高抽象層面的“意圖”(語義,意圖如何可計算地定義出來是一個人工智能問題),在這個層面上可能性便大大縮減了,我們關心的抽象層面越高,可能性越小。單詞集合和句子的對應是多對一的,句子和語義的對應又是多對一的,語義和意圖的對應還是多對一的,這是個層級體系。神經科學的發現也表明大腦的皮層大致有一種層級結構,對應著越來越抽象的各個層面,至于如何具體實現一個可放在計算機內的大腦皮層,仍然是一個未解決問題,以上只是一個原則(principle)上的認識,只有當 computational 的 cortex 模型被建立起來了之后才可能將其放入電腦。

5.2 為什么樸素貝葉斯方法令人詫異地好——一個理論解釋

樸素貝葉斯方法的條件獨立假設看上去很傻很天真,為什么結果卻很好很強大呢?就拿一個句子來說,我們怎么能魯莽地聲稱其中任意一個單詞出現的概率只受到它前面的 3 個或 4 個單詞的影響呢?別說 3 個,有時候一個單詞的概率受到上一句話的影響都是絕對可能的。那么為什么這個假設在實際中的表現卻不比決策樹差呢?有人對此提出了一個理論解釋,并且建立了什么時候樸素貝葉斯的效果能夠等價于非樸素貝葉斯的充要條件,這個解釋的核心就是:有些獨立假設在各個分類之間的分布都是均勻的所以對于似然的相對大小不產生影響;即便不是如此,也有很大的可能性各個獨立假設所產生的消極影響或積極影響互相抵消,最終導致結果受到的影響不大。具體的數學公式請參考 這篇 paper

6. 層級貝葉斯模型

層級貝葉斯模型 是現代貝葉斯方法的標志性建筑之一。前面講的貝葉斯,都是在同一個事物層次上的各個因素之間進行統計推理,然而層次貝葉斯模型在哲學上更深入了一層,將這些因素背后的因素(原因的原因,原因的原因,以此類推)囊括進來。一個教科書例子是:如果你手頭有 N 枚硬幣,它們是同一個工廠鑄出來的,你把每一枚硬幣擲出一個結果,然后基于這 N 個結果對這 N 個硬幣的 θ (出現正面的比例)進行推理。如果根據最大似然,每個硬幣的 θ 不是 1 就是 0 (這個前面提到過的),然而我們又知道每個硬幣的 p(θ) 是有一個先驗概率的,也許是一個 beta 分布。也就是說,每個硬幣的實際投擲結果 Xi 服從以 θ 為中心的正態分布,而 θ 又服從另一個以 Ψ 為中心的 beta 分布。層層因果關系就體現出來了。進而 Ψ 還可能依賴于因果鏈上更上層的因素,以此類推。

6.1 隱馬可夫模型(HMM)

?

吳軍在數學之美系列里面介紹的 隱馬可夫模型 (HMM)就是一個簡單的層級貝葉斯模型:

那么怎么根據接收到的信息來推測說話者想表達的意思呢?我們可以利用叫做“隱含馬爾可夫模型”(Hidden Markov Model)來解決這些問題。以語音識別為例,當我們觀測到語音信號 o1,o2,o3 時,我們要根據這組信號推測出發送的句子 s1,s2,s3。顯然,我們應該在所有可能的句子中找最有可能性的一個。用數學語言來描述,就是在已知 o1,o2,o3,...的情況下,求使得條件概率 P (s1,s2,s3,...|o1,o2,o3....) 達到最大值的那個句子 s1,s2,s3,...

吳軍的文章中這里省掉沒說的是,s1, s2, s3, .. 這個句子的生成概率同時又取決于一組參數,這組參數決定了 s1, s2, s3, .. 這個馬可夫鏈的先驗生成概率。如果我們將這組參數記為 λ ,我們實際上要求的是:P(S|O, λ) (其中 O 表示 o1,o2,o3,.. ,S表示 s1,s2,s3,..)

當然,上面的概率不容易直接求出,于是我們可以間接地計算它。利用貝葉斯公式并且省掉一個常數項,可以把上述公式等價變換成

P(o1,o2,o3,...|s1,s2,s3....) * P(s1,s2,s3,...)

其中

P(o1,o2,o3,...|s1,s2,s3....) 表示某句話 s1,s2,s3...被讀成 o1,o2,o3,...的可能性, 而 P(s1,s2,s3,...) 表示字串 s1,s2,s3,...本身能夠成為一個合乎情理的句子的可能性,所以這個公式的意義是用發送信號為 s1,s2,s3...這個數列的可能性乘以 s1,s2,s3.. 本身可以一個句子的可能性,得出概率。

這里,s1,s2,s3...本身可以一個句子的可能性其實就取決于參數 λ ,也就是語言模型。所以簡而言之就是發出的語音信號取決于背后實際想發出的句子,而背后實際想發出的句子本身的獨立先驗概率又取決于語言模型。

7. 貝葉斯網絡

吳軍已經對貝葉斯網絡作了科普,請直接跳轉到 這里 。更詳細的理論參考所有機器學習的書上都有。

參考資料:

一堆機器學習,一堆概率統計,一堆 Google ,和一堆 Wikipedia 條目,一堆 paper

轉載于:https://my.oschina.net/dake/blog/196748

總結

以上是生活随笔為你收集整理的转:数学之美番外篇:平凡而又神奇的贝叶斯方法 收藏的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。

在线看欧美 | 日韩区在线观看 | 日产乱码一二三区别免费 | 久久久久久久国产精品影院 | 国产一级二级三级视频 | av高清一区二区三区 | 日韩女同一区二区三区在线观看 | 久久久久9999亚洲精品 | 黄色毛片大全 | 99久高清在线观看视频99精品热在线观看视频 | 婷婷色在线播放 | 成人在线视频网 | 精品国产欧美一区二区三区不卡 | 久久久久伊人 | 在线天堂中文在线资源网 | 麻豆视频免费在线 | 色偷偷88888欧美精品久久 | 中文字幕国产在线 | 奇米影视8888在线观看大全免费 | 人人模人人爽 | av电影免费在线看 | 婷婷久草| 国产精品久久久久久一区二区三区 | 久久精品最新 | 在线三级播放 | 成人久久电影 | 91自拍视频在线 | 少妇高潮流白浆在线观看 | 国产精品久久久久久久久蜜臀 | 美女久久久久久久 | 999久久| 日韩有码中文字幕在线 | 人人爽人人干 | 国产精品久久久久亚洲影视 | 久久精品1区 | 韩日三级av| 人人要人人澡人人爽人人dvd | 久久亚洲免费视频 | 精品久久影院 | 欧美激情视频一区二区三区免费 | 国产麻豆剧传媒免费观看 | 天天射天天舔天天干 | 8x8x在线观看视频 | 在线导航av | 一区二区影院 | 97视频在线 | 天天拍天天色 | 亚洲天天做 | 天天干天天操天天操 | 99久久婷婷国产精品综合 | 久久国产手机看片 | 日韩女同一区二区三区在线观看 | 天天摸夜夜添 | 欧美视频一区二 | 亚洲精品成人网 | 激情综合色播五月 | 日韩电影中文,亚洲精品乱码 | 日日爱网址 | 欧美精品视 | 国产成人av免费在线观看 | 国产精品无av码在线观看 | 丁香激情综合国产 | 国产精品久久一区二区三区, | 午夜美女网站 | 激情五月***国产精品 | 久操97| 色天天中文| 国产精品一区二区吃奶在线观看 | 男女日麻批 | 一区二区不卡在线观看 | 91精品在线看 | 国产在线中文字幕 | 天天操夜操视频 | 成人av电影在线 | 在线影视 一区 二区 三区 | 欧美性天天 | 国产精品美女视频网站 | 国产午夜精品理论片在线 | 999色视频| 超碰人人超 | 国产xx视频 | 丁香5月婷婷 | 欧美一性一交一乱 | 午夜 免费 | 欧美日韩国产一区二区在线观看 | 91免费在线看片 | 一区 二区电影免费在线观看 | 久久香蕉电影 | 91视频免费看 | 久久久久激情视频 | 综合亚洲视频 | 超碰个人在线 | 99久久精品免费看国产一区二区三区 | 精品国产成人av在线免 | 日日夜夜av | 国产99久久久欧美黑人 | 97电影手机版 | 国产精品网在线观看 | 亚洲综合视频在线播放 | 国产免费黄色 | 中中文字幕av | 欧美韩国日本在线观看 | 亚洲国产网址 | 在线观看成人av | 国产成在线观看免费视频 | 亚洲成人在线免费 | 毛片区 | 97视频精品 | 国产精品久久久久久久午夜片 | av日韩精品 | 久久久久国产一区二区三区四区 | 天天色天天操天天爽 | 天天操天天色天天射 | 精品一区二区av | 久草在线国产 | 成年人免费电影在线观看 | 三上悠亚一区二区在线观看 | 久久久人人人 | 日韩特级片 | 亚洲永久av | 久青草电影 | 中文字幕精品视频 | 欧美精品午夜 | 欧美最新大片在线看 | 国产精品99爱 | 国产精品高潮久久av | 十八岁以下禁止观看的1000个网站 | 国产这里只有精品 | 怡红院av久久久久久久 | 在线 影视 一区 | 亚洲综合在线五月天 | 久久观看免费视频 | 97碰碰视频| 亚洲传媒在线 | 久久99精品国产麻豆婷婷 | 在线观看黄色大片 | 国内视频在线 | 99热 精品在线 | 亚洲精品在线电影 | 久久情侣偷拍 | 国产精品九九久久久久久久 | www国产一区 | 青青河边草免费 | 日韩午夜视频在线观看 | 久久久久久久久久网站 | 久草在线在线精品观看 | 国产成人精品午夜在线播放 | 成人小视频在线播放 | 午夜精品一区二区三区在线 | 久草资源在线观看 | 干干日日 | a级黄色片视频 | 国产综合激情 | 91精品国产99久久久久久久 | 在线观看网站黄 | 在线观看国产中文字幕 | 日韩理论电影在线观看 | 免费在线观看毛片网站 | 久久人人爽人人爽 | 有码视频在线观看 | 91精品在线播放 | 国产在线精品一区二区三区 | 中文字幕乱码视频 | 日韩av视屏| 四季av综合网站 | 人人超碰在线 | 久久久一本精品99久久精品 | 天海翼一区二区三区免费 | 国产精品成人在线 | 九九热免费观看 | 亚洲特级毛片 | 国产一区视频导航 | 狠狠操在线 | 国产麻豆剧传媒免费观看 | 91久久国产露脸精品国产闺蜜 | 精品人妖videos欧美人妖 | 成人免费一区二区三区在线观看 | 日韩a在线播放 | 天天操天天操天天操天天操天天操天天操 | 97视频在线免费播放 | 日韩av一区在线观看 | 一级做a爱片性色毛片www | 99国产在线 | 中文字幕在线看视频国产 | 午夜精品一区二区三区免费 | 99久久国产免费免费 | 亚洲精选国产 | 国产免费黄视频在线观看 | 久久午夜电影 | 91网址在线观看 | 国产精品亚洲片在线播放 | 亚洲综合少妇 | 久久久久久久久网站 | 欧美一区二区三区在线播放 | 日韩中文在线播放 | 人人爱爱人人 | 久久午夜国产精品 | 成年人免费在线播放 | 成年人视频在线观看免费 | 国产精品一区二区久久久久 | 久久高清| 午夜91在线 | 丝袜美女视频网站 | 国产精品女人网站 | 国产日韩精品在线 | 国产99久久九九精品免费 | 狠狠干夜夜操 | 午夜精品久久久久99热app | 中文字幕国语官网在线视频 | av+在线播放在线播放 | 久久亚洲精品国产亚洲老地址 | 中文字幕在线国产 | 久久综合狠狠综合 | 91麻豆免费看 | 丁香在线观看完整电影视频 | 96久久 | 日韩精品中文字幕久久臀 | 一级特黄av| 亚洲无人区小视频 | 亚洲h在线播放在线观看h | 中文字幕婷婷 | 手机在线中文字幕 | 午夜性福利 | 亚洲黄色小说网址 | 免费高清在线观看电视网站 | 91麻豆操 | 69热国产视频 | 国产一级免费观看 | 色视频在线免费观看 | 久久你懂的 | 久久精品免费电影 | 一本色道久久精品 | 男女激情麻豆 | 日韩电影久久 | 国产精品综合久久久久 | 色视频一区 | 亚洲精品一区二区三区在线观看 | 色噜噜日韩精品一区二区三区视频 | 日本九九视频 | 日韩久久久久久 | a天堂免费 | 麻豆精品91 | 九九精品无码 | 又黄又爽的视频在线观看网站 | 91亚洲精品国偷拍 | 国产精品久久久久久久99 | 色婷婷伊人 | 久久国产精品久久国产精品 | 欧美精品久久99 | 国产精品免费在线视频 | 91亚洲精品在线观看 | 久久久久激情电影 | 99精品热视频| 麻豆免费在线视频 | 国产婷婷在线观看 | 天天色天天艹 | 亚洲视频,欧洲视频 | 天天干天天做 | 国产一级片久久 | av免费在线网 | 国产美女精品久久久 | 六月丁香激情综合 | 久久综合久久综合这里只有精品 | 国产成人精品一区二区三区在线观看 | 日本黄色大片免费看 | 黄色毛片视频免费观看中文 | 激情 婷婷| 中文字幕在线第一页 | 97在线观看视频国产 | 激情综合网天天干 | 97超碰人人爱| 国产91粉嫩白浆在线观看 | 日韩xxxxxxxxx | 国产精品国产三级国产 | 六月婷色 | 国产亚洲精品美女 | 日韩啪视频 | 在线观看视频99 | 免费十分钟 | 中文字幕在线影视资源 | 日韩视频1| 国产日韩在线视频 | 欧美精品久久久久久久免费 | 亚洲国产精品成人综合 | 99精品视频免费在线观看 | 亚洲精品午夜久久久久久久久久久 | 亚洲理论在线观看电影 | 中文字幕在线乱 | 日韩a级黄色 | 亚洲涩综合 | 久久综合精品国产一区二区三区 | 婷色在线 | 中文乱码视频在线观看 | 美女网站在线观看 | 久久麻豆精品 | 国产精品一区二区久久精品 | 欧美在线你懂的 | 成人黄色大片在线免费观看 | 免费看毛片在线 | 国产91免费观看 | 狠狠狠色丁香婷婷综合激情 | 欧美另类xxxx | 18网站在线观看 | 97国产大学生情侣酒店的特点 | 久久a国产 | 天天色天天综合 | 国产黄色电影 | 日韩电影在线观看一区 | 欧美日韩在线视频一区二区 | 超碰免费成人 | 欧美一级电影在线观看 | 国产在线看一区 | 婷婷成人亚洲综合国产xv88 | 欧美午夜a| 天天射天天操天天 | 在线视频a | 天天射天天爱天天干 | 成人污视频在线观看 | 五月天激情视频 | 欧美少妇xx| 久久国产精品一区二区三区 | 91精品导航 | 日本公妇在线观看 | 国产精品黑丝在线观看 | 国产高清一区二区 | 一二区av | 韩国av电影在线观看 | 日韩xxxx视频 | 精品夜夜嗨av一区二区三区 | 久久精品视频在线看 | 成人一级在线观看 | 精品久久久久久久久久久久久久久久久久 | 国产破处视频在线播放 | 国产粉嫩在线 | 成人av资源网站 | 日韩精品视频在线观看免费 | 婷婷激情站 | 天天搞天天干天天色 | 成人精品一区二区三区电影免费 | 久久久国产一区二区 | 久久精精品 | 中文字幕在线视频一区二区三区 | 成人一级视频在线观看 | 婷婷在线视频 | 亚洲精品在线看 | 色www精品视频在线观看 | 日韩在线高清免费视频 | 日韩精品播放 | 在线视频黄 | 久久久五月婷婷 | 欧美一区在线看 | 亚洲一区精品二人人爽久久 | 婷婷开心久久网 | 亚洲视频在线观看免费 | 91视频在线观看下载 | 国产成人精品日本亚洲999 | 久久国产精品久久精品国产演员表 | 亚洲 欧洲av| 日本视频久久久 | 99久久国产免费,99久久国产免费大片 | 国产精品久久久久aaaa | 中文字幕有码在线观看 | 亚洲丝袜中文 | 日日干天天操 | 成年人免费看片 | 337p日本大胆噜噜噜噜 | 久久高清免费视频 | 中文字幕色在线 | jizzjizzjizz亚洲| 免费日韩视| 91精品国产成人 | 奇米四色影狠狠爱7777 | 中文字幕丝袜一区二区 | 69国产精品成人在线播放 | 在线观看自拍 | 特级aaa毛片 | 99免在线观看免费视频高清 | 久久精品久久精品久久39 | 亚洲九九影院 | 久久久国产电影 | 天天色视频 | 日韩精品一卡 | 99精品免费视频 | 久久免费观看少妇a级毛片 久久久久成人免费 | 欧美一区二区三区免费看 | 在线影院av| 一区二区三区日韩在线观看 | 亚洲视频一区二区三区在线观看 | 国产91在线 | 美洲 | 精品国产乱码久久久久久三级人 | 国产亚洲精品久久19p | 9ⅰ精品久久久久久久久中文字幕 | 香蕉成人在线视频 | 久草在线 | 国产破处视频在线播放 | 精品久久一区 | 久久国产91 | 亚洲四虎在线 | av888av.com | 色播99 | 中国一级片在线 | 高清一区二区三区 | 91看片成人 | 国产手机av | 欧美韩国日本在线 | 日韩一级黄色大片 | 国产视频综合在线 | 一级片色播影院 | 日韩理论影院 | 免费在线观看成人小视频 | 九九国产视频 | 国产精品色婷婷视频 | 天天想夜夜操 | 黄色成人av网址 | 99免费在线观看 | 99精品免费在线观看 | 黄色av免费电影 | 精品国内自产拍在线观看视频 | 免费黄色特级片 | 1024手机基地在线观看 | 夜夜夜影院 | 午夜精品福利一区二区三区蜜桃 | 日韩二区三区在线 | 嫩模bbw搡bbbb搡bbbb | 国产一区电影在线观看 | 久草男人天堂 | 少妇av片 | 欧美aa在线观看 | 在线看日韩 | 色播五月激情综合网 | 欧美一级片播放 | 中文字幕一区二区三区在线视频 | 国产黄色免费看 | 色婷婷综合久久久 | 久久精品视频免费播放 | 久久成人人人人精品欧 | 日本精品久久久久久 | 欧美福利视频一区 | 婷婷色网 | 国产精品短视频 | 91激情小视频 | 9在线观看免费高清完整版在线观看明 | 国产高清视频免费最新在线 | 欧美激情视频一区二区三区免费 | 国产免费高清视频 | 日韩欧美一区二区三区在线 | 中文字幕一区二区三区四区 | 久久免费精品国产 | 91亚洲精品久久久久图片蜜桃 | 在线久热 | 欧洲精品一区二区 | 日韩理论在线播放 | 狠狠地操 | 国产福利在线不卡 | 精品成人免费 | 人人看黄色 | 亚洲精品国产精品乱码不99热 | 美女视频黄是免费的 | 久久久久日本精品一区二区三区 | 日本免费久久高清视频 | 手机成人在线电影 | 精品国产精品久久一区免费式 | 成人av一二三区 | 亚洲免费一级 | 狠狠干狠狠久久 | 91看片在线播放 | 91九色在线视频观看 | 欧美日韩免费视频 | 日日日日 | av在线永久免费观看 | 欧美成亚洲 | 久久精品—区二区三区 | 午夜 在线 | 久久久国产在线视频 | 在线视频一区观看 | 精品国产一区二区在线 | 九九久久精品视频 | 亚洲午夜久久久影院 | 黄色电影网站在线观看 | 免费看片亚洲 | 国产剧在线观看片 | 亚洲高清在线观看视频 | 久操操 | 亚洲精品欧美成人 | 日韩天天干 | 国产精品久久久久久久毛片 | 成人午夜电影免费在线观看 | 黄www在线观看 | 久久99精品国产麻豆宅宅 | 亚洲视屏一区 | 久久香蕉国产 | 韩国一区二区在线观看 | 丝袜美腿亚洲 | 日韩欧美在线一区 | www国产亚洲精品久久网站 | 在线免费观看一区二区三区 | av在线永久免费观看 | 五月婷婷丁香色 | 国产精品资源在线观看 | 国产大片黄色 | 欧洲激情综合 | 91色蜜桃| 久久黄色小说视频 | 国产拍揄自揄精品视频麻豆 | 激情图片区 | 97热视频 | 久久亚洲国产精品 | 久久不卡国产精品一区二区 | 精品视频久久 | 九色视频网站 | 69人人| 久久99精品久久只有精品 | 国产中文字幕视频在线观看 | 国产精品永久免费观看 | 天天干,天天射,天天操,天天摸 | 国产色啪| 一区中文字幕在线观看 | 亚洲资源网 | 国产成人av电影在线观看 | 黄色毛片在线 | 久久精品1区2区 | 日韩一区精品 | 久草精品视频在线播放 | 天天插综合网 | 欧美成人亚洲成人 | 欧美日韩首页 | 国产视频1 | 91黄色免费看 | 久久蜜臀一区二区三区av | 亚洲自拍av在线 | 91中文字幕在线观看 | 国产精品一区二区av日韩在线 | 六月激情 | 色综合国产 | 欧美另类高潮 | 在线播放你懂 | 国产精品免费一区二区 | 黄色大片国产 | 91在线免费观看网站 | 国产在线观看免费观看 | 国产精品不卡视频 | 三级黄色理论片 | 久久久99精品免费观看乱色 | 日韩电影一区二区在线观看 | 91久久黄色| 在线国产日韩 | 伊人中文字幕在线 | 国色综合| 亚洲成人频道 | 婷婷六月综合网 | 欧美伦理一区二区三区 | 亚洲成人av在线电影 | 久久久久久久久久久久av | 久久综合色婷婷 | 九月婷婷综合网 | 国产视频一二区 | av电影一区二区三区 | 免费观看一级特黄欧美大片 | 色www免费视频 | 欧美一级片在线播放 | 黄色小网站在线 | 国产偷v国产偷∨精品视频 在线草 | 欧美精品亚洲精品日韩精品 | 激情影院在线 | 黄色特级一级片 | 亚洲综合欧美日韩狠狠色 | 亚洲成人精品影院 | 日韩激情精品 | 999久久久精品视频 日韩高清www | 国产一区二区三区视频在线 | 97精品免费视频 | 国产精品18videosex性欧美 | 亚洲欧美一区二区三区孕妇写真 | 国产午夜精品免费一区二区三区视频 | 久久色视频 | 久久亚洲国产精品 | 麻豆视频免费播放 | 欧美一区二区三区四区夜夜大片 | 2000xxx影视| 亚洲精品视频中文字幕 | 青青河边草观看完整版高清 | 天天综合天天综合 | 九色福利视频 | 久久美女免费视频 | 日日操天天爽 | 久在线观看 | 国产99re| 国产精品男女 | 操操色 | 国产精品一区欧美 | 色偷偷人人澡久久超碰69 | 97操碰| 国产精品永久免费观看 | 天天干天天干天天操 | 国产在线一线 | 在线观看视频在线观看 | 九月婷婷人人澡人人添人人爽 | 久久久久国产成人精品亚洲午夜 | 亚洲最新av在线网址 | 国产精品岛国久久久久久久久红粉 | 久久理论视频 | 狠狠色丁香婷婷综合橹88 | 日韩在线视频观看免费 | 在线国产黄色 | 国产精品成人a免费观看 | 一级免费看 | 中文字幕国产一区 | 一区二区三区免费播放 | 一区二区视频在线看 | 午夜黄色大片 | 国产视频精品免费播放 | 日韩激情视频在线 | 毛片网在线| 免费三级a | 欧美成天堂网地址 | 国产精品久久久久久久av大片 | 综合天天| 久久超级碰视频 | 91在线视频免费 | 欧美激情亚洲综合 | 国产成人av片 | 国产一区免费在线观看 | 久久高清国产 | 久久精品综合网 | 97人人艹| 日本一区二区三区视频在线播放 | 久久精品99国产精品日本 | 欧美大片第1页 | 三级黄色网络 | 国产精品毛片一区二区在线 | 国产xx在线 | 亚洲,播放 | av网站在线免费观看 | 色视频网站在线 | 免费在线观看一级片 | 91亚洲精品乱码久久久久久蜜桃 | 久久国产精品一国产精品 | 日韩欧美视频免费观看 | 国产高清在线视频 | 亚洲欧洲精品视频 | 日韩国产精品毛片 | 亚洲欧洲国产日韩精品 | 中文字幕在线播放一区二区 | 天天草天天 | 国产精品久久久久久久久久久久午夜 | 国产精品一二三 | 一级成人网 | 毛片网站免费在线观看 | 久久爱www. | 51久久夜色精品国产麻豆 | 婷婷综合av| 丁香亚洲 | 国产特级毛片 | 天堂av网址 | 国产一级一级国产 | 亚洲视频 中文字幕 | 久久久久久久av | 国产91精品久久久久久 | 中文字幕日本特黄aa毛片 | 在线观看av片 | 中文字幕精品www乱入免费视频 | 91久久久久久国产精品 | 久久久久亚洲精品成人网小说 | 91九色在线观看 | 国产精品专区在线观看 | 日韩字幕 | 国产精品久久久久aaaa九色 | 成人免费色 | 天天玩夜夜操 | 狠狠色丁香久久婷婷综合五月 | 福利二区视频 | www黄色| 国产亚洲精品久久 | 18做爰免费视频网站 | 99综合电影在线视频 | 午夜视频一区二区三区 | 国产精品亚洲视频 | 成人网中文字幕 | 精品亚洲成a人在线观看 | 国内精品久久久久久久97牛牛 | 国产一区二区电影在线观看 | 日本精品小视频 | 91麻豆精品国产91久久久无需广告 | 国产成人在线观看免费 | 久久精品国产v日韩v亚洲 | 亚洲精品欧美成人 | www日韩在线 | 一区二区三区在线影院 | 日韩久久精品一区二区三区 | 久久99精品久久久久久秒播蜜臀 | 日韩小视频网站 | 久久久久国产精品一区二区 | 99精品亚洲 | 视频二区在线 | 亚洲美女在线国产 | 99热免费在线 | 国产成人免费在线 | 天天干天天干天天干 | 亚洲永久精品一区 | 99欧美 | 久久久国产精华液 | 天堂av中文字幕 | 亚洲精品小区久久久久久 | 18性欧美xxxⅹ性满足 | 亚洲精品一区二区久 | 在线视频a| 麻豆播放| 日本精品二区 | 在线看成人 | 激情丁香综合 | 日韩在线一二三区 | 久久艹免费 | 国产精品区在线观看 | 成人一区二区三区在线 | 黄色片免费电影 | 色视频在线看 | 日本中文字幕在线看 | 欧美日韩大片在线观看 | 国产亚洲综合精品 | 中文在线资源 | 91精品久久久久久综合乱菊 | 免费在线色 | av在线永久免费观看 | 深夜免费福利视频 | 日韩电影中文字幕在线观看 | 91精品久久香蕉国产线看观看 | 黄色h在线观看 | 叶爱av在线| 深夜福利视频一区二区 | 国产精品久久久久久麻豆一区 | 国产精品美女久久久久久久网站 | 91精品国产自产在线观看 | 久久精品一区二区三 | 亚洲综合成人专区片 | 国产精品18久久久久久不卡孕妇 | 久久精品一二区 | 国产精品国产三级国产不产一地 | 狠狠干美女 | 精品国产一区二区三区四 | av在线免费在线观看 | 最近字幕在线观看第一季 | 天天爽天天爽 | 婷婷激情综合五月天 | 黄色免费观看网址 | 波多野结衣视频一区二区 | 91精品综合在线观看 | 狠狠色婷婷丁香六月 | 久久综合色8888 | 国产区免费在线 | h视频在线看 | 特级a老妇做爰全过程 | 夜夜操天天 | 在线观看视频亚洲 | 天天干.com | 免费视频久久久久久久 | 久久久久久激情 | 国产精品免费一区二区三区 | 激情网色| 99久久久久久久久久 | 亚洲人在线7777777精品 | 日韩午夜av电影 | 亚洲人成网站精品片在线观看 | 久久激情片 | 99视频在线免费看 | 成人免费一级片 | 三上悠亚一区二区在线观看 | 在线一二区 | 蜜臀av性久久久久蜜臀aⅴ涩爱 | 欧美日韩性 | 少妇bbw搡bbbb搡bbbb | 色婷婷www | 欧美日韩不卡在线观看 | 99精品国产99久久久久久97 | 嫩嫩影院理论片 | 最新色视频| 精品国内自产拍在线观看视频 | 日韩中文字幕视频在线 | 欧美特一级片 | 国产一区二区不卡在线 | 激情五月婷婷激情 | 欧美激情综合五月色丁香小说 | 日日摸日日碰 | 日韩欧美专区 | 91九色视频国产 | 国产日韩一区在线 | 国产视频 亚洲视频 | 中文字幕一区二区三区视频 | 亚洲va韩国va欧美va精四季 | 在线观看一区视频 | 激情综合亚洲 | 伊人黄 | 国产五十路毛片 | 四虎8848免费高清在线观看 | 日韩激情视频在线观看 | 国产欧美综合在线观看 | 久99久精品 | 久色 网| 久久久久久影视 | 国产精品免费视频观看 | 久久精品国产成人精品 | 婷婷电影在线观看 | 日韩有码欧美 | 国产精品第10页 | 精品国产乱码久久久久久浪潮 | 久久一区精品 | 精品视频不卡 | 日韩欧美91| 五月天激情在线 | 亚洲精品视频国产 | 国产1区2区| 日韩欧美v | 欧美一级艳片视频免费观看 | japanesexxxxfreehd乱熟| 国产视频一区精品 | 久久黄色精品视频 | 成年人看片| 麻豆国产精品视频 | 在线观看v片 | 成人国产精品电影 | 日韩欧美成人网 | 色诱亚洲精品久久久久久 | 国产毛片久久 | 色婷婷电影 | 人人干人人模 | 一级免费看 | 免费国产在线精品 | 性色av一区二区三区在线观看 | 伊人久操| 国产精品国产毛片 | 午夜久久福利影院 | 亚洲精品国产成人 | 在线看成人av | 97国产在线观看 | 日日爽视频 | 蜜臀aⅴ国产精品久久久国产 | 免费开视频 | 久草久草视频 | 国产成人三级三级三级97 | 中文字幕精品一区久久久久 | 伊人久久在线观看 | 免费av网址大全 | 天天干天天拍天天操天天拍 | 免费情缘| 欧美日韩精品在线一区二区 | 一区二区三区四区在线免费观看 | 国产免费又爽又刺激在线观看 | 日韩久久精品一区二区三区 | 久久这里有| 草久在线 | 久草视频在线免费看 | 毛片网站免费在线观看 | 网站在线观看你们懂的 | 久久99国产一区二区三区 | 91丨九色丨蝌蚪丨老版 | 亚洲一区在线看 | 国产高清不卡av | 手机看片午夜 | 国产精品麻豆一区二区三区 | 精品不卡av| 日韩一区二区三免费高清在线观看 | 手机在线观看国产精品 | 深爱五月网 | 国产成人精品一区二区三区福利 | 欧美日韩国产一区二区三区 | 正在播放国产一区 | 国产精品理论片在线播放 | 久久久国产精品一区二区中文 | 亚洲精品高清视频在线观看 | 99视频这里有精品 | 亚洲日本欧美在线 | 免费激情网 | 国产永久免费高清在线观看视频 | 99精品国产在热久久下载 | 国产精品精 | 99久久久国产精品免费99 | 中文字幕第一 | 亚洲妇女av | 91亚洲精品国产 | 日韩免费视频播放 | 天天操网站 | 欧美一级裸体视频 | 国产精品精品国产 | 久久黄色免费视频 | 狠狠色狠狠综合久久 | 亚洲婷婷在线视频 | 夜夜骑天天操 | 制服丝袜在线 | 午夜精品中文字幕 | 成人夜晚看av | 91精品视频免费 | 久久久久久久久久伊人 | 久久高视频| 午夜av色| 国产一区网 | 国产成人精品三级 | 一级黄色免费 | 日韩免费在线一区 | av在线小说 | 亚洲九九九在线观看 | 97超碰人人澡人人爱 | 久久毛片网站 | 成人一级免费电影 | 久久久久久久久免费 | 欧美作爱视频 | 天天操天天干天天爽 | 欧美精品亚洲精品 | 国产美女在线观看 | 99国产精品免费网站 | 水蜜桃亚洲一二三四在线 | 在线视频观看亚洲 | 国产亚洲视频在线观看 | 丝袜少妇在线 | 亚洲精品日韩av | 天天爱天天操天天爽 | 亚洲国内精品视频 | 午夜18视频在线观看 | 欧美色精品天天在线观看视频 | 日韩高清在线一区二区三区 | 九九视频在线 | 久久精品网站视频 | 99久久精品国产免费看不卡 | 人人狠狠综合久久亚洲婷 | 成人在线免费观看视视频 | 麻豆视频在线观看免费 | 日韩免费视频线观看 | 国产精品96久久久久久吹潮 | 欧美一级电影在线观看 | 97视频在线观看网址 | 久久久www成人免费精品 | 久久婷婷国产 | 成人毛片在线观看 | 中文在线8新资源库 | av片一区二区 | 久久久精品二区 | 久久伦理视频 | 成年人在线电影 | 成人小视频免费在线观看 | 在线欧美最极品的av | 高清不卡一区二区三区 | 丁香久久激情 | av网站播放| 日本精品午夜 | 开心激情网五月天 | 亚洲黄色一级大片 | 亚洲va韩国va欧美va精四季 | 九九九在线 | 国产欧美综合视频 | 天天色官网 | 日韩h在线观看 | 国产成人精品久久二区二区 | 不卡av在线 | 中文字幕二区在线观看 | 久草剧场 | 美女免费视频一区二区 | 欧美日韩精品综合 | 日韩欧美一区二区三区在线观看 | 天天五月天色 | 国产黄色精品在线观看 | 日韩久久精品一区二区 | 免费一级片在线观看 | 精品久久九九 | 中文字幕在线电影 | 一级片免费视频 | 亚洲第一久久久 | 久久99国产精品视频 | 午夜精品久久久久久久99婷婷 | 夜夜躁天天躁很躁波 | 激情在线免费视频 | 成年人免费在线观看网站 | 久久国产精品久久精品国产演员表 | 丁香花五月 | 成人一级免费电影 | 国模精品在线 | 黄色片免费在线 | 九色91在线视频 | 国产福利电影网址 | 欧美aa级 | 久久精品欧美日韩精品 | 五月天伊人网 | 亚洲一区精品二人人爽久久 | 免费观看mv大片高清 | 国产成人精品av在线 | 国产小视频免费在线网址 | 午夜一级免费电影 | 日韩欧美精品在线观看视频 | 国产亚洲精品日韩在线tv黄 | 最新的av网站| 久久久久成| 国产精品成人久久 | 国产91aaa| 91精品在线免费视频 | 欧美成人999 | 精品亚洲男同gayvideo网站 | 黄色av网站在线观看 | 九九视频免费在线观看 | 在线草| 日韩精品首页 | 国产色小视频 |