This C program was written for inference of the age of a 
rare, selectively-neutral mutation in the paper
<ul>
    Graham and Thompson (2000). A coalescent model of ancestry for
    a rare allele. Genetics. 156:375-384.
</ul>

and for calculations in my 
<a href="http://www.stat.sfu.ca/~jgraham/Papers/papers.html">PhD thesis</a>, 
pages 52-57.
The program is for age inference given the 
current copy number of the allele and 
demographic information on the overall population. 
Subtree coalescent methods of Griffiths and Tavare' (1998) are used, in which 
age is viewed as a random quantity with a probability distribution.  

<p>
Note that developed software is available from Griffiths and Tavare', and 
from Rannala and Slatkin.  Rannala and Slatkin view age as a parameter to be 
estimated, rather than a random quantity, and their software can be found at 
http://allele.bio.sunysb.edu/software.html.

<p>
The gzipped tarfile <a href="ageinf.tar.gz">ageinf.tar.gz</a> 
contains the following:
<ul> README (this file),  makefile,  ran2.c,  ageinf.c.
</ul>

<p>
To unpack the gzipped tarfile, type "gzip -d ageinf.tar.gz" or 
"gunzip ageinf.tar.gz". Then type "tar xvf ageinf.tar". 

<p>
To make the executable, called ageinf.out, type "make".
Currently, parameters are defined as constants in the program. 
(Awkward, yes, but to be improved when I get some time ...)

<p>
If you find this code useful, please email me at jgraham@stat.sfu.ca
so that I can update you when there are improvements.
Any bug reports or fixes, please email me as well.

<p>
Thanks for your interest in this work.

<p>
Jinko Graham 07/09/99


<p>
P.S. There are also C programs available for disequilibrium mapping of a 
rare allele 
<a href="http://www.stat.sfu.ca/~jgraham/Papers/Programs/DisequilibriumMapping/">here</a>.