soource code to estimate selection coefficients

This source code implements a Bayesian model to estimate selection coefficients. The software was written in C++ and requires GSL to compile. The software takes a genotype and a phenotype file (the phenotypes are just 0 or 1 and indicate whether an individual died = 0 or lived = 1). The input files contain the following,

genotype file:
line 1: number of loci, number of populations
line 2: pop1 sample size, pop2 sample size, ...
line 3: block for pop1, block for pop2, ....
line 4: treatment for pop1, treatment for pop2, ....
line 5-N: population number, locus number, individual number, Pr(AA), Pr(Aa), Pr(aa)

phenotype file:
line 1: phenotypes for pop1 (in the same order as in the genotype file)
line 2: phenotypes for pop2

Here are the command line options,

./wigs version 0.1 -- 21 August 2012

Usage: wigs -g genotypefile -s survivalfile [options]
-g Infile with genotype posterior probabilities
-s Infile with suvival data for each individual
-x Number of Markov chain Monte Carlo steps
-b Number of MCMC steps to discard as a burnin
-t Specify chain thinning, write samples every nth step
-d Prior probability a SNP experienced selection [0.5]
-f Outfile prefix