「統計学関連なんでもあり」の過去ログ---033

★ 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 も同じなのか同なのか，わかりません。

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