LinearMixedModel

R2026a

線形混合効果モデル

説明

LinearMixedModel オブジェクトは、固定効果と変量効果をもつ応答変数を表します。このオブジェクトは、モデルの説明、近似係数、共分散パラメーター、計画行列、残差、残差プロットおよび線形混合効果モデルクラスのその他の診断情報で構成されます。predict 関数を使用してモデルの応答を予測でき、random 関数を使用して新しい計画点に乱数データを生成できます。

作成

LinearMixedModel モデルの作成には fitlme または fitlmematrix を使用します。データが table またはデータセット配列に格納されている場合は、fitlme(tbl,formula) を使用して線形混合効果モデルを当てはめることができます。あるいは、式で表現することが困難なモデルの場合は、固定効果と変量効果を定義する行列を作成し、fitlmematrix(X,y,Z,G) を使用してモデルを当てはめることができます。

プロパティ

すべて展開する

係数推定値

`Coefficients` — 固定効果係数の推定値
データセット配列

固定効果の係数推定値とそれに関連する統計。以下のフィールドを含むデータセット配列として格納されます。

`Name`	項の名前。
`Estimate`	係数の推定値。
`SE`	係数の標準誤差。
`tStat`	係数が 0 になるという帰無仮説を検定するための t 統計。
`DF`	t 検定に対する自由度。`DF` の計算方法は `'DFMethod'` 名前と値のペア引数で指定します。`Coefficients` は常に `'Residual'` を `'DFMethod'` メソッドに使用します。
`pValue`	t 検定の p 値
`Lower`	係数の信頼区間の下限。`Coefficients` は常に 95% の信頼水準 (つまり `'alpha'` が 0.05) を使用します。
`Upper`	係数の信頼区間の上限。`Coefficients` は常に 95% の信頼水準 (つまり `'alpha'` が 0.05) を使用します。

'DFMethod' と 'alpha' は、信頼区間の計算中または固定効果と変量効果を含む帰無仮説の検定中に、coefCI メソッドと coefTest メソッドを使用して変更できます。

`CoefficientCovariance` — 推定された固定効果係数の共分散
p 行 p 列の行列

線形混合効果モデルの固定効果係数の推定値の共分散。p 行 p 列の行列として格納されます。ここで p は固定効果係数の数です。

変量効果に関連する共分散パラメーターは、covarianceParameters メソッドを使用して表示できます。

データ型: double

`CoefficientNames` — 固定効果係数の名前
1 行 p 列の文字ベクトルの cell 配列

線形混合効果モデルの固定効果係数の名前。1 行 p 列の文字ベクトルの cell 配列として格納されます。

データ型: cell

`NumCoefficients` — 固定効果係数の数
正の整数値

近似線形混合効果モデルの固定効果係数の数。正の整数値として格納されます。

データ型: double

`NumEstimatedCoefficients` — 推定固定効果係数の数
正の整数値

近似線形混合効果モデルの推定固定効果係数の数。正の整数値として格納されます。

データ型: double

近似法

`FitMethod` — 線形混合効果モデルの当てはめ手法
`ML` | `REML`

線形混合効果モデルを当てはめる方法。以下のいずれかとして格納されます。

当てはめの方法が最尤法の場合は ML
当てはめの方法が制限付き最尤法の場合は REML

データ型: char

入力データ

`Formula` — 固定効果項、変量効果項、グループ化変数の仕様
オブジェクト

線形混合効果モデルを定義する固定効果項、変量効果項およびグループ化変数の仕様。オブジェクトとして格納されます。

式を使用して近似するモデルを指定する方法の詳細については、式を参照してください。

`NumObservations` — 観測値の数
正の整数値

近似に使用される観測の数。正の整数値として格納されます。これは、table またはデータセット配列または計画行列の行の数から、除外された行または NaN のある行を差し引いたものです。

データ型: double

`NumPredictors` — 予測子の数
正の整数値

線形混合効果モデルで予測子として使用される変数の数。正の整数値として格納されます。

データ型: double

`NumVariables` — 変数の総数
正の整数値

応答や予測子を含む変数の総数。正の整数値として格納されます。

