请高手帮忙~~ 我从事卫生统计,要学习的东西很多~~
SAS help 中 catmod process 例子28.8 中:
title 'Growth Curve Analysis';
data growth2;
input Diagnosis $ Treatment $ week1 $ week2 $ week4 $ count @@;
datalines;
mild std n n n 16 severe std n n n 2
mild std n n a 13 severe std n n a 2
mild new n n n 31 severe new n n n 7
mild new n n a 0 severe new n n a 2
;
\\Diagnosis (mild or severe)
Treatment (std or new).
The times used in the design matrix (0, 1, 2) correspond to the logarithms (base 2) of the actual times (1, 2, 4). \\
proc catmod data=growth2 order=data;
title2 'Reduced Logistic Model';
weight count;
population Diagnosis Treatment;
response logit;
model week1*week2*week4=(1 0 0 0, /* mild, std */
1 0 1 0,
1 0 2 0,
1 0 0 0, /* mild, new */
1 0 0 1,
1 0 0 2,
0 1 0 0, /* severe, std */
0 1 1 0,
0 1 2 0,
0 1 0 0, /* severe, new */
0 1 0 1,
0 1 0 2)
(1='Mild diagnosis, week 1',
2='Severe diagnosis, week 1',
3='Time effect for std trt',
4='Time effect for new trt')
/ freq design;
contrast 'Diagnosis effect, week 1' all_parms 1 -1 0 0;
contrast 'Equal time effects' all_parms 0 0 1 -1;
quit;



雷达卡



京公网安备 11010802022788号







