그 전 포스팅은 두 모집단의 평균 비교하기 위해 R을 이용했다면, 이번엔 SAS로 하는 방법을 올립니다.

 

데이터는 그 전 포스팅과 똑같습니다. 조금 더 설명해보자면, 데이타에는 총 7개 judge그룹이 있는데요~ 하나는 spock's judge라는 group이고, 나머지 6개는 A, B, C, D, E, F,라는 judge groups 입니다. 이 6개의 groups을 합쳐서, Spocks과 평균값이 차이가 나는지 알아보려고 합니다. 

 

(F까지 있는데 이미지가 좀 잘렸어요;;)

 

 

 

 

두 그룹의 평균값을 비교하는거라서 t-test를 이용하고요~ 물론 glm으로도 알아볼 수 있지만, 그건 다음편에 정리하겠습니다.  

 

결과물 포함한 자세한 내용의 포스팅은

http://statnmath.blogspot.ca/2015/05/case-study-two-sample-t-test-in-sas.html  기 참조하세요~

 

먼저 data를 all이라고 지정했어요. 6개 그룹을 하나로 합쳐야하니까 데이타 정리가 필요하겠죠. if statement 이용해서 두 개의 그룹으로 나눴습니다. 자세한 코드는 이미지 참조하세요~ 

 

Box Plot SAS Code

proc sgplot data=all;

  title 'Percent of women on juries by judge';

  vbox pcwomen / category=judge1;

run;

 

T-test SAS code

proc ttest data=all;

  title "T-test: means of Spock's judge versus all other judges";

  class judge1;

  var pcwomen;

run;

 

 

결과값은 R로 돌린것과 같습니다.  

반응형

+ Recent posts