標本データが table またはデータセット配列 tbl に格納されている場合、NumVariables は応答変数を含む tbl の変数の総数です。
近似が行列入力に基づいている場合、NumVariables は予測子行列と応答ベクトルの列の総数です。

NumVariables には、予測子または応答として使用されない変数があればそれらが含まれます。

データ型: double

`ObservationInfo` — 観測についての情報
table

近似に使用される観測についての情報。table として格納します。

ObservationInfo には、観測ごとに 1 つの行と以下の 4 つの列があります。

`Weights`	対象の観測の重み付き変数の値。既定値は 1 です。
`Excluded`	名前と値のペア引数 `'Exclude'` を使用して観測を近似から除外した場合は `true`、それ以外の場合は `false`。1 は `true` を表し、0 は `false` を表します。
`Missing`	応答または予測子の値が欠落しているために観測を近似から除外した場合は `true`、それ以外の場合は `false`。欠損値は、数値の場合は `NaN`、cell 配列の場合は空のセル、文字配列の場合は空行、categorical 配列の場合は `<undefined>` 値です。
`Subset`	観測値が近似に使用された場合は `true`、欠落していたか除外されたため観測値に使用されなかった場合は `false`。

データ型: table

`ObservationNames` — 観測の名前
文字ベクトルの cell 配列

当てはめに使用される観測値の名前。文字ベクトルの cell 配列として格納されます。

データが、観測名を含む table またはデータセット配列 tbl に含まれている場合、ObservationNames にはそれらの名前が格納されます。
データが、観測名を含まない行列、table またはデータセットにある場合、ObservationNames は空の cell 配列です。

データ型: cell

`PredictorNames` — 予測子の名前
文字ベクトルの cell 配列

当てはめで予測子として使用する変数の名前。NumPredictors と同じ長さの文字ベクトルの cell 配列として格納されます。

データ型: cell

`ResponseName` — 応答変数の名前
文字ベクトル

当てはめで応答変数として使用される変数の名前。文字ベクトルとして格納されます。

データ型: char

`Variables` — 変数
table

変数。table として格納します。

近似が table またはデータセット配列 tbl に基づいている場合、Variables は tbl と同じです。
近似が行列入力に基づいている場合、Variables は予測子行列のすべての変数と応答変数が格納された table です。

データ型: table

`VariableInfo` — 変数の情報
table

近似に使用される変数についての情報。table として格納します。

VariableInfo には、観測ごとに 1 つの変数と以下の 4 つの列があります。

`Class`	変数のクラス (`'double'`、`'cell'`、`'nominal'` など)。
`Range`	変数の値の範囲。数値変数の場合は、`[min,max]` 形式の 2 要素ベクトルです。セルまたはカテゴリカル変数の場合は、変数の一意の値をすべて含むセルまたは categorical 配列です。
`InModel`	変数が近似モデルの予測子である場合は `true`。変数が近似モデルに含まれていない場合は `false`。
`IsCategorical`	変数の型が、セル、論理、カテゴリカルなどのカテゴリカル予測子として処理される型である場合や、変数が `fit` メソッドの名前と値のペア引数 `'Categorical'` でカテゴリカルとして指定された場合は `true`。変数が連続予測子である場合は `false`。

データ型: table

`VariableNames` — 変数の名前
文字ベクトルの cell 配列

当てはめで使用する変数の名前。文字ベクトルの cell 配列として格納されます。

標本データが table またはデータセット配列 tbl に含まれている場合、VariableNames には tbl の変数の名前が格納されます。
標本データが行列形式である場合、VariableInfo にはモデルの近似中に提供した変数名が格納されます。変数名を指定しなかった場合、VariableInfo には既定の名前が格納されます。

データ型: cell

要約統計量

`DFE` — 残差自由度
正の整数値

残差の自由度。正の整数値として格納されます。DFE = n - p で、n は観測値の数、p は固定効果係数の数です。

これは fixedEffects および randomEffects メソッドでの自由度の計算の 'Residual' メソッドに対応します。

データ型: double

`LogLikelihood` — 最大化された対数尤度または制限付き対数尤度
スカラー値

選択した近似法に応じた近似線形混合効果モデルの最大化された対数尤度または最大化された制限付き対数尤度。スカラー値として格納されます。

データ型: double

`ModelCriterion` — モデル基準
データセット配列

