贝叶斯统计茆诗松版本大部分课后练习习题答案.docxVIP

习题讲解 一、 1,3,5,6,10,11,12,15 记样本为 x. p( x 0.1) C82 *0.1 2 *0.9 6 0.1488 p( x 0.2) C82 *0.2 2 *0.8 6 0.2936 后验分布： 0.1 x 0.1488*0.7 0.5418 0.1488*0.7 0.2936*0.3 0.2 x 0.2936*0.3 0.4582 0.1488*0.7 0.2936*0.3 m x p(x | ) d C83 3(1 )5 *2(1 )d 112 3 (1 )6 d 2 1 1 1 0 0 0 15 x p(x | ) 840 3 (1 )6 ,0 1 m x 由题意设 x 表示等候汽车的时间，则其服从均匀分布 U (0, ) 1 , 0 x p( x) 0, 其它 因为抽取 3 个样本，即 X ( x1, x2 , x3 ) ，所以样本联合分布为 p( X ) 13 , 0 x1 , x2, x3 0, 其它 又因为 192/ 4 , 4 ( ) 0, 4 所以，利用样本信息得 h( X , ) p( X ) ( ) 1 192 192 8,0 x1, x2 , x3 ) 3 4 7( 于是 m( X ) 8 h( X , )d 8 1927d 的后验分布为 h( X , ) 192 / 7 6 86 ( X ) 192 7 m( X ) d 8 7 6 86 8 ( X ) 7 , 0, 8 样本联合分布为： p( x ) 1n ,0 x ( ) 0 / 1 , 0 0, 0 ( x) p( x ) ( ) 0 / n 1 1/ n 1, 1 max 0 , x1,L , xn 因此 的后验分布的核为 1/ n 1 ，仍表现为 Pareto 分布密度函数的核 即 ( n) 1 n / n 1 , 1 ( x) 0, 1 即得证。 n 1 样本的似然函数： p( x ) n e xi ne nx i 1 ( ) 1e 参数 的后验分布 ( x) p( x ) ( ) 服从伽马分布 Ga n, nx . 0.0002 (2) 4, 20000. 0.00012  n 1e ( nx ) 二、 1,2,3,5,6,7,8,10,11,12 解： 由题意，变量 t 服从指数分布： p(t ) e t 样本联合分布 ( ) n t i p T e 且 ~ Ga( , ) ( ) 1e , 0，E() 0.2 Var ( ) 1 由伽玛分布性质知： 0.2 0.04, 0.2 1 又已知 n=20, t 3.8 n n ti 20 3.8 76 ，所以 n 20.04, ti 76.2 i 1 i 1 由于伽玛分布是指数分布参数的共轭先验分布，而且后验分布 ( t) p(T ) ( ) ne ti 1e n 1e ( ti ) 即后验分布为 Ga( n, ti ) Ga(20.04,76.2) E |T( ) n ti 20.04 0.263 76.2 1 服从倒伽玛分布 IGa( n, ti ) IGa (20.04,76.2) E |T( ) E |T( 1) ti n 4.002 1 可以算出 的后验分布为 Ga (11,4) ， 的后验期望估计的后验方差为 11 . n 36 16 . 的先验分布为： ( ) 0 / 1, 0 0, 0 令 1 max 0 , x1 ,L , xn 可得后验分布为： ( x) ( n) 1 n / n 1, 1 0, 1 则 的后验期望估计为： E( x) ( n) 1 ， n 12 1 后验方差为： Var ( x) ( n) . 1)2 (n (n 2) 由 x ~ Ga( n , 1 ), ~ IGa( , ) 可以得出 2 2 ( 1 )n2 n 1 p(x ) 2 x 2 1 x 0 e 2 , x n ( ) ( ) ( 1)e, 0 ( ) （1） 的后验分布为： ( x) p(x ) ( ) 即为倒伽玛分布 IGa ( n , x ) 的核。 2 2 所以 的后验分布为 IGa( n , x ) 2 2 x x 2 （2）后验均值为 E( x) 2 n 1 n 2 2 2  x 2 (1) 2 e 2 ( x ) 2 后验方差为 Var ( x) 2 ( n 1)2 ( n 2) 2 2 （3）样本分布函数为： n n 1 n nn n (2 ) 2 1 2 xi p( x ) p(xi ) xi2 e i 1 (n / 2) i 1 i 1 所以 的后验分布为