data example1;
input hp sex path counts n @@;
ln=log(n);
cards;
1 1 1 1 4 1 1 2 13 102 1 1 3 53 1296 1 1 4 42 2902 1 1 5 3 145
2 1 1 1 6 2 1 2 9 75 2 1 3 28 612 2 1 4 21 429 2 1 5 0 1
3 1 1 0 1 3 1 2 7 51 3 1 3 16 165 3 1 4 2 31 3 1 5 0 0
4 1 1 0 1 4 1 2 6 33 4 1 3 9 35 4 1 4 1 2 4 1 5 0 0
5 1 1 1 3 5 1 2 5 18 5 1 3 3 10 5 1 4 0 1 5 1 5 0 0
1 2 1 0 1 1 2 2 0 121 1 2 3 1 934 1 2 4 1 199 1 2 5 0 1
2 2 1 0 2 2 2 2 1 158 2 2 3 0 186 2 2 4 0 24 2 2 5 0 0
3 2 1 0 2 3 2 2 0 352 3 2 3 0 18 3 2 4 0 0 3 2 5 0 0
4 2 1 0 33 4 2 2 1 218 4 2 3 1 4 4 2 4 0 0 4 2 5 0 0
5 2 1 0 19 5 2 2 0 25 5 2 3 0 6 5 2 4 0 0 5 2 5 0 0
;
proc genmod data=example1;
class hp(param=ref ref='3' '4') sex path
model counts=hp sex path/dist=poisson
link=log offset=ln type1 type3;
run;
我在做APC模型(以poisson回归模型为基础),现在对照组问题,我想设置下hp这个因素,如果以第三组为参照的时候ref='3',但是我现在想设置第三组和第四组的参数都为0,改如何操作呢?