公開日: 04/01/2021

分割実験の例

変量効果とは、「その水準が、母集団から無作為抽出された標本である」と見なされる効果です。モデルとしては、変量効果は平均が0の正規分布に従うと仮定されます。変量効果の分散は、分散成分と呼ばれることもあります。

この定義から考えると、どんな線形回帰モデルにも少なくとも1つの変量効果があると言えます。それは残差誤差(residual error)と呼ばれる効果です。この残差誤差の項は、通常、母集団において平均が0で正の分散(s2)を持つ正規分布に従うと仮定されます。

典型的な変量効果モデルとして、分割実験モデルが挙げられます。この分割実験モデルは、誤差が複合対称の反復測定モデルと等価です。表4.2は、分割実験モデルの効果と、それに等価な反復測定モデルの効果をまとめたものです。これらのモデルでは、実験は2層に分かれています。まず、処理(処置)のうちのいくつかが、1次単位または被験者(個体)に施されます。次に、それらが2次単位に分割され(反復測定実験では、被験者ごとに複数回、測定が行われ)、2次単位に別の処理が施されます。このモデルでは、1次単位(または被験者)と、2次単位の効果が変量効果です。通常、2次単位の効果はモデル効果としては含めず、残差誤差として扱います。

表4.2 分割実験モデルに見られる効果の種類

分割実験モデル

効果の種類

反復測定モデル

1次単位への処理

固定効果

被験者(個体)に対する処置

1次単位のID

変量効果

被験者のID

2次単位への処理

固定効果

被験者内の処置

2次単位のID

変量効果

測定時点(反復測定のID)

これらのモデルを、多層構造のモデル(layered model)として扱う計算方法があります。そこでは、モデルは次のような2つの異なる実験で構成されているとみなされます。

1. 1次単位への処理。区画または被験者を実験単位として、誤差が計算されます。

2. 2次単位への処理。個々の測定値を実験単位として、誤差が計算されます。この誤差は、残差誤差に相当します。

従来の古典的方法では、1次単位への処理(第1層における効果)を検定する際、次の方法のいずれかを用いていました。

1次単位(区画、被験者)ごとに平均を計算し、それをデータとして、1次単位に対する処理の一元配置分散分析を実行する。

1次単位に対する処理の平均平方を、1次単位IDの平均平方で割ったものを、F値として計算する。

まず、2次単位(時点)ごとの測定値を複数の列で保持している形式でデータを準備する。そして、そのデータにMANOVAモデル(多変量分散分析モデル)をあてはめ、単変量の統計量を計算する。

これらの計算方法は、モデルが単純で、欠測値がない釣合い型データ(バランスが取れたデータ)でないと適用できません。また、変量効果を含むモデルには、分割実験モデルや反復測定モデルより複雑なモデルも存在します。固定効果と変量効果の両方を含むモデルは、混合モデルと呼ばれています。

釣合い型分割実験は、複数の層から構成される計画のなかで最もよく使われています。なお、反復測定モデルは、分割実験モデルと見なすこともできます。分割実験では、ある効果を1次単位に割り付けた後、その1次単位を分割し(反復測定の場合は、時点によって分割し)、分割された2次単位に他の効果を割り付けます。

サンプルデータのフォルダにある「Animals.jmp」データを見てみましょう(これは架空のデータです)。このデータは、キツネとコヨーテの餌の捕獲習性について、季節ごとの差を調べたものです。3頭ずつのキツネとコヨーテに印をつけ、それを1年間にわたって定期的に(季節ごとに)観察し、穴からどれぐらい離れて徘徊したかを平均距離(単位はマイルで、四捨五入した数値)として記録しました。モデルは次の要素で定義されます。

「距離(マイル)」。連続尺度の応答変数。

「種別」。固定効果。値は「キツネ」または「コヨーテ」

「季節」。固定効果。値は「秋」「冬」「春」「夏」

「個体」。各個体を識別する番号。名義尺度。キツネとコヨーテで重複する番号1、2、3が使われています。

モデルには2つの層があります。

1. 上の層は個体間変動(被験者間変動)を扱う層です。キツネとコヨーテを水準とする因子(「種別」効果)について、個体間変動を考慮して検定します。

2. 下の層は、個体内変動(被験者内変動)を扱う層です。4つの季節を水準とする因子(「季節」効果)について、個体内変動を考慮して検定します。なお、個体内変動は、残差誤差として表されます。

「季節」効果に対しては、残差誤差をそのままF統計量の分母として使用することができます。一方、個体間変動は、残差誤差ではなく、「種別」から枝分かれしている「個体」の効果(個体種別])によって把握しなければなりません。個体間の効果である「種別」F統計量を計算するときは、残差誤差ではなく、この枝分かれ効果を分母にします。

注: JMP Proのユーザは、手法として[混合モデル]を使っても、以下と同じ分析を行えます。

このデータの分割実験モデルを指定するには、次の手順に従います。

1. [ヘルプ]>[サンプルデータライブラリ]を選択し、「Animals.jmp」を開きます。

2. [分析]>[モデルのあてはめ]を選択します。

3. 「距離(マイル)」を選択し、[Y]をクリックします。

4. 「種別」「個体」を選択し、[追加]をクリックします。

5. 「列の選択」リストで、「種別」を選択します。

6. 「モデル効果の構成」リストで「個体」を選択します。

7. [枝分かれ]ボタンをクリックします。

「種別」から枝分かれしている「個体」の効果(個体[種別])がモデルに追加されます。

8. 枝分かれ効果の「個体[種別]」を選択します。

9. [属性]>[変量効果]を選択します。

これで、枝分かれ効果が「個体[種別]&変量効果」のように表示され、「種別」効果の誤差項として指定されたことになります。

10. 「列の選択」リストで「季節」を選択し、[追加]をクリックします。

[属性]ポップアップメニューを使って効果を変量効果として指定すると、ダイアログボックスの右上の方に、「方法」オプション([REML][EMS])が表示されます。デフォルトでは[REML]が選択されています。図4.19に指定が完了した起動ウィンドウを示します。

図4.19 「モデルのあてはめ」ダイアログボックス 

Image shown here

11. [実行]をクリックします。

図4.20のようなレポートが作成されます。固定効果である「種別」「季節」は、両方とも有意です。「REML法による分散成分推定値」レポートには、「個体」と残差誤差に対して、分散の推定値が表示されています。

図4.20 REML分析のレポート(一部) 

Image shown here

より詳細な情報が必要な場合や、質問があるときは、JMPユーザーコミュニティで答えを見つけましょう (community.jmp.com).