逻辑回归与最大熵模型

本文为《统计学习方法》第6章笔记。

概论

逻辑回归与最大熵模型都属于对数线性模型,逻辑回归求解似然函数的极大值,得到参数w,最大熵模型先转对偶问题,求得条件概率模型,也是通过极大值求解得到w。涉及到最优化算法的部分都比较晦涩,由于本人理解得也不深刻,这里就不详细展开了。

逻辑回归模型

逻辑斯蒂分布

设X为连续随机变量,服从逻辑斯蒂分布的分布函数和密度函数分别为:
F(x)=P(X \leq x)=\frac{1}{1+e^{-(x-u)/\gamma}}\\ f(x)=F'(x)=\frac{e^{-(x-u)/\gamma}}{\gamma(1+e^{-(x-u)/\gamma})^2}

二项逻辑斯蒂回归模型

条件概率分布为:
P(Y=1|x)=\frac{exp(w \odot x+b)}{1+exp(w \odot x+b)}\\ P(Y=0|x)=\frac{1}{1+exp(w \odot x+b)}
将b扩充到w向量内,1扩充到x向量内,则方程简化为
P(Y=1|x)=\frac{exp(w \odot x)}{1+exp(w \odot x)}\\ P(Y=0|x)=\frac{1}{1+exp(w \odot x)}
几率:事情发生与不发生的概率之比。
对数几率logit(p)=log\frac{p}{1-p}
对逻辑回归而言,对数几率为w\odot x

模型参数估计

P(Y=1|x)=\pi(x)\\ P(Y=0|x)=1-\pi(x)
似然函数为:
\prod_{i=1}^N \pi(x_i)^{y_i}[1-\pi(x_i)]^{1-y_i}
对数似然函数为:
L(w)=\sum_{i=1}^N [y_ilog\pi(x_i)+(1-y_i)log(1-\pi(x_i))]\\ =\sum_{i=1}^N [y_i log\frac{\pi(x_i)}{1-\pi(x_i)}+log(1-\pi(x_i))]\\ =\sum_{i=1}^N [y_i (w \odot x_i)-log(1+exp(w \odot x_i))]
对L(w)求极大值,得到w的估计值。一般使用梯度下降或拟牛顿法(代码中常见的BFGS算法)

多项逻辑斯蒂回归

P(Y=k|x)=\frac{exp(w_k \odot x)}{1+\sum_{k=1}^{K-1}exp(w_k \odot x)},k=1,2,...K-1\\ P(Y=K|x)=\frac{1}{1+\sum_{k=1}^{K-1}exp(w_k \odot x)}

最大熵模型

假设离散随机变量X的概率分布为P(X),则其熵为:
H(P)=-\sum_x P(x)logP(x)
熵满足不等式:
0 \leq H(p) \leq log|X|
|X|为X的取值个数,当X服从均匀分布时,熵最大。
熵最大原理:最好的模型是熵最大的模型。在没有更多信息时,不确定的部分都是等概率的。
给定训练数据集,联合经验分布和边缘经验分布为:
\widetilde P(X=x,Y=y)=\frac {v(X=x,Y=y)}{N}\\ \widetilde P(X=x)=\frac {v(X=x)}{N}\\ v(X=x,Y=y)表示训练数据中(x,y)出现的频次,v(X=x)表示x出现的频次
用特征函数f描述x和y之间的某一个事实:
f(x,y)= \left\{\begin{array}\\ 1,x与y满足某一个事实\\ 0,其他 \end{array}\right.
f关于联合经验分布的期望值为:
E_{\widetilde P}(f)=\sum_{x,y}\widetilde P(x,y)f(x,y)
f关于边缘经验分布及模型P(Y|X)的期望为:
E_P(f)=\sum_{x,y}\widetilde P(x)P(y|x)f(x,y)
如果模型能够获取数据中的信息,那么2个期望值应该相等
E_{\widetilde P}(f)=E_P(f)
定义在P(Y|X)上的条件熵为:
H(P)=-\sum_{x,y}\widetilde P(x)P(y|x)logP(y|x)
满足所有约束的模型集合C为:
C \equiv {P \in \rho|E_{\widetilde P}(f_i)=E_P(f_i),i=1,2,...n}
模型集合C中条件熵最大的模型成为最大熵模型,公式中对数为自然对数。
在约束条件下的求解,需要用到拉格朗日函数,然后转对偶问题,这里略去推导过程(可自行查询资料),对偶问题为:
max_w min_{P \in C} L(P,w)
先求得min_{P \in C} L(P,w)的解为
P_w(y|x)=\frac{1}{Z_w(x)}exp(\sum_{i=1}^n w_if_i(x,y))\\ Z_w(x)=\sum_y exp(\sum_{i=1}^n w_if_i(x,y)) ,称为规范化因子
再求极大化问题,得到参数w的解。这需要用到最优化算法,常见的有改进的迭代尺度法(IIS)、梯度下降法、牛顿法、拟牛顿法(典型如BFGS)
这里稍微做点介绍:
IIS思路:假设w+\delta能使得模型的似然函数增大,则更新w\rightarrow w+\delta
IIS最后求解的方程为:
\sum_ {x,y}\widetilde P(x)P_w(y|x) f_i(x,y)exp(\delta_i f^{\#}(x,y))=E_{\widetilde P}(f_i) -----------(1)\\ f^{\#}(x,y)=\sum_i f_i(x,y),表示特征在(x,y)出现的次数
如果f^{\#}(x,y)是常数M(对任何x,y),则:
\delta_i=\frac{1}{M}log\frac{E_{\widetilde P}(f_i)}{E_{P}(f_i)}
如果不是常数,可以通过牛顿法求解:
g(\delta_i)=0表示以上(1)式,则:
\delta_i^{k+1}=\delta_i^k-\frac{g(\delta_i^k)}{g'(\delta_i^k)}
拟牛顿法思路:牛顿法需要求解黑塞矩阵H的逆,较为麻烦,拟牛顿法考虑用一个n阶矩阵逼近H(BFGS算法)或H的逆(DFP算法)
BFGS算法更新逼近矩阵:
B_{k+1}=B_k+\frac{y_ky_k^T}{y_k^T\delta_k}-\frac{B_k\delta_k\delta_k^TB_k}{\delta_k^TB_k\delta_k}\\ y_k=g_{k+1}-g_k , \delta_k=w^{k+1}-w^k \\ 其中B_k用于逼近H

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

推荐阅读更多精彩内容