data;
infile '\turbo\stat2\e00\00-1.dat' firstobs=2;
input PRODUKT$ PROD_NR KAMP MAALING EXPOSURE F_R$  TYPE$;
w=125*(MAALING/100)*(1-MAALING/100);
maaling=log(maaling)-log(100-maaling);


proc glm;
   class f_r type;
   weight w;
   model  MAALING = EXPOSURE F_R  TYPE / SOLUTION;
   output out=fitted predicted=fit;

data;
set fitted;
prob=exp(fit)/(1+exp(fit));
maaling=100*exp(maaling)/(1+exp(maaling));


proc print;


proc gplot;
plot maaling*prob;
run;