近似線形混合効果モデルを比較するためのモデル基準。以下の列をもつデータセット配列として格納します。

`AIC`	赤池情報量基準
`BIC`	ベイズ情報量基準
`Loglikelihood`	モデルの対数尤度値
`Deviance`	モデルの対数尤度の –2 倍

n がモデルによる近似に使用される観測の数で、p が固定効果係数の数である場合、AIC と BIC の計算に使用される値は以下のようになります。

パラメーターの総数は nc + p + 1 。ここで、nc は残差分散を除く変量効果の共分散のパラメーターの総数
有効な観測の数は以下のとおりです。
- 当てはめの方法が最尤法 (ML) の場合は n
- 当てはめの方法が制限付き最尤法 (REML) の場合は n – p

`MSE` — ML または REML の推定値
正のスカラー値

σ² の計算に使用される近似方法に基づく ML または REML の推定値。正のスカラー値として格納されます。σ² は線形混合効果モデルの観測誤差項の残差分散または分散です。

データ型: double

`Rsquared` — 当てはめたモデルで説明される応答の変動性の比率
構造体

近似モデルで説明される応答の変動性の比率。構造体として格納されます。これは複数の相関係数または決定係数です。Rsquared には 2 つのフィールドがあります。

Ordinary 決定係数値。構造体のスカラー値として格納されます。Rsquared.Ordinary = 1 – SSE./SST

Adjusted

固定効果係数の数に対する自由度調整済み決定係数値。構造体のスカラー値として格納されます。

Rsquared.Adjusted = 1 – (SSE./SST)*(DFT./DFE),

ここで DFE = n – p および DFT = n – 1 です。n は観測値の総数、p は固定効果係数の数です。

データ型: struct

`SSE` — 誤差の二乗和
正のスカラー

二乗和誤差。正のスカラーを指定します。SSE は二乗した条件付き残差に等しい値で、次のようになります。

SSE = sum((y – F).^2),

ここで y は応答ベクトルで、F は線形混合効果モデルの近似条件付き応答です。条件付きモデルには、固定効果と変量効果の両方からの寄与が含まれます。

観測値の重みを使用して学習させたモデルの場合、SSE の計算の二乗和は重み付き二乗和です。

データ型: double

`SSR` — 回帰二乗和
正のスカラー

回帰二乗和。正のスカラーを指定します。SSR は線形混合効果回帰で説明される二乗和で、当てはめた値と応答の平均の間における偏差の二乗和に等しい値です。

SSR = sum((F – mean(y)).^2),

ここで F は線形混合効果モデルの近似条件付き応答で、y は応答ベクトルです。条件付きモデルには、固定効果と変量効果の両方からの寄与が含まれます。

観測値の重みを使用して学習させたモデルの場合、SSR の計算の二乗和は重み付き二乗和です。

データ型: double

`SST` — 二乗の総和
正のスカラー

二乗総和。正のスカラーを指定します。

切片をもつ線形混合効果モデルでは、SST は次のように計算されます。

SST = SSE + SSR,

ここで、SST は二乗総和、SSE は残差平方和、SSR は回帰二乗和です。

切片のない線形混合効果モデルでは、SST は観察された応答値のその平均からの偏差の二乗和として計算され、次のようになります。

SST = sum((y – mean(y)).^2),

ここで、y は応答ベクトルです。

観測値の重みを使用して学習させたモデルの場合、SST の計算の二乗和は重み付き二乗和です。

データ型: double

オブジェクト関数

`anova`	線形混合効果モデルの分散分析
`coefCI`	線形混合効果モデルの係数の信頼区間
`coefTest`	線形混合効果モデルの固定効果と変量効果についての仮説検定
`compare`	線形混合効果モデルの比較
`covarianceParameters`	線形混合効果モデルの共分散パラメーターの抽出
`designMatrix`	固定効果と変量効果の計画行列
`fitted`	線形混合効果モデルからの近似応答
`fixedEffects`	固定効果と関連する統計の推定
`partialDependence`	部分依存の計算
`plotPartialDependence`	部分依存プロット (PDP) および個別条件付き期待値 (ICE) プロットの作成
`plotResiduals`	線形混合効果モデルの残差のプロット
`predict`	線形混合効果モデルの応答予測
`random`	近似線形混合効果モデルからのランダム応答の生成
`randomEffects`	変量効果と関連する統計の推定
`residuals`	近似線形混合効果モデルの残差
`response`	線形混合効果モデルの応答ベクトル

