読者です 読者をやめる 読者になる 読者になる

要約データを用いたt検定

SASの便利機能。
各群の平均、標準偏差、サンプルサイズをデータセットに入れると検定結果が出る。

/* ダミーデータを作る */
data ttest;
  do group = 1 to 2;
  do i = 1 to 10;
    Value = group * 2 + rand('normal');
    output;
  end; end;
/* 要約データを用いたt検定 */
proc summary data = ttest;
  var Value;
  by group;
  output out = summary;
proc ttest data = summary;
  class group;
  var Value;
run;

/* 手打ちする場合 */
data summary2;
input group _STAT_ $ Value;
cards;
1 N 10
1 MEAN 1.5945
1 STD 0.7567
2 N 10
2 MEAN 3.6526
2 STD 1.5461
;
proc ttest data = summary2;
  class group;
  var Value;
run;