大家看下我传的图片，里面是分析在不同年龄和时间，发病的人数和总人数。用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;