例

すべて折りたたむ

カテゴリカル予測子をもつランダム切片モデル

ライブスクリプトを開く

標本データを読み込み、table の形式に変換します。

load flu
flu = dataset2table(flu)

flu=52×11 table
         Date          NE      MidAtl    ENCentral    WNCentral    SAtl     ESCentral    WSCentral     Mtn      Pac     WtdILI
    ______________    _____    ______    _________    _________    _____    _________    _________    _____    _____    ______

    {'10/9/2005' }     0.97    1.025       1.232        1.286      1.082      1.457          1.1      0.981    0.971    1.182 
    {'10/16/2005'}    1.136     1.06       1.228        1.286      1.146      1.644        1.123      0.976    0.917     1.22 
    {'10/23/2005'}    1.135    1.172       1.278        1.536      1.274      1.556        1.236      1.102    0.895     1.31 
    {'10/30/2005'}     1.52    1.489       1.576        1.794       1.59      2.252        1.612      1.321    1.082    1.343 
    {'11/6/2005' }    1.365    1.394        1.53        1.825       1.62      2.059        1.471      1.453    1.118    1.586 
    {'11/13/2005'}     1.39    1.477       1.506          1.9      1.683      1.813        1.464      1.388    1.204     1.47 
    {'11/20/2005'}    1.212    1.231       1.295        1.495      1.347      1.794        1.303      1.371    1.137    1.611 
    {'11/27/2005'}    1.477    1.546       1.557        1.855      1.678      2.159        1.739      1.628    1.443    1.827 
    {'12/4/2005' }    1.285     1.43       1.482        1.635      1.577      1.903         1.53      1.701    1.516    1.776 
    {'12/11/2005'}    1.354     1.45        1.46        1.794      1.583      1.894        1.831      2.364    2.094    1.941 
    {'12/18/2005'}    1.502    1.622       1.638        1.988      1.947       2.22        2.577       3.89     2.66     2.34 
    {'12/25/2005'}     1.86    1.915       1.955         2.38      2.343      3.027        3.219      4.862    2.595    3.086 
    {'1/1/2006'  }    2.114    2.174       2.065        2.557      2.275      2.498        2.644      3.352    2.181     3.26 
    {'1/8/2006'  }    1.815    1.932       1.822        2.046      1.969      1.805        2.189      2.132    1.717    2.613 
    {'1/15/2006' }    1.541    1.695       1.581        2.008      1.718      1.662        2.156      1.694    1.351    2.247 
    {'1/22/2006' }    1.632    1.758       1.711        2.217      1.866      2.194        2.268      1.826    1.384    2.352 
      ⋮

table flu には、変数 Date と、インフルエンザ推定罹患率 (Google® 検索から推定される 9 地域の値と疾病対策センター (CDC) による全国の推定値) が格納されている 10 個の変数が含まれています。

線形混合効果モデルを当てはめるには、データが適切な形式の table になっていなければなりません。インフルエンザ罹患率を応答として、地域を予測子変数として線形混合効果モデルを当てはめるため、地域に対応する 9 個の列を 1 つの配列にまとめます。新しい table flu2 には、応答変数 FluRate、各推定の元になっている地域を示すノミナル変数 Region、およびグループ化変数 Date が含まれなければなりません。

flu2 = stack(flu,2:10,NewDataVariableName="FluRate",...
    IndexVariableName="Region");
flu2.Date = nominal(flu2.Date);

地域に対する固定効果と、Date で変化するランダム切片で、線形混合効果モデルを当てはめます。

地域はノミナル変数であるため、fitlme は最初の地域 NE を参照として受け取り、他の 8 つの地域を表す 8 つのダミー変数を作成します。たとえば、 $I [M i d A t l]$ は地域 MidAtl を表すダミー変数です。詳細については、ダミー変数を参照してください。

対応するモデルは以下のとおりです。

