Post date: Jun 08, 2017 8:58:32 PM
7vi17.
1. split size and coord files by sub-group, this splits the gentype and
phenotype files
perl mkSubs.pl resid-size-6vi17-subgroups.csv pheno_size.txt
mod_geno_size.txt
perl mkSubs.pl resid-coord-centroidsize-procrustes-6vi17-subgroups.csv
pheno_coord.txt mod_geno_coord.txt
This generates 8 geno and pheno files for size and coordinates. These
are pre-pended with the sub-group ID, and are all in
king:/uufs/chpc.utah.edu/common/home/u6000989/projects/lycaeides_wings/Gemma/
2. Run gemma for all 8 sub-groups with the 17 size traits and 52
position traits, here are the perl commands (the MCMC details are
identical to the setting for all). This runs 552 jobs (552 nodes), each
of which forks 5 chains/processes to run in parallel.
perl wrap_qsub_slurm_gemma_sub.pl *_mod_geno_size.txt
perl wrap_qsub_slurm_gemma_sub.pl *_mod_geno_coord.txt