★ SASでのカイ二乗検定 ★

6774. SASでのカイ二乗検定 ウェンディ 2005/05/24 (火) 21:32
├6786. Re: SASでのカイ二乗検定 TY 2005/05/25 (水) 11:29
└6776. Re: SASでのカイ二乗検定 青木繁伸 2005/05/24 (火) 22:07


6774. SASでのカイ二乗検定 ウェンディ  2005/05/24 (火) 21:32
 1標本2分類のカイ二乗検定をSASで行いたいのですが,プログラムの書き方が分かりません。
 (もし,分かる方がいましたらS-Plusでの操作方法も教えてください。)

     好き 嫌い 計
 実側値  70 30  100
 期待値  50 50 100

 結果はカイ2乗値…16.0になるはずです。一応,私が書いたプログラムです。
 tablesの後に行と列を指定するはずなんですが,どれが行でどれが列なのか分かりません。というよりも,freqでクロス集計をだす必要があるのでしょうか。

 PROC FREQ DATA=KAI;
TABLES (行*列)/CHISQ;
RUN;

ということは分かっています。

 data kai;
input suki kirai kei;
cards;
70 30 100

;
run;
proc freq data=kai;
tables
run;

     [このページのトップへ]


6786. Re: SASでのカイ二乗検定 TY  2005/05/25 (水) 11:29
> TABLES (行*列)/CHISQ;

の部分が分かっていないようです。次のコードを試してみてください。それから,SASマニュアルは使う前に良く読みましょう。
* code1;
data kai;
input s_k $ n;
cards;
kirai 30
suki 70
;
proc freq data=kai;
tables s_k / chisq;
weight n;
run;
* code2;
data kai;
do i = 1 to 70;
s_k = 'suki ';
output;
end;
do i = 71 to 100;
s_k = 'kirai';
output;
end;
run;
proc freq data=kai;
tables s_k / chisq;
run;

     [このページのトップへ]


6776. Re: SASでのカイ二乗検定 青木繁伸  2005/05/24 (火) 22:07
>  (もし,分かる方がいましたらS-Plusでの操作方法も教えてください。)

SAS は分かりませんし,S-Plus も R と違うのかどうかも分かりませんが,R だと
> chisq.test(c(70,30))

Chi-squared test for given probabilities

data: c(70, 30)
X-squared = 16, df = 1, p-value = 6.334e-05

となるのでしょう。S も同じだと思いますが,S-Plus も同じなのか同なのか,わかりません。

     [このページのトップへ]


● 「統計学関連なんでもあり」の過去ログ--- 033 の目次へジャンプ
● 「統計学関連なんでもあり」の目次へジャンプ
● 直前のページへ戻る