Data set:  argconc.dat (the contents are described in the book on p. 136)

1. Estimation based on individual estimates.

argmixsetup.R:  R Program to calculate individual estimates as referenced
                on the bottom of page 125 - top of page 126 and bottom of
                page 136 - top of page 127

argmix.dat      The output of argmixsetup.R

argmixmod.sas   The program to fit the second stage model using proc mixed
                as described on the the top third of page 137.

argmixmod.lst   Output

2. First-order method using the nlinmix macro.

argnlinmix.sas         Program
argnlinmix_zero.lst    Output

3. First-order conditional methods using the nlinimix macro.

Give that the program here is identical to argnlinmix.sas above except
for one change described top of page 138, I don't give the program, just
the output:

argnlinmix_eblup.lst

4. First-order conditional method using R function nlme.

argnlme.R     Program
argnlme.Rout  Output

5.  Maximum likelihood via adpative Gaussian quadrature to perform the
integrations using SAS proc nlmixed.

argnlmixed.sas   Program
argnlmixed.lst   Output