view evaluate_population_numbers.bash @ 36:51cd0307fb70

Phylip's extra ouputs are now stored in the job working directory
author Richard Burhans <burhans@bx.psu.edu>
date Wed, 20 Nov 2013 16:32:01 -0500
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