data example1;
input hp path c n @@;
ln=log(n);
cards;
1 1 23 158 1 2 41 655 1 3 31 1048 1 4 8 308 1 5 8 483 1 6 7 1003 1 7 1 605 1 8 1 158
2 1 14 78 2 2 32 338 2 3 19 462 2 4 4 85 2 5 1 157 2 6 2 262 2 7 7 438 2 8 2 178
3 1 7 23 3 2 13 86 3 3 15 143 3 4 0 35 3 5 0 204 3 6 0 40 3 7 0 47 3 8 1 67
4 1 4 4 4 2 7 14 4 3 4 30 4 4 5 16 4 5 0 239 4 6 1 44 4 7 1 15 4 8 1 30
5 1 1 3 5 2 2 3 5 3 4 8 5 4 1 4 5 5 1 52 5 6 0 68 5 7 1 7 5 8 2 18
;
proc genmod data=example1 order=data;
class hp path;
model c=hp path/dist=poisson
link=log offset=ln type1 type3;
run;
proc genmod data=example1 order=data;
class hp path;
make 'obstats' out=d1;
model c=path/dist=poi link=log offset=ln type1 type3 obstats;
run;
set example1;
if hp=1;
rate=exp(xbeta);
run;
proc print;
run;
data g1;
set d1;
rr=pred/exp(ln);
proc print data=g1;
run;
这是老师给的程序,我加的自己的数,但是总觉得有点运行不出来。我QQ:469691061 方便的话QQ聊下,万分感谢