LDA漫游系列(二)-数学基础

要想了解LDA背后的算法原理,我们首先要有一定的数学基础,这一节,我们一起来学习LDA算法中涉及的一些基础的数学知识。

1、gamma函数

所谓的gamma函数其实就是阶乘的函数形式,中学时候我们都学过n的阶乘,比如我们知道n!=123....n,如果我问你3的阶乘是多少,你能够立即回答出12*3=6,但是咱们学过的阶乘知识仅限于整数范围,如果我问你0.5的阶乘是多少,你可能会哑口无言,不过没关系,学过gamma 函数之后你就能够回答了,欧拉经过不懈努力,终于发现阶乘的更一般的函数形式gamma函数,gamma函数如下所示:

不难发现:



接下来,我们共同推导一下0.5的阶乘:



在上面的推导过程中,大家可能对雅可比行列式不是十分了解,我也一样,所以经过百度,终于明白了其中的道理,雅可比行列式在求解积分的过程中往往可以发挥重要的作用:

所以,在极坐标变换中:



为什么gamma函数可以表示阶乘呢,这里我们只需证明下面的式子即可:

证明过程也很简单:

上面提到了分部积分法,说实话是真忘了,如果大家对分部积分法有点遗忘的话,希望下面的图片可以帮你回忆起来,如果你没有忘,下面的部分可以直接跳过:

所以,其实gamma函数表示的是x-1的阶乘。即:

2、二项分布

二项分布相信大家都接触过,二项分布即重复n次独立的伯努利试验,在每次实验中只有两种可能的结果,成功或者失败,每次成功的概率为p,失败的概率就为1-p,假设在n次实验中,成功了k次,则二项分布的概率密度函数是:


3、beta分布

在概率论中,beta分布是指一组定义在区间(0,1)的连续概率分布,有两个参数α和 β,且 α ,β > 0
Beta分布的概率密度函数是:



随机变量服从参数为α和 β的beta分布通常写作,X~Beta(α,β),式中的分布的函数称为B函数
可以看到beta分布中用到了gamma分布,所以我们首先来验证一个公式,这个公式也就是B函数和Gamma函数的关系:



这里给出一种证明方法

证明过程太复杂了,没有点数学功底是搞不出来的,不过我们就先记住结论吧。

接下来我们来看一下Beta分布的期望:


讲了这么多Beta分布的公式,那么我们该如何通俗的理解Beta分布呢?用一句话来说,beta分布可以看作一个概率的概率分布,当你不知道一个东西的具体概率是多少时,它可以给出了所有概率出现的可能性大小。
举一个简单的例子,熟悉棒球运动的都知道有一个指标就是棒球击球率(batting average),就是用一个运动员击中的球数除以击球的总数,我们一般认为0.266是正常水平的击球率,而如果击球率高达0.3就被认为是非常优秀的。现在有一个棒球运动员,我们希望能够预测他在这一赛季中的棒球击球率是多少。你可能就会直接计算棒球击球率,用击中的数除以击球数,但是如果这个棒球运动员只打了一次,而且还命中了,那么他就击球率就是100%了,这显然是不合理的,因为根据棒球的历史信息,我们知道这个击球率应该是0.215到0.36之间才对啊。对于这个问题,我们可以用一个二项分布表示(一系列成功或失败),一个最好的方法来表示这些经验(在统计中称为先验信息)就是用beta分布,这表示在我们没有看到这个运动员打球之前,我们就有了一个大概的范围。beta分布的定义域是(0,1)这就跟概率的范围是一样的。接下来我们将这些先验信息转换为beta分布的参数,我们知道一个击球率应该是平均0.27左右,而他的范围是0.21到0.35,那么根据这个信息,我们可以取α=81,β=219



之所以取这两个参数是因为:

beta分布的均值是

从图中可以看到这个分布主要落在了(0.2,0.35)间,这是从经验中得出的合理的范围。

在这个例子里,我们的x轴就表示各个击球率的取值,x对应的y值就是这个击球率所对应的概率。也就是说beta分布可以看作一个概率的概率分布。

那么有了先验信息后,现在我们考虑一个运动员只打一次球,那么他现在的数据就是”1中;1击”。这时候我们就可以更新我们的分布了,让这个曲线做一些移动去适应我们的新信息。结果很简单:
其中α0和β0是一开始的参数,在这里是81和219。所以在这一例子里,α增加了1(击中了一次)。β没有增加(没有漏球)。这就是我们的新的beta分布Beta(81+1,219),我们跟原来的比较一下:

可以看到这个分布其实没多大变化,这是因为只打了1次球并不能说明什么问题。但是如果我们得到了更多的数据,假设一共打了300次,其中击中了100次,200次没击中,那么这一新分布就是:



注意到这个曲线变得更加尖,并且平移到了一个右边的位置,表示比平均水平要高。
一个有趣的事情是,根据这个新的beta分布,我们可以得出他的数学期望为:
,这一结果要比直接的估计要小
你可能已经意识到,我们事实上就是在这个运动员在击球之前可以理解为他已经成功了81次,失败了219次这样一个先验信息。
因此,对于一个我们不知道概率是什么,而又有一些合理的猜测时,beta分布能很好的作为一个表示概率的概率分布。

4、多项分布

多项分布是二项分布的扩展,在n次独立试验中每次只输出k种结果中的一个,且每种结果都有一个确定的概率p。多项分布给出了在多种输出状态的情况下,关于成功次数的各种组合的概率。
举一个掷骰子的例子,掷n次骰子,每次可能有六种输出,1点出现概率为p1,2点出现概率为p2,依次类推,在n次试验中,骰子1点出现次数为x1次,骰子2点出现次数为x2次,依次类推,这个组合出现的概率为:



也可以用gamma分布来表示:


5、狄利克雷分布(dirichlet distribution)

狄利克雷分布是beta分布在多项情况下的推广,也是多项分布的共轭先验分布(共轭先验分布我们在接下来会讲到),dirichlet分布的概率密度函数如下:


还有一个重要的公式,实在是太难推导了,大家记住结论就好:

跟前面beta分布类似,我们可以得到参数p的期望:

参考文献:https://www.zhihu.com/question/30269898/answer/123261564

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,456评论 5 477
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,370评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,337评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,583评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,596评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,572评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,936评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,595评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,850评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,601评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,685评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,371评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,951评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,934评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,167评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 43,636评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,411评论 2 342

推荐阅读更多精彩内容