机器学习的许多公式推导都涉及了数理统计的内容,特别是参数估计对理解机器学习很重要。这里三篇文章就对三种参数估计方法进行简单介绍。
对一些数理统计的基本概念的介绍,可参考之前的文章“数理统计学的基本概念”。
参数的点估计问题
设有一个统计总体,以f(x;θ1,⋯,θk)f(x;θ1,⋯,θk)记其概率密度函数(若总体分布为连续型的)或其概率函数(若其总体分布为离散型的)。避免重复交代这两种情况,我们约定称f(x;θ1,⋯,θk)f(x;θ1,⋯,θk)为“总体分布”,其具体含义视其为连续型或离散型而定。这个分布包含k个位置参数θ1,⋯,θkθ1,⋯,θk。例如,对正态总体分布N(μ,δ2)N(μ,δ2),有θ1=μθ1=μ,θ2=δ2θ2=δ2,而
若总体有二项分布B(n,p)B(n,p),则θ1=pθ1=p,而
当k=1,即只有一个参数时,就用θθ代替θ1θ1
参数估计问题的一般提法是:设有了从总体中抽出样本X1,⋯,XnX1,⋯,Xn(独立同分布),要依据这些样本去对参数θ1,⋯,θkθ1,⋯,θk的未知值做出估计。当然我们也可以只要求估计θ1,⋯,θkθ1,⋯,θk中的一部分,或估计他们的某个已知函数g(θ1,⋯,θk)g(θ1,⋯,θk)。例如,为要估计θ1θ1,我们需要构造出适当的统计量θ1^=θ1^(X1,⋯,Xn)θ1^=θ1^(X1,⋯,Xn)。每当有了样本X1,⋯,XnX1,⋯,Xn,就代入函数θ1^(X1,⋯,Xn)θ1^(X1,⋯,Xn)中计算出一个值,用来作为θ1θ1的估计值。为着这样的特定目的而构造的统计量θ1^θ1^叫做θ1θ1的估计量。由于未知参数θ1θ1是数轴上的一个点,用θ1^θ1^去估计θ1θ1,等于用一个点估计另一个点,所以这样的估计叫做点估计,以别与区间估计。
矩估计法
矩估计法的思想比较简单:设总体分布为f(x;θ1,⋯,θk)f(x;θ1,⋯,θk),则它的矩(原点矩和中心矩都可以,此处以原点矩为例)
依赖于θ1,⋯,θkθ1,⋯,θk。另一方面,至少在样本大小n较大时,样本原点矩amam应该接近于αmαm。于是
取m=1,⋯,km=1,⋯,k,并将上面的近似式改成等式,就得到一个方程组:
解此方程组,得起根θ1^(X1,⋯,Xn)(i=1,⋯,k)θ1^(X1,⋯,Xn)(i=1,⋯,k),就以θ̂ iθ^i作为θiθi的估计(i=1,⋯,k)(i=1,⋯,k)。如果要估计的是θ1,⋯,θkθ1,⋯,θk的某函数g(θ1,⋯,θk)g(θ1,⋯,θk),则用ĝ (X1,⋯,Xn)=g(θ̂ 1,⋯,θ̂ k)g^(X1,⋯,Xn)=g(θ^1,⋯,θ^k)去估计它。这样定出的估计量就叫做矩估计。
例子
设X1,⋯,XnX1,⋯,Xn是从正态总体N(μ,δ2)N(μ,δ2)中抽取的样本,要估计μμ和δ2δ2。μμ是总体的一阶原点矩,按矩估计,用样本的一阶原点矩,即样本均值X⎯⎯⎯⎯X¯去估计。δ2δ2是总体方差,即总体的二阶中心矩,可用样本的二阶中心矩m2m2去估计。一般地,在估计方差时,常用样本方差S2S2而不用m2m2,即对矩估计做了一定的修正。
参考书目
《概率论与数理统计》——陈希孺