$\begin{array}{l} y_{i m} = β_{0} + β_{1} I {[M i d A t l]}_{i} + β_{2} I {[E N C e n t r a l]}_{i} + β_{3} I {[W N C e n t r a l]}_{i} + β_{4} I {[S A t l]}_{i} \\ + β_{5} I {[E S C e n t r a l]}_{i} + β_{6} I {[W S C e n t r a l]}_{i} + β_{7} I {[M t n]}_{i} + β_{8} I {[P a c]}_{i} + b_{0 m} + ε_{i m}, m = 1, 2, . . ., 52, \end{array}$

ここで、 $y_{i m}$ はグループ化変数 Date の水準 $m$ に対する観測値 $i$ 、 $β_{j}$ は固定効果係数 ( $j$ = 0、1、...、8)、 $b_{0 m}$ はグループ化変数 Date の水準 $m$ に対する変量効果、 $ε_{i m}$ は観測値 $i$ の観測誤差です。変量効果の事前分布は $b_{0 m} \sim N (0, σ_{b}^{2})$ 、誤差項の分布は $ε_{i m} \sim N (0, σ^{2})$ です。

lme = fitlme(flu2,"FluRate ~ 1 + Region + (1|Date)")

lme = 
Linear mixed-effects model fit by ML

Model information:
    Number of observations             468
    Fixed effects coefficients           9
    Random effects coefficients         52
    Covariance parameters                2

Formula:
    FluRate ~ 1 + Region + (1 | Date)

Model fit statistics:
    AIC       BIC       LogLikelihood    Deviance
    318.71    364.35    -148.36          296.71  

Fixed effects coefficients (95% CIs):
    Name                        Estimate    SE          tStat      DF     pValue        Lower        Upper    
    {'(Intercept)'     }          1.2233    0.096678     12.654    459     1.085e-31       1.0334       1.4133
    {'Region_MidAtl'   }        0.010192    0.052221    0.19518    459       0.84534    -0.092429      0.11281
    {'Region_ENCentral'}        0.051923    0.052221     0.9943    459        0.3206    -0.050698      0.15454
    {'Region_WNCentral'}         0.23687    0.052221     4.5359    459    7.3324e-06      0.13424      0.33949
    {'Region_SAtl'     }        0.075481    0.052221     1.4454    459       0.14902     -0.02714       0.1781
    {'Region_ESCentral'}         0.33917    0.052221      6.495    459    2.1623e-10      0.23655      0.44179
    {'Region_WSCentral'}           0.069    0.052221     1.3213    459       0.18705    -0.033621      0.17162
    {'Region_Mtn'      }        0.046673    0.052221    0.89377    459       0.37191    -0.055948      0.14929
    {'Region_Pac'      }        -0.16013    0.052221    -3.0665    459     0.0022936     -0.26276    -0.057514

Random effects covariance parameters (95% CIs):
Group: Date (52 Levels)
    Name1                  Name2                  Type           Estimate    Lower     Upper  
    {'(Intercept)'}        {'(Intercept)'}        {'std'}        0.6443      0.5297    0.78368

Group: Error
    Name               Estimate    Lower      Upper
    {'Res Std'}        0.26627     0.24878    0.285

$p$ 値 7.3324e-06 および 2.1623e-10 は、地域 WNCentral および ESCentral のインフルエンザ罹患率の固定効果が地域 NE のインフルエンザ罹患率に対して有意に異なることをそれぞれ示しています。

変量効果項の標準偏差 $σ_{b}$ の信頼限界 (0.5297, 0.78368) には、変量効果項が有意であることを示す 0 が含まれていません。compare メソッドを使用して、変量効果の項の有意性をテストすることもできます。

観測値の推定値は、その観測値に対応するグループ化変数水準の固定効果および変量効果の値の合計です。たとえば、2005 年 10 月 9 日の週における地域 WNCentral の罹患率に対する最良線形不偏予測量 (BLUP) の推定値は次のようになります。

$\begin{aligned} {y_{}^{ˆ}}_{W N C e n t r a l, 10 / 9 / 2005} & = {β_{}^{ˆ}}_{0} + {β_{}^{ˆ}}_{3} I [W N C e n t r a l] + {b_{}^{ˆ}}_{10 / 9 / 2005} \\ = 1.2233 + 0.23687 - 0.1718 \\ = 1.28837 . \end{aligned}$

