第3章先验分布的确定

合集下载

第3章先验分布的确定

第3章先验分布的确定

m(
x)
p( x
|
)
(
)d, 为连续时
p( x | ) ( ), 为离散时
(3.1)
当先验分布含有未知参数时,譬如π(θ)=π(θ|λ),那么
边缘分布m(x)依赖于λ,可记为m(x|λ).
例 设总体X~N(, σ2) 其中σ2已知。
取另一正态分布N(μπ,τπ2)作为正态均值的先验分布 , 则可以算得X的边缘分布为N(μπ,τπ2+σ2)
2.若还有σ2 (θ)= σ2(与θ无关的常数 ),则
m2 2 2
其中
2
E
2
例3.14 设X
~
N
(
,1),的先验分布类为{N
(
,
2
)}
主观经验预测X的均值为1, 方差为3,即m
10和
2 m
3,
试确定其先验分布.
解 :已知m
1;
2 m
3,由X
~
N ( ,1)知 2
1,利用推论1,
( ) m 10
2.注意的问题
(1)所咨询专家应是声誉良好的和富有经验的
(2)这两个方法相比,决策者更愿意使用变分度法.
§3.3利用边缘分布m(x)确定先验密度
一、边缘分布m(x)
设总体X的密度函数为p(x|θ),它含有未知参数θ,若θ 的先验分布选用形式已知的密度函数π(θ),则可算得 X的边缘分布(即无条件分布)

d
d
2
[ln
(
2
)]
2(
2
n
2
)
ns2
2(
2
2
)2
0

2

第3篇先验分布的确定

第3篇先验分布的确定

3.假如有历史数据,要尽量利用,帮助形成初步概念, 然后再做一些对比修正,再形成个人信念.
注意:1.利用先验信息确定主观概率没有固定模式; 2.主观概率必须满足概率的3条公理.
总结 1.理解主观概率的定义 2.了解主观概率确定的常用方法
§3.2 利用先验信息确定先验分布 在贝叶斯统计方法中关键的一步是确定先验分布。1.
要进一步分析先验信息.先验信息很分散;柯西分布
先验信息较为集中:正态分布
3.两个先验分布都满足给定的先验信息。
(1)如果两个先验分布差别不大,对后验分布的影响 也不大,那可任选一个。
(2)假如面临两个差异较大的先验分布可供选择时, 应慎重选择。因不同的选择对后验分布的影响也会很大.
三、定分度法与变分度法 两种方法的共同点:通过咨询专家获得各种主观概 率,然后经过整理加工可得到累积概率分布曲线.
例3.12 设总体X~N(, σ2) 其中σ2已知。
取另一正态分布N(μπ,τπ2)作为正态均值的先验分布
则可以算得X的边缘分布为N(μπ,τπ2+σ2)
设X 在给定时条件分布为N ( ,
2 ),
(
)
~
N (
,
2
),
则边缘分布m( x
)
~
N (
,
2
2)
n
由 m( X ) m(xi )
i 1
第三章 先验分布的确定
§3.1 主观概率 §3.2 利用先验信息确定先验分布 §3.3 利用边缘分布m(x)确定先验密度 §3.4 无信息先验分布 §3.5 多层先验
概率的公理化定义
定义:设Ω 为一个样本空间, F为Ω的某些子集组成 的一个事件域,如果对任一事件A∈F,定义在F上 一个实值函数P(A)满足下列条件:

先验分布的确定

先验分布的确定

