options ls=79 ps=200 nocenter;
/* See Section 7.3 of Venables and Ripley (1999) MASS */

/* table of means */
data dat;
   infile 'brandx.dat' firstobs=2;
   input Brand $ Temp $ Muser $ Soft $ Fr;
proc freq;
   table Brand*Temp*Muser*Soft/ norow nocol nopercent;
   weight Fr;
proc genmod;       
   class Brand Temp Muser Soft;
   model Fr = Brand Temp | Muser | Soft / rorder=data type1 type3
       dist = poisson;