固定効果と変量効果の両方による推定に対する寄与が含まれているので、これは近似された条件付き応答です。この値を以下のように計算することもできます。

beta = fixedEffects(lme);
[~,~,STATS] = randomEffects(lme); % Compute the random-effects statistics (STATS)
STATS.Level = nominal(STATS.Level);
y_hat = beta(1) + beta(4) + STATS.Estimate(STATS.Level=="10/9/2005")

y_hat = 
1.2884

fitted メソッドを使用して、近似値を簡単に表示できます。

F = fitted(lme);
F(flu2.Date == "10/9/2005" & flu2.Region == "WNCentral")

ans = 
1.2884

2005 年 10 月 9 日の週における地域 WNCentral の近似限界応答を計算します。

F = fitted(lme,Conditional=false);
F(flu2.Date == "10/9/2005" & flu2.Region == "WNCentral")

ans = 
1.4602

ランダムな傾きをもつ線形混合効果モデル

ライブスクリプトを開く

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

load carbig

ガロンあたりの走行マイル数 (MPG) の線形混合効果モデルを当てはめます。ここで、加速度、馬力、気筒数は固定効果で、モデル年度によってグループ化される切片と加速度は無相関の変量効果です。このモデルは以下に対応します。

${M P G}_{i m} = β_{0} + β_{1} {A c c}_{i} + β_{2} H P + b_{0 m} + {b_{1}}_{m} {A c c}_{i m} + ε_{i m}, m = 1, 2, 3,$

変量効果項の事前分布は次のようになります。

$b_{m} = (\begin{array}{l} b_{0 m} \\ b_{1 m} \end{array}) \sim N (0, (\begin{array}{cccccccccccccccccccc} σ_{0}^{2} & σ_{0, 1} \\ σ_{0, 1} & σ_{1}^{2} \end{array})),$

ここで、 $m$ はモデル年を表します。

最初に、線形混合効果モデルを当てはめるための計画行列を準備します。

X = [ones(406,1) Acceleration Horsepower];
Z = [ones(406,1) Acceleration];
Model_Year = nominal(Model_Year);
G = Model_Year;

次に、定義した計画行列とグループ化変数で fitlmematrix を使用してモデルを当てはめます。'fminunc' 最適化アルゴリズムを使用します。

lme = fitlmematrix(X,MPG,Z,G,'FixedEffectPredictors',....
{'Intercept','Acceleration','Horsepower'},'RandomEffectPredictors',...
{{'Intercept','Acceleration'}},'RandomEffectGroups',{'Model_Year'},...
'FitMethod','REML')

lme = 
Linear mixed-effects model fit by REML

Model information:
    Number of observations             392
    Fixed effects coefficients           3
    Random effects coefficients         26
    Covariance parameters                4

Formula:
    y ~ Intercept + Acceleration + Horsepower + (Intercept + Acceleration | Model_Year)

Model fit statistics:
    AIC       BIC       LogLikelihood    Deviance
    2202.9    2230.7    -1094.5          2188.9  

Fixed effects coefficients (95% CIs):
    Name                    Estimate    SE           tStat      DF     pValue        Lower       Upper   
    {'Intercept'   }          50.064       2.3176     21.602    389    1.4185e-68      45.507       54.62
    {'Acceleration'}        -0.57897      0.13843    -4.1825    389    3.5654e-05    -0.85112    -0.30681
    {'Horsepower'  }        -0.16958    0.0073242    -23.153    389    3.5289e-75    -0.18398    -0.15518

Random effects covariance parameters (95% CIs):
Group: Model_Year (13 Levels)
    Name1                   Name2                   Type            Estimate    Lower       Upper   
    {'Intercept'   }        {'Intercept'   }        {'std' }           3.72       1.5215      9.0954
    {'Acceleration'}        {'Intercept'   }        {'corr'}        -0.8769     -0.98275    -0.33845
    {'Acceleration'}        {'Acceleration'}        {'std' }         0.3593      0.19418     0.66483

Group: Error
    Name               Estimate    Lower     Upper 
    {'Res Std'}        3.6913      3.4331    3.9688

