options ls=79 ps=50 nocenter;
/* See Section 7.2 of Venables and Ripley (1999) MASS */
    
/* table of means */
data dat;
   infile 'budworm.dat' firstobs=2;
   input ldose dead sex $;
   total = 20;
   pctdead = dead / total;
proc freq;
   table sex*ldose/ norow nocol nopercent;
   weight dead;
proc genmod;       
   class sex pctdead;
   model dead/total = ldose | sex / type1 type3
       dist = binomial;
   output out = plot p = pdead;

proc plot;
    plot ldose * pdead = '*' ldose * pctdead = sex / overlay;
    
proc genmod;       
   class sex;
   model dead/total = ldose sex / noint
       dist = binomial;