Perfect Simulation from Nonneutral Population Genetics Models with Populationg Growth and Population Subdivision
This is code used in the paper:
Fearnhead, P. (2006) Perfect Simulation from Non-neutral Population Genetic Models: Variable Population Size and Population sub-division. Submitted
Download a tar archive of the code here
This code is run within R.
To run on a Linux machine:
- Save in a new directory and type tar -xvf NonNeutralSimulation.tar
- Compile the C code with R SHLIB MultiPop_R.c and R SHLIB SinglePop_R.c
Some versions of R use R CMD SHLIB MultiPop_R.c and R CMD SHLIB SinglePop_R.c
- The files MultiPop_PSc.R and SinglePop_PSc.R are the R functions that perform simulation. Details of the input and output of these functions can be found in the headers of these files. Examples of their use are found in the remaining file.
Any queries email: p.fearnhead@lancs.ac.uk.
Last update:12 May 2006.