pytable_filter=tables.Filters(complevel=9, complib='zlib')
gen=[]
for i,j in enumerate(args.genotype):
gen.append(Reader('genotype'))
gen[i].start(j,hdf5=True, study_name=args.study_name[i], ID=False)
RSID=[]
SUB_ID=[]
for i in gen:
SUB_ID.append(i.folder._data.get_id())
mapper.cluster=args.cluster
mapper.node=args.node
while True:
if args.cluster=='n':
SNPs_index, keys=mapper.get_next()
else:
chunk=mapper.chunk_pop()
if chunk is None:
SNPs_index=None
break
print chunk
SNPs_index, keys=mapper.get_chunk(chunk)
if SNPs_index is None:
break
RSID.append(keys)
data=merge_genotype(gen, SNPs_index) #TODO (high) add mapper
print data.shape
if args.cluster=='n':
