ドキュメンテーション

最新のリリースでは、このページがまだ翻訳されていません。 このページの最新版は英語でご覧になれます。

covarianceParameters

クラス: GeneralizedLinearMixedModel

一般化線形混合効果モデルの共分散パラメーターの抽出

説明

psi = covarianceParameters(glme) は一般化線形混合効果モデル glme における変量効果予測子に推定される前の共分散パラメーターを返します。

[psi,dispersion] = covarianceParameters(glme) は、分散パラメーターの推定も返します。

[psi,dispersion,stats] = covarianceParameters(glme) は、共分散パラメーターの推定および関連する統計を含む cell 配列 stats も返します。

[___] = covarianceParameters(glme,Name,Value) は、1 つ以上の Name,Value のペアの引数により指定される追加オプションを使用し、上記から任意の出力引数を返します。たとえば、共分散パラメーターの信頼限界に対する信頼度を指定できます。

入力引数

すべて展開する

一般化線形混合効果モデル。GeneralizedLinearMixedModel オブジェクトとして指定します。このオブジェクトのプロパティとメソッドについては、GeneralizedLinearMixedModel を参照してください。

名前と値のペアの引数

オプションの Name,Value 引数のコンマ区切りペアを指定します。Name は引数名で、Value は対応する値です。Name は引用符で囲まなければなりません。Name1,Value1,...,NameN,ValueN のように、複数の名前と値のペアの引数を、任意の順番で指定できます。

有意水準。'Alpha' と [0,1] の範囲にあるスカラー値から構成されるコンマ区切りのペアとして指定します。値が α の場合、信頼度は 100 × (1 – α)% です。

たとえば、99% の信頼区間の場合は、次のように信頼度を指定できます。

例: 'Alpha',0.01

データ型: single | double

出力引数

すべて展開する

変量効果予測子に推定される前の共分散パラメーター。長さ R の cell 配列として返されます。ここで、R はモデル内で使用されたグループ化変数の数です。psi{r} には、グループ化変数 gr に関連付けられた変量効果の共分散行列が含まれ、r = 1, 2, ..., R になります。psi のグループ化変数の順序は、モデルを近似したときに入った順序と同じです。グループ化変数の詳細は、グループ化変数を参照してください。

分散パラメーター。スカラー値として返されます。

共分散パラメーターの推定値および関連する統計量。長さが (R + 1) の cell 配列として返されます。R は、モデルで使用するグループ化変数の個数です。stats の最初の R セルにはそれぞれ以下の列をもつデータセット配列が含まれます。

列名説明
Groupグループ化変数名
Name1最初の予測子変数の名前
Name22 番目の予測子変数の名前
Type

Name1Name2 が同じ場合、Typestd (標準偏差) です。

Name1Name2 が異なる場合、Typecorr (相関) です。

Estimate

Name1Name2 が同じ場合、Estimate は予測子 Name1 または Name2 に関連付けされた変量効果の標準偏差です。

Name1Name2 が異なる場合、Estimate は予測子 Name1Name2 に関連付けされた変量効果の間の相関です。

Lower共分散パラメーターの信頼区間の下限
Upper共分散パラメーターの信頼区間の上限

セル R + 1 は分散パラメーターの関連する統計を含みます。

glme の共分散パラメーターの有無には、尤度比検定を使用する compare メソッドを使用した検定が推奨されます。

fitglme および最尤近似メソッド ('Laplace' または 'ApproximateLaplace') のいずれかを使用して GLME モデルをあてはめる場合、covarianceParametersstats の信頼区間を導出するため、一般化線形混合効果モデルの対数尤度へのラプラス近似を基にします。

fitglme および疑似尤度近似メソッド ('MPL' または 'REMPL') のうちいずれかを使用して GLME モデルを近似させる場合、covarianceParametersstats の信頼区間を派生するため、疑似尤度の最後の反復から近似された線形混合効果モデルに基づきます。

すべて展開する

標本データを読み込みます。

load mfr

