Script for this run: DEF2a.rTo 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 .