无失效数据的贝叶斯分析.pdfVIP

无失效数据的贝叶斯分析 翟艳敏 (山东经济学院 统计与数学学 院，济南 250014) 摘 要 ：对于无失效数据 问题 ，文章利用多层先验分布给 出了失效概率的多层Ar叶斯估计 ，并证 明其 “保序性” 关键词 ：贝叶斯 分析 ；无失效数据 ；失效概率 中图分类号 ：0212 文献标识码 ：A 文章编号 ：1002—6487(2008)11-0016—03 0 引言 选取 Beta分布为第一层先验分布 ，即： a-t b-t 在可靠性试验中，由于受试验时 问、费用等条件的限制 ， 订t(pila,b)= (1) 常采用截尾方式获得数据。在定时截尾试验中，有 时会遇到 无失效数据 (即在规定 的时间内无一样品发生失效)。随着产 其中B(a，b)=J。x(1-x)dx为贝塔函数，a，bl为两个 品质量的 日益提高 。在实际中无失效情况更易产生 。向文献 超参数。 『11发表 以来 ，对无失效数据的研究逐渐引起重视 ，并取得一 由于在无失效情况下。失效概率大的可能性小，小的可 定的成果 。 能性大 。故先验分布应选取减 函数 。随着 a，b的变化 ，由 MatlablSl~61软件得到Beta密度 曲线p(x)有如下几种类型 (见图 1 无 失 效 数 据 问题 的统 计 模 型 1) 从一批产品中随机地抽取 n个样品．分成 k组 。各组样 k 品数分别为n。，n，…，nk，且 ni=n，分别对各组样品做定时 i= 1 截尾试验 ，截尾时间分别为tl，t2，… ，tk，且 tt2…tk，结果无 一 失效 。若记 si=n．+n¨+…+n，表示在时刻 t，前未失效的样品 (a)al，bl 总数 ，则根据试验获得一组无失效数据ft…S)；i=l，2，…，k。 问题是利用上述信息估计分布 数中未知参数或产品 可靠性指标 。若产 品寿命分布类型 已知 ，如指数分布 、正态分 布、威布尔分布等 ，则可根据文献 2【]提 出的思路首先估计失 效概率pi=F(t，)=P{Tt，)，然后利用配合 曲线法得到分布函数中 参数 的估计 ，从而再得到可靠性指标的估计。对关键性的的 (c)a1，0bl (d)Oa1b1 估计多采用贝叶斯估计 ，其不 同是先验分布 的选取 。文献 2【] 图 1 几种典型的密度曲线 选取均匀分布为先验分布 ，文献 3【]选取不完全 Beta分布为 从 图形上看 。当al，bl时，密度函数递减。如何具体确 先验分布。文献 4【]选取 Garoma分布为先验分布。另外，我们 定 a。b的值?一方面 ，由于两参数 a，b是无信息 ，根据贝叶斯 应注意到失效概率必须满足 0p。≤p2≤…≤p1，为 了满足 假设 ，取 (0，1)，(1，c)上 的均匀分布分别 a，b作为 的先验分 布 。于是第二层先验分布密度为： “保序性”，文献 3【]对范围加以限制 ，取pp，。 本文提出以Beta分布为先验分布 。利用多层贝叶斯方 1T2(a)=1(0