なんかの測定値の平均の95%信頼区間を求めたいときは往々にしてあると思います.
あちこちに載っていますが,そんな時はproc meansで出せますよ,という話です.
proc univariateでも出せますが,結果をデータセット化するのはproc meansのほうが簡単だと思います.
meansでlclm(信頼区間の下限),uclm(信頼区間の上限)を指定するだけです.簡単ですね.
/*----- test data -----*/
data aa ;
a = 1 ; b = 1 ; output ;
a = 2 ; b = 1 ; output ;
a = 3 ; b = 1 ; output ;
a = 4 ; b = 1 ; output ;
a = 5 ; b = 1 ; output ;
a = 6 ; b = 1 ; output ;
a = 9 ; b = 2 ; output ;
a = 8 ; b = 2 ; output ;
a = 7 ; b = 2 ; output ;
a = 5 ; b = 2 ; output ;
a = 6 ; b = 2 ; output ;
run ;
proc means data = AA noprint nway ;
var A ;
class B ;
output out = HOGE n = lclm = uclm = / autoname ;
run ;
0 件のコメント:
コメントを投稿