AggregateAssociation.R fails if all variants included in the aggregation units are monomorphic in the sample. The workflow should be robust to this and return an empty file/array