%macro ml(data=,bin=,norm=,x=,id=,print=yes,method=ISAMP,seed=,qpoints=200,
tech=NRRIDG,it_his=no,maxit=);

%let p=0;
%do %while(%scan(&x,&p+1)^=); %let p=%eval(&p+1); %end;

data dataz(keep=&id &bin &norm &x r_bin r_norm zzz);
  set &data;
  zzz=1;
    r_bin=1;
      if (&bin=.) then do;
	    r_bin=0; &bin=0;
	  end;
	  r_norm=1;
	      if (&norm=.) then do;
		    r_norm=0; &norm=0;
	  end;
  %do k=1 %to &p;
     %let covk=%scan(&x,&k);
        if (&covk=.) then do;
	       r_bin=0; r_norm=0; &covk=0;
		end;
  %end;

run;


/* CC estimate of beta */


ods listing close;

 proc genmod data=dataz DESCENDING;
   model &bin =&x / dist=bin link=logit ;
 ods output parameterestimates=glim_est;
 where r_bin=1;
 run;


data glim_est(keep=estimate);
  set glim_est;
    if parameter ne 'Scale';
run;

data names;
length Parameter $ 12 ;
   %do m=1 %to &p+1 ;
    Parameter = "b&m";
	output;
  %end;
run;

data glim_est;
 merge names glim_est;
 run;
 
  proc genmod data=dataz  ;
    model &norm =&x  ;
  ods output parameterestimates=glim_est2;
  where r_norm=1;
  run;
 
 
 data glim_est2(keep=estimate);
   set glim_est2;
  run;
 
 data names;
 length Parameter $ 12 ;
    %do m=1 %to &p+1 ;
     Parameter = "a&m";
 	output;
   %end;
      Parameter = "sig_nor";
 	output;  
 run;
 
 data glim_est2;
  merge names glim_est2;
 run;

data add;
length Parameter $ 12 ;
 Parameter = "std_re_bri";
 estimate =  1;
 output;
  Parameter = "std_re_nor";
  estimate =  1;
 output;
  Parameter = "rho_bri_nor";
  estimate =  .01;
 output;
run;

data glim_est;
 set glim_est glim_est2 add;
 run;

proc datasets;
 delete names glim_est2 add ;
 run;

proc sort data=dataz;
 by &id;
 run;



ods listing  close;

proc nlmixed    data=dataz tech=&tech method=&method NOAD
%if &seed^= %then %do;
 seed=&seed 
%end;
%if &maxit^= %then %do;
 MAXIT=&maxit 
%end;
 /* fd */ qpoints=&qpoints ;
 parms / data= glim_est ;

pi = constant('pi');

/* binary */
uni = probnorm(z1);
phi = 1.0/sqrt(1 + 3/pi/pi*std_re_bri*std_re_bri);
B_bin = 1/phi*log(sin(pi*uni*phi)/sin( phi*pi* (1-uni) )   );

 xb_bin= B_bin + b1
  %do j=1 %to &p ;
      %let ji=%eval(1+&j);
     %let covj=%scan(&x,&j);
        + b&ji * &covj
  %end;
  ;

  p  = exp(xb_bin)/(1 + exp(xb_bin) ); 
 
   llik1 =  r_bin*( &bin * log( p ) + (1 - &bin  ) * log(1 - p ) ) ;

/* normal */

xb_nor=std_re_nor*z2 + a1
  %do j=1 %to &p ;
      %let ji=%eval(1+&j);
     %let covj=%scan(&x,&j);
        + a&ji * &covj
  %end;
  ;

llik2 = r_norm*( (-0.5)*log(2*pi*sig_nor*sig_nor ) - 
        (( &norm - xb_nor)**2)/(2*sig_nor*sig_nor)
) ;

llik = llik1 + llik2 ;

   model zzz ~ general(llik);

random z1 z2 ~ normal( [0,0], [1, rho_bri_nor, 1] )  subject= &id ;


ESTIMATE 'b1' phi*b1  ;
  %do j=1 %to &p ;
      %let ji=%eval(1+&j);
ESTIMATE "b&ji"  b&ji * phi ;
  %end;
ESTIMATE 'phi' phi ;

 ods output  ParameterEstimates=pars_j(keep= Estimate StandardError tValue
 Probt) ;
  ods output  IterHistory=it_his(drop= Rho Flag1) ;
ods output  AdditionalEstimates =marg_est(drop=Label df Alpha Lower Upper);
 

run;

ods listing;

%if &it_his=yes %then %do;
title 'Iteration history for ML';
    proc print noobs data=it_his ;
	run;
%end;