固定効果係数の表示には、推定値、標準誤差 (SE) および 95% 信頼区間の限界 (Lower および Upper) が含まれます。 $p$ 値 (pValue) は、3 つの固定効果係数がすべて有意であることを示しています。

標準偏差の信頼区間と、切片に対する変量効果と加速度の間の相関には 0 は含まれていません。そのため、これらは有意と考えられます。変量効果を検定するには、compare メソッドを使用します。

推定固定効果係数の共分散行列を表示します。

lme.CoefficientCovariance

ans = 3×3

    5.3711   -0.2809   -0.0126
   -0.2809    0.0192    0.0005
   -0.0126    0.0005    0.0001

対角要素は、固定効果係数の推定値の分散を示しています。たとえば、切片の推定値の分散は 5.3711 です。推定値の標準誤差は分散の平方根です。たとえば、切片の標準誤差は 2.3176 (sqrt(5.3711)) です。

非対角要素は、固定効果係数の推定値間の相関を示しています。たとえば、切片と加速度間の相関は –0.2809 で、加速度と馬力間の相関は 0.0005 です。

モデルの決定係数を表示します。

lme.Rsquared

ans = struct with fields:
    Ordinary: 0.7866
    Adjusted: 0.7855

調整値は、モデルの予測子の数に対する自由度調整済み決定係数値です。

詳細

すべて展開する

式

一般に、モデル仕様の式は 'y ~ terms' という形式の文字ベクトルまたは string スカラーです。線形混合効果モデルでは、この式は 'y ~ fixed + (random1|grouping1) + ... + (randomR|groupingR)' の形式になります。ここで、fixed および random には固定効果および変量効果の項が含まれます。

table tbl に以下のものが格納されていると仮定します。

応答変数 y
連続変数またはグループ化変数である予測子変数 X_j
グループ化変数 g₁、g₂、...、g_R

ここで、X_j および g_r のグループ化変数は、categorical 配列、logical 配列、文字配列、string 配列、または文字ベクトルの cell 配列が可能です。

この場合、'y ~ fixed + (random₁|g₁) + ... + (random_R|g_R)' の形式の式において、項 fixed は固定効果の計画行列 X の仕様に対応し、random₁ はグループ化変数 g₁ に対応する変量効果の計画行列 Z₁ の仕様であり、同様に random_R はグループ化変数 g_R に対応する変量効果の計画行列 Z_R の仕様です。fixed 項および random 項はウィルキンソンの表記法で表現できます。

ウィルキンソンの表記法は、モデルに存在する因子を記述します。この表記法は、モデルに存在する因子に関係するものであり、それらの因子の乗数 (係数) に関係するものではありません。

ウィルキンソンの表記法	標準表記の因子
`1`	定数 (切片) 項
`X^k`、`k` は正の整数	`X`, `X²`, ..., `X^k`
`X1 + X2`	`X1`, `X2`
`X1*X2`	`X1`, `X2`, `X1.*X2 (elementwise multiplication of X1 and X2)`
`X1:X2`	`X1.*X2` のみ
`- X2`	`X2` は含めない
`X1*X2 + X3`	`X1`, `X2`, `X3`, `X1*X2`
`X1 + X2 + X3 + X1:X2`	`X1`, `X2`, `X3`, `X1*X2`
`X1X2X3 - X1:X2:X3`	`X1`, `X2`, `X3`, `X1X2`, `X1X3`, `X2*X3`
`X1*(X2 + X3)`	`X1`, `X2`, `X3`, `X1X2`, `X1X3`

Statistics and Machine Learning Toolbox™ 表記は、-1 を使用して項を明示的に削除しない限り、常に定数項を含みます。線形混合効果モデルの仕様例を次にいくつか挙げます。

次に例を示します。

