This archive provides computer code for the marginal
data augmentation algorithm using a full matrix working
parameter for posterior sampling in a two-level Gaussian
model. The algorithms are full described in
van Dyk and Meng (2000).
This code uses routines from Numerical Recipes in C.
It is violation of copyright to down-load these routines
unless you have purchased them from the publisher.
Click
here to purchase the Numerical Recipes in C.
This archive contains the source code. Once you
have down-loaded the code, you must uncompress
(gunzip) and un-tar the archive. The code can
be complied by typing make in the subdirectory
FullMarDAsampler which is created when the archived
is un-tarred:
~: gunzip marDA.tar.gz
~: tar xvf marDA.tar
~: cd FullMarDAsampler
~: make
Click here to donwload the code.
The README file tells how to format your data file
and other inputs, such as starting values. A sample
data file and inputs file are included (the mare
data set of Pierson and Ginther, Animal Reproduction
Science, 1987).
Good Luck!
Return to David van Dyk's homepage.