data names;
length Outcome $ 15 Parameter $ 20 ;

   Outcome="&bin";
   Parameter="Intercept";
    output;
  %do k=1 %to &p;
     %let covk=%scan(&x,&k);
    Parameter="&covk";
    output;
  %end;
   Outcome="&norm";
   Parameter="Intercept";
    output;
  %do k=1 %to &p;
     %let covk=%scan(&x,&k);
    Parameter="&covk";
    output;
  %end;  
    Parameter = "sigma norm";
 	output; 
   Outcome="&bin";
 Parameter = "sigma re bridge";
   output;
      Outcome="&norm";
  Parameter = "sigma re nor";
    output;
	      Outcome=" ";
  Parameter = "rho bri nor";
 output;
run;

data pars_j;
 merge names pars_j;
 run;

proc datasets;
 delete names;
 run;
 
title "ML estimate,  conditional model for binary";

proc print noobs;
run;

title '  ';

data names;
length Outcome $ 15 Parameter $ 20 ;

   Outcome="&bin";
   Parameter="Intercept";
    output;
  %do k=1 %to &p;
     %let covk=%scan(&x,&k);
    Parameter="&covk";
    output;
  %end;
   Parameter="phi";
    output;
run;

data marg_est;
 merge names marg_est;
 run;

title "ML estimate,  marginal model for binary";

proc print noobs;
run;

proc datasets;
 delete DATAZ GLIM_EST pars_j names it_his marg_est;
run;

%mend;


data EG_data;

