for i in pard ]))) while True: if mapper.cluster == 'n': SNPs_index, keys = mapper.get() else: ch = mapper.chunk_pop() if ch is None: SNPs_index = None break SNPs_index, keys = mapper.get(chunk_number=ch) if isinstance(SNPs_index, type(None)): break Analyser.rsid = keys if np.sum(PD) == 0: genotype = np.array([]) with Timer() as t_g: genotype = merge_genotype(gen, SNPs_index, mapper) genotype = genotype[:, row_index[0]] print "Time to get G {}s".format(t_g.secs) #TODO (low) add interaction a_test = np.array([]) b_cov = np.array([]) C = np.array([]) a_cov = np.array([]) b4 = np.array([]) if args.protocol is not None:
while True: if mapper.cluster=='n': SNPs_index, keys=mapper.get() else: ch=mapper.chunk_pop() if ch is None: SNPs_index=None break print ch SNPs_index, keys=mapper.get(chunk_number=ch) if isinstance(SNPs_index, type(None)): break Analyser.rsid=keys if np.sum(PD)==0: genotype=np.array([]) genotype=merge_genotype(gen, SNPs_index, mapper, flip_flag=False) genotype=genotype[:,row_index[0]] #TODO (low) add interaction a_test=np.array([]) b_cov=np.array([]) C=np.array([]) a_cov=np.array([]) b4=np.array([]) if args.protocol is not None: if protocol.enable: