catmod1.sas

SAS for Windows统计分析系统教程,里面有源码

options nodate nonumber;
data catmod1;
  input hospital treat$ y count@@;
cards;
1 a 0 23  1 a 0.5  7  1 a 1 2
1 b 0 23  1 b 0.5 10  1 b 1 5
1 c 0 20  1 c 0.5 13  1 c 1 5
1 d 0 24  1 d 0.5 10  1 d 1 6
2 a 0 18  2 a 0.5  6  2 a 1 1
2 b 0 18  2 b 0.5  6  2 b 1 2
2 c 0 13  2 c 0.5 13  2 c 1 2
2 d 0  9  2 d 0.5 15  2 d 1 2
3 a 0  8  3 a 0.5  6  3 a 1 3
3 b 0 12  3 b 0.5  4  3 b 1 4
3 c 0 11  3 c 0.5  6  3 c 1 2
3 d 0  7  3 d 0.5  7  3 d 1 4
4 a 0 12  4 a 0.5  9  4 a 1 1
4 b 0 15  4 b 0.5  3  4 b 1 2
4 c 0 14  4 c 0.5  8  4 c 1 3
4 d 0 13  4 d 0.5  6  4 d 1 4
5 a 0  8  5 a 0.5  7  5 a 1 6
5 b 0 11  5 b 0.5  8  5 b 1 4
5 c 0  2  5 c 0.5  7  5 c 1 5
5 d 0  1  5 d 0.5  3  5 d 1 4
;
proc catmod order=data;
  model y=treat hospital/freq oneway;
  response 0 0.5 1;
  weight count;
run;