options ls=79 ps=200 nocenter;

/* table of means */
data dat;
   infile 'growth.dat' firstobs=2;
   input trt BA TDZ code $ advplt;
   /*
   if advplt=0 then delete;
   */
   if trt=20 then delete;
   nonzero = ( advplt > 0 );
run;
proc freq;
    table BA*TDZ;
    table nonzero*TDZ / chisq cellchi2 nopercent;
    table BA*nonzero*TDZ / chisq cellchi2 nocol norow nopercent;
options ls=79 ps=200 nocenter;
/* See Growth Example of Yandell (1997) PDA */

proc genmod;       
   class BA TDZ;
   model nonzero = BA | TDZ / rorder=data type1 type3 dist = binomial;

proc genmod;       
   class BA TDZ;
   model advplt = BA | TDZ / rorder=data type1 type3 dist = poisson;

proc genmod;       
   class trt;
   model advplt = trt / noint dist = poisson;
   lsmeans trt / diff;