Yc@ sdZddlmZddlZddlZddlZddlZddlmZddl Z ddl Z ddl j Zddl Z ddlZddlZddlZddZdZded Zdd Zd Zed Zd ZdZdZdZdZdZdZdZ dZ!dZ"dZ#dZ$dZ%dZ&dZ'dZ(dS(s: Module containing important auxilary screening functions i(tdivisionN(t defaultdicticC sO|dkrdSttjtjt|}|d|}t||S(sd Function rounds number in a reasonable manner. Default is to three significant digits. ii(tinttmathtfloortlog10tabstround(txtnumtorder_of_magnitudetdigits((su/mnt/lab_data/kundaje/users/oursu/bassik/pgCTCF/results/2017_10_09/results/genes/weissman/castle/Scripts/screenFun.pytsigDigs  $cC std}t|}xQt|D]C}|dksU|dksU|dksU|dkr|d7}qqWt||S(NitCtGtctgi(tlentlisttfloat(tguidetnumGCttotaltnuc((su/mnt/lab_data/kundaje/users/oursu/bassik/pgCTCF/results/2017_10_09/results/genes/weissman/castle/Scripts/screenFun.pytgetGC0s  0tGenRefc C satjj|d}tjj|d}tjj|d}td}td}td}td}yZt|dE} tj| d d } x&| D]} | d j|| d MscS sdS(NsN/A((((su/mnt/lab_data/kundaje/users/oursu/bassik/pgCTCF/results/2017_10_09/results/genes/weissman/castle/Scripts/screenFun.pyRNscS sdS(NsN/A((((su/mnt/lab_data/kundaje/users/oursu/bassik/pgCTCF/results/2017_10_09/results/genes/weissman/castle/Scripts/screenFun.pyROscS sdS(NsN/A((((su/mnt/lab_data/kundaje/users/oursu/bassik/pgCTCF/results/2017_10_09/results/genes/weissman/castle/Scripts/screenFun.pyRPstrt delimitert,iiNs$Ensembl information file not found. sUse -r to change file locations iis"Mouse information file not found. s"Human information file not found. ( tostpathtjoinRtopentcsvtreadertuppertIOErrortnext(tref_basetmousetorg_file_humantorg_file_mousetens_filet geneID2Namet geneID2Infot geneName2IDt geneEns2Nametens_opentens_csvtlinetorg_opentorg_csv((su/mnt/lab_data/kundaje/users/oursu/bassik/pgCTCF/results/2017_10_09/results/genes/weissman/castle/Scripts/screenFun.pyt retrieveInfo?sN &    &    &  c C stjj|d}tt}tt}tt}tjj|s]dGH|||fSt|d}tj|dd}x|D]}t |dkr|ddkrqn|d d kr||d c|d d 7 sF            K 2 2  ^ K J n U   S P   # 6 5