求助：全新产品零失效可靠性验证方案的置信区间如何计算？

xxex

liaojenyi 發表於 2014-4-23 20:22
？那如果我再增加一些约束条件呢？比如相似产品的数据分析或经验推断，甚至组织人拍脑袋？
May be you can  ...

请问,廖老师，bayesian和贝叶斯是什么关系？和条件概率又是什么关系？

hlperng

xxex 發表於 2014-4-23 21:14
请问,廖老师，bayesian和贝叶斯是什么关系？和条件概率又是什么关系？

貝氏(Bayes' or Bayesian)就是中國大陸稱的貝葉斯，貝葉斯定理(Bayes' theorem or Bayesian theorem)，都是以條件機率為基礎的統計推論技法。

一般經由測試資料計算得到的統計推論值稱為概似函數，過去的測試或經驗值所推論得到的稱為先驗機率(a prior)，後驗機率(a posterior)為先驗機率 x 概似函數。

雖然貝氏理論的推導與驗證過程相當複雜，不過根據貝氏定理進行二項分布與波桑分布參數的推定，數學模型很直覺。

二項分布的參數為 p，定義 p' 為先驗機率參數、p" 為後驗機率參數，p 為貝他分布。假設試驗樣本數為 n，不符數(或不良數或失效數)為 r ，則根據試驗結果的概似函數推定 p 為：



若先驗樣本數為 n'、不符數為 r' ，則先驗參數 p' 為：



則根據貝氏定理，後驗參數 p" 的計算式為：



不符率(或不良率或失效機率)為不符數(不良數或失效數)除以樣本數，這是基本概念，學理上是可以以貝氏定理印證的。已知試前的樣本數為 n'、不符數為 r' ，當下試驗的樣本數為 n、不符數為 r。假如可以證明兩次試驗的樣本取自同一群體，則兩個樣本可混合在一起、當作同一個新的樣本，混合後的樣本數為 n" = n'+n、不符數為 r" = r'+r，有了樣本數(n")與不符數(r")，自然可算出混合後的不符率(p")。


在波桑分布(或指數分布)也有類似的程序，波桑分布參數為 λ，λ' 為先驗機率參數、λ" 為後驗機率參數，λ 為伽瑪分布。假設試驗總時間為 T、失效數為 r，則試驗結果的概似函數推定 λ 為：



若先驗試驗時間為 T'、失效數為 r'，則先驗參數 λ' 為：


則根據貝氏定理，後驗參數λ" 的計算式為：



其道理同二項分布。

應用貝氏定理最大的困擾在於，先驗參數 p' 或 λ' 是已知的，但是計算 p' 的 r' 與 n' ，或 計算 λ' 的 r' 與 T' ，卻有各種不同的組合，結果是完全不同的。因此，有人評論貝氏方法相當主觀，沒有一致的結果，這也是貝氏定理奧妙的地方。重點在於，雖然當下看到的是一樣，但是來源不同，其表現應該會有不同。統計推論的目的在於推論未來的可能性。

舉例如下：
試驗樣本 n =10、不符數 r=1，所以不符率 p = 0.1。假設先驗參數 p'=0.2，

r	1	1	1	1	1
n	10	10	10	10	10
p	0.1	0.1	0.1	0.1	0.1
p'	0.2	0.2	0.2	0.2	0.2
r'	2	20	200	2000	20000
n'	10	100	1000	10000	100000
r"	1+2=3	1+20=21	1+200=201	1+2000=2001	1+20000=20001
n"	10+10=20	10+100=110	10+1000=1010	10+10000=10010	10+100000=100010
p"	0.15	0.1909091	0.1990099	0.1999001	0.1999900

xxex

hlperng 發表於 2014-4-24 08:51
Bayesian就是中國大陸稱的貝葉斯，貝葉斯定理(Bayes' theorem)，都是以條件機率為基礎的統計推論技法。

...

请问，廖老师，Bayes和Bayesian是否同一人?

hlperng

xxex 發表於 2014-4-24 09:40
请问，廖老师，Bayes和Bayesian是否同一人?

是！是同一人的！
貝氏定理或貝葉斯定理：http://zh.wikipedia.org/wiki/%E8 ... F%E5%AE%9A%E7%90%86
貝氏統計或貝葉斯統計：http://wiki.mbalib.com/zh-tw/%E8 ... F%E7%BB%9F%E8%AE%A1
貝氏法則或貝葉斯法則：http://wiki.mbalib.com/zh-tw/%E8 ... F%E6%B3%95%E5%88%99
貝氏公式或貝葉斯公式：http://baike.baidu.com/view/541856.htm

hlperng

你的問題應該有兩個可能：

• 參數抽樣檢定問題：0 收 1 退，抽樣計畫的應用，必須知道兩個風險，而不是信心水準。類似IC 產品鑑定試驗中的壽命試驗(life test)。例如：22個樣本，執行1000小時壽命試驗，假設型 I 風險 (α) 為 5 %、型 II 風險 (β) 為 10%，假如通過試驗，代表產品的失效率優於2,300 fit (λ < λ₀ = 2,300 fit）、若不通過試驗，代表產品的失效率比 100,000 fit 還要差 (λ > λ₁ = 100,000 fit)。（參閱：http://redi.org.tw/forum.php?mod=viewthread&tid=85&extra=page%3D1）
• 參數推定問題：根據試驗結果，計算失效率的 95%心水準的信心上限。
  

在產品壽命(或失效發生時間) T 符合指數分布的假設下，求得失效率後，平均壽命(或平均失效發生時間)、可靠度的推定值(或估計值)，可根據公式計算得到。

根據壽命試驗資料，計算失效率的 95% 信心水準的信心上限， λ_U， (失效率為望小，悲觀為上限)推定值為：



平均壽命或平均失效發生時間的 95% 信心水準的信心下限推定值，MTBF_L ，(平均壽命為望大，悲觀為下限)推定值為：



任務時間為 t 時的可靠度的 95% 信心水準的信心下限，R_L ， (可靠度為望大，悲觀為下限)推定值為：

xxex

hlperng 發表於 2014-4-24 16:30
你的問題應該有兩個可能：
1. 參數抽樣檢定問題：0 收 1 退，抽樣計畫的應用，必須知道兩個風險，而不是信 ...

请问大师，两个风险怎么算？

hlperng

xxex 發表於 2014-4-24 16:50
请问大师，两个风险怎么算？

這兩個風險不是算出來，是必須在決定抽樣計畫之前，事先決定的。就像所舉例，型 I 風險或稱生產者風險 (α) 取 5 %、型II 風險或稱消費者風險 (β) 取 10%。

xxex

hlperng 發表於 2014-4-24 16:59
兩個風險不是算出來，是必須在決定抽樣計畫之前，事先決定的。就像所舉例，型 I 風險或稱生產者風險 (α)  ...

请问大师，两类风险和第3，第n类 风险的关系