上哪个网站做forma和co写软文能赚钱吗
最大熵模型中的对数似然函数的解释
最近在学习最大熵模型,看到极大似然估计这部分,没有看明白条件概率分布p(y∣x)p(y|x)p(y∣x)的对数似然函数。上网查了很多资料都没有一个合理的解释。基本直接给出对数似然函数的一般形式:
Lp‾=∏xp(x)p‾(x).L_{\overline{p}}=\prod_{x} p(x)^{\overline{p}(x)}.Lp=x∏p(x)p(x).
其实并没有解决问题。为了方便以后其他人的学习和理解,我结合自己的理解给出完整的解释。
其实第一眼之所以不理解,因为这是最大似然函数的另外一种形式。一般书上描述的最大似然函数的一般形式是各个样本集XXX中各个样本的联合概率:
L(x1,x2,...,xn;θ)=∏i=1np(xi;θ).L(x_1,x_2,...,x_n;\theta)=\prod_{i=1}^{n} p(x_i;\theta).L(x1,x2,...,xn;θ)=i=1∏np(xi;θ).
其实这个公式和上式是等价的。x1,x2,...,xnx_1,x_2,...,x_nx1,x2,...,xn是样本具体观测值。随机变量XXX是离散的,所以它的取值范围是一个集合,假设样本集的大小为nnn,XXX的取值有kkk个,分别是v1,v2,...,vkv_1,v_2,...,v_kv1,v2,...,vk。用C(X=viC(X=v_iC(X=vi)表示在观测值中样本viv_ivi出现的频数。所以L(x1,x2,...,xn;θ)L(x_1,x_2,...,x_n;\theta)L(x1,x2,...,xn;θ)可以表示为:
L(x1,x2,...,xn;θ)=∏i=1kp(vi;θ)C(X=vi).L(x_1,x_2,...,x_n;\theta)=\prod_{i=1}^{k} p(v_i;\theta)^{C(X=v_i)}.L(x1,x2,...,xn;θ)=i=1∏kp(vi;θ)C(X=vi).
对等式两边同时开nnn次方,可得
L(x1,x2,...,xn;θ)1n=∏i=1kp(vi;θ)C(X=vi)n.L(x_1,x_2,...,x_n;\theta)^{\frac{1}{n}}=\prod_{i=1}^{k} p(v_i;\theta)^{\frac{C(X=v_i)}{n}}.L(x1,x2,...,xn;θ)n1=i=1∏kp(vi;θ)nC(X=vi).
因为经验概率p‾(x)=C(X=vi)n\overline{p}(x)=\frac{C(X=v_i)}{n}p(x)=nC(X=vi),所以简写得到:
L(x1,x2,...,xn;θ)1n=∏xp(x;θ)p‾(x).L(x_1,x_2,...,x_n;\theta)^{\frac{1}{n}}=\prod_{x} p(x;\theta)^{\overline{p}(x)}.L(x1,x2,...,xn;θ)n1=x∏p(x;θ)p(x).
很明显对L(x1,x2,...,xn;θ)L(x_1,x_2,...,x_n;\theta)L(x1,x2,...,xn;θ)求最大值和对L(x1,x2,...,xn;θ)1nL(x_1,x_2,...,x_n;\theta)^{\frac{1}{n}}L(x1,x2,...,xn;θ)n1求最大值的优化的结果是一样的。整理上式所以最终的最大似然函数可以表示为:
L(x;θ)=∏xp(x:θ)p‾(x).L(x;\theta)=\prod_{x} p(x:\theta)^{\overline{p}(x)}.L(x;θ)=x∏p(x:θ)p(x).
忽略θ\thetaθ,更一般的公式就是本文的第一个公式。结合公式一,参考v_JULY_v博客中的最大熵模型中的数学推导(http://m.blog.csdn.net/v_july_v/article/details/40508465),可得到联合概率密度的似然函数,即最大熵中的对数似然函数:
Lp‾=log∏x,yp(x,y)p‾(x,y)=∑x,yp‾(x,y)logp(x,y)=∑x,yp‾(x,y)log[p‾(x)p(y∣x)]=∑x,yp‾(x,y)logp(y∣x)+∑x,yp‾(x,y)logp‾(x)\begin{aligned} L_{\overline{p}} &=\log\prod_{x,y} p(x,y)^{\overline{p}(x,y)} \\ &=\sum_{x,y}{\overline{p}(x,y)}\log p(x,y)\\ &=\sum_{x,y}{\overline{p}(x,y)}\log [{\overline{p}(x)}p(y|x)] \\ &=\sum_{x,y}{\overline{p}(x,y)}\log p(y|x)+\sum_{x,y}{\overline{p}(x,y)}\log {\overline{p}(x)} \end{aligned} Lp=logx,y∏p(x,y)p(x,y)=x,y∑p(x,y)logp(x,y)=x,y∑p(x,y)log[p(x)p(y∣x)]=x,y∑p(x,y)logp(y∣x)+x,y∑p(x,y)logp(x)
上述公式第二项是一个常数项(都是样本的经验概率),一旦样本集确定,就是个常数,可以忽略。所以最终的对数似然函数为:
Lp‾=∑x,yp‾(x,y)logp(y∣x).L_{\overline{p}}=\sum_{x,y}{\overline{p}(x,y)}\log p(y|x).Lp=x,y∑p(x,y)logp(y∣x).
上式就是最大熵模型中用到的对数似然函数。