先验分布的确定幻灯⽚67 其步骤如下:(1)写出样本的对数似然函数∑∏====ni i n i i x p x p x l 11)|(ln )|(ln )|(θθθ(2)求样本的信息阵pj i l E I j i x ,...,2,1,,)(2|=-=θθθθ2|2(),x l I Eθθθ=-??在单参数(p=1)场合,(3)Θ的⽆信息先验密度为2/1)]([det )(θθπI =1/2()[()]I πθθ=在单参数(p=1)场合,幻灯⽚682122(,,...,)(,),(,).n X x x x N Jeffreys µσθµσ==设是来⾃正态分布的⼀组样本试求的先验2211:()ln[]2i x ni l x e µσθπσ--==∑写出样本的对数似然函数22111(,)ln(2)ln ().22ni i l n x µσπσµσ=?=---∑22222222()()0:(,);20()()ll nE E Fisher I n ll E E µµσσµσσµσσ-- ??==?? -- ??其信息阵42),(det -=?σσµn I22,(,):(,)2.Jeffreys n µσπµσσσ--=∝所以的先验为幻灯⽚6911:,(),:()1;,()2,();,(,);nI I n σµπµσµσσπσσµσπµσσ---=∝=∝∝注当已知当已知当和独⽴幻灯⽚70 例3.22关于成功概率的⽆信息先验分布⾄今已有4种π1(θ)=1 ——正常π2(θ)=θ-1(1-θ)-1 ——不正常π3(θ)=θ-1/2(1-θ)-1/2 ——正则化后可成为正常π4(θ)=θθ(1-θ)(1-θ) ——正则化后可成为正常注意:1.⼀般说来,⽆信息先验不是唯⼀的.但它们对贝叶斯统计推断的影响都很⼩,很少对结果产⽣较⼤的影响2.任何⽆信息先验都可以采⽤。

一先验分布和后验分布

一先验分布和后验分布

2[ E( | x)][E( | x) d( x)]h( | x)d
又因为 E( | x) h( | x)d 则 [ E( | x)][E( | x) d( x)]h( | x)d
[E( | x) d( x)][ E( | x)]h( | x)d
[E( | x) d( x)][E( | x) E( | x)] 0
由第一小节内容可知,给定损失函数以后,风险函数定 义为
R( , d ) E (L( , d( X ))
L( , d( x))q( x | )dx
此积分仍为的函数,在给定的先验分布()时,定义
R(d ) E (R( , d ))
R( , d )π( )d
为决策函数d在给定先验分布()下的贝叶斯风险,简 称为d的贝叶斯风险.
0,
0,
x 0,
设Y X 1,则Y的密度函数为
f
(
y;
,
)
(
)
(
1 y
)
1
e
y,y
0,
0,
y 0,
此分布密度为倒分布的密度函数, 设 ²的先验分布为倒
分布,即
(
2
)
(
)
(1
2
) 1
e 2,y
0,
0,
y 0,
则 ²的后验分布为
h( 2 | x) q( x | 2 )π( 2 )
min R(d ) min m( x){ [ d( x)]2 h( | x)d }dx
min a.s [ d( x)]2 h( | x)d
又因为
[ d( x)]2 h( | x)d
[ E( | x) E( | x) d( x)]2 h( | x)d

贝叶斯统计3.4,3.5教材

贝叶斯统计3.4,3.5教材

27
例3.22
关于成功概率的无信息先验分布至今已有4种 π1(θ)=1 π2(θ)=θ-1(1-θ)-1 π3(θ)=θ-1/2(1-θ)-1/2 ——正常 ——不正常 ——正则化后可成为正常
π4(θ)=θθ(1-θ)(1-θ) ——正则化后可成为正常
注意:1.一般说来,无信息先验不是唯一的.
但它们对贝叶斯统计推断的影响都很小,很少对结 果产生较大的影响
2.任何无信息先验都可以采用。
28
总结
1. 掌握贝叶斯假设
2.掌握位置参数和尺度参数的无信息先验分布
3.会用Fisher信息阵确定无信息先验
29
§3.5 多层先验
一、多层先验 二、多层模型
30
一、多层先验
1.定义
当所给先验分布中超参数难于确定时,可以对超参数 再给出一个先验,第二个先验称为超先验。由先验和 超先验决定的一个新先验称为多层先验。
试求分布参数 与的无信息先验.
取为位置参数, 为尺度参数, 令 1, ln( ), w ln( x), 则有
p( w; , )
1
w
d * 由随机变量函数知, ( ) ( ) 1 , 2 ( ) 1 , d
浙江财经学院本科教学课程经济数学三概率统计精品文档贝叶斯统计34352第三章先验分布的确定31主观概率32利用先验信息确定先验分布33利用边缘分布mx确定先验密度34无信息先验分布35多层先验334无信息先验分布一贝叶斯假设二位置尺度参数族的无信息先验三用fisher信息阵确定无信息先验4所谓参数??的无信息先验分布是指除参数??的取值范围和??在总体分布中的地位之外再也不包含??的任何信息的先验分布
例3.23 设对某产品的不合格品率了解甚少,只知道 它比较小。现需要确定θ的先验分布。决策人经过 反复思考,最后把他引导到多层先验上去,他的思 路是这样的: (1)开始他用(0,1)上的均匀分布U(0,1)作为θ的先 验分布。

第三章 先验分布的确定

第三章 先验分布的确定


则称该分布族{p(x | , )} 为Cramer-Rao正则分布族,简称 C—R正则族。

在C—R正则族前提下,记分向量S (x) 的方差协方差阵

' 称为该分布族中参数 (1 ,..., p ) 的Fisher信息阵,简称θ 的信息阵。


二、Jeffreys先验 设总体密度函数为p(x | , ) ,又设参数θ 的无信息先验 ( ) ,由于一一 为π (θ )。若对参数θ 作一一对应变换; 对应变换不会增加或减少信息,故新参数η 的无信息先验 * ( ) 与 ( ) 在结构上应完全相同,即 ( ) * ( ) 。另一方 面,按随机变量函数的运算规则,θ 与η 的密度函数间应满 足如下关系式

假如混合样本 x (x1 , x 2 ,..., x n ) 所涉及的先验密度函数的形式 已知,未知的仅是其中的超参数,即先验密度函数族可表示 如下:
{ ( | ), }

ˆ 使得 这时寻求ML-II先验是较为简单的事,只要寻求这样的

这可用最大化似然函数方法来实现。
第一节 主观概率 第二节 利用先验信息确定先验分布 第三节 利用边际分布m(x) 确定先验密度 第四节 无信息先验分布 第五节 多层先验


3.1.1 主观概率 贝叶斯统计中要使用先验信息,而先验信息主要是指经验和 历史资料。因此如何用人们的经验和过去的历史资料确定概 率和先验分布是贝叶斯学派要研究的问题。 贝叶斯学派认为:一个事件的概率是人们根据经验对该事件 发生可能性所给出个人信念。这样给出的概率称为主观概率。


此外σ 的参数空间与η 的参数空间都为 R ,可见(X,σ )问题 ( ) 与(y,η )问题的统计结构完全相同,故σ 的无信息先验 与η 的无信息先验 * ( ) 应相同,即 ( ) * ( ) 另一方面,由变换 c 可以得η 的无信息先验

贝叶斯讲义先验分布的确定解析

贝叶斯讲义先验分布的确定解析

13
14
15

说明:如果有两个甚至多个先验分布都满足给定的先 验信息,则要看情况选择:假如这两个先验分布差异 不大,对后验分布影响也不大,则可任选一个;如果 我们面临着两个差异极大的先验分布可供选择时,一 16 定要根据实际情况慎重选择。
三、定分度法与变分度法
基本概念: (1)定分度法:把参数可能取值的区间逐次分 为长度相等的小区间,每次在每个小区间上 请专家给出主观概率. (2)变分度法:该法是把参数可能取值的区间 逐次分为机会相等的两个小区间,这里的分 点由专家确定. 例3.9(自学)
17
§3.3 利用边缘分布m(x)确定先验密度
一、边缘分布m(x)
二、混合分布 三、先验选择的ML-II方法
四、先验选择的矩方法
18
一、边缘分布m(x) 设总体X的密度函数为p(x|θ ),它含有未 知参数θ ,若θ 的先验分布选用形式已知的 密度函数π (θ ),则可算得X的边缘分布(即 无条件分布): 当为连续时 p( x | ) ( )d , m( x) p ( x | ) ( ) , 当 为离散时
P( Ai ) P( Ai )
i 1 i 1


(2)如果发现所确定的主观概率与上述三个公理 及其推出的性质相悖,必须立即修正。直到两者一 致为止。(例3.5)
9
10
§3.2 利用先验信息确定先验分布
一、直方图法 二、选定先验密度函数形式再估计其超参数 三、定分度法与变分度法
11
一、直方图法
23
24
三、先验选择的ML-Ⅱ方法
定义:设 { ( | ), } 为所考虑的先 验类,且x=(x1,x2,…,xn)是来自Г 中某一 ˆ ) 满足(对 ˆ ( 分布的样本,若存在 n 观测数据x): ˆ ) sup m(x | m( xi | )

一先验分布和后验分布

一先验分布和后验分布

i 1
i 1
第二种方法 设总体X的分布密度为p(x|),统计量
T ( X ) T ( X1X2 , , X n )是参数的充分统计量,则有
定理3.1 设f ( )为任一固定的函数,满足条件
(1) f ( ) 0, ,
(2) 0 gn(t | ) f ( )d

D f
{
gn (t | ) f ( ) gn (t | ) f ( )d
由第一小节内容可知,给定损失函数以后,风险函数定 义为
R( , d ) E (L( , d( X ))
L( , d( x))q( x | )dx
此积分仍为的函数,在给定的先验分布()时,定义
R(d ) E (R( , d ))
R( , d )π( )d
为决策函数d在给定先验分布()下的贝叶斯风险,简 称为d的贝叶斯风险.
X N (9.80, 0.12 )
这个信息就是重力加速度的先验信息.
在统计学中,先验信息可以更好的帮助人们解决 统计决策问题. 贝叶斯将此思想应用于统计决策中,形成了 完整的贝叶斯统计方法.
2、先验分布
对未知参数的先验信息用一个分布形式()来 表示,此分布()称为未知参数的先验分布.
例如 例1中重力加速度的先验分布为
则的贝叶斯估计为
d*( x) E(( ) | x) E(( ) | x)
证明略,此证明定理3.2的证明类似.
定理3.4 设参数为随机向量,先验分布为()
和损失函数为二次损失函数
L( , d ) (d )T Q(d )
1:改进生产设备后,高质量产品可占90%,
:改进生产设备后,高质量产品可占70%,
2
经理根据以往的经验,两个顾问建议可信度分别为
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1.把参数空间分成一些小区间.
2.在每个小区间上决定主观概率或依据历史数据确定其 频率.
3.绘制直方图
4.在直方图上做一条光滑的曲线,此曲线就是( ).
注意:这样得到的先验密度常常仅限于有限区间上,有 时使用也不方便。第二种方法更为适用.
二、选定先验密度函数形式再估计其超参数
1.要点
(1).根据先验信息选定的先验密度函数的形式( ) 。
F(x)称为F(x|θ1)和F(x|θ2) 的混合分布。 这里的π和1-π可以看作一个新的随机变量θ的分布,即
P( 1) , P( 2 ) (1 )
F(x) F(x |1) (1 )F(x |2 )
P( 1) , P( 2 ) (1 )
从混合分布F(x)中抽取一个样品x1,相当于如下两个抽样: 第一次,从π(θ) 中抽取一个样品θ。 若θ= θ1,则从F(x|θ1)中再抽一个样品,这个样品就 是x1,若θ=θ2 ,则从F(x|θ2)中再抽一个样品,这个 样品就是x1. 若从混合分布抽取一个容量为n的样本x1, x2,…,xn,那么 其中约有nπ(θ1) 个来自F(x|θ1),约有nπ(θ2)个来自 F(x|θ2),这样的样本有时也称为混合样本.
2
}求其最大值则可 )

d
d
2
[ln
(
2
)]
2(
2
n
2
)
ns2
2(
2
2
)2
0

2
S2
2
注意到 : 若S 2
2时
2
S2
2为最大点
而若S 2
2时
2
S2
2
0则取
2
0为最大点
ML2
ˆ ~ N (ˆ ,ˆ2 )
, 其中ˆ
1 n
n i 1
xi
,ˆ2 max{0, S 2 2}.
四、先验选择的矩方法
矩方法用于先验密度函数形式π(θ|λ)已知,利用先验矩 与边缘分布矩之间的关系寻求超参数的估计.
定理 设总体分布若p(x|θ)的期望μ(θ)和方差σ2 (θ),假设
分别为Xm ,的m边2 缘分布m(x)的均值与方差,假设以上值都存 在,则
m E ( ) ,
2 m
E
2( )
E ( ) m 2
例3.11
三、先验选择的ML—Ⅱ方法
m(x)
p( x | ) ( )d,为连续时
p( x | ) ( ),为离散时
(3.1)
若p(x|θ)已知,则m(x)大小反映π(θ)的合理程度,这里
把m(x)记为m(x|π) 或mπ(x) 是由无限个不可数的密度
函数混合而成.
设Γ为所考虑的先验类, ˆ 满足(对观察值x)
例如,许多经济现象都是不能重复或不能大量重复 的随机现象。
在经典统计中有一种习惯,对所得到概率都要给出 频率解释。这些在有些场合是难于作出的。
例如:天气预报。
§3.1 主观概率
一、主观概率
1.定义:一个事件的概率是人们根据经验对该事件发 生的可能性大小所给出的个人信念,这样给出的概率 称为主观概率
素(Ω, F ,P)为概率空间(Probability space) .
概率是定义在σ-域F上的一个非负的、正则的、可列可加的集函数
在经典统计中,概率用非负性、正则性和可加性三 条公理定义的。
确定概率的方法主要有两种。一是古典方法(包括 几何方法),另一种是频率方法。实际中大量使用的 是频率方法,所以经典统计的研究对象是能大量重复 的随机现象,不是这类随机现象就不能用频率的方法 去确定其有关事件的概率。这无疑把统计学的应用和 研究领域缩小了。
推论1 若μ(θ)=θ,则 m E (先验均值).
2.若还有σ2 (θ)= σ2(与θ无关的常数 ),则
m2 2 2 其中 2 E 2
先验选择的矩方法的步骤如下(当先验分布中有两个超 参数时) (1, 2 )
1.定义
定分度法是把参数可能取值的区间逐次分为长度相等的 小区间,每次在每个小区间上请专家给出主观概率.
变分度法是把参数可能取值的区间逐次分为机会相等的 两个小区间,这里分点由专家确定.
2.注意的问题
(1)所咨询专家应是声誉良好的和富有经验的
(2)这两个方法相比,决策者更愿意使用变分度法.
例3.9 一开发商希望获知一个新建仓库的租金可能 达到的水平是什么?为此向一位推租经纪人咨询。
2.用专家意见来确定主观概率的方法(最常用的).
注意:(1).向专家提的问题要设计好,既要使专家 易懂又要使专家回答不是模棱两可。 (2).要对专家本人比较了解,以便做出修正,形成 决策者自己的主观概率.
(3).通过向多位专家咨询后,经修正和综合获得主 观概率,关键在于把问题设计好,便于往后综合,即 在提出问题时,就要想到如何综合。
当总体参数是离散时,即参数空间Θ只含有限个或可
数个点时,可对Θ中每个点确定一个主观概率。
2.当总体参数是连续时,即参数空间Θ是实数轴或其 上某个区间时,要构造一个先验密度( ),就有些困
难了.
当的先验信息足够多时,下面有三个方法可供使用.
直方图法
选定先验密度函数形式再估计其超参数
定分度法与变分度法
一、直方图法
则可以算得X的边缘分布为N(μπ,τπ2+σ2)
二、混合分布
设随机变量X以概率π在总体F1中取值,以概率1-π在 总体F2中取值.若F(x|θ1)和F(x|θ2)分别是这两个总体 的分布函数,则X的分布函数为
F(x) F(x |1) (1 )F(x |2 )
或用密度函数或概率函数表示
p(x) p( x |1) (1 ) p( x |2)
第三章 先验分布的确定
§3.1 主观概率 §3.2 利用先验信息确定先验分布 §3.3 利用边缘分布m(x)确定先验密度 §3.4 无信息先验分布 §3.5 多层先验
概率的公理化定义
定义:设Ω 为一个样本空间, F为Ω的某些子集组成 的一个事件域,如果对任一事件A∈F,定义在F上 一个实值函数P(A)满足下列条件:
m(x)
p( x | ) ( )d,为连续时
p( x | ) ( ),为离散时
(3.1)
当先验分布含有未知参数时,譬如π(θ)=π(θ|λ),那么
边缘分布m(x)依赖于λ,可记为m(x|λ). 例3.10 设总体X~N(, σ2) 其中σ2已知。
取另一正态分布N(μπ,τπ2)作为正态均值的先验分布
1 n
n i 1
xi
, s2
1 n
n
( xi
i 1
x )2
m(x
)
[2
(
2
2
n
)] 2
exp{
ns2
2(
2
2
}exp{ )
n(x
2(
2
)2
2
} )
由exp{ 故只需令
(2((x2 )2[2)2(2)}2知, 不2 )]论 n2 exp2如 { 何 2(xn2 s2
可使m(x )达最大
要进一步分析先验信息.先验信息很分散;柯西分布
先验信息较为集中:正态分布
3.两个先验分布都满足给定的先验信息。
(1)如果两个先验分布差别不大,对后验分布的影响 也不大,那可任选一个。
(2)假如面临两个差异较大的先验分布可供选择时, 应慎重选择。因不同的选择对后验分布的影响也会很大.
三、定分度法与变分度法 两种方法的共同点:通过咨询专家获得各种主观概 率,然后经过整理加工可得到累积概率分布曲线.
•用主观方法确定经验的例子
(1)明天下雨的概率为60%
(2)某新产品在未来市场上畅销的概率为80% (3)我班研究生考取大概为25%
说明:1.主观概率不是随意决定的,而是要求当事人 对所考察的事件有较透彻的了解和丰富的经验,甚至 是这方面的专家。并能对周围信息和历史信息进行仔 细分析,在这个基础上确定的主观概率就能符合实际。 所以应把主观概率与主观臆造,瞎说一通区别开来。
其中,
m(xi
)
N
(
,
2
2)
于是
m(x )
n i 1
1
1
[2
(
2
2)2
n
exp{ (xi )2 }
n
2(
2
2)
(xi )2
[2
(
2
2)2
exp{
i 1
2(
2
2)
}
[2
(
2
2
)]
n 2
exp{
n
n
(xi x )2
i 1
n
2(
2
2)
}exp{
n(x
2(
2
)2
2
} )
取x
例3.12 设总体X~N(, σ2) 其中σ2已知。
取另一正态分布N(μπ,τπ2)作为正态均值的先验分布
则可以算得X的边缘分布为N(μπ,τπ2+σ2)
设X 在给定时条件分布为N ( ,
2 ),
(
)
~
N (
,
2
),
则边缘分布m( x
)
~
N (
,
2
2)
n
由 m( X ) m(xi )
i 1
(1)非负性公理:对于每一事件A,有P(A)≥0;
(2)正则性(规范性)公理:P(Ω)=1;
(3)可列可加性(完全可加性)公理:设A1,A2,…是互
不相容的事件,即对于i≠j,AiAj=,i,j=1,2,…,则有
P(
Ai
)
P(
相关文档
最新文档