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