Mercurial > repos > glogobyte > isoread

comparison mirbase.py @ 29:0e7bd3f72c2c 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:19 +0000
parents 6e365f2eb24a
children 5d232330e81f
comparison
equal deleted inserted replaced
28:537dcf5b0133 29:0e7bd3f72c2c
212 control_group=list(con_list) 212 control_group=list(con_list)
213 treated_group=list(tre_list) 213 treated_group=list(tre_list)
214 214
215 # Filters low count mirnas (otpional) 215 # Filters low count mirnas (otpional)
216 if int(args.per)!=-1: 216 if int(args.per)!=-1:
217 217 if int(args.per)>0 and float(args.per)<=100 and int(args.count)>0:
218 fil_con_group=manager.list() 218
219 fil_tre_group=manager.list() 219 fil_con_group=manager.list()
220 220 fil_tre_group=manager.list()
221 ps_low_counts = Process(target=filter_low_counts,args=(control_group,treated_group,fil_con_group,fil_tre_group,args.per,args.count)) 221
222 ps_low_counts.start() 222 ps_low_counts = Process(target=filter_low_counts,args=(control_group,treated_group,fil_con_group,fil_tre_group,args.per,args.count))
223 ps_low_counts.join() 223 ps_low_counts.start()
224 224 ps_low_counts.join()
225 fil_con_group=list(fil_con_group) 225
226 fil_tre_group=list(fil_tre_group) 226 fil_con_group=list(fil_con_group)
227 fil_tre_group=list(fil_tre_group)
228 else:
229 sys.exit("Not acceptable values for filter")
230
227 231
228 if "fil_con_group" not in locals() or "fil_con_group" not in globals(): 232 if "fil_con_group" not in locals() or "fil_con_group" not in globals():
229 fil_con_group=control_group 233 fil_con_group=control_group
230 fil_tre_group=treated_group 234 fil_tre_group=treated_group
231 235
385 n_treated_group=list(n_tre_list) 389 n_treated_group=list(n_tre_list)
386 390
387 391
388 # Filters low count mirnas (otpional) 392 # Filters low count mirnas (otpional)
389 if int(args.per)!=-1: 393 if int(args.per)!=-1:
390 394 if int(args.per)>0 and float(args.per)<=100 and int(args.count)>0:
391 n_fil_con_group=manager.list() 395
392 n_fil_tre_group=manager.list() 396 n_fil_con_group=manager.list()
393 397 n_fil_tre_group=manager.list()
394 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)) 398
395 ps_low_counts.start() 399 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))
396 ps_low_counts.join() 400 ps_low_counts.start()
397 401 ps_low_counts.join()
398 n_fil_con_group=list(n_fil_con_group) 402
399 n_fil_tre_group=list(n_fil_tre_group) 403 n_fil_con_group=list(n_fil_con_group)
400 404 n_fil_tre_group=list(n_fil_tre_group)
405
406 else:
407 sys.exit("Not acceptable values for filter")
408
401 if "n_fil_con_group" not in locals() or "n_fil_con_group" not in globals(): 409 if "n_fil_con_group" not in locals() or "n_fil_con_group" not in globals():
402 n_fil_con_group=n_control_group 410 n_fil_con_group=n_control_group
403 n_fil_tre_group=n_treated_group 411 n_fil_tre_group=n_treated_group
404 412
405 ps_write = Process(target=write_main,args=(n_control_group, n_treated_group,n_fil_con_group, n_fil_tre_group, n_con_file_order, n_tre_file_order,2,args.group1,args.group2,args.per)) 413 ps_write = Process(target=write_main,args=(n_control_group, n_treated_group,n_fil_con_group, n_fil_tre_group, n_con_file_order, n_tre_file_order,2,args.group1,args.group2,args.per))