Examples

Small example

We provide a small test example in the data/ folder.

HapMap example

We can run terastructure on the HapMap dataset available here: http://pngu.mgh.harvard.edu/~purcell/plink/res.shtml

In the terastructure directory (assuming you have the program built), run the commands:

wget  http://pngu.mgh.harvard.edu/~purcell/plink/dist/hapmap_r23a.zip
unzip hapmap_r23a.zip
./src/terastructure -file hapmap_r23a.bed -n 270 -l 4098136 -k 3 -rfreq 25000 -label hapmap -nthreads 4

(adjust the -nthreads option as appropriate for your machine). Note that for HapMap, the structure is so "strong" that the algorithm converges very quickly, and we choose a low -rfreq value to reflect this.