Sas で lsmeans ステートメントを使用する方法 (例付き)

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

ANOVA 表全体のp 値が一定の有意性レベルを下回っている場合、グループ平均の少なくとも 1 つが他の平均とは異なると言える十分な証拠があります。

どのグループの平均が異なるかを正確に調べるには、 事後テストを実行する必要があります。

SAS でLSMEANSステートメントを使用すると、さまざまな事後テストを実行できます。

次の例は、 LSMEANSステートメントを実際に使用する方法を示しています。

例: SAS で LSMEANS ステートメントを使用する方法

研究者が研究に参加する 30 人の学生を募集するとします。学生は、試験の準備のために 3 つの学習方法のいずれかを使用するようにランダムに割り当てられます。

各生徒の試験結果は以下のとおりです。

次のコードを使用して、SAS でこのデータセットを作成できます。

 /*create dataset*/
data my_data;
    input Method $Score;
    datalines ;
At 78
At 81
At 82
At 82
At 85
At 88
At 88
At 90
B 81
B 83
B 83
B85
B 86
B 88
B90
B91
C 84
C 88
C 88
C 89
C 90
C 93
C 95
C 98
;
run ;

次に、 proc ANOVA を使用して一元配置分散分析を実行します。

 /*perform one-way ANOVA*/
proc ANOVA data =my_data;
classMethod ;
modelScore = Method;
run ;

これにより、次の ANOVA テーブルが生成されます。

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

• 全体のF値： 5.26
• 対応する p 値: 0.0140

一元配置分散分析では次の帰無仮説と対立仮説が使用されることを思い出してください。

• H ₀ :すべてのグループ平均が等しい。
• H _A :少なくとも 1 つのグループの平均が異なります 休む。

ANOVA テーブルの p 値 ( 0.0140 ) は α = 0.05 より小さいため、帰無仮説を棄却します。

これは、試験の平均点が 3 つの勉強方法で同じではないことを示しています。

どのグループ平均が異なるかを正確に判断するには、 PROC GLIMMIXステートメントとLSMEANSステートメントおよびADJUST=TUKEYオプションを使用して、Tukey の事後テストを実行できます。

 /*perform Tukey post-hoc comparisons*/
proc glimmix data =my_data;
    classMethod ;
    modelScore = Method;
    lsmeans Method / adjust =tukey alpha = .05 ;
run ;

最後の結果表は、Tukey の事後比較の結果を示しています。

Adj P列を見ると、グループ平均の差について調整された p 値が表示されます。

この列には、調整された p 値が 0.05 未満の行が 1 行だけあることがわかります。これは、グループ A とグループ C の平均差を比較する行です。

これは、グループ A とグループ C の間に試験の平均点に統計的に有意な差があることがわかります。

具体的には次のことがわかります。

• グループ A の学生とグループ B の学生の試験の平均点の差は – 6.375でした。 (つまり、グループ A の学生は、グループ C の学生よりも試験の平均点が 6.375 ポイント低かった)
• 平均の差の調整された p 値は0.0137です。
• これら 2 つのグループ間の平均試験得点の真の差の調整済み 95% 信頼区間は[-11.5219, -1.2281]です。

他のグループの平均値の間に統計的に有意な差はありません。

注: この例では、 ADJUST=TUKEYを使用して Tukey のポストホック比較を実行しましたが、 BON 、 BUNNET 、 NELSON 、 SCHEFFE 、 SIDAK 、およびSMMを指定して、他のタイプのポストホック比較を実行することもできます。

関連: テューキー vs.ボンフェローニ vs.シェッフェ: どのテストを使用する必要がありますか?

追加リソース

次のチュートリアルでは、ANOVA モデルに関する追加情報を提供します。

ANOVA で事後テストを使用するためのガイド
SAS で一元配置分散分析を実行する方法
SAS で二元配置分散分析を実行する方法