统计推断 @ Probability & Statistics Chapter 7
这是 DeGroot, M. H., & Schervish, M. J. (2012). Probability and statistics. Pearson Education. 一书的笔记。
7.1 Statistical Inference 统计推断
这一章首首先介绍了灯泡寿命问题。公司想要估计电子元件的failure rate,用未知的参数表示,它服从指数分布(gamma分布)。根据大数定律6.2 我们知道独立同分布的样本均值 (依概率收敛)。根据Theorem 6.2.5 我们知道给这个式子套一层连续函数也是成立的,即 . 所以,我们可以用样本均值的倒数来近似failure rate。
因为我们不可能测量所有元件的寿命,参数只能是hypothetically observable的。
接下来定义了统计模型和统计推断(statistical inference)。统计推断的例子:
- 构造一个随机变量 使得 。 这个RV 可以是观测到的一系列样本的函数。
- 在观测到一系列样本后,我们有多大把握认为。
后面提了下模型参数和参数空间。接着定义了统计量(statistic):是样本的实值函数,例如样本均值和样本最大值.
7.2 先验 & 后验
先验分布就是在做出观测之前,我们已知或假定参数在参数空间上的可能分布。
e.g. 灯泡寿命为指数分布,其损坏率假定其先验分布为Gamma分布。
灯泡的损坏率服从Gamma分布,它的概率密度为
注意这里会被绕晕:先验概率分布是theta的函数,刻画了theta取各个值的可能性的密度。一般好像RV会用大写,而p.d.f.会用小写?
然后介绍了后文model样本分布的一个notation,即。这里,它反过来model了给定条件下的分布,每个观察值都是在给定theta时条件独立同分布(conditional i.i.d.) 的,这时我们可以说构成了概率密度函数的一组随机样本(form a random sample)。
Sensitivity Analysis是分析不同的Prior产生的Posterior有多大程度的不同。这里有一个fun fact是如果实验观测数据足够大,那么不同的先验对试验结果的影响是微乎其微的。这提示我们不必过于纠结怎么挑选先验 (experimenters might be less inclined to spend time specifying a prior distribution if it is not going to matter much),但如果不选择先验的话,就没法计算后验概率了。
后验分布是给定一组样本下的条件概率。The conditional distribution of θ given is called the posterior distribution of θ, typically denoted as .
与后验分布相关的最重要的东西就是贝叶斯定理了,在随机样本的形势下,贝叶斯定义可以表述为
其中为样本的联合概率分布(可由分数线上面部分对积分,求得样本的边际分布得到),整个式子即后验=先验x似然 / 证据 (,似然刻画了由先验到后验结果的可能性,证据则是各种结果的出现的可能性)
这节例子为对之前先验分布为Gamma的灯泡损坏率求后验。
- 似然为给定theta下每个样本的条件概率密度,由于灯泡寿命为指数分布,
- evidence是样本的联合概率分布,
- 最后的后验概率为 ,是一个关于的函数,其中样本通过观察得到,视为已知的量。
- 图像上的改变:
