このページの翻訳は最新ではありません。ここをクリックして、英語の最新版を参照してください。
この例では、打ち切りを使用した寿命データを分析する方法を示します。この分析の生物学的または医学的な応用は、生存時間分析として知られており、生物の生存時間、または病気が治癒するまでの時間を表す場合もあります。この分析の工学的な応用は信頼性分析として知られており、装置の故障までの時間を表すことがあります。
この例では,自動車の燃料噴射システムのスロットルが故障するまでの時間をモデル化しています。
寿命データには,他のタイプのデータにはない,顕著な特徴があります。第1に寿命データの値は常に正であり,通常は時間を表します。第2に寿命データの中には,正確に観測されず,一部の値よりも大きいことのみがわかるようなデータがあります。第3の特徴は,よく使用される分布と分析方法が寿命データに特化したものであるということです。
100個のスロットルを故障するまでテストした結果のシミュレーションを実行してみましょう。ほとんどのスロットルは寿命がかなり長く、少数のスロットルは非常に早く故障する傾向があった場合に観測されるデータを生成します。
rng(2,“龙卷风”); 寿命=[wblrnd(15000,3,90,1);wblrnd(1500,3,10,1)];
この例では、現場で実際に 100時間使用した状況を各テスト時間で再現できるように、厳しい条件下でスロットルをテストすると仮定します。多くの場合、実践的な理由のため、信頼性試験は一定の時間が経過すると停止されます。この例では、140時間を採用します。これは、実際の稼働時間が 14,000 時間に達した状況に相当します。試験中に故障するものもあれば、140時間故障しないものもあります。実際の試験では、後者の時間は 14,000 として記録されます。これを、シミュレートされたデータで再現します。また、故障時間を並べ替えることもよく行われます。
T = 14000;obstime = sort(min(T, lifetime));
試験で故障しなかったスロットルもいずれは故障することは理解できますが、実際の故障時間を観測できるほど試験期間は長くありません。この試験から判断できるのは、それらのスロットルの寿命が 14,000 時間以上であることのみです。このような値を、打ち切られた値と呼びます。次のプロットは、データの約 40% が 14,000 で打ち切られていることを示しています。
失败= obstime (obstime < T);nfailed =长度(失败);幸存= obstime (obstime = = T);nsurvived =长度(存活);censored = (obstime >= T);情节([0(大小(obstime)), obstime]”,repmat(1:长度(obstime), 2, - 1),...“颜色”,“b”,“线型”,'-')行([T;3e4],repmat(n失败+(1:n被接受),2,1),“颜色”,“b”,“线型”,“:”); 行([T;T],[0;n失败+n无效],“颜色”,“k”,“线型”,'-')文本(T, 30岁,“<--这里的未知生存时间”)包含(的生存时间); 伊拉贝尔(“观察数量”)
データの分布を観測する前に、確率分布を調べるさまざまな方法を検討してみましょう。
確率密度関数 (PDF)は、さまざまな時間における故障の相対的確率を示します。
生存時間関数は,生存の確率を時間の関数として与えます。1から累積分布関数を引くだけで得られます(1-CDF)。
ハザード率は、特定の時間まで生存している場合の、故障の瞬間的な確率を与えます。PDFを生存時間関数で除算することによって得られます。この例では、ハザード率が上昇していることがわかります。これは、時間が経過するにつれて故障しやすくなること (経年劣化) を意味しています。
確率プロットは再スケーリングされた CDFであり、データを近似分布と比較するために使用します。
ワイブル分布を使用した、これら 4.つのプロット タイプの例を以下に示しています。ワイブル分布は、寿命データのモデル化によく使用されています。
x=linspace(130000);子地块(2,2,1);地块(x,wblpdf(x,14000,2),x,wblpdf(x,18000,2),x,wblpdf(x,14000,1.1))标题(的概率。密度Fcn”次要情节(2,2,2);情节(x, 1-wblcdf (x, 14000, 2), x, 1-wblcdf (x, 18000, 2), x, 1-wblcdf (1.1 x 14000))标题(“幸存者Fcn”)子地块(2,2,3);wblhaz=@(x,a,b)(wblpdf(x,a,b)/(1-wblcdf(x,a,b));地块(x,wblhaz(x,14000,2),x,wblhaz(x,18000,2),x,wblhaz(x,14000,1.1))标题(“风险率Fcn”)子批次(2,2,4);probplot(“威布尔”wblrnd(14000 2, 40岁,1))标题(“概率图”)
ワイブル分布は、指数分布を一般化したものです。寿命が指数分布に従う場合、寿命のハザード率は一定です。これは、寿命が経年変化しないことを意味します。つまり、ある区間において故障が観測される確率は、その区間の開始までの生存が与えられていれば、区間がどこで始まろうが関係ない、ということを意味します。ワイブル分布では、ハザード率が増減することがあります。
寿命データのモデル化に使用される他の分布として,対数正規分布,ガンマ分布,およびバーンバウム・サンダース分布が挙げられます。
可能な各生存時間までに故障する比率を示す、データの経験累積分布関数をプロットします。点線の曲線によって、これらの確率の 95% 信頼区間が与えられます。
次要情节(1 1 1);[empF x, empFlo, empFup] = ecdf (obstime,“审查”、审查);楼梯(x,empF);保持在;楼梯(x、empFlo、,“:”);楼梯(x, empFup,“:”);持有从包含(“时间”); 伊拉贝尔(“比例失败”);标题(“经验CDF”)
このプロットからわかるように,たとえば,4000時間までに故障する比率は約12%であり,この時間までの故障確率の95%信頼限界は6% ~ 18%です。14000年試験は時間しか実施されていないため,経験累積分布関数を使用して計算できるのは,その限度までの故障確率にすぎません。14000年データの約半分がで打ち切られているため,経験累積分布関数は,1.0ではなく約0.53まで上がるだけです。
ワイブル分布は通常、装置故障に適したモデルです。関数wblfit
は、ワイブル分布をデータ (打ち切りを伴うデータを含む) にあてはめます。パラメーター推定を計算した後、それらの推定を使用して近似ワイブル分布の累積分布関数を評価します。累積分布関数の値は推定されたパラメーターに基づくため、それらの信頼限界を計算します。
paramEsts = wblfit (obstime,“审查”,截取);[nlogl,paramCov]=wbllike(参数,obstime,截取);xx=linspace(1,2*T,500);[wblF,wblFlo,wblFup]=wblcdf(xx,参数(1),参数(2),paramCov);
経験累積分布関数と近似累積分布関数のプロットを重ね合わせることができます。こうすることで、スロットルの信頼性データがワイブル分布によって、どの程度適切にモデル化されるかを判断できます。
楼梯(x,empF);保持在处理=情节(xx, wblF“r-”,xx,wblFlo,“:”xx wblFup,“:”);持有从包含(“时间”); 伊拉贝尔(“安装失败概率”);标题(“威布尔模型与经验模型”)
ワイブルモデルを使用すると,試験が終了した後の時間における故障確率を予想し,計算することができます。しかし,近似曲線がデータと十分に一致しないように思われます。2000年時間より前の初期故障が,ワイブルモデルで予想されるよりもはるかに多く,その結果,約7000 ~ 13000年約時間における初期故障が非常に少なくなっています。これは驚くことではありません。この種の動作を伴うデータを生成したことを思い出してください。
统计和机器学习工具箱™に用意されている事前定義済み関数には,このような過剰な初期故障がある分布は含まれていません。その代わり,関数ksdensity
を使用して,経験累積分布関数を介した滑らかなノンパラメトリック曲線を描画することができます。ワイブル累積分布関数の信頼帯を取り除き,2つの曲線を追加します。一方の曲線は既定の平滑化パラメーターを使用したもので,他方の曲線は既定値の1/3の平滑化パラメーターを使用したものです。平滑化パラメーターが小さいほど,曲線がより密接にデータに追従します。
delete(handles(2:end)) [npF,ignore,u] = ksdensity(obstime,xx,)“岑”,审查,“函数”,“cdf”); 生产线(xx,npF,“颜色”,“g”);npF3 = ksdensity (obstime, xx,“岑”,审查,“函数”,“cdf”,“宽度”,u/3);生产线(xx,npF3,“颜色”,“米”);xlim (1.3 * T[0])标题(“威布尔和非参数模型与经验模型”)传奇(“经验的”,“拟合威布尔”,'非参数,默认',“非参数,1/3违约”,...“位置”,“西北”);
より小さい平滑化パラメーターを使用したノンパラメトリック推定は,十分にデータと一致します。しかし,経験累分布関数と同様に,テストが終了した後にノンパラメトリックモデルを外挿することは不可能です。つまり,推定された累積分布関数が,最後の観測の後は平坦になります。
このノンパラメトリック近似のハザード率を計算し、データの範囲全体にわたってプロットしてみましょう。
hazrate = ksdensity (obstime, xx,“岑”,审查,“宽度”u / 3) / (1-npF3);情节(xx, hazrate)标题(“非参数模型的风险率”) xlim ([0, T])
この曲線は "バスタブ" のような形状をしており、ハザード率が 2,000 付近で高く、その後、低い値になり、次に再び上昇しています。これは、寿命の早期に故障しやすく (初期故障)、寿命の後期に再び故障しやすい (経年劣化) コンポーネントに典型的なハザード率です。
また,ノンパラメトリックモデルにおける打ち切られていない最大の観測を超えてハザード率を推定することはできず,グラフがゼロまで下がることに注意してください。
この例で使用した、シミュレートされたデータの場合、ワイブル分布は適切な近似ではないことがわかりました。ノンパラメトリック近似を使用すると、十分にデータを近似できましたが、そのモデルは、そのデータの範囲内でのみ役に立ちました。
代替法は,異なるパラメトリック分布を使用するというものです。统计和机器学习工具箱には,対数正規分布,ガンマ分布,バーンバウム・サンダース分布など他の一般的な寿命分布向けの関数と,寿命モデルではあまり使用されない他の分布向けの関数が含まれています。また,カスタム一変量分布の近似 (第 2.部)の例で説明されているように、カスタム パラメトリック モデルを定義して、寿命データにあてはめることもできます。
さらに別の方法は、2.つのパラメトリック分布 (一方は初期故障を表し、他方は残りの分布を表す) を組み合わせて使用するというものです。分布の混合のあてはめについては、カスタム一変量分布の近似の例で説明します。