covariateSurvivalModel

残存耐用期間を推定する比例ハザード生存モデル

説明

covariateSurvivalModel は、コンポーネントの残存耐用期間 (RUL) を比例ハザード生存モデルを使用して推定するために使用します。このモデルは、コンポーネントの寿命に関する履歴情報と関連する共変量を使用してテスト コンポーネントの生存確率を示します。共変量は環境変数または説明変数で、コンポーネントの製造元や動作状態などがあります。共変量生存モデルは、同じ仕様で製造された複数のマシンなどの類似コンポーネントのアンサンブルに関する故障時間と関連する共変量のデータしかない場合に便利です。生存モデルの詳細については、比例ハザード生存モデルを参照してください。

特定のタイプのコンポーネントの covariateSurvivalModel オブジェクトを構成するには、fit を使用します。これにより、故障時間データと関連する共変量の集合を使用してモデルの係数が推定されます。共変量生存モデルのパラメーターを構成した後、類似コンポーネントの残存耐用期間を predictRUL を使用して予測できます。RUL 予測の基本的な例については、データの着信に応じた RUL 予測の更新を参照してください。

寿命の測定値のみがあり、共変量の情報がない場合は、reliabilitySurvivalModel を使用します。

残存耐用期間の予測に関する一般的な情報については、残存耐用期間を予測するモデルを参照してください。

作成

構文

mdl = covariateSurvivalModel

mdl = covariateSurvivalModel(initModel)

mdl = covariateSurvivalModel(___,Name,Value)

説明

mdl = covariateSurvivalModel は、RUL を推定するための共変量生存モデルを作成し、既定の設定でモデルを初期化します。

mdl = covariateSurvivalModel(initModel) は、共変量生存モデルを作成し、既存の covariateSurvivalModel オブジェクト initModel を使用してモデルのパラメーターを初期化します。

mdl = covariateSurvivalModel(___,Name,Value) は、ユーザー設定可能なモデル プロパティを、名前と値のペアで指定します。たとえば、covariateSurvivalModel('LifeTimeUnit',"days") は、日数をライフタイムの単位として使用する共変量生存モデルを作成します。複数の名前と値のペアを指定できます。各プロパティ名を引用符で囲みます。

入力引数

initModel — 共変量生存モデル

covariateSurvivalModel オブジェクト

共変量生存モデル。covariateSurvivalModel オブジェクトとして指定します。

プロパティ

BaselineCumulativeHazard — ベースラインのハザード率関数

2 列配列

この プロパティ は読み取り専用です。

生存モデルのベースラインのハザード率。2 列の配列として指定し、関数 fit で推定します。2 列目にベースラインの生存時間関数の値が含まれ、1 列目に対応するライフタイム値が含まれます。

生存モデルの詳細については、比例ハザード生存モデルを参照してください。

EncodingMethod — 符号化手法

"dummy" (既定値) | "binary"

EncodedVariables のカテゴリカル特徴の符号化手法。次のいずれかとして指定します。

• "dummy" — N 個のカテゴリがあるカテゴリカル特徴について、(N - 1) ビットを使用して変数を符号化します。

• "binary" — バイナリ符号化

次の場合に、EncodingMethod を指定できます。

• モデル作成時に名前と値のペアを使用する

• モデル作成後にドット表記を使用する

Standardize — 共変量の特徴を標準化するためのフラグ

false (既定値) | true

コックス回帰のパラメーターの計算時に共変量の特徴を標準化するためのフラグ。logical 値として指定します。Standardize が true の場合、共変量 X が (X-mean(X))/std(X) になるように数値の共変量変数が標準化されます。

標準化は、符号化されたカテゴリカル変数には影響しません。

次の場合に、Standardize を指定できます。

• モデル作成時に名前と値のペアを使用する

• モデル作成後にドット表記を使用する

Ties — 同順位の故障時間を扱う方法

"breslow" (既定値) | "efron"

同順位の故障時間を扱う方法。"breslow" または "efron" のいずれかとして指定します。これらの方法の詳細については、コックス比例ハザード モデルを参照してください。

次の場合に、Ties を指定できます。

• モデル作成時に名前と値のペアを使用する

• モデル作成後にドット表記を使用する

Options — 数値と表示の設定

構造体

コックス回帰の数値と表示の設定。statset('coxphfit') を使用して作成された構造体として指定します。構造体に含まれるオプションはドット表記を使用して変更できます。