input Litter_ID  Dose  Fetal_Wt  Fetal_Mal ;
dose=dose/1000;
cards;
          60      0      0.903      0 
          60      0      0.828      0 
          60      0      0.953      0 
          60      0      0.954      0 
          60      0      1.070      0 
          60      0      1.065      0 
          60      0      0.959      0 
          61      0      1.025      0 
          61      0      1.006      0 
          61      0      1.069      0 
          61      0      0.825      0 
          61      0      1.013      0 
          61      0      0.994      0 
          61      0      1.059      0 
          61      0      1.041      0 
          61      0      1.020      0 
          61      0      1.041      0 
          61      0      0.911      0 
          61      0      0.988      0 
          61      0      0.947      0 
          61      0      0.976      0 
          68      0      1.058      0 
          68      0      0.950      0 
          68      0      0.860      0 
          68      0      1.015      0 
          68      0      0.907      0 
          68      0      0.856      0 
          68      0      1.075      0 
          68      0      0.958      0 
          68      0      1.082      0 
          68      0      0.916      0 
          68      0      1.054      0 
          68      0      0.921      0 
          68      0      0.874      0 
          69      0      0.999      0 
          69      0      0.950      0 
          69      0      1.065      0 
          69      0      0.944      0 
          69      0      0.949      0 
          69      0      1.044      0 
          69      0      1.025      0 
          69      0      0.987      0 
          69      0      0.879      0 
          69      0      0.908      0 
          76      0      1.071      1 
          76      0      1.056      0 
          76      0      1.084      0 
          76      0      0.851      0 
          76      0      0.925      0 
          76      0      0.832      0 
          76      0      0.979      0 
          76      0      1.096      0 
          76      0      0.954      0 
          76      0      0.875      0 
          76      0      1.057      0 
          76      0      1.245      0 
          76      0      1.018      0 
          76      0      1.123      0 
          76      0      1.030      0 
          76      0      0.968      0 
          77      0      0.852      0 
          77      0      0.940      0 
          77      0      0.905      0 
          77      0      0.787      0 
          77      0      0.918      0 
          77      0      0.745      0 
          77      0      0.907      0 
          77      0      0.818      0 
          77      0      0.896      0 
          77      0      0.858      0 
          77      0      0.935      0 
          77      0      0.669      0 
          77      0      0.930      0 
          77      0      0.890      0 
          84      0      0.983      0 
          84      0      0.974      0 
          84      0      1.005      0 
          84      0      0.897      0 
          84      0      1.110      0 
          84      0      1.035      0 
          84      0      0.997      0 
          84      0      1.059      0 
          84      0      1.067      0 
          84      0      0.911      0 
          85      0      1.122      0 
          85      0      1.134      0 
          85      0      1.117      0 
          85      0      0.999      0 
          85      0      1.200      0 
          85      0      1.154      0 
          85      0      1.077      0 
          85      0      0.980      0 
          85      0      0.989      0 
          85      0      1.161      0 
          85      0      1.101      0 
          85      0      1.153      0 
          93      0      1.051      0 
          93      0      1.082      0 
          93      0      1.062      0 
          93      0      0.955      0 
          93      0      1.095      0 
          93      0      0.994      0 
          93      0      1.093      0 
          93      0      1.002      0 
          93      0      1.038      0 
          93      0      1.084      0 
          93      0      1.026      0 
         100      0      0.975      0 
         100      0      1.110      0 
         100      0      1.120      0 
         100      0      1.170      0 
         100      0      1.130      0 
         100      0      1.070      0 
         100      0      1.100      0 
         100      0      1.020      0 
         101      0      1.190      0 
         101      0      1.040      0 
         101      0      0.821      0 
         101      0      1.070      0 
         101      0      0.665      0 
         101      0      1.120      0 
         106      0      1.044      0 
         106      0      0.933      0 
         106      0      0.978      0 
         106      0      0.981      0 
         106      0      0.958      0 
         106      0      1.005      0 
         106      0      1.053      0 
         106      0      1.059      0 
         106      0      0.917      0 
         106      0      1.079      0 
         106      0      1.077      0 
         106      0      1.060      0 
         106      0      1.031      0 
         106      0      1.006      0 
         106      0      1.007      0 
         107      0      0.960      0 
         107      0      0.844      0 
         107      0      0.833      0 
         107      0      0.870      0 
         107      0      0.887      0 
         107      0      0.835      0 
         107      0      0.716      0 
         107      0      0.964      0 
         107      0      0.933      0 
         107      0      0.753      0 
         107      0      0.803      0 
         107      0      0.742      0 
         114      0      0.885      0 
         114      0      0.903      0 
         114      0      0.815      0 
         114      0      0.937      0 
         114      0      0.936      0 
         114      0      0.984      0 
         114      0      0.835      0 
         114      0      0.954      0 
         114      0      0.890      0 
         114      0      0.930      0 
         114      0      0.924      0 
         114      0      0.873      0 
         115      0      1.045      0 
         115      0      1.087      0 
         115      0      1.021      0 
         115      0      0.999      0 
         115      0      0.908      0 
         115      0      1.027      0 
         115      0      1.041      0 
         115      0      0.968      0 
         115      0      0.971      0 
         115      0      1.029      0 
         115      0      0.923      0 
         115      0      0.944      0 
         115      0      0.953      0 
         122      0      1.032      0 
         122      0      1.045      0 
         122      0      0.967      0 
         122      0      1.065      0 
         122      0      0.981      0 
         122      0      0.964      0 
         122      0      1.080      0 
         122      0      1.100      0 
         122      0      0.994      0 
         122      0      0.949      0 
         123      0      1.080      0 
         123      0      1.081      0 
         123      0      1.012      0 
         123      0      0.856      0 
         123      0      0.915      0 
         123      0      0.952      0 
         123      0      0.983      0 
         123      0      0.872      0 
         123      0      0.950      0 
         123      0      0.939      0 
         130      0      1.014      0 
         130      0      0.907      0 
         130      0      0.841      0 
         130      0      0.892      0 
         130      0      0.868      0 
         130      0      0.837      0 
         130      0      0.931      0 
         130      0      0.871      0 
         130      0      0.901      0 
         130      0      0.987      0 
         130      0      0.986      0 
         131      0      0.933      0 
         131      0      0.935      0 
         131      0      1.033      0 
         131      0      1.081      0 
         131      0      0.773      0 
         131      0      1.054      0 
         131      0      0.946      0 
         131      0      0.994      0 
         131      0      0.724      0 
         131      0      0.883      0 
         131      0      1.047      0 
         131      0      0.951      0 
         138      0      0.806      0 
         138      0      0.858      0 
         138      0      0.815      0 
         138      0      0.914      0 
         138      0      0.841      0 
         138      0      0.866      0 
         138      0      0.792      0 
         138      0      0.773      0 
         138      0      0.882      0 
         138      0      0.887      0 
         138      0      0.847      0 
         138      0      1.035      0 
         138      0      0.918      0 
         139      0      0.886      0 
         139      0      0.923      0 
         139      0      1.027      0 
         139      0      0.991      0 
         139      0      0.909      0 
         139      0      1.061      0 
         139      0      1.054      0 
         139      0      1.018      0 
         139      0      1.029      0 
         139      0      1.100      0 
         139      0      0.932      0 
         139      0      0.886      0 
         139      0      1.024      0 
         139      0      0.930      0 
         146      0      1.025      0 
         146      0      1.014      0 
         146      0      0.954      0 
         146      0      1.081      0 
         146      0      0.967      0 
         146      0      0.912      0 
         146      0      1.030      0 
         146      0      0.871      0 
         146      0      0.940      0 
         146      0      0.996      0 
         146      0      1.028      0 
         146      0      0.794      0 
         146      0      0.879      0 
         147      0      1.058      0 
         147      0      1.100      0 
         147      0      1.121      0 
         147      0      1.093      0 
         147      0      1.017      0 
         147      0      1.058      0 
         147      0      1.090      0 
         147      0      1.052      0 
         147      0      1.046      0 
         147      0      1.103      0 
         147      0      1.026      0 
         147      0      1.076      0 
         147      0      1.095      0 
         153      0      1.037      0 
         153      0      0.987      0 
         153      0      1.025      0 
         153      0      0.967      0 
         153      0      0.785      0 
         153      0      0.861      0 
         153      0      0.893      0 
         153      0      0.941      0 
         153      0      0.762      0 
         153      0      0.976      0 
         153      0      1.019      0 
         153       0     0.889      0 
         153       0     0.938      0 
         153       0     0.961      0 
         154       0     0.999      0 
         154       0     0.880      0 
         154       0     0.971      0 
         154       0     0.924      0 
         154       0     0.939      0 
         154       0     0.966      0 
         154       0     0.981      0 
         154       0     0.974      0 
         154       0     0.965      0 
         154       0     0.938      0 
         154       0     0.903      0 
         154       0     1.094      0 
         154       0     0.867      0 
         154       0     1.007      0 
          54     750     0.768      0 
          54     750     0.731      0 
          54     750     0.716      0 
          54     750     0.693      1 
          54     750     0.769      0 
          54     750     0.731      0 
          54     750     0.791      0 
          54     750     0.765      1 
          54     750     0.611      1 
          54     750     0.745      0 
          59     750     0.747      0 
          59     750     0.755      0 
          59     750     0.600      0 
          59     750     0.796      0 
          59     750     0.719      0 
          59     750     0.458      0 
          59     750     0.599      1 
          59     750     0.635      0 
          59     750     0.775      1 
          59     750     0.660      0 
          59     750     0.736      0 
          59     750     0.745      0 
          59     750     0.732      0 
          59     750     0.653      1 
          62     750     0.869      0 
          62     750     0.834      1 
          62     750     0.892      0 
          62     750     0.542      0 
          62     750     0.940      0 
          62     750     0.814      1 
          62     750     0.872      0 
          62     750     0.772      0 
          62     750     0.847      0 
          62     750     0.946      0 
          62     750     0.742      0 
          67     750     0.959      0 
          67     750     0.931      0 
          67     750     0.967      0 
          67     750     0.851      0 
          67     750     1.061      0 
          67     750     0.866      0 
          67     750     0.838      0 
          67     750     0.911      0 
          67     750     0.885      0 
          67     750     1.041      0 
          67     750     0.940      0 
          67     750     0.951      0 
          70     750     0.836      0 
          70     750     0.897      0 
          70     750     0.738      0 
          70     750     0.793      1 
          70     750     0.764      0 
          70     750     0.759      0 
          70     750     0.872      0 
          70     750     0.802      0 
          70     750     0.857      0 
          70     750     0.875      0 
          70     750     0.830      0 
          70     750     0.778      0 
          75     750     0.860      0 
          75     750     0.871      0 
          75     750     0.794      1 
          75     750     0.879      0 
          75     750     0.674      1 
          75     750     0.817      0 
          75     750     0.913      0 
          75     750     0.846      0 
          75     750     0.931      1 
          75     750     0.953      0 
          75     750     0.809      0 
          75     750     0.918      0 
          75     750     0.850      0 
          78     750     0.940      0 
          78     750     0.871      0 
          78     750     0.834      0 
          78     750     0.963      0 
          78     750     0.945      0 
          78     750     0.966      0 
          78     750     0.911      0 
          78     750     0.934      0 
          78     750     1.018      0 
          78     750     0.948      0 
          78     750     0.948      0 
          78     750     0.947      0 
          78     750     0.963      0 
          78     750     0.837      0 
          78     750     0.855      0 
          86     750     0.815      0 
          86     750     0.961      0 
          86     750     0.920      0 
          86     750     0.963      0 
          86     750     0.812      0 
          86     750     0.960      0 
          86     750     0.982      0 
          86     750     0.966      0 
          86     750     0.874      0 
          86     750     0.854      0 
          86     750     0.900      0 
          91     750     0.943      0 
          91     750     0.898      0 
          91     750     0.895      0 
          91     750     0.938      0 
          91     750     0.870      0 
          91     750     0.878      0 
          91     750     0.849      0 
          91     750     0.802      0 
          94     750     1.014      0 
          94     750     1.043      0 
          94     750     0.967      0 
          94     750     0.967      0 
          94     750     0.936      0 
          94     750     0.996      1 
          94     750     0.909      0 
          94     750     0.962      0 
          94     750     0.714      0 
          94     750     0.979      0 
          94     750     0.921      0 
          99     750     1.010      0 
          99     750     1.053      0 
          99     750     1.006      0 
          99     750     0.961      0 
          99     750     1.043      0 
          99     750     0.984      0 
          99     750     1.034      0 
          99     750     1.031      0 
          99     750     0.894      0 
          99     750     1.004      0 
         102     750     0.970      0 
         102     750     0.952      0 
         102     750     0.932      0 
         102     750     0.954      0 
         102     750     1.030      0 
         102     750     0.979      0 
         102     750     0.890      0 
         102     750     1.100      0 
         102     750     0.965      0 
         102     750     1.110      0 
         102     750     0.965      0 
         102     750     1.010      0 
         102     750     1.080      1 
         102     750     1.020      0 
         108     750     0.812      0 
         108     750     0.813      0 
         108     750     0.833      0 
         108     750     0.739      0 
         108     750     0.819      0 
         108     750     0.874      0 
         108     750     0.800      0 
         108     750     0.728      0 
         108     750     0.774      1 
         108     750     0.877      0 
         116     750     0.564      0 
         116     750     0.835      0 
         116     750     0.722      0 
         116     750     0.778      0 
         116     750     0.810      0 
         116     750     0.816      0 
         116     750     0.838      0 
         116     750     0.691      0 
         116     750     0.695      0 
         116     750     0.903      0 
         116     750     0.731      0 
         116     750     0.811      0 
         116     750     0.770      0 
         116     750     0.791      0 
         121     750     0.928      0 
         121     750     0.879      0 
         121     750     0.906      1 
         121     750     0.841      0 
         121     750     0.906      0 
         121     750     0.899      0 
         121     750     0.888      0 
         121     750     0.865      0 
         121     750     0.914      0 
         121     750     0.859      0 
         121     750     0.927      0 
         121     750     0.905      0 
         124     750     0.936      0 
         124     750     0.930      0 
         124     750     0.845      0 
         124     750     0.757      0 
         124     750     0.868      0 
         124     750     0.759      0 
         124     750     0.931      0 
         124     750     0.920      0 
         124     750     0.921      1 
         124     750     0.848      0 
         129     750     1.016      0 
         129     750     0.988      0 
         129     750     0.900      0 
         129     750     0.820      0 
         129     750     0.949      0 
         129     750     0.826      1 
         129     750     0.828      0 
         129     750     1.003      0 
         129     750     0.924      0 
         129     750     0.879      0 
         129     750     0.913      0 
         132     750     0.835      1 
         132     750     0.876      0 
         132     750     0.894      0 
         132     750     0.841      0 
         132     750     0.853      0 
         132     750     0.733      0 
         132     750     0.860      0 
         132     750     0.814      0 
         132     750     0.835      0 
         132     750     0.995      1 
         132     750     0.864      0 
         132     750     0.870      0 
         132     750     0.802      1 
         137     750     0.875      0 
         137     750     0.925      0 
         137     750     0.813      0 
         137     750     0.819      0 
         137     750     0.961      0 
         137     750     0.957      0 
         137     750     0.903      0 
         137     750     0.895      0 
         137     750     0.876      0 
         137     750     0.943      0 
         137     750     0.882      0 
         137     750     0.831      0 
         137     750     0.911      0 
         137     750     0.857      0 
         137     750     0.912      0 
         140     750     0.874      0 
         140     750     0.867      0 
         140     750     0.815      0 
         140     750     0.832      0 
         140     750     0.875      0 
         140     750     0.876      0 
         140     750     0.820      0 
         140     750     0.857      0 
         140     750     0.861      0 
         140     750     0.858      0 
         140     750     0.878      0 
         140     750     0.836      0 
         140     750     0.828      0 
         140     750     0.903      0 
         140     750     0.874      0 
         145     750     0.957      1 
         145     750     1.074      0 
         145     750     0.993      0 
         145     750     0.854      0 
         145     750     1.043      1 
         145     750     0.912      0 
         145     750     0.952      0 
         148     750     1.008      0 
         148     750     1.048      0 
         148     750     0.911      0 
         148     750     0.998      0 
         148     750     0.851      0 
         148     750     0.868      1 
         148     750     0.842      0 
         148     750     0.744      0 
         148     750     1.037      0 
         148     750     1.027      0 
         148     750     1.015      0 
         148     750     1.033      0 
         152     750     1.045      0 
         152     750     1.026      0 
         152     750     1.116      0 
         152     750     0.794      1 
         152     750     1.038      0 
         152     750     1.105      0 
         152     750     0.939      0 
         155     750     0.809      0 
         155     750     0.914      0 
         155     750     0.772      1 
         155     750     0.837      0 
         155     750     0.960      0 
         155     750     0.913      0 
         155     750     0.883      0 
         155     750     0.892      0 
         155     750     0.831      0 
          55    1500     0.673      1 
          55    1500     0.750      1 
          55    1500     0.680      0 
          55    1500     0.501      0 
          55    1500     0.712      1 
          55    1500     0.662      1 
          55    1500     0.631      1 
          55    1500     0.647      1 
          55    1500     0.708      1 
          55    1500     0.657      1 
          58    1500     0.670      1 
          58    1500     0.716      0 
          58    1500     0.749      0 
          58    1500     0.672      1 
          58    1500     0.676      1 
          58    1500     0.682      1 
          58    1500     0.735      1 
          58    1500     0.717      1 
          58    1500     0.727      0 
          58    1500     0.679      0 
          58    1500     0.674      0 
          63    1500     0.673      0 
          63    1500     0.644      0 
          63    1500     0.605      0 
          63    1500     0.627      0 
          63    1500     0.457      1 
          63    1500     0.656      1 
          63    1500     0.665      0 
          63    1500     0.748      0 
          63    1500     0.755      1 
          63    1500     0.743      1 
          63    1500     0.744      1 
          63    1500     0.645      0 
          66    1500     0.574      0 
          66    1500     0.619      1 
          66    1500     0.607      1 
          66    1500     0.637      1 
          66    1500     0.573      1 
          66    1500     0.600      1 
          66    1500     0.623      1 
          66    1500     0.554      1 
          66    1500     0.658      1 
          66    1500     0.624      0 
          66    1500     0.582      1 
          66    1500     0.526      1 
          66    1500     0.467      1 
          71    1500     0.722      1 
          71    1500     0.419      1 
          71    1500     0.704      1 
          71    1500     0.718      1 
          71    1500     0.695      0 
          71    1500     0.867      1 
          71    1500     0.800      1 
          71    1500     0.827      0 
          71    1500     0.756      0 
          74    1500     0.677      1 
          74    1500     0.749      0 
          74    1500     0.686      1 
          74    1500     0.788      0 
          74    1500     0.719      0 
          74    1500     0.652      1 
          74    1500     0.771      0 
          74    1500     0.740      1 
          74    1500     0.570      1 
          74    1500     0.778      0 
          74    1500     0.809      0 
          74    1500     0.658      1 
          74    1500     0.758      1 
          79    1500     1.154      0 
          82    1500     0.760      0 
          82    1500     0.786      1 
          82    1500     0.828      0 
          82    1500     0.791      0 
          82    1500     0.644      1 
          82    1500     0.815      0 
          82    1500     0.730      0 
          82    1500     0.872      0 
          82    1500     0.943      0 
          82    1500     0.833      0 
          82    1500     0.809      1 
          87    1500     0.925      0 
          87    1500     0.866      0 
          87    1500     0.789      1 
          87    1500     0.857      0 
          87    1500     0.904      1 
          87    1500     0.793      1 
          87    1500     0.800      1 
          87    1500     0.822      1 
          87    1500     0.789      1 
          87    1500     0.711      1 
          87    1500     0.703      1 
          90    1500     0.654      0 
          90    1500     0.707      0 
          90    1500     0.674      0 
          90    1500     0.629      0 
          90    1500     0.772      0 
          90    1500     0.784      0 
          90    1500     0.861      0 
          90    1500     0.700      0 
          90    1500     0.763      1 
          90    1500     0.834      0 
          90    1500     0.645      1 
          90    1500     0.716      0 
          90    1500     0.712      0 
          90    1500     0.714      0 
          95    1500     0.942      0 
          95    1500     0.871      0 
          95    1500     0.973      0 
          95    1500     0.893      0 
          95    1500     0.788      0 
          95    1500     0.835      0 
          95    1500     0.917      0 
          95    1500     0.821      1 
          95    1500     0.870      0 
          95    1500     0.816      0 
          95    1500     0.848      0 
          95    1500     0.951      0 
          95    1500     0.874      0 
          98    1500     0.881      0 
          98    1500     0.804      0 
          98    1500     0.795      1 
          98    1500     0.733      1 
          98    1500     0.656      1 
          98    1500     0.850      1 
          98    1500     0.737      1 
          98    1500     0.851      0 
          98    1500     0.883      1 
          98    1500     0.814      1 
          98    1500     0.715      1 
          98    1500     0.728      1 
          98    1500     0.831      0 
          98    1500     0.792      0 
          98    1500     0.711      1 
         103    1500     0.895      0 
         103    1500     0.931      1 
         103    1500     0.916      1 
         103    1500     0.848      0 
         103    1500     0.654      0 
         103    1500     0.613      1 
         103    1500     0.842      0 
         103    1500     0.669      1 
         103    1500     0.831      0 
         103    1500     0.788      1 
         103    1500     0.835      0 
         105    1500     0.877      0 
         105    1500     0.954      0 
         105    1500     0.886      0 
         105    1500     0.895      0 
         105    1500     0.840      0 
         105    1500     0.723      0 
         105    1500     0.797      0 
         105    1500     0.836      0 
         105    1500     0.726      0 
         105    1500     0.727      1 
         105    1500     0.837      0 
         105    1500     0.715      0 
         112    1500     0.810      0 
         112    1500     0.812      0 
         112    1500     0.782      1 
         112    1500     0.757      0 
         112    1500     0.749      0 
         112    1500     0.779      0 
         112    1500     0.686      1 
         112    1500     0.760      1 
         112    1500     0.795      0 
         112    1500     0.785      0 
         112    1500     0.708      0 
         112    1500     0.660      0 
         112    1500     0.644      0 
         117    1500     0.785      0 
         117    1500     0.762      0 
         117    1500     0.744      0 
         117    1500     0.709      0 
         117    1500     0.718      0 
         117    1500     0.785      0 
         117    1500     0.820      0 
         117    1500     0.674      0 
         117    1500     0.811      0 
         117    1500     0.845      0 
         117    1500     0.693      0 
         117    1500     0.743      0 
         117    1500     0.723      0 
         125    1500     0.705      1 
         125    1500     0.607      1 
         125    1500     0.709      1 
         125    1500     0.715      0 
         125    1500     0.713      1 
         125    1500     0.765      1 
         125    1500     0.723      1 
         128    1500     0.820      0 
         128    1500     0.909      0 
         128    1500     0.925      0 
         128    1500     0.884      0 
         128    1500     0.788      1 
         128    1500     0.822      0 
         128    1500     0.816      0 
         128    1500     0.789      1 
         136    1500     0.877      0 
         136    1500     1.000      0 
         136    1500     0.817      1 
         141    1500     0.808      0 
         141    1500     0.879      0 
         141    1500     0.816      0 
         141    1500     0.939      0 
         141    1500     0.824      0 
         141    1500     0.922      0 
         141    1500     0.947      0 
         144    1500     0.792      0 
         144    1500     0.868      0 
         144    1500     0.526      1 
         144    1500     0.816      0 
         144    1500     0.770      1 
         144    1500     0.760      1 
         144    1500     0.771      1 
         144    1500     0.854      0 
         144    1500     0.921      0 
         144    1500     0.840      0 
         149    1500     0.844      0 
         149    1500     0.957      0 
         149    1500     0.929      0 
         149    1500     0.809      0 
         149    1500     0.883      0 
         149    1500     0.856      0 
         149    1500     0.755      0 
         149    1500     0.693      0 
         149    1500     0.914      0 
         149    1500     0.905      0 
         149    1500     0.860      0 
         149    1500     0.828      0 
          56    3000     0.678      1 
          56    3000     0.765      0 
          56    3000     0.518      1 
          56    3000     0.629      1 
          56    3000     0.763      0 
          56    3000     0.658      1 
          56    3000     0.732      0 
          57    3000     0.550      1 
          57    3000     0.633      1 
          57    3000     0.498      1 
          57    3000     0.633      1 
          57    3000     0.619      0 
          57    3000     0.521      1 
          57    3000     0.507      1 
          57    3000     0.696      1 
          57    3000     0.599      1 
          57    3000     0.605      1 
          64    3000     0.718      0 
          64    3000     0.615      0 
          64    3000     0.634      0 
          64    3000     0.514      0 
          64    3000     0.716      1 
          64    3000     0.708      0 
          64    3000     0.668      1 
          64    3000     0.693      0 
          64    3000     0.627      1 
          64    3000     0.623      0 
          64    3000     0.734      1 
          64    3000     0.599      0 
          65    3000     0.702      1 
          65    3000     0.678      1 
          65    3000     0.733      1 
          65    3000     0.656      1 
          65    3000     0.647      1 
          65    3000     0.719      1 
          65    3000     0.644      1 
          65    3000     0.628      1 
          65    3000     0.653      1 
          65    3000     0.597      1 
          65    3000     0.531      1 
          72    3000     0.624      1 
          72    3000     0.827      1 
          72    3000     0.846      1 
          72    3000     0.778      1 
          72    3000     0.764      0 
          72    3000     0.838      1 
          72    3000     0.841      1 
          72    3000     0.919      1 
          72    3000     0.852      0 
          72    3000     0.833      0 
          73    3000     0.665      1 
          73    3000     0.665      1 
          73    3000     0.563      0 
          73    3000     0.649      1 
          73    3000     0.617      1 
          73    3000     0.617      1 
          73    3000     0.550      1 
          73    3000     0.621      1 
          73    3000     0.607      1 
          73    3000     0.599      1 
          80    3000     1.025      1 
          80    3000     0.490      1 
          80    3000     1.004      1 
          80    3000     0.968      0 
          81    3000     0.687      1 
          81    3000     0.593      1 
          81    3000     0.631      1 
          81    3000     0.713      1 
          81    3000     0.747      1 
          81    3000     0.685      1 
          81    3000     0.791      1 
          89    3000     0.660      0 
          89    3000     0.640      0 
          89    3000     0.756      0 
          89    3000     0.727      1 
          96    3000     0.619      1 
          96    3000     0.601      1 
          96    3000     0.665      1 
          96    3000     0.572      1 
          96    3000     0.472      1 
          96    3000     0.561      1 
          96    3000     0.637      1 
          96    3000     0.627      1 
          96    3000     0.621      1 
          96    3000     0.637      1 
          96    3000     0.596      1 
          96    3000     0.591      1 
          97    3000     0.686      1 
          97    3000     0.685      1 
          97    3000     0.464      1 
          97    3000     0.649      1 
          97    3000     0.595      1 
          97    3000     0.523      1 
          97    3000     0.640      1 
          97    3000     0.618      1 
          97    3000     0.566      1 
          97    3000     0.539      1 
          97    3000     0.606      1 
          97    3000     0.585      1 
         104    3000     0.684      1 
         104    3000     0.748      1 
         104    3000     0.701      1 
         104    3000     0.484      0 
         104    3000     0.779      1 
         104    3000     0.734      1 
         104    3000     0.777      1 
         104    3000     0.893      1 
         104    3000     0.812      0 
         104    3000     0.852      1 
         104    3000     0.770      1 
         104    3000     0.734      1 
         104    3000     0.733      0 
         104    3000     0.751      1 
         110    3000     0.686      1 
         110    3000     0.762      0 
         110    3000     0.798      0 
         110    3000     0.804      0 
         110    3000     0.779      0 
         110    3000     0.731      0 
         110    3000     0.830      1 
         110    3000     0.644      1 
         110    3000     0.764      0 
         110    3000     0.750      1 
         110    3000     0.721      1 
         119    3000     0.669      0 
         119    3000     0.716      0 
         119    3000     0.717      0 
         119    3000     0.629      0 
         119    3000     0.690      0 
         119    3000     0.755      1 
         119    3000     0.663      0 
         119    3000     0.618      1 
         119    3000     0.717      1 
         119    3000     0.690      1 
         119    3000     0.675      0 
         119    3000     0.677      0 
         126    3000     0.727      0 
         126    3000     0.677      0 
         126    3000     0.623      1 
         126    3000     0.618      1 
         126    3000     0.641      1 
         126    3000     0.592      0 
         126    3000     0.768      1 
         126    3000     0.816      0 
         126    3000     0.498      0 
         126    3000     0.644      1 
         126    3000     0.743      0 
         126    3000     0.662      0 
         127    3000     0.713      0 
         127    3000     0.411      1 
         127    3000     0.768      0 
         127    3000     0.665      0 
         127    3000     0.711      0 
         127    3000     0.719      0 
         127    3000     0.710      0 
         127    3000     0.722      0 
         127    3000     0.603      1 
         127    3000     0.662      0 
         127    3000     0.649      1 
         127    3000     0.614      0 
         134    3000     0.402      1 
         134    3000     0.507      1 
         134    3000     0.451      1 
         134    3000     0.517      1 
         134    3000     0.475      1 
         134    3000     0.576      1 
         134    3000     0.651      0 
         134    3000     0.618      1 
         134    3000     0.539      1 
         134    3000     0.705      1 
         135    3000     0.687      0 
         135    3000     0.847      0 
         135    3000     0.817      0 
         135    3000     0.842      0 
         135    3000     0.699      0 
         135    3000     0.838      0 
         135    3000     0.801      0 
         135    3000     0.635      0 
         135    3000     0.761      0 
         142    3000     0.696      0 
         142    3000     0.573      1 
         142    3000     0.656      1 
         142    3000     0.688      1 
         142    3000     0.689      1 
         142    3000     0.587      0 
         142    3000     0.608      0 
         142    3000     0.643      1 
         142    3000     0.618      0 
         143    3000     1.198      0 
         143    3000     0.758      0 
         143    3000     0.929      0 
         143    3000     0.930      0 
         143    3000     1.005      1 
         143    3000     0.917      1 
         143    3000     0.910      0 
         150    3000     0.733      0 
         150    3000     0.731      1 
         150    3000     0.810      1 
         150    3000     0.934      0 
         150    3000     0.785      0 
         150    3000     0.799      0 
         150    3000     0.742      0 
         150    3000     0.921      0 
         150    3000     0.957      0 
         150    3000     0.894      0 
         150    3000     0.870      0 
         150    3000     0.791      1 
         151    3000     0.883      0 
         151    3000     0.900      0 
         151    3000     0.756      1 
         151    3000     0.776      0 
         151    3000     0.857      0 
         151    3000     0.769      0 
         151    3000     0.890      1 
         156    3000     0.811      1 
         156    3000     1.002      0 
         156    3000     0.877      0 
         156    3000     0.790      0 
         156    3000     0.902      0 
         156    3000     0.864      0 
         156    3000     0.856      0 
         156    3000     0.797      0 
         156    3000     0.839      0 
         156    3000     0.867      0 
         156    3000     0.724      0 
         156    3000     0.829      0 
;
run;


options nocenter mprint;

%ml(data=EG_data,
     bin=Fetal_Mal,norm=Fetal_Wt,
     x=Dose,
    id=Litter_ID, 
method=ISAMP,
qpoints=200,maxit=100,
tech=NRRIDG,it_his=yes);