式	説明
`'y ~ X1 + X2'`	切片の固定効果: `X1` および `X2`。これは、`'y ~ 1 + X1 + X2'` と等価です。
`'y ~ -1 + X1 + X2'`	切片のない `X1` と `X2` の固定効果。`-1` を含めることによって暗黙的な切片の項は抑制されます。
`'y ~ 1 + (1 \| g1)'`	グループ化変数 `g1` の水準ごとの切片の固定効果と切片の変量効果の和。
`'y ~ X1 + (1 \| g1)'`	固定傾きのランダム切片モデル。
`'y ~ X1 + (X1 \| g1)'`	相関があり得るランダムな切片と傾き。これは、`'y ~ 1 + X1 + (1 + X1\|g1)'` と等価です。
`'y ~ X1 + (1 \| g1) + (-1 + X1 \| g1)'`	切片と傾きの独立した変量効果項。
`'y ~ 1 + (1 \| g1) + (1 \| g2) + (1 \| g1:g2)'`	`g1` と `g2` に対する独立したメイン効果のあるランダムな切片モデル + 独立した交互作用効果。

バージョン履歴

R2013b で導入

参考

fitlme | fitlmematrix

LinearMixedModel

説明

作成

プロパティ

係数推定値

Coefficients — 固定効果係数の推定値 データセット配列

CoefficientCovariance — 推定された固定効果係数の共分散 p 行 p 列の行列

CoefficientNames — 固定効果係数の名前 1 行 p 列の文字ベクトルの cell 配列

NumCoefficients — 固定効果係数の数 正の整数値

NumEstimatedCoefficients — 推定固定効果係数の数 正の整数値

近似法

FitMethod — 線形混合効果モデルの当てはめ手法 ML | REML

入力データ

Formula — 固定効果項、変量効果項、グループ化変数の仕様 オブジェクト

NumObservations — 観測値の数 正の整数値

NumPredictors — 予測子の数 正の整数値

NumVariables — 変数の総数 正の整数値

ObservationInfo — 観測についての情報 table

ObservationNames — 観測の名前 文字ベクトルの cell 配列

PredictorNames — 予測子の名前 文字ベクトルの cell 配列

ResponseName — 応答変数の名前 文字ベクトル

Variables — 変数 table

VariableInfo — 変数の情報 table

VariableNames — 変数の名前 文字ベクトルの cell 配列

要約統計量

DFE — 残差自由度 正の整数値

LogLikelihood — 最大化された対数尤度または制限付き対数尤度 スカラー値

ModelCriterion — モデル基準 データセット配列

MSE — ML または REML の推定値 正のスカラー値

Rsquared — 当てはめたモデルで説明される応答の変動性の比率 構造体

SSE — 誤差の二乗和 正のスカラー

SSR — 回帰二乗和 正のスカラー

SST — 二乗の総和 正のスカラー

オブジェクト関数

例

カテゴリカル予測子をもつランダム切片モデル

ランダムな傾きをもつ線形混合効果モデル

詳細

式

バージョン履歴

参考

`Coefficients` — 固定効果係数の推定値
データセット配列

`CoefficientCovariance` — 推定された固定効果係数の共分散
p 行 p 列の行列

`CoefficientNames` — 固定効果係数の名前
1 行 p 列の文字ベクトルの cell 配列

`NumCoefficients` — 固定効果係数の数
正の整数値

`NumEstimatedCoefficients` — 推定固定効果係数の数
正の整数値

`FitMethod` — 線形混合効果モデルの当てはめ手法
`ML` | `REML`

`Formula` — 固定効果項、変量効果項、グループ化変数の仕様
オブジェクト

`NumObservations` — 観測値の数
正の整数値

`NumPredictors` — 予測子の数
正の整数値

`NumVariables` — 変数の総数
正の整数値

`ObservationInfo` — 観測についての情報
table

`ObservationNames` — 観測の名前
文字ベクトルの cell 配列

`PredictorNames` — 予測子の名前
文字ベクトルの cell 配列

`ResponseName` — 応答変数の名前
文字ベクトル

`Variables` — 変数
table

`VariableInfo` — 変数の情報
table

`VariableNames` — 変数の名前
文字ベクトルの cell 配列

`DFE` — 残差自由度
正の整数値

`LogLikelihood` — 最大化された対数尤度または制限付き対数尤度
スカラー値

`ModelCriterion` — モデル基準
データセット配列

`MSE` — ML または REML の推定値
正のスカラー値

`Rsquared` — 当てはめたモデルで説明される応答の変動性の比率
構造体

`SSE` — 誤差の二乗和
正のスカラー

`SSR` — 回帰二乗和
正のスカラー

`SST` — 二乗の総和
正のスカラー