次の場合に、Options を指定できます。

• モデル作成時に名前と値のペアを使用する

• モデル作成後にドット表記を使用する

ParameterValues — 共変量乗算係数

ベクトル

この プロパティ は読み取り専用です。

生存モデルの共変量乗算係数。スカラーとして指定し、関数 fit で推定します。生存モデルの詳細については、比例ハザード生存モデルを参照してください。

ParameterCovariance — 共変量乗算係数の共分散

配列

この プロパティ は読み取り専用です。

共変量乗算係数の共分散。係数の数と等しいサイズの正の配列として指定し、関数 fit で推定します。

ParameterNames — 共変量乗算係数の名前

string 配列

この プロパティ は読み取り専用です。

共変量乗算係数の名前。string 配列として指定し、モデルに学習させるときに関数 fit を使用して割り当てます。

数値の共変量に対応する係数は、DataVariables の対応するデータ変数と同じ名前になります。符号化された変数の場合は、係数の名前に EncodedVariables の対応する符号化変数の名前と符号化ビット順序の表現が含まれます。

CensorVariable — 打ち切り変数

"" (既定値) | string

打ち切り変数。有効な MATLAB^® 変数名を含む string として指定します。打ち切り変数は、data のどの寿命測定値が耐用寿命末期の値でないかを示す 2 値変数です。

CensorVariable は、DataVariables または LifeTimeVariable のいずれの string とも一致してはなりません。

次の場合に、CensorVariable を指定できます。

• モデル作成時に名前と値のペアを使用する

• 関数 fit を呼び出す際の引数として使用する

• モデル作成後にドット表記を使用する

LifeTimeVariable — ライフタイム変数

"" (既定値) | string

ライフタイム変数。有効な MATLAB 変数名を含む string として指定します。生存モデルの場合、コンポーネントの過去の寿命測定値がライフタイム変数に格納されます。

次の場合に、LifeTimeVariable を指定できます。

• モデル作成時に名前と値のペアを使用する

• 関数 fit を呼び出す際の引数として使用する

• 手動でドット表記を使用する

LifeTimeUnit — ライフタイム変数の単位

"" (既定値) | string

ライフタイム変数の単位。string として指定します。

ライフタイム変数の単位は時間ベースとする必要はありません。テスト コンポーネントの寿命は、走行距離 (マイル) や消費燃料 (ガロン) などの使用量変数で測定できます。

DataVariables — 共変量データ変数

"" (既定値) | string | string 配列

共変量データ変数。string または string 配列として指定します。DataVariables の string は、有効な MATLAB 変数名でなければなりません。共変量は環境変数または説明変数とも呼ばれます。

次の場合に、DataVariables を指定できます。

• モデル作成時に名前と値のペアを使用する

• 関数 fit を呼び出す際の引数として使用する

• モデル作成後にドット表記を使用する

EncodedVariables — 符号化された共変量変数

"" (既定値) | string | string 配列

符号化された共変量変数。string または string 配列として指定します。EncodedVariables の string は、有効な MATLAB 変数名でなければなりません。符号化された変数は、通常は非数値のカテゴリカル特徴で、fit で当てはめる前に数値ベクトルに変換されます。少数のセットから値を取る logical 値または数値を符号化の対象として指定することもできます。

符号化の方法を指定するには、EncodingMethod を使用します。

次の場合に、EncodedVariables を指定できます。

• モデル作成時に名前と値のペアを使用する

• 関数 fit を呼び出す際の引数として使用する

• モデル作成後にドット表記を使用する

EncodedVariables の string は、DataVariables の string のサブセットでなければなりません。

UserData — その他のモデル情報

[] (既定値) | 任意のデータ型またはデータ形式

記録用のその他のモデル情報。任意のデータ型またはデータ形式として指定します。モデルでは、この情報を使用しません。

次の場合に、UserData を指定できます。

• モデル作成時に名前と値のペアを使用する

• モデル作成後にドット表記を使用する

オブジェクト関数

predictRUL	テスト コンポーネントの残存耐用期間を推定
fit	履歴データを使用して残存耐用期間モデルのパラメーターを推定
plot	共変量残存耐用期間モデルの生存時間関数のプロット

例

共変量生存モデルの学習

学習データを読み込みます。

load('covariateData.mat')

