ワイブルフィット

概要

ワイブルフィットは、生存関数と故障時間の間の関係を分析するパラメータ法の1つです。解析によって、ワイブル分布のハザード関数と生存関数を決めるパラメータの推定値を得ることができます。

ワイブル分布:

f(x)=\frac{c}{\sigma }({\frac{x-\theta }{\sigma }})^{c-1}\exp (-\left({\frac{x-\theta }{\sigma }}\right)^{c}) ここで c,\sigma > 0 のとき、x >\theta

生存関数:

S(x)=\exp(-\left({\frac{x-\theta }{\sigma }}\right)^{c})

ハザード関数:

h(x)=\frac{c}{\sigma }\left({\frac{x-\theta }{\sigma }}\right)^{c-1}

ここで、cは形状パラメータ、\sigma はスケールパラメータ、\theta は位置パラメータです。Originでは、OriginではWeibull fitはcと \sigma についてのみ議論し、\theta = 0と仮定します。

もしc>1ならハザードは増加し、c=1の場合にはハザードは一定(exponentioal model)になり、c<1なら減少します。

必要なOriginのバージョン: OriginPro 9.1 SR0以降

学習する項目

このチュートリアルでは、以下の項目について解説します。

  • ワイブルフィットを実行
  • 分析レポートの説明方法

ワイブルフィットを実行

  1. 新しいワークブックを用意します。ヘルプ: フォルダを開く: サンプルフォルダを選択して、サンプルフォルダを開きます。このフォルダ内のStatisticsサブフォルダにあるweibull fit.dat ファイルを探します。空のワークシートにファイルをドラッグアンドドロップしてインポートします。
  2. 統計:生存分析:ワイブルフィット選択してダイアログを開きます。
  3. 時間範囲A(X)列を入力します。同様に、打切り範囲B(Y)列を入力します。
  4. 打切り値のドロップダウンリストから打切り値として0を選択します。
  5. プロットの項目で、生存プロットハザードプロットにチェックを付けます。
    W01.png
  6. OK ボタンをクリックして、ワイブルフィットを実行します。
W04.png

結果の解釈

分析レポートのワークシートWeibullFit1を開きます。

  • 「イベントと打ち切り変数のサマリー」表から、打切りは19で打切りパーセントは0.2111である事がわかります。
W02.png
  • 「パラメータ推定の解析」の表では、Weibull 分布の全てのパラメータの推定値を入手できます。 切片=\theta=4.1959( \theta は小さな極値分布の切片、\theta = ln(Weibull Scale))
    Weibullスケール= \sigma =66.4153, Weibull 形状=c=2.0204
    スケール=0.495 (scale=1 / c)
W03.png
c>1なので、時間とともにハザードが増加すると結論付けられます。
  • さらに、生存関数とハザード関数を得ることができます。
W05.png W06.png