我们从数据中能得到以下信息:
总体信息。总体所属分布或者所属的分布族带来的信息;
样本信息。从总体中抽样得来的样本给我们提供的信息;
- 以上两种信息进行的统计推断称为经典统计学。它的观点是把样本看成来自具有一定概率分布的总体。
先验信息。在抽样之前,对总体的基本认知,一般来自经验或历史资料。
- 利用以上三种信息进行的统计推断称为贝叶斯统计。它的观点是:任一未知量θ\thetaθ都可看做一个随机变量,应用一个概率分布去描述对θ\thetaθ的未知状况。这个概率分布是在抽样前就有的关于θ\thetaθ的先验信息的概率陈述。这个分布被称之为先验(
Prior
)分布。
关于未知量θ\thetaθ的一些讨论:
- 依赖于参数θ\thetaθ的密度函数在经典统计中记为p(x;θ)p(x;\theta)p(x;θ)或pθ(x)p_{\theta}(x)pθ(x),它表示在参数空间Θ={θ}\Theta=\{\theta\}Θ={θ}中不同的θ\thetaθ对应不同的分布。可以在贝叶斯统计中记为p(x∣θ)p(x|\theta)p(x∣θ),他表示在随机变量θ\thetaθ给定某个值时,总体指标XXX的条件分布。
- 根据参数θ\thetaθ的先验信息确定先验分布π(θ)\pi(\theta)π(θ)。
- 从贝叶斯的观点看,样本x=(x1,⋅⋅⋅Xn,⋅⋅⋅)x=(x_1 ,···X_n,···)x=(x1,⋅⋅⋅Xn,⋅⋅⋅)的产生分两步进行。首先设想从先验分布π(θ)\pi(\theta)π(θ)产生一个样本θ\thetaθ,这一步是“老天爷”做的,人们是看不到的,故用“设想”二字。第二步是从总体分布p(x∣θ)p(x|\theta)p(x∣θ)产生一个样本x=(x1,⋅⋅⋅xn,⋅⋅⋅)x=(x_1,···x_n,···)x=(x1,⋅⋅⋅xn,⋅⋅⋅),这个样本是具体的,人们能看得到的,此样本xxx发生的概率是与如下联合密函数成正比。p(x∣θi)=∏i=1np(xi∣θi)p(x|\theta^i)=\prod_{i=1}^n{p(x_i|\theta^i)}p(x∣θi)=∏i=1np(xi∣θi)这个联合密度函数是综合了总体信息和样本信息,常称为似然函数,记为L(θi)L(\theta^i)L(θi)。频率学派和贝叶斯学派都承认似然函数,二派认位:在有了样本观察值x=(x1,⋅⋅⋅xn,⋅⋅⋅)x=(x_1,···x_n,···)x=(x1,⋅⋅⋅xn,⋅⋅⋅)后,总体和样本所含θ\thetaθ的信息都被包含在似然函数L(θi)L(\theta^i)L(θi)之中,可在使用似然函数做统计推断时,两派还是有差异的。
- 由于θ\thetaθ是设想出来的,他仍然是未知的,他是按先验分布π(θ)\pi(\theta)π(θ)而产生的,要把先验信息进行综合,不能只考虑θ\thetaθ,而应对θ\thetaθ的一切可能加以考虑。故要用π(θ)\pi(\theta)π(θ)参与进一步综合。这样一来,样本xxx和参数θ\thetaθ的联合分布h(x,θ)=p(x∣θ)π(θ)h(x,\theta)=p(x|\theta)\pi(\theta)h(x,θ)=p(x∣θ)π(θ)把三种可用的信息都综合进去了。
- 我们的任务是要对未知数θ\thetaθ做出统计推断。在没有样本信息时,人们只能根据先验分布对θ\thetaθ做出判断。在有样本观察值x=(x1,⋅⋅⋅xn,⋅⋅⋅)x=(x_1,···x_n,···)x=(x1,⋅⋅⋅xn,⋅⋅⋅)后,我们应该依据h(x,θ)h(x,\theta)h(x,θ)对θ\thetaθ作出推断。为此我们需要把h(x,θ)h(x,\theta)h(x,θ)作如下分解:h(x,θ)=π(θ∣x)m(x)h(x,\theta)=\pi(\theta|x)m(x)h(x,θ)=π(θ∣x)m(x)其中m(x)m(x)m(x)是xxx的边缘密度函数。m(x)=∫θh(x,∣θ)dθ=∫θp(x∣θ)π(θ)m(x)=\int_\theta{h(x,|\theta)d\theta}=\int_\theta{p(x|\theta)\pi(\theta)}m(x)=∫θh(x,∣θ)dθ=∫θp(x∣θ)π(θ)他与θ\thetaθ无关,或者说是,m(x)m(x)m(x)中不含θ\thetaθ的任何信息。因此能用来对θ\thetaθ做出推断的仅是条件分布π(θ∣x)\pi(\theta|x)π(θ∣x)。他的计算公式为π(θ∣x)=h(x∣θ)m(x)=p(x∣θ)π(θ)∫θp(x∣θ)π(θ)dθ\pi(\theta|x)=\frac{h(x|\theta)}{m(x)}=\frac{p(x|\theta)\pi(\theta)}{\int_\theta{p(x|\theta)\pi(\theta)}d\theta}π(θ∣x)=m(x)h(x∣θ)=∫θp(x∣θ)π(θ)dθp(x∣θ)π(θ),这就是贝叶斯公式的密度函数形式。这个在样本xxx给定下,θ\thetaθ的条件分布被称为θ\thetaθ的后验分布。他是集中了总体、样本和先验三种信息中包含有θ\thetaθ的一切信息,而又是排除一切与θ\thetaθ无关的信息之后所得到的结果。故基于后验分布π(θ∣x)\pi(\theta|x)π(θ∣x)对θ\thetaθ进行统计推断是更为有效,也是合理的。
- 在θ\thetaθ是离散随机变量时,先验分布可用先验分布列pi(θi),i=1,2⋅⋅⋅\\pi(\theta_i),i=1,2···pi(θi),i=1,2⋅⋅⋅,表示。这时后验分布也是离散形式。π(θi∣x)=p(x∣θi)π(θi)∑ip(x∣θi)π(θi),i=1,2,⋅⋅⋅\pi(\theta_i|x)=\frac{p(x|\theta_i)\pi(\theta_i)}{\sum_i{p(x|\theta_i)\pi(\theta_i)}},i=1,2,···π(θi∣x)=∑ip(x∣θi)π(θi)p(x∣θi)π(θi),i=1,2,⋅⋅⋅假如总体XXX也是离散的,那么只要把密度安徽省农户p(x∣θ)p(x|\theta)p(x∣θ)看作是概率函数P(X=x∣θ)P(X=x|\theta)P(X=x∣θ)即可。
一般来说,先验分布π(θ)\pi(\theta)π(θ)是反映人们在抽样分布前对θ\thetaθ的认识,后验分布π(θ∣x)\pi(\theta|x)π(θ∣x)是反映人们在抽样后θ\thetaθ的认识。之间的差异是由于样本xxx出现后人们对θ\thetaθ认识的一种调整。所以后验分布π(θ∣x)\pi(\theta|x)π(θ∣x)可以看作是人们用总体信息和样本信息对先验分布π(θ)\pi(\theta)π(θ)做调整的结果。