data example1;

input hp path c n @@;

ln=log(n);

cards;

1 1 2 8 1 2 3 24 1 3 2 30 1 4 1 11 1 5 0 6 1 6 0 18 17 0 26 1 8 0 27

2 1 0 3 2 2 3 16 2 3 0 16 2 4 1 8 2 5 1 13 2 6 1 19 27 1 50 2 8 0 28

3 1 0 0 3 2 1 2 3 3 3 10 3 4 0 2 3 5 0 12 3 6 0 4 3 7 011 3 8 0 14

4 1 0 1 4 2 0 1 4 3 1 3 4 4 0 1 4 5 0 73 4 6 1 40 4 71 8 4 8 1 16

5 1 0 1 5 2 0 1 5 3 0 1 5 4 0 1 5 5 0 28 5 6 0 63 5 71 6 5 8 2 16

;

proc genmod data=example1;

class hp  path;

model c=hp path/dist=zinb offset=ln;

zeromodel hp path;

output out=zinb predicted=pred pzero=pzero;

ods output parameterestimates=zinbparms;

ods output modelfit=fit;

run;