case $# in 0|1|2|3|4|5|6) echo 'Usage: dotests_seq n_sequences n_sites file n_permutations n_localities n1 n2 ...' 1>&2; exit 1 esac nsam=$1; ns=$2; file=$3; shift; shift; shift echo 'Data file: ' $file seqtomatrix $nsam $ns <$file | permtest K 0 y $* seqtomatrix $nsam $ns <$file | logtrans | permtest 'K*' 2 n $* seqtomatrix $nsam $ns <$file | haplotrans | permtest H 2 n $* seqtomatrix $nsam $ns <$file | ranktrans | permtest Z 2 n $* seqtomatrix $nsam $ns <$file | ranktrans | logtrans | permtest 'Z*' 2 n $* seqtomatrix $nsam $ns <$file | permchi n $*