***msstats README****
Copyright 2002-2007 Kevin Thornton, University of Chicago, UC Irvine

current version 0.3.0 (Sept. 22, 2011).  All instances of outputting "not a number" are now calls to strtod("NAN",NULL), for consistency purposes. 

version 0.3.0 (Sept. 5, 2011).  Version 0.2.9 incorrectly output the ZnS statistic (a.k.a. the mean r^2).  The statistic wasn't divided by the number of site pairs in the analysis.  This bug is now fixed.

version 0.2.9 (March 20, 2010).  Rewrote the R2 and R2E calculations, and they are now orders of magnitude more efficient.  Also upgraded to more efficient calculation of ZnS (the mean r^2 LD statistic).  REQUIRES libsequence 1.6.8 OR GREATER

version 0.2.8 (Nov. 2, 2009).  In previous versions, Fu and Li's D and F were  swapped in the output.  This is now fixed.

version 0.2.7 (16 April 2007) Handling of multiple populations now has a command-line like ms:
ms 20 1000 -t 10 -I 2 10 10 0.1 | msstats -I 2 10 10 &
Also, a bug was found where Myers' and Griffiths Rm was mis-calculated for multiple populations, as were Rozas' statistics.
Kelly's Zns (mean r^2 in the sample) is now reported.  By default, all variable sites are included.  To change the minimum minor allele frequency used in the calculation, use the -m option.  For example, -m 2 requires that the minor allele be present at least twice in the sample (i.e. singletons are excluded).

version 0.2.6 (16 May 2006).  Fixed output order of Fu and Li's D* and F*.  The program now can process multiple populations, like this:
ms 20 1000 -t 10 -I 2 10 10 0.1 | msstats -I 10 10 &

version 0.2.4 (22 Sept 2005).  No longer uses deprecated libsequence functions.  Outputs more stats and output has colummn headers.x

version 0.2.3 (07 Nov 2003): should now compile with g++ 2.95

version 0.2.2 (4 Nov 2003)--compatible w/libsequence >= 1.3.3

version 0.2.1: g++3.1 compatible on OS X

version: 0.2.0 -- updated to libsequence 1.0.4, resulting in MUCH faster code