/*  Fit the DYME data using PROC GENMOD to get starting values
    for the NLINMIX macro */

options ps=59 ls=80 nodate; run;

data dyme; 
  infile "dymedat.dat";
  input mouse dose birthwt;
run;

/*  We can trick proc genmod into fitting the model we want even though
    this is not a "generalized linear model" situation.  We want to 
    fit an exponential model, which is equivalent to using the
    "link=log" option.  The power model with theta = 1 for variance
    can be obtained by specifying the gamma distribution */

proc genmod data=dyme;
  class mouse;
  model birthwt = dose / link=log dist=gamma;
  repeated subject=mouse / type=cs corrw modelse v6corr;
run;