贝叶斯统计ch贝叶斯决策.pptx

《贝叶斯统计ch贝叶斯决策.pptx》由会员分享，可在线阅读，更多相关《贝叶斯统计ch贝叶斯决策.pptx（41页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1第五章 贝叶斯决策第一节 贝叶斯决策问题第二节 后验风险决策第三节 常用损失函数下的贝叶斯估计第四节 抽样信息期望值第五节 最佳样本量的确定第六节 二行动线性决策问题的EVPI第1页/共41页2第一节 贝叶斯决策问题一、决策问题分类 (1)仅使用先验信息的决策问题称为无数据(或无样本信息)的决策问题;(2)仅使用抽样信息的决策问题称为统计决策问题;(3)先验信息和抽样信息都使用的决策问题称为贝叶斯决策问题.第2页/共41页3二、贝叶斯决策问题 先验信息和抽样信息都用的决策问题称为贝叶斯决策问题。若以下条件已知，则我们认为一个贝叶斯决策问题给定了。(4)定义在 上的二元函数 称为损失函数第3页

2、/共41页4三、贝叶斯决策的优缺点1.优点主要表现在:(1)贝叶斯决策充分利用各种信息,使决策结果更加科学化;(2)能对调查结果的可能性加以数量化的评价;(3)贝叶斯决策巧妙地将调查结果和先验知识有机地结合起来;(4)贝叶斯决策过程可以不断地使用,使决策结果逐步完善.2.缺点:(1)贝叶斯决策所需要的数据多，分析计算也比较复杂,如果解决的问题比较复杂时,这个矛盾就更加突出;(2)在决策的过程中，有些数据必须要使用主观概率,有些人不是很相信,这也妨碍了贝叶斯决策方法的推广和使用.O第4页/共41页5第二节 后验风险决策1.后验风险函数 我们把损失函数 对后验分布 的期望称为后验风险,记 ,即 后

3、验风险就是用后验分布计算的平均损失.11第5页/共41页62.决策函数定义5.1 在给定的贝叶斯决策问题中,从样本空间 到行动集A上的一个映照 称为该决策问题的一个决策函数,表示所有从样本空间到A上的决策函数组成的类,称为决策函数类.在贝叶斯决策中我们面临的是决策函数类D,要在D中选择决策函数 ,使其风险最小.例题分析第6页/共41页73.后验风险准则 定义 在给定的贝叶斯决策问题中 是其决策函数类,则称为决策函数 的后验风险.假如在决策函数类中存在这样的决策函数 ,它在D中有最小的风险,即则称 为后验风险准则下的最优决策函数,或称贝叶斯决策,或贝叶斯解,或贝叶斯估计贝叶斯估计。注:(1)定义

4、中的条件:给定的贝叶斯决策问题 (2)定义中的先验分布允许是广义的.第7页/共41页8例1 设 是来自正态分布N(,1,1)的一个样本。又设参数的先验分布为共轭先验分布N(0,2),其中2已知.而损失函数为0-1损失函数试求参数的贝叶斯估计。解:分三步求解:(1)求参数的后验分布(2)对于任意一个决策函数 计算后验风险函数:(3)求出使得上述风险函数达到最小时的决策函数:第8页/共41页9例2 在市场占有率的估计问题中，已知损失函数为：药厂厂长对市场占有率无任何先验信息，另外在市场调查中，在n个购买止痛剂的顾客中有x人买了新药，试在后验风险准则下对作出贝叶斯估计。解:(1)求参数的后验分布:结

5、果为 Be(x+1,n-x+1)(2),计算风险函数(3)求最优行动使上述风险函数达到最小.令:则得:(4)数值计算:第9页/共41页10例3 如何判断一个样本是来自密度函数为p0(x)的总体还是来自密度函数为p1(x)的总体。解:两个假设:问题:接受H0还是H1?(1)把假设检验问题转化为贝叶斯决策问题:参数空间=0,1行动空间A=0,1先验分布:P(=0)=0,P(=1)=1损失函数:决策正确无损失,决策错误的损失为1.则(2)求后验分布:(3)计算每个行动下的后验风险:R(a=0|x)=P(=1|x)R(a=1|x)=P(=0|x)(4)找出最佳行动,即确定拒绝域.第10页/共41页11

