大家看下我传的图片,里面是分析在不同年龄和时间,发病的人数和总人数。用poisson回归模型,但是现在我要加入一个队列因素,队列是用时间减去年龄得到。
年龄和时间两个因素的模型程序如下,其中hp代表时间,path代表年龄,c是发病人数,n是总人数
data exapmle1; input hp path counts n@@;
ln=log(n);
cards;
1 1 2 8 1 2 0 3 1 3 0 0 1 4 0 0 1 5 0 1
2 1 3 21 2 2 3 13 2 3 1 1 2 4 0 0 2 5 0 0
...
8 1 0 27 8 2 0 28 8 3 0 14 8 4 1 16 8 5 2 16
;
proc genmod data=example2;
class hp(param=ref ref=first) path(param=ref ref=first);
model counts=hp path/dist=poisson link=log
offset=ln type1 type3;
run;