Sas에서 lsmeans 문을 사용하는 방법(예제 포함)


일원 분산 분석은 3개 이상의 독립 그룹의 평균 간에 통계적으로 유의미한 차이가 있는지 여부를 확인하는 데 사용됩니다.

ANOVA 테이블의 전체 p-값이 특정 유의 수준보다 낮으면 그룹 평균 중 하나 이상이 다른 평균과 다르다고 말할 수 있는 충분한 증거가 있는 것입니다.

정확히 어떤 그룹의 평균이 다른지 확인하려면 사후 테스트를 수행해야 합니다.

SAS에서 LSMEANS 문을 사용하여 다양한 사후 테스트를 수행할 수 있습니다.

다음 예에서는 LSMEANS 문을 실제로 사용하는 방법을 보여줍니다.

예: SAS에서 LSMEANS 문을 사용하는 방법

한 연구자가 연구에 참여할 학생 30명을 모집한다고 가정해 보겠습니다. 학생들은 시험 준비를 위해 세 가지 학습 방법 중 하나를 사용하도록 무작위로 배정 됩니다.

각 학생의 시험 결과는 다음과 같습니다.

다음 코드를 사용하여 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 ;

그러면 다음과 같은 분산 분석 테이블이 생성됩니다.

SAS의 일원 분산 분석

이 표에서 우리는 다음을 볼 수 있습니다:

  • 전체 F값: 5.26
  • 해당 p-값: 0.0140

일원 분산 분석에서는 다음과 같은 귀무 가설과 대립 가설을 사용합니다.

  • H 0 : 모든 그룹 평균이 동일합니다.
  • H A : 적어도 하나의 그룹 평균이 다릅니다.   나머지.

ANOVA 테이블의 p-값( 0.0140 )이 α = 0.05보다 작으므로 귀무 가설을 기각합니다.

이는 세 가지 학습 방법의 평균 시험 점수가 동일하지 않음을 나타냅니다.

정확히 어떤 그룹 평균이 다른지 확인하려면 LSMEANS 문 및 ADJUST=TUKEY 옵션과 함께 PROC GLIMMIX 문을 사용하여 Tukey의 사후 테스트를 수행할 수 있습니다.

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

마지막 결과 테이블은 Tukey의 사후 비교 결과를 보여줍니다.

SAS의 LSMEANS 문

Adj P 열을 보면 그룹 평균의 차이에 대해 조정된 p-값을 볼 수 있습니다.

이 열에는 조정된 p-값이 0.05보다 작은 행이 하나만 있음을 알 수 있습니다. 행은 그룹 A와 그룹 C 간의 평균 차이를 비교하는 행입니다.

이는 그룹 A와 그룹 C 사이의 평균 시험 점수에 통계적으로 유의미한 차이가 있음을 말해줍니다.

구체적으로 다음을 볼 수 있습니다.

  • A그룹 학생과 B그룹 학생의 평균 시험 점수 차이는 – 6.375 입니다. (즉, 그룹 A의 학생들은 그룹 C의 학생들보다 평균 시험 점수가 6.375점 낮았습니다)
  • 평균의 차이에 대한 수정된 p-값은 0.0137 입니다.
  • 이 두 그룹 간의 평균 시험 점수의 실제 차이에 대한 수정된 95% 신뢰 구간은 [-11.5219, -1.2281] 입니다.

다른 그룹의 평균 간에는 통계적으로 유의미한 차이가 없습니다.

참고 : 이 예에서는 ADJUST=TUKEY를 사용하여 Tukey 사후 비교를 수행했지만 BON , BUNNET , NELSON , SCHEFFE , SIDAKSMM 을 지정하여 다른 유형의 사후 비교를 수행할 수도 있습니다.

관련 항목: Tukey 대. 본페로니 vs. Scheffe: 어떤 테스트를 사용해야 합니까?

추가 리소스

다음 자습서에서는 분산 분석 모델에 대한 추가 정보를 제공합니다.

ANOVA를 통한 사후 테스트 사용 가이드
SAS에서 일원 분산 분석을 수행하는 방법
SAS에서 양방향 ANOVA를 수행하는 방법

의견을 추가하다

이메일 주소는 공개되지 않습니다. 필수 필드는 *로 표시됩니다