data e;
input year month$ n rem@@;
ln=log(n);
cards;
1983  jan  0    0  1983  feb  10   0  1983  mar  26   0
1983  apr  192  0  1983  may  238  0  1983  jun  374  1
1983  jul  356  0  1983  aug  358  4  1983  sep  594  0
1983  oct  786  1  1983  nov  885  2  1983  dec  981  2
1984  jan  1044  0  1984  feb  1266  1  1984  mar  1533  1
1984  apr  1510  2  1984  may  1818  2  1984  jun  2210  7
1984  jul  2003  2  1984  aug  2489  4  1984  sep  2841  2
1984  oct  3236  3  1984  nov  3104  5  1984  dec  2919  3
1985  jan  2849  3  1985  feb  3119  2  1985  mar  3596  7
1985  apr  4003  5  1985  may  4067  5  1985  jun  3690  2
1985  jul  3509  2  1985  aug  3653  9  1985  sep  4186  3
1985  oct  4562  2  1985  nov  4277  2  1985  dec  3838  2

1986  jan  4253  3  1986  feb  4242  5  1986  mar  5119  5
1986  apr  5281  7  1986  may  5163  4  1986  jun  4977  2
1986  jul  4663  4  1986  aug  5465  4  1986  sep  5314  5
1986  oct  5485  4  1986  nov  5688  6  1986  dec  5403  0
1987  jan  5749  2  1987  feb  5682  4  1987  mar  6460  3
1987  apr  6681  3  1987  may  7215  3  1987  jun  6650  8
1987  jul  7094  2  1987  aug  6600  6  1987  sep  7649  3
1987  oct  7615  9  1987  nov  6733  4  1987  dec  6540  10
1988  jan  6680  4  1988  feb  7646  6  1988  mar  8139  4
1988  apr  7829  4  1988  may  8220  3  1988  jun  7671  5
1988  jul  7120  3  1988  aug  7293  4  1988  sep  8045  5
1988  oct  8567  3  1988  nov  7682  6  1988  dec  7048  3
1989  jan  7369  2  1989  feb  7270  6  1989  mar  8124  1
1989  apr  7636  5  1989  may  7512  5  1989  jun  7049  4
1989  jul  7286  2  1989  aug  7624  2  1989  sep  7623  2
1989  oct  7970  5  1989  nov  7569  1  1989  dec  7156  10

1990  jan  7404  3  1990  feb  7447  8  1990  mar  7951  12
1990  apr  8056  7  1990  may  7742  3  1990  jun  7109  2
1990  jul  7229  4  1990  aug  7249  3  1990  sep  7366  0
1990  oct  7955  6  1990  nov  7044  6  1990  dec  3929  3
;
proc genmod data=e order=data;
class month year;
model rem=month year/dist=poi link=log offset=ln type1 type3;
run;
proc genmod data=e order=data;
class month year;
make 'obstats' out=d;
model rem=year/dist=poi link=log offset=ln type1 type3 obstats;
run;

data rates;
keep year rate month rem ln;
set d;
set e;
if month='dec';rate=exp(xbeta);
run;
proc print;
run;
data g; set d;
rr=pred/exp(ln);
proc print data=g;
run;


proc genmod data=e order=data;
class month year;
model rem/n=month year/dist=poi link=log  type1 type3;
run;