Mercurial > repos > miller-lab > genome_diversity
view evaluate_population_numbers.bash @ 31:a631c2f6d913
Update to Miller Lab devshed revision 3c4110ffacc3
author | Richard Burhans <burhans@bx.psu.edu> |
---|---|
date | Fri, 20 Sep 2013 13:25:27 -0400 |
parents | 8ae67e9fb6ff |
children |
line wrap: on
line source
#!/usr/bin/env bash if [ $# -ne 3 ]; then echo "usage" exit 1 fi input_ped_file="$1" output_file="$2" max_populations="$3" ADMIXTURE=admixture for (( i=1; $i <= $max_populations; i++ )); do $ADMIXTURE --cv "$input_ped_file" $i 2>&1 | grep CV | perl -ne 's/CV error/CVE/; print;' >> "$output_file" done