このデータには、電池の放電時間や関連する共変量情報が含まれています。共変量変数は以下のとおりです。

• 温度

• 負荷

• 製造元

製造元情報は、符号化する必要があるカテゴリカル変数です。

共変量生存モデルを作成します。

mdl = covariateSurvivalModel;

学習データを使用して、ライフタイム変数、データ変数、符号化された変数を指定し、生存モデルに学習させます。この学習データには打ち切り変数がありません。

fit(mdl,covariateData,"DischargeTime",["Temperature","Load","Manufacturer"],[],"Manufacturer")

Successful convergence: Norm of gradient less than OPTIONS.TolFun

モデルのベースラインの生存時間関数をプロットします。

plot(mdl)

共変量生存モデルを使用した RUL の予測

学習データを読み込みます。

load('covariateData.mat')

このデータには、電池の放電時間や関連する共変量情報が含まれています。共変量変数は以下のとおりです。

• 温度

• 負荷

• 製造元

製造元情報は、符号化する必要があるカテゴリカル変数です。

共変量生存モデルを作成し、学習データを使用して学習させます。

mdl = covariateSurvivalModel('LifeTimeVariable',"DischargeTime",'LifeTimeUnit',"hours",...
   'DataVariables',["Temperature","Load","Manufacturer"],'EncodedVariables',"Manufacturer");
fit(mdl,covariateData)

Successful convergence: Norm of gradient less than OPTIONS.TolFun

たとえば、製造元が B 社で、30 時間稼働したバッテリー パックがあるとします。使用時間 DischargeTime、測定された周囲温度 TestAmbientTemperature、電流 TestBatteryLoad を含むテスト データ table を作成します。

TestBatteryLoad = 25;
TestAmbientTemperature = 60; 
DischargeTime = hours(30);
TestData = timetable(TestAmbientTemperature,TestBatteryLoad,"B",'RowTimes',hours(30));
TestData.Properties.VariableNames = {'Temperature','Load','Manufacturer'};
TestData.Properties.DimensionNames{1} = 'DischargeTime';

バッテリーの RUL を予測します。

estRUL = predictRUL(mdl,TestData)

estRUL = duration
   38.332 hr

バッテリーの共変量データの生存時間関数をプロットします。

plot(mdl,TestData)

アルゴリズム

比例ハザード生存モデル

covariateSurvivalModel オブジェクトは、次の比例ハザード生存モデルを実装します。

$$h\left(X,t\right)=h_0\left(t\right)e^{b^{T}X}$$

ここで、

• X は共変量値のベクトルです。

• b は共変量乗算係数のベクトルです。これらの係数は、モデルの ParameterValues プロパティに対応します。

• h₀(t) はベースラインのハザード率関数です。これは、モデルの BaselineCumulativeHazard プロパティに対応します。

• h(X,t) は t の時点の X のハザード率です。

このモデルのパラメーターを見つけるために、関数 fit では関数 coxphfit を使用します。

比例ハザード モデルの詳細については、コックス比例ハザード モデルを参照してください。

拡張機能

C/C++ コード生成

MATLAB® Coder™ を使用して C および C++ コードを生成します。

使用に関するメモと制限:

• predictRUL コマンドは、この RUL モデル タイプの MATLAB Coder™ を使用したコード生成に対応しています。このモデルを使用するコードを生成する前に、saveRULModelForCoder を使用してモデルを保存する必要があります。例については、残存耐用期間の予測用コードの生成を参照してください。

• covariateSurvivalModel については、table または timetable の形式で string または文字ベクトルの値として指定された共変量データを使用して predictRUL のコードを生成することはできません。代わりに、string または文字ベクトルのデータを別のデータ型に変換するか、データを table から削除してください。

• 読み取り専用のプロパティに加え、生存モデルの次のプロパティは実行時に変更できません。

  • LifeTimeVariable

  • LifeTimeUnit

  • DataVariables

  • CensorVariable

  • EncodedVariables (covariateSurvivalModel のみ)

  • EncodingMethod (covariateSurvivalModel のみ)

• コード生成に関するその他の制限については、predictRUL を参照してください。

バージョン履歴

R2018a で導入

参考

関数

• coxphfit | reliabilitySurvivalModel | fit | predictRUL

トピック

• データの着信に応じた RUL 予測の更新
• RUL 推定器モデルを使用した RUL の推定
• コックス比例ハザード モデル