data; infile 'c:\turbo\stat2\e00\00-1.dat' firstobs=2; input PRODUKT$ PROD_NR KAMP MAALING EXPOSURE F_R$ TYPE$; n=125; maaling=125*maaling/100; proc genmod; class f_r type; model MAALING/n = EXPOSURE F_R TYPE / link=logit dist=bin scale=p; output out=fitted p=prob; proc print; run; proc gplot; plot maaling*prob; run;