コード例 #1
0
ファイル: hase.py プロジェクト: urmovosa/hase
                            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:
コード例 #2
0
ファイル: hase.py プロジェクト: roshchupkin/hase
		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: