title1 'Beetle mortality and pesticide dose';
        data beetle;
           input y n dose;
        label
           y    = 'number killed in group'
           n    = 'number in dose group'
           dose = 'insecticide dose' ;

        datalines;
           15  50  1.08
           24  49  1.16
           26  50  1.21
           24  50  1.26
           29  50  1.31
           29  49  1.35
        run;

        proc print;
        run;

        title2 'Fit a linear dose effect to the binomial data';
        proc genmod;
          model y/n=dose / dist=binomial link=identity obstats;
        run;

        title2 'Logistic regression in GENMOD';
        proc genmod;
          model y/n=dose / dist=binomial obstats;
        run;

        title2 'Probit regression in GENMOD';
        proc genmod;
           model y/n=dose / dist=binomial link=probit obstats;
        run;