Stata で二元配置分散分析を実行する方法

二元配置分散分析は、 2 つの因子に分割された 3 つ以上の独立したグループの平均間に統計的に有意な差があるかどうかを判断するために使用されます。

二元配置分散分析の目的は、応答変数に対する 2 つの因子の影響を判断し、応答変数に対する 2 つの因子間に交互作用があるかどうかを判断することです。

このチュートリアルでは、Stata で二元配置分散分析を実行する方法について説明します。

例: Stata での二元配置分散分析

この例では、 systolicという組み込みの Stata データセットを使用して、二元配置分散分析を実行します。このデータセットには、58 人の異なる個人に関する次の 3 つの変数が含まれています。

• 使用した薬剤
• 患者の病気
• 最高血圧の変化

次の手順を使用して二元配置分散分析を実行し、使用される薬剤の種類と患者の疾患の種類が収縮期血圧の変化に重大な影響を与えるかどうかを判断します。

ステップ 1: データをロードします。

まず、コマンド ボックスにwebuse systolicと入力し、Enter をクリックしてデータをロードします。

ステップ2: 生データを表示します。

二元配置分散分析を実行する前に、まず生データを見てみましょう。上部のメニュー バーから、 [データ] > [データ エディター] > [データ エディター (参照)]に移動します。これにより、58 人の患者全員の実際のデータが表示されます。

ステップ 3: 二元配置分散分析を実行します。

上部のメニュー バーから、 [統計] > [線形モデルおよび関連モデル] > [ANOVA/MANOVA] > [分散および共分散の分析]に移動します。

[従属変数] には、 systolicを選択します。次に、 [モデル]の下のドロップダウン矢印の横にある 3 つの点をクリックします。

新しい画面が表示されます。 [因子] 変数は選択したままにしておきます。二元配置分散分析を実行しているため、[仕様] で [二元配置完全階乗]を選択します。変数 1 では、薬物を選択し、ベースとしてデフォルトを使用します。変数 2 では、病気を選択し、ベースとしてデフォルトを使用します。次に、下部にある「変数リストに追加」をクリックします。次に、 「OK」をクリックします。

元のウィンドウが表示され、 Model の下に Drug##disease が表示されます。ここではもう何もすることはありません。 「OK」をクリックするだけです。

二元配置分散分析の結果が自動的に表示されます。

結果から次のことがわかります。

• 薬物と病気の間に有意な相互作用はない (p 値 = 0.3958)
• この病気は収縮期血圧に有意な影響を与えません (p 値 = 0.1637)
• この薬は収縮期血圧に対して統計的に有意な効果があります (0.0001)

ステップ 4: 結果を報告します。

最後に、二元配置分散分析の結果を報告します。これを行う方法の例を次に示します。

収縮期血圧に対する薬と病気の影響を調べるために、58 人の個人に対して二元配置分散分析が実行されました。

収縮期血圧に対する薬物と疾患の影響の間には有意な相互作用はありませんでした (p = 0.3958)。収縮期血圧に対する疾患の有意な影響はありませんでした (p = 0.1637)。この薬は収縮期血圧に重大な影響を及ぼしました (0.0001)。