6、1.1.平方损失函数下的贝叶斯估计平方损失函数下的贝叶斯估计定理5.1 在平方损失函数 下,的贝叶斯估计为后验均值,即Pr 在平方损失函数下,任何一个决策函数 的后验风险为第三节 常用损失函数下的贝叶斯估计第11页/共41页12定理5.2 在加权平方损失函数 下,的贝叶斯估计为:其中()为参数空间上的正值函数.定理5.3 在参数向量 的场合下,对多元二次损失函数 ,Q为正定阵,的贝叶斯估计为后验均值向量:第12页/共41页13例4 设 是来自泊松分布的一个样本.若的先验分布用其共轭先验分布G(G(,),),即其中参数与已知.求平方损失函数下的贝叶斯估计.解:解题过程分为以下三步:(1)根据题意

7、求出的后验分布(2)写出后验均值(3)结论:由定理5.1知的贝叶斯估计为:第13页/共41页14例5 设 是来自均匀分布U(0,U(0,)的一个样本.又设的先验分布为ParetoPareto分布.在损失函数分别为绝对值损失函数和平方损失函数下求的贝叶斯估计.解题步骤:第一步:求的后验分布:第二步:在绝对值损失函数下的贝叶斯估计:恰为后验分布的中位数.第三步:平方损失函数下的贝叶斯估计:Pareto分布的分布函数:密度函数为:期望:方差:中位数:第14页/共41页15例例6 6 贝叶斯决策在可靠性统计中的应用贝叶斯决策在可靠性统计中的应用 问题的描述:假设某产品的寿命T服从指数分布,其分布函数为

8、:把指定时间t0后该产品才失效的概率称为产品在t0时刻的可靠度。在平方损失函数下怎样估计可靠度R(t0)？第15页/共41页162.2.线性损失函数下的贝叶斯估计线性损失函数下的贝叶斯估计证明:(1)证明的思路:设m为(|x x)的中位数,要证明定理成立,即要证:,都有 即当m时,2-(m+)m的情形,此时:因此(3)同理可证m的情形。第16页/共41页17定理5.5 在线性损失函数:下,的贝叶斯估计n(x)为后验分布(|x x)的 分位数.证明:(1)计算任一决策函数=(x)的后验风险:(3)结论:的贝叶斯估计是后验分布的 分位数.例9(p191例题5.13)自学(2)求R(|x)的驻点:令

9、第17页/共41页183.3.有限个行动问题的假设检验有限个行动问题的假设检验(1)一般问题:设A=a1,a2,ar,在ai下的损失为L(,ai),如何从这些行动当中选择一个最优行动?(使后验期望损失 达到最小的行动)(2)特例:r=2的情形,即二个行动的假设检验问题(3)特例:r=3的情形(三个行动的假设检验问题)(儿童智商检验的实例P193)第18页/共41页19二个行动的假设检验问题二个行动的假设检验问题1.问题的描述:有两个假设:两个行动:a0:表示接受H0的行动,a1:表示接受H1的行动.决策方法:如果 ,则认为a0为最优行动;如果 ,则认为a1为最优行动.2.损失函数的确定:确定原

10、则:决策正确无损失,决策错误损失ki个单位.则得到损失矩阵:第19页/共41页203.计算每个行动的后验期望损失：假设后验分布已经求出,则4.按照后验风险准则,确定最优行动:如果 ,即:则应选择a1,拒绝a0.如果 ,即:则应选择a0,拒绝a1.5.确定拒绝域:若 ,则则由上一步可算得:即拒绝域第20页/共41页21第四节第四节 抽样信息期望值抽样信息期望值一、基本概念 1.完全信息：对需要作决策的问题，假如决策者所获得的信息足以肯定那一个状态即将发生，则该信息就称为(该状态的)完全信息。2.完全信息期望值(EVPI)：设某决策问题有n种状态1,2,n,且各种状态的先验概率(i)已知,又有m种

11、行动a1,a2,am。设Qij为出现i采取行动aj的收益，a为使 取得最大时的行动，则称为完全信息期望值，记为EVPI。第21页/共41页223.先验EVPIEVPI：在一个决策问题中()是状态集=上的先验分布。a是先验期望准则下的最优行动，则在a下的损失函数L(,a)的先验期望 称为完全信息先验期望值，记为先验EVPI。4.两者的关系:5.例题:对给定的Q或L怎样计算EVPI和先验EVPI？如：第22页/共41页23二、抽样信息期望值二、抽样信息期望值1.定义：在一个贝叶斯决策问题中,a是先验期望准则下的最优行动，是后验风险准则下的最优决策函数。则先验EVPI与后验EVPI期望值的差称为抽样

12、信息期望值，记为：2.计算一个EVSI的基本步骤：第一步：计算先验EVPI；第二步：计算的后验分布；第三步：计算每个行动的后验期望损失 ；第四步：确定最优决策函数；第五步：计算后验EVPI；第六步：计算后验EVPI的期望值；第七步：计算抽样信息期望值。第23页/共41页24案例分析 甲厂的某一零件由乙厂生产，每批1000只，其次品率的概率分布如下表所示：甲厂在整机装配时，如发现零件是次品，必须更换，每换一只，乙厂赔偿2.20元的损失费，但也可以在送装前采取全部检查的办法，使每批零件的次品率降为1%，但乙厂必须支付每只0.10元的检查费。乙厂面临如下两种选择：a1:一批中一件都不检查 a2:一批

13、中每件都检查若乙厂厂长想从每批中任取三只零件进行抽查，根据不合格品个数来决定是采取行动a1还是行动a2，并想知道这样能否带来更大的收益？0.02 0.05 0.10()0.45 0.39 0.16第24页/共41页25分析过程一、计算先验EVPI：支付函数：由此的支付矩阵和损失矩阵：计算每个行动下的先验期望损失：由此得在先验期望准则下，a1是最优行动，则：先验EVPI=15.68第25页/共41页262.计算的后验分布3.计算各行动的后验期望损失x 0 1 2 3m(x)0.8745 0.1176 0.0076 0.0002x 0 1 2 31=0.020.4843 0.2202 0.0658

14、 0.00282=0.050.3824 0.4490 0.3684 0.01903=0.100.1333 0.3308 0.5658 0.9782x 0 1 2 313.0634 32.4148 55.4484 95.863642.3642 22.5636 9.5532 0.4464第26页/共41页274.确定最优决策函数：5.计算后验EVPI：x=0时，后验EVPI=13.0634 x=1时，后验EVPI=22.5636 x=2时，后验EVPI=9.5532 x=3时，后验EVPI=0.44946.计算后验EVPI的期望值:7.计算抽样信息期望值：EVSI=15.68-14.15思考：该厂

15、长所确定的抽取三件产品检查，是否是最好？第27页/共41页28第五节 最佳样本量的确定一、抽样净益1.抽样成本:抽样费用称为抽样成本,记为 其中Cf是固定成本,Cv是可变成本.2.抽样净益:从抽样信息期望值中扣除抽样成本所剩下的净收益.记为ENGS,即二、最佳样本量及其上界1.最佳样本量:使得抽样净益达到最大的样本量n*称为最佳样本量即2.上界的确定:第28页/共41页29三、最佳样本量的求法第29页/共41页30第六节 二行动线性决策问题的EVPI一、正态分布下二行动线性决策问题的先验EVPI二、贝塔分布下二行动线性决策问题的先验EVPI三、伽玛分布下二行动线性决策问题的先验EVPI第30页

16、/共41页31案例分析 某厂的产品每100件装成一箱运交顾客.在向顾客交货前面临如下二中选择:a1:一箱中逐一检查;a2一箱中一件也不检查.若工厂选择a1,则可保证交货时每件产品都是合格品,但因每件产品的检查费为0.8元,为此工厂要支付检查费80元/箱.若工厂选择a2,虽可免去每箱80元的检查费,但一旦顾客发现不合格品时,不仅要免费更换,而且每件还要支付12.5元的赔偿费.假设该厂产品的不合格率没有超过0.12的记录.根据下列不同的情形,作出相应的最优决策.(1)不作任何抽样,怎样做决策?（解答）(2)每箱中抽取两件进行检查,结果如何?（解答）(3)每箱中抽取三件进行检查,结果又如何?B 第3

17、1页/共41页32x012决策函数第32页/共41页33x01222.203015.04832.798622.203015.048336.892422.203055.97832.798622.203055.978336.892415.615915.04832.798615.615915.048336.892415.615955.97832.798615.615955.978336.8924各决策函数的后验风险第33页/共41页34（1）第一步：由支付函数确定损失函数：第34页/共41页35第二步：计算每一个行动的先验期望损失B第35页/共41页36（2）第一步：由支付函数确定损失函数：第36页/共41页37第二步：计算后验分布：第37页/共41页38第38页/共41页39第三步：利用后验风险准则做决策 计算每一个决策行动的后验风险第39页/共41页40B第40页/共41页41感谢您的观看。第41页/共41页