Script for this run: DEF2a.r
To reproduce results, run this script as follows:  Rscript DEF2a.r outfilename aafilename v227 v203,v204,V670,V7,V53,V21,V205 bio.5,v149,v665,WY v227 extra=v872,v157 ignore true true false true false SCCS v227 out.html .