Mercurial > repos > glogobyte > isoread

comparison mirgene.py @ 30:d119424f67b6 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author glogobyte
date Wed, 20 Oct 2021 16:24:32 +0000
parents 537dcf5b0133
children 23a88a3ec37d
comparison
equal deleted inserted replaced
29:0e7bd3f72c2c 30:d119424f67b6
204 control_group=list(con_list) 204 control_group=list(con_list)
205 treated_group=list(tre_list) 205 treated_group=list(tre_list)
206 206
207 # Filters low count mirnas (otpional) 207 # Filters low count mirnas (otpional)
208 if int(args.per)!=-1: 208 if int(args.per)!=-1:
209 209 if int(args.per)>0 and float(args.per)<=100 and int(args.count)>0:
210 fil_con_group=manager.list() 210 fil_con_group=manager.list()
211 fil_tre_group=manager.list() 211 fil_tre_group=manager.list()
212 212
213 ps_low_counts = Process(target=filter_low_counts,args=(control_group,treated_group,fil_con_group,fil_tre_group,args.per,args.count)) 213 ps_low_counts = Process(target=filter_low_counts,args=(control_group,treated_group,fil_con_group,fil_tre_group,args.per,args.count))
214 ps_low_counts.start() 214 ps_low_counts.start()
215 ps_low_counts.join() 215 ps_low_counts.join()
216 216
217 fil_con_group=list(fil_con_group) 217 fil_con_group=list(fil_con_group)
218 fil_tre_group=list(fil_tre_group) 218 fil_tre_group=list(fil_tre_group)
219
220 else:
221 sys.exit("Not acceptable values for filter")
219 222
220 if "fil_con_group" not in locals() or "fil_con_group" not in globals(): 223 if "fil_con_group" not in locals() or "fil_con_group" not in globals():
221 fil_con_group=control_group 224 fil_con_group=control_group
222 fil_tre_group=treated_group 225 fil_tre_group=treated_group
223 226
377 n_control_group=list(n_con_list) 380 n_control_group=list(n_con_list)
378 n_treated_group=list(n_tre_list) 381 n_treated_group=list(n_tre_list)
379 382
380 # Filters low count mirnas (otpional) 383 # Filters low count mirnas (otpional)
381 if int(args.per)!=-1: 384 if int(args.per)!=-1:
382 385 if int(args.per)>0 and float(args.per)<=100 and int(args.count)>0:
383 n_fil_con_group=manager.list() 386
384 n_fil_tre_group=manager.list() 387 n_fil_con_group=manager.list()
385 388 n_fil_tre_group=manager.list()
386 ps_low_counts = Process(target=filter_low_counts,args=(n_control_group,n_treated_group,n_fil_con_group,n_fil_tre_group,args.per,args.count)) 389
387 ps_low_counts.start() 390 ps_low_counts = Process(target=filter_low_counts,args=(n_control_group,n_treated_group,n_fil_con_group,n_fil_tre_group,args.per,args.count))
388 ps_low_counts.join() 391 ps_low_counts.start()
389 392 ps_low_counts.join()
390 n_fil_con_group=list(n_fil_con_group) 393
391 n_fil_tre_group=list(n_fil_tre_group) 394 n_fil_con_group=list(n_fil_con_group)
395 n_fil_tre_group=list(n_fil_tre_group)
396 else:
397 sys.exit("Not acceptable values for filter")
392 398
393 if "n_fil_con_group" not in locals() or "n_fil_con_group"not in globals(): 399 if "n_fil_con_group" not in locals() or "n_fil_con_group"not in globals():
394 n_fil_con_group=n_control_group 400 n_fil_con_group=n_control_group
395 n_fil_tre_group=n_treated_group 401 n_fil_tre_group=n_treated_group
396 402