このシミュレーションされたデータは、世界中で 50 の工場を操業している製造企業から取得しており、各工場が完成品の生産のためにバッチ処理を実行しています。同社は各バッチの欠陥数を減少させるために新たな製造プロセスを開発しました。新しいプロセスの効果を検定する実験のため、参加させる 20 件の工場を無作為に選びました。10 工場では新プロセスを実施しますが、残りの 10 工場では旧プロセスの実行を続けます。各 20 工場で、同社は 5 つのバッチ (合計 100 バッチ) を実行し以下のデータを記録しました。

  • 新しいプロセスがバッチに使用されたかどうかを示すフラグ (newprocess)

  • 各バッチの処理時間。時間単位 (time)

  • バッチの温度。摂氏 (temp)

  • バッチで使用する化学薬品の供給業者 (AB または C) を示すカテゴリカル変数 (supplier)

  • バッチ内の欠陥数 (defects)

またデータに含まれる time_devtemp_dev は、摂氏 20 度で 3 時間の標準プロセスから得られる時間と温度の絶対偏差をそれぞれ表します。

固定効果予測子として newprocesstime_devtemp_dev および supplier を使用して一般化線形混合効果モデルを近似します。工場特有の変動に起因して品質に差がある可能性を考慮するために、factory 別にグループ化された切片の変量効果項を含めます。応答変数 defects はポアソン分布であり、このモデルの適切なリンク関数は対数です。係数の予測にラプラス近似メソッドを使用します。ダミー変数エンコードを 'effects' として指定すると、ダミー変数の係数の合計が 0 になります。

欠陥数はポアソン分布を使用してモデル化できます

defectsijPoisson(μij)

これは一般化線形混合効果モデルに対応します

log(μij)=β0+β1newprocessij+β2time_devij+β3temp_devij+β4supplier_Cij+β5supplier_Bij+bi,

ここで、

  • defectsij は、バッチ j 処理中の工場 i で実行されたバッチで観測された欠陥数です。

  • μij は、バッチ j (j=1,2,...,5) 処理中の工場 i (i=1,2,...,20) に対応する欠陥の平均数です。

  • newprocessijtime_devij および temp_devij は、バッチ j 処理中の工場 i に対応する各変数の測定値です。たとえば newprocessij は、工場 i で実行されたバッチ j 処理中に新プロセスが使用されたかどうかを示します。

  • supplier_Cij および supplier_Bij は効果 (合計はゼロ) の符号化を使用するダミー変数であり、バッチ j 処理中に工場 i で実行されたバッチに対して、それぞれ会社 C または B が加工化学薬品を供給したかどうかを示します。

  • biN(0,σb2) は、工場特有の品質変動に相当する、各工場 i の変量効果の切片です。

glme = fitglme(mfr,'defects ~ 1 + newprocess + time_dev + temp_dev + supplier + (1|factory)','Distribution','Poisson','Link','log','FitMethod','Laplace','DummyVarCoding','effects');

変量効果の予測子について、前の共分散パラメーターの推定を計算して表示します。

[psi,dispersion,stats] = covarianceParameters(glme);
psi{1}
ans = 0.0985

psi{1} は、最初のグループ化変数の前にある共分散行列を推定しています。この例ではグループ化変数が 1 つだけ (factory) なので、psi{1}σb2 の推定値です。

分散パラメーターを表示します。

dispersion
dispersion = 1

予測子に関連付けられた変量効果の推定標準偏差を表示します。stats の最初のセルは factory に関する統計を含み、2 番目のセルは分散パラメーターに関する統計を含みます。

stats{1}
ans = 
    Covariance Type: Isotropic

    Group      Name1                Name2                Type         Estimate
    factory    '(Intercept)'        '(Intercept)'        'std'        0.31381 


    Lower      Upper  
    0.19253    0.51148

予測子に関連付けられた変量効果の推定標準偏差は 0.31381 です。95% 信頼区間は [0.19253 , 0.51148] です。信頼区間が 0 を含まないため、無作為の切片は 5% の有意水準において有意です。