From cffef5f66a704858fc0a3a40c2d9339d61b9b473 Mon Sep 17 00:00:00 2001 From: Venkat Malladi <Venkat.Malladi@utsouthwestern.edu> Date: Mon, 3 Apr 2017 11:18:30 -0500 Subject: [PATCH] Add in overlap of enhacners. --- intervene_test/Enhancer_percentages.png | Bin 0 -> 17440 bytes intervene_test/enhancer_percentages.py | 62 ++++++++++++++++++++++++ 2 files changed, 62 insertions(+) create mode 100644 intervene_test/Enhancer_percentages.png create mode 100644 intervene_test/enhancer_percentages.py diff --git a/intervene_test/Enhancer_percentages.png b/intervene_test/Enhancer_percentages.png new file mode 100644 index 0000000000000000000000000000000000000000..a90fdab321f3941a01141c6c713e2377a307f3a0 GIT binary patch literal 17440 zcmdsf2UJwqnr<nrbhm*vH!5I2zytyUlBEq;l2jA~2`WjlP~=pGR_rzx0+LaZWF><H zBT7akNJhyyha%PcHhOR0nRV~XdvDEJGw-ohF9fPi?X&kEzVM&^<Cv=AvY*%gOrcPg zQI8xvL7}j9QYZ^Xe_V{;BnYtY#<w46zfn*Ah%dJv4R7G*pDrHJqERRu=g5Bxp2$a; zQYf1!)Po8qodSCsTpv{1UzVTi*tB$$)YW^%mzE#-bBo@{$~(>x!Xjcqt=HuxhEDWv zueev(8rSJ^Nb-5(pS-%cp^;iX+G&5XAAhE?|4k8NZS3L8Uzcuw`D5XwJ$dp{CH(@W zk;~m0Z{2+-IY4vD(;FxnxjGgdHBb~4mzkL<i8bw~+V8N(cglVRb_(UvZ^z0PP$;34 z)pq!m0A&#uzWw~|PYeG){!?i4+Fhsb3AuMEaqc<0Kg+(yhs)vfW<@UhlerT)1F06Y z%2S$cxzqhS&mE;ShwU)8va%Ai|B-UGCD`ff9~+%p?{_uj(GTBV@7>+qy-O>Jm)2ff zP$ZNW=&R?lotji4Ue~ImKXoEfB|xM01?7cjRc6Hk%B2a%HgC<yq$H7%_7bJCVA(Sz zyU(P>9gL=2QhB+%pK|Q1P~ME~_-BK=6AwmsyFNv!F~Z~GcExHYmd2+y`X?nNJ@EE! z8~voF#+d0+qZdl#4&CletF!8A%yD&@=nwKzlar*+wD2_1vRPA!timszJi+}9&rRNO zITo+c)VUw69u*lE*DG5mKiPa$a<u%wr!d947SH?AB-<Yv49|>|rBA$LzjSwvX8D59 z7mJtvdaowgq$53V)^fa_R(9KpHr<!r5vxCICO6xq*3p|>rpoBe9x0_VoJZQ*{@@RO z|Nj2g<d!?q?RV)j23A>JksEE^tZ@EvZC%n>)fs-Op8HrO*|_aHcNS&Xw0y*8&rIjB zT(cbdvObL5UN5%ewI@fJShbE<dZWuui3U6M<&v%UWJZcbYfPhp#7p!SQ|vAVE3jW$ z7$9tL_>SCECvQQN{A|}om(EAbUK9C%t#@TabhGSCYx-w{73f&(rySP_4FjAE)zZ?^ z3is*Owo{+V2y|?8X3kCIg+zsTNKN-Yw`LeJhub*!UO46>q#u0f=31kuAcw-3Q?V!1 z)xF4`N;J9Aii9_;U!PQG^}LN|?hB9hZ8i#BK0SkN|79-wisOE=t`51AZ-v@3JA#Ye z{IN>ZZE~>p>FedpHX%Bb-p|S_Vj{wZ^EX%<`szARaq=TL#q?s85INb@)YRQ=GE|qb zZf|up+0@D*<Kl<+*X-12eW0#d(36<gm}HS!slh2_cY)h|^n=4>1!Zukb*&>sAcKuV zX3s=DO|2?HH=rZfRdAxw)jHvMf0kR9l73-wnOwl;1J|<$a;@`b9r%<V-oxp)S8D2) zuySpz+6ujQ)8{^g%sr^Bt)1x2D-wL3qD>!qyUvUE_4Q?Soy;9OcF2{Ruu<kbuPmpW zxiLIEywR!ovgF_kw&b3K=k2#dOW!v*HU%jOX12OdNMaSD_T3e?I7RtzJJE4*%34}l zhwp6kTed;mt$8WEu4BB;TFu^m%(Y$Gah$%J()@Zj?=j`+BB|y}Yz*GSoN;YNUT%TA zfya1qh*8Cp<IHsWbh^(oG4m?Jw|}mQPuqYTEwYYbSr2ee*A;b4jdpouWo2=&NKOqd zm7=Wv9ClaE|8`E<*kyL^BMg!hR0g-KOLmkUe}Tps%dz0uJfyYJ?9^C%{&fzXU=rCi zCLxxbXL8+~`!ibZOju|hKYl+pA!|oZRlI8%b8e<Xo)r|bl;R{YpAMz66ZP4MW67a6 zt3;g!o8xq!i>t{D9rPIX3Na+{!f6tCpSOn2njIN@&CyY7o~V-NF)MPh`NdK$slD&G zT`jEgX2(xzXc)i0x0%EXx5wm{QDkSOvtUN^m6=XO`NHFhT;8-Y+4MG(s;P;AkhUV< zbe%_uPoBJ2k!X03h~TZkS&G=UzPG@ow6dLS=33<EwuUf^O@?O&==5Zk?1KG@%_dcp z&=-+vA;HfFa;3=<^_<(h$#Q3>C(1r3@&w%!%o1>)8l^IJ>16CN32}F5FA4BXZ}N~1 zJAAt^#q?9s`L}=G;*u7hXrj}r;<UWV+}cGmzg*uSf~Y1b9q-}OhrKIt?g;Yb4{;Nz zv(9rre3v^2``wCqpyclE9upIzo#x*A?49jcrKU+~prompS>bZ&-g^f)WcPLi(p1Nr z=&X|2>K`avR}f&?iNXf?t)B8T=QQHAz2nnrV=8PYp$A5AviDjO<qO91a)$CZWX_J) z(cW+N6zHHl<#bxiMLEK#sH$qevGvaT{^tY1*RNk!_7mB5`}S>4c@MYfumvyPk!$l( za|f>pwVKFJYgy$o+fhM{uzN)hwyWEIKIiEzwVvX(i(Cv{?}aEDI^`JUXDNJJeaWO; z-dNn?^SuZqU*kH<bQQ#ObPN8n?PGm?AZmgrqDd)AHKAec<tuZ-$#z{ID=I1;NC6HU z`$frTSHT_W{?pyR9#mj`u5qZd^qLvV`+d!tHQQ9#D3k|B4Gav-%*?c<7Dhyc(u-&M z(#HS%^Jmh02RqF#JC);O{)@WsXi4kPJDxI|cU-<6vr`f1Cld{h`mxOR_I4e&p|@Ri z>DF0}1HnkA?8&b;<I}W|5UR3c6(>)g3^yqYmh6d3PBN>EDPam@NuLj{9xX8(>#oH4 zyNj5948L8f7PwDXqS48KPc1kwlhzTbc$b?~&dmvd^C_t$Ix5(|BZMv|YFzw2(PJ!D zW^S@Ya=JHJqFNDGJG<(8U%hcR0L`u5bs$Ib%kQfgGh^M(NS@(%YxkW^bK~h=ygM}$ zluw;Hg?8$5f2(4<W8Spx^`)Hijpq{E@*@mo6jJ@kSL=ogcP57G=XubpGdFD5psuCm z=e<)i^6}&IpWOBX<P6$f{c5u)hq89J_JgUZsfx<V)LW;JlAYy`NJ=sc$9rqKP#TmZ z)fm_iyY7nZw2o5uCIbmEv505So>eXcY8tW=w5Krec3vu$6gz}i_sPVMqiEQbZXKmT zpU|*<EX9pq{r*$t=treGNt<csFJ5)NQVY*K+)SUOMQ?q2ZAkVtdz<TZj>3thJP!sq z8}D9ee0h6=grLE@JC<7378c4XDk?v%N==*7yz4DRK6P`(+Z;+vY^)dWvSrKS^U;j7 z+wo(vMNN{%$UydpDtcFAQc1pDU#%FTEPgHa_fJIyp-Y$9*dw#Dq(n_hEjZ0BEQ;$h zY3gUr1n@@3$Cv!MW~T^!ZpO!w8-3a0^V7`-JeN20`T-AK(6_Q_QYP65Kof1`{^9W< z=kXd7a%-alP4Wi$w+aZNP`J9-sQJZZl5kUFJ=MdJh^;2s5Y|k`NNS?{ZZCedQZzo> zQppa(*7x`66C3Im$lZ}S%0{u+jzDTfcog1O<gQ7YOYeB42ry;(DYEU#62bS~m9aE@ zoC%6QFu(y3lq);bWup-3{*0}pe#lMNltrpXsd5URKWLRblwUHZg&V#=_h6qjh`f+z zpV54+Ly7D1?qa<I;{x_Sr0vNynGAI+8&q(<Da~{8c;i>S6e}b+LB7BIKY)BkuQ-pC zNFa+IF8y`SriSN^WiQ!z0)e|yizOr^I=<c%Qf<n0CwK=j-S&2!{;qT9u^+WpujuX3 zpY7HZ#evm(O!sL$J*yx$ky#dzkRX_tJ7vZ&{QBmPOw?9U=i%1kAZZC%m(llwUtT90 zm-w&b2R&#GJ667=wJ*J?w2VGCqKsVk<Z?LeP-EmfQ<OKOME+54_esUn%Bpm%T0exe zmUldKxA8{TfqE1@Ex<9<=G+tS$hEeTYXvzr79dtW#h*#nXd!oE?J*LtOFKniY-}t# zHrCuPHZHEX{LukX*YVyGu61B24&kBe-lL-j0ITSDFx$vQs}l7?cxF2lN!seU_N7)- zSMyUz*Xh$db*fTi!V39gYGGlKid{G+^RvrKa`$AFqzNVHK5qvWE8ZRz$WRTEtCx@Y zX%$~MDv(ihbhK&58@wlleNoO?heqd)Vt=unu^RD(jBQH8tzK%z4O#ZgR(^Ruz<yrp zfddC(-Pe;;Y_zP;u+hr2H7ZInEMTb2apD@2&8tVRbbkm_oEq;7$eU>sa_)TOX*)CY z&IBnFosdvk>?e9A?XR4C{-NdHWN#ELnytJR9rFQdHGwqFT@k*+K?835SJ^pA3Esex z3Fgw&u{pGO+4=`ab>o@Ip|0vgebs$mem|l3w@o(&-s6<8(jiD-w7YWW*2A~ndu!_X z3F>C?QVINs)60<}NI)a*^!)?NC3!(T=w$&-NRQ7+55y2+zEVE~sG(==*ts*=wyp5D z`u^O1m0>x91%30WcTavMvtk8fb>fC*&UteA!#8p{-)09_#IY(3BJNP^en|*LN8G<a z;NLV~bZ!b|QRsjBpY%5MR40~A4z-8^%GnJx%1RGp(**NohWAQHq+5If61XXpyU%@k zTxFuaL0wbRcb&f5v)`jh3FY|N<+2W~L@cu_GJpUy`Ml|d0oxL>YLax!#);hOq!N$? zTFY$-<Kfl<KWRrRm(kAqUtW8fl=zEnQ^pRwP^<)Toa7-e?vDK2<VIStFqIMHP`h1z zZs_g^0+Nr}Q#k^fAh%A}?l^e~Z-HfcgAY=n-m%G@Q_|M3APL($_#xwp-X?Y)w~NN5 zfqtyH0hTI52C$u3V(sAt>?}o1y71vKvM4GcA;Bo5u&|K8)?|}1L8>thqSkX~=OSBB zxyWbFgh})PQ+uP}cc9-ElP1p_otUV`IC<(+YHcDu+o-?6-`}6$p`^yopELgWNS1J% zECGZ~jRr#6S}HwI{7Fqd=!|ua7-6ga_{Z1*9vZqP2!w?e&{3^0!5n>u6ol>#xOGTp z4!&9<3iR8XTCK+{)XZg!l*!Al;!{4tIH{@0PbGQSS6yCSj_hyAzg<-5z4MIEHYFd6 zS`M7X3(v=>BKQn%)gW!w)*7?e_@r~Muk#kXxXh-ms~hAeYLZqgAuev*U!RHA*w88b z+@a6^&p-dvSb{aYZQ&eUC<&fmO9(ujHCJPBS9Z8qw5%0Czr@cZ*pHXc+K0jhwr|&H z8kUDDaPGT!HvcY92ryUM@2gY=Y+P^|xt|y1>?ity5z>^IFOhrSsK_S*c!e27pNb-I zY<Hc`#g~h;a$M|7<)(Y%PCwh?xzY9rLGxA(S>Ahye#3K0i-xc015cTI#l?x50CcqL zO<(^KJ*V7?6`x=AiW%NLp?G#l%U}3H?#W}fveohbE)oClRGv_jCG4k9uMbmqKF@C- zI2v^6vF-nra&)OQV)=`Ab-PEN`w)1v!SjI>4$ETw|9(wyx!jxbNfo>AS|)<@;VtkN zGbjH5x00^2NLD#~>(~3B1}me0%R^-aIv0spH)dy#*P8DxLhWSmg2C7_S|NUv%}z9h z5YB|=SNJhPiT07o?=2!iA}}1cJSd;m<T0($qp9z1Pa4BU@HbNMBbEIas0pIr%}NaA zASp3i!nM8F&+F=PsVA%NlS06s5W23SHr0G2jAvHBB-q&)=|yl?+?jMidwcucbxXNr z_7Qg0Ta(;0vG3otE(Jg;cdcjV52?z}j`C8$t_q8aJ|Q(V>VQ8Za&zb0?~0y#^*gv_ z+Ac9Mu>cWq5+#f)ddJiwl?w>tLlC-$x7-50M5tJ3`RvmAsQMFGnOd-dP$)@B=gBXg zpq>*cw@ztln)KDCk`IiH?F86Fr}Ahrum`TxeeTc@?u_?McIqftza}ja`JWOe&usPL z6g5^R{7S{8k;qJHA#y}i(#|JD56o!WJpMeBaq-+*A-zB}aut6xr%NhoxDIQAHA{CM z{h$bn#5ovtM{0+f9Mf=)8N!l;T2jL&*>Qc@1`(n!KqCs0W|b{l?(OYOh>3_p?`N5* z&M-fDx<~20mRq2ME?1?|OxQz1)4veqibz%bRFA1oA;!47;#F!g_X&jj6sc0j+u??r zPIeml8tB~aH{b{5q{kMN)%eYg6_QYIPUz}3bc;J-FVl2_u}lD2?cgSyy64DSArIjK zeR?CQM7YjOeTSScdi3p5`<|-ZUVJKr;Z!g2f=uj`0GDQY-BZwiLsWR?SRy?L8`U}z z!g6C6c72TK0s#8s6JA$-U!c=6qut+p_%+Ym9>t9|2UR889zbK0qcZg8hJAOZ=aMGr zdhzD;^aa2eetBlWF1_r11o;CQJom}OOQliX?0=fS9LoKioPR(-GL)V0+}zw#U^Jv{ z07!2?PGVYW%c5C#LU67wURJolI{SFH0DkMXlq_SaJ9dP?jJ(-CYbCBlIH%>MqMere z`cgg(SDmpo3UME&Iy*aqHa41hKswuY6uGruF${Y?)srY>{{yZh$G3V^90wLtu>wC$ zJ{(*+#ro3bt6zF(oi3Cu$!$eGyOkI+5F(xDm>o)iOYtdp4q}_*w33sm^<0fmhoPpq z67I=cfY9`X))fAx3m5A=f=lZ-eYuwT_0OFmK(-HjebaAoY!KfCsc-wa`MX=cC}-J{ zwY%O!mYfV8l5?MG1w0i}*n+o9?ZDq;UB`>fp*sZ+As_jqd!HWoX(c&^<J&tee?`;G zMXQd?&6P7UGOBrc?!erLJWG^NQCh8Ol(!T&F2i->O5=5+7oi=hB4r6zg5cH)*{_Iz zb3oxI6+>jX`)8KBY(UP^igu^B0;-o@S+s=QmRq-aNMWeL)uUX?$QB<~U`eB4hDSz5 zR*=5e`J(=7<MkC=54ZZ8Cpec)5~wg-!0)?>Y$%sEJ{?F~hOpj7Q4eHK2Wb93+#x{g zQD<kSYyNY)E=GQS{z_iZ$k4D82uEyh8wnuS{DG{Vgr%G!;XnLAKJZC4Z^E*ccX!LE z{6P<V^*to#Yg%!jq$oifG;llLeMBhX`U4*hp;wS72?<%<F*A@~m4JVeAwhOPbuPt! zdfg;->^A+3bAM=Xr$G4KM)mGr$*&;dvZRS1m6$zz6y>1atY$S4hi9N>84w+tXx2EJ zLg=BDB2E7kV%c~Z;<(O^951gM#Fn;~AK)MWO!VxFOG`KI+r!4byjWGbk09c-X(sX~ zuT?aJshowa00!D4ls)*{Twl}NZpb>iMdGXYRf}&+G(Ifw;x|Fqbpf&yjS%W%iYqy4 zah~*y^e;qr2fH7E0=|RbJuR&oKG`OsR|D#x@TVH}(pfhz5nY8oJC@n`_$I;n{8S=5 z0~irAp$H1OgsB|4Bo~6POm5nuy(}c;<yCfTpOIOtvKZWp9I^@IlbUm$T*J4m-L~pa z*H$xlZ?4^a6ty3`Lwg(m={rTTAYJ$DlJ>O|;y6d%#JssFW9cTIxj|1t%ZR4U4SC#c zwnx8=$aMs4QBUaWhww(n#F(Y(qk$CFKS$H}OWpWhhiqRiTG#sQ9d{9G!?S{${V2gn z=znc7u?dAya%0hf?&CGg+Ui^%Y#P1OK=t9C%_z#dm5>S2e5BRp;6hJgV~YYMX~Q*T z9^$Ax1z?1}eQ9-q$Ucj)uW92S)2*9G_*H%eNo*X>pZq(kviv8k%ApK-@<OqLb5xuZ z{eUqyB7G)Z{PIe}fHm913R=Lnn3vhJqy@T&C`v&uIx&@mJ|VXHtf^`Vmts)ANb~B1 z0P%e|1Rg)c@oNrwx8&df0su(bl8j~za<J|ZWs?a1q^a>zq5b9G;`G{?kg<gnOn^pF zkD1BQo@ycfz>BYpOW+FbpFUB}S`HAaL0Vj@d3D<QCMa5UV>0aAPG|dbCK^ctMn}_H zynN$B0@<U2UCfHa*XdkT&`8t^A_aV>MjQmrjSNUhQD@Sv`us`pEm1;QB$`T~O+Kgd z;7fK1DX3#p4)Ps1p8x>Sarcg*!a`ec@wTKQLGuY3TE#Iby_MD1A*&JrolpxR^5L=_ zDBvY1I>SZ!b4lOp6KUo!PZ;_$U)M=Ok^hJePILz1hQSCB0k|N<BjV1FxJ&jgKf7!o z>qSxJrIOw=!z`o23PLY*p6ZHH;{0CCtY0JYtA>+Jp!cbHOupuscOO~{yguR7S4Lrt zp<8Wkx^Bn?2WiKfl^8m<1>2a=-RD)Fwzf7uaZ7E}*Q6{zLNHCfqtCqaksruo%`u5{ z(ZrgYASRsL&>w0K_7K9&b^FTP0*jonDx#Ag<-#YxQpLSNZV<JsO(BI0U8e}4tX(2M z{C?{{^sRuV8DpX@8r+fObM4PG0YW><^%J?kQCvXRvB5iJZZwMdYF+NJ`9^7X(Ze`? z8P<qR*ZuZUcNZX~o2IJLZC@{w2#48PaU5&vY_9)0C&(XUcOrI7k=vP<LvHdcarl6V zA}Z@)DkMvWkvqJ+y>s9MN{ZR+n^kJNj<3X9<y2RXI{Hc3n?-NMPlLj(-&iySAA_%X z21FEHD3pZFXcParT#xh&&??!R>x4?dTFG8R^7fV@7vy??cyF3%mzW14?h_5^e_vwb zb$|c$d6x<0G=B80kC$Crr}vN){)jt7fMoALjd5cZ>6a1`1r<k<pbb>ikN@K^fa%1@ zEK#R$h}yf+o~``!lZ{n~MFFGKB)~GQ;KqurAHfv+37`5r130nk^wUlJuy$SZr{OJC z>(kx60%^#s_W%i1N!YSHQ$K0e%m|G&JJEq+RZ3tA?mf`%BQ<~*P?i$l>-#445fpRh zzSL^s`)IV3$&4sFz+_OdEiDa_j>^iKl5GJ1FDBl@9{wF~|27NI5#(t6`RO?-19@kR zo?3>QX<M$qZpo0rg)O)fLetMdK4tbh<^@1#CKd)dGqa2}Q|2~Z*Hj97G>uxyL?4VK zT$0FBXq{pR0aZqEh`c<MQzczp-J&!CRakSgL|Ri}%$F8LME5C8U2n9YsK`h@qIBKc ze1JG3ZOvCX2w1D=@7MtIkt7ExACUo8g^v&2w1t$_ULw(?{10?w8eQ%i9jSp2g1Y;i zIwNj?&^mn1SQS8-3|zEitzl)%DP~2KoHz7iOU4KnjLqG2NQOD0cAe#<@(@Zw)JzWh zVn}nXAC-hG6K{~c*z)sK(<3aD@SNFfg5uzK`V-Tg6f5F{!a7aK`m?dT{OJdiDwZ*1 zIZkN@iy{|@S5!H|<t=G9T`$LLr|uwTJCWwV`{dyp7Z5qFeVw4r!?!orMWv;UyWjnQ zDh?kcMH_Du;MRnCS=Th@0ADAAEK<o=uw<LB4g-GGj;KJ}qtM1^Raz!Sy`P_TRmN(f za@zK1c0{J7RYCR!?nL75CkkAkZQ*ujh(1a1vpVh_SYCkRCjaKA68Dy1a_DnM)Ya7q z2;rG&X5XfYlDglOsHb>^87&vB!EdkW2|gA76h8u6pJobl4h-wF#`NhEY1qaIK4qlr zi+6uqnG+-iHxdD!RrUvFUI}UeaX|wl=YgCzUI1XR;~QkZj(S6;A->m}@K6?)9$+@% z8R(goVsqd<&koK*#S6%pnJ*l40<a(AL{nS@b6_n|Wv-%x?0?>MvYT>f9NLm7L1IZz zMPZ$GFgrrZTJFjRvn}7Z3AuGTx@ErdWuPh%KNe0v0Y9oa*?}M_`$EzP$?3yWED|k~ zQST-ZK;**obGU<fPmMQZ@(n-9LteW3{rXqemVEC%`D;eSdqE$ZaR7hwV1e?K!CFW^ zcY|rW%bp~c)Gw%}L^>lKjW@G5sfbX(c{n-Hm_w`9b@cuV@<N$QpK(S9FBqV+rg!P* z%FN@(0CKvavQCS0mJxE_=r&TKVUXi8R!9mZ7M_d&4>W=s*s^+;o_J?WCNzo8!4dC3 z_uQ%QjgzEjSRX`)g0c}!k+_Rcn$7*r<axMr!+9s>((i<_cVg#ncVfG6rYAHse9)T4 zWJf?#h}EQYZ=OaSRcN}Q$m71Zq@;vX+~Tyy+(<~2W}@D18u}GEU-U!GSHF|Z-23y? zVIpKu85nzrKwbq&IiKV~+fyN!<I(~);sp~|l}MKiUf14dilvyyfYjjgeOr`XTRL?| zM<tDvN_T)v6b-i(CA%?)b)d=D0CAEz2jA~}X(uQoO6WX$POrCBl5(8u(XWJd;zg`M zA|P<nXNFetgFc4BcBmVyP*+Rc(LwzPTal<@1ezRuIDd!vyZbGlz!?Af_GXi4E);h% z#j@=P-nm~f0__Xn#<o<dCvFC{umCy*5#{Uct96;(ntA?UpAp}Y59yZl*=joRcEJ}? z5Ngs?24wNp9(|TeLhe-ecc22EaZFDyc!O2OW+;yD$)pc2xxZw@weP-f;Nwhv-Nkl4 z6H*q7fwF>f#%p(C8#L1zU3%sP_F(6BKG_M(U<7>=?jwyXAoOy|y4VrcW0_XNM-u<Q zf!BN`IO!Mdq!AOi27`Szc!)Uk=?xBbeAqGPk$If}sqPp%L}XCn)%Sy85W9lZB`~;@ z0$UvMQ>H-Yh7X@VKO{yd2EeRyna1lkY&hdHZ%8y>X8#za7NTO?e=1g^4f|cPJH4I( zh@N2N2G77x#OP=OZeuC9_4Gg<JyJEuA6;?_z;z9;7Owu;P3CIu3Y1|o@PUX-!f;5g zNiHe@N^%#U<V5F(Fd)Zw;qaS?OHnjA!yHM5lsRv`F8fRNqLFqt>O`b32*1Put>-2B z$s`GPC$L+R70_$jFqPqqS(!=^O&kXRKXH5eJQ`6rcdE(H+F_JoBAB?b=r~bm#_{>F zF)<)KtlW2jd(OS$oAw2$qonXGa^JdKerAxj0KM#E9-W1rYB@|3#2^>JN5n28eq+NO zX{-7S%tE9Y<v{pG)9G_?!~1JJP0H)`fNPeVDD}!~e3c6qwioBq1#?x(u|H5N*|-e` zlQQU{`l(ai=px2sjAh*X!i6K4f|=${M?H;x^5nuddkXU48_Xxv0EC2AekR5a+HuG8 zrdccO1(Zl%f?P>x%uj5PwvKtG?OT9v%kRKq`V3l?45=h-eGix`La9N-<Dd4Sad-V+ z099W298hSwqT|?L-&GN<PV`Yo4oi9DWaKAbavbGCkG$j%8p;aN1C1rri5MUlt9h8~ zK5aEjlLLB<uX>JwIjR?*(Z_IVz)U^1#6l09`3pH6l6So*1zPD=x@*x*@v=6|3klfo z#ml6Aw6<xwfnZD3&O^Z>4n}(ce+ZOx884!}0@AQ_9uv37?jUMYF9G}LIv<i6Q3Q5j zSk5Tt(WA|%V-Yuq=EUB%cH*`Zun19nl>@}*#kXY=_4hD$RiuU5386NK^b-hvJMERp z0)D~KJ3z&m0i96vz_eTk=5^-DF+mYRIin|tw@C01x@^Q9gv);RuYfS5Yx_x9>ji@z zkSQ1<3ynknXxm8ktK;XXX_RU9VVs&EAv|*+q_RXwChiC^SP1b02h;o3#Wp2`TE9up z9g8QNj$C%GD`ajvG#6(QmNK)WpSr5z&wyvEwSkrH`fKR~hk;LS!s%&}xgFF7S?USo z0!k7wbkv9*voea*(5=VIB_3lOLo4ONaRv(DwgKmMH%uQ9^soEp1+vG-u*b@v?igUw zMpy6!8I~c7h~A3aD}U)stv?*-(j7fd%H+}SiJ;&z-mBGpguMGAZn_f+CPGekUfOX} z&Uj03#eeDiXwH5zOV|tkjEp&Q=+GfEbMp~^)3oA5tz<q-5vrn7Cu8E$5W=nLg$n}T ztnZsJn*h`dND(I>SWR~9B$%HM0Ek?fJoMLpsrl~{!!7hZn%DdZ)`RA+=WIk40ffj; zuKB;K`D057gfK$ABI6eSQ2Q*%3>lfEAQs!X6ZpW|tFYLT&%V3@5gYEiSdDRjr7fTw znHM9Y7wS0UU|H9++T56!9q`HrT(gn3d&R_3<6v7^CH3MOZcZpA#h~3jFmuI03`?r> z=EjI@4#n%0r(EOjKdmApGYL9liUC2jwb-UT-)U?TYi`;E_8Xb)ndh(=8$q2SQzF`y zP37hHpw3LlHez-&2^O%A5h5yJYMz>IhOLs0>WtMGfgq{*3hIiWNBCr3Te8l0=<8eD zynE9u47J+fWJ~%Pe69d@8SFS*$`-Xf3Wge)Z=it)wG+*Wc;4tHf%a8e1iX{kNH}qU z$UIzWY_TKZ<@@0Z!aZw)`BMav=7ZOh7@7CjD~>B}`_I9&Vc~=AmZgsWAla)EEGE9L z?7!W0{y#LoT`rR(CNfhR00OWgKUAsTP<%g|1&|ra(asOVX-Gp%-TGB{3)F*cM}edE zyB;rJVhE{`KoDYC0fQNm0SDq|fWf3&(4ftC%uNrJkzr4CH6y4Nl3n3mM<1@o8H|_D zo6{H;uuheOz(PKFq@%11!le)3hB5Xx6M#w-lSRA*Xf?^;B4ij)vNQDR(G?p=-E4+O z9<P%rLj8!zBHn_3l_Au_Bvy$Dor*9cGY|99ADSK?wY#V1@41ooEjKpEs)9a2mesDp z7?^}EsdaLD!QafV6U;tqdH^wdwsMhZ|GU||eKTIwM};9qE0yHWo2pD4?hLzY46oXP zGZE7oFz}LrhBfLW_U;YO%5n*Nw;m%E8VNczNf-b(^fD&X7O%0R(0(Mzj1rOjWFRvU zD>#<<k1R8v^u=#1(;@6#FRl!-I4|I=zH5Rb!79ndjoj3q*#d#hvWV6O-bTi}{}Yax zb`kqYlwT@i>z#d%D<`q*G>CHll3E5xZYOuGf@Fel&7uVWn~|mtdil|q;OQ6bk3-qY zA2o>NnUlG-Oq(Oot0;puZQE?{zQ2!3{!f|K_-A?Ja(RC05a#{IMEyJfJ)*>$i~Vz7 z;n_h`biNB}*hT``Nhc>P`vdcjRs8caL(nfp9sBEXLR-nWI?QHq+I;yw({X~mnZ6NH zh9tyk38`uT*8gZOqqY4O5sjs3i`<ty`-j<!I2igE!z&?Puwx!6@#5lX-vkJ}^7!Ah zdoevS=aeo?<{!ZRiZH95Or?SK4}KBo<bvw)SCNjmL(cKAK<Akof7jQCV2S%=-tPK0 zp?4+s+OLS!@AA3Nv+vS5M!L5t6SU<zxOQTpR;t-4VEYK4F;K=f^a_<ky&M9e|AR@S zm65~|Am;lpI7Zw1t#>W7cu?3GhM4sxD(}CTK~hV%YB;(E?CO#Onci5TxbZ%sbRtLy zvlh@@ePC`{h(_Y3&(g>gI2qdgpB#zbSD3rVAL0$cpx_w0$ywtegKv<^4WM^zJG?n+ z{vECWIG((O)S&x^5PS~tdgnZfmRpEWph0JBH`kxmCeAl8f5@~Y;hpf)=kziEj96HK zL%&Dp$J`8)4A*h?pdcU!M4$k!<O6|d4(rzGrp(L|tRJVJCkI1lo-z$jH?cgzrqR$N zFcW9w4l+1zrKqeMWjj4R9uRTo#o9;#R*g~7mwz{9GP7F|(<y&=4NCm%!Y%>pKnKAx zz?m2wy<K7bWT(fZT7)Xa&Mw#Fl1luACQ-siQPh;Uu41VV)ZD(VB4RcXXU(#K^eZbX zH#C_1O~jibE)>*-8CG8dUAeN7=iSh{EqSdmn*`UiDjrci{>!36TYZapCJml+{;=%r z!*dbgf?~GYA8uEFc|V}0m}f{)DsZz_*xpU{7v_FkW29B8yiNU+cY7I=t#smUZoNI5 z|N7=iS)MwR=M#Y=k{xcP4rR01105rh)kQ^0%a<>|U8^sX1w>R(U9E`;u$$VcW@8!8 z*_yY%kzorJ7Jetgc}qCUV>T;0Cx^#TbCCW+r`j(jCMMTBJ!8Ye4+5$BjB38nJ#peO z;Fr{djg8Iyn>TN&_<MRPhE1$l76w-27QNuY$zoB;`c<@Bvn_?0?wj}1Uu84f1N>#2 zleTQxQW2%blLMHR<vbF}@bcnutVYmND69nwrdIK7-mKs&q_5%Pk|n5@{p9J>J(f+m zIaRMnp{88oUxJqi04+A2sU5PrXk&By@ZrPp?^aPlAMGilTu?3+ei3G~CoY<Uy!wSf zb*~Lj7Ww*xfFWeqx*y53ECdC%A3Di*%I5_O?LwDBL`{OSuDjE{IziX=;>B{{2K%pX zembP67>jukTYLLcXU~Sh<xF&)=+`O|(#be>?%X+X`v*Qgg|A;907LpHUwxstsL1Q; z)dNA664x;oY6i%QCr%hwuU>_19uTfeN%OZY)q3s_kH;c})YR0RU0lpCV<qJ@7+NBS zgR!@Q{BT>8Lh-LaR~4hRD=%8KsGzd)6b4z;Po8`rH#2wuxJHSw`j=nM&i1BX`0(52 z&6~ART@3mq2S;11t1c<Fwzh_chw}oeD{*skOSz0bIdS5|NnPFe;^JbdiS4^~RiLJh zLbcl|E8D;~04m@wYH|k)YFgdk%iECctPYHQ<!eieH-qufuy3K0(=xFRw00#_=yFIj zLYkVIc-l+(!1bjcpw<ZS@bFMS0v}s+GrMCG_Zk(hz73thhGXUauxSXT<VgVz4h|z@ z;~l)b2L%NMW1c)Ihp5Mk_ubwod4iF0;lq}Q@bFI1#M&gIn^<gXOG_Bgk!*VyC0I{K zft~v}b{<hu8K%g=oSrs^g>Y0^`J}0-X;7PlWo<aj!3&t)NWC%w)U=I{uNxo_hjljI zoo(U56sMu{+qP{hIRET@Q<HuOeJ)c_KR51p*kSxeVDH|O=sDukeKpC7m|OA<2ncX< zcNa6rUy6PB<S%dfkvAeL${a_uY|j?Tdi4V(k{cwfD)4N{>a}YlA+KG)X~ZNXXi%wt zq$*F3_gSV{)}?KilBzZKII5(i4&9{k@u4M)7B4n?d1aBf=E$}#Oj#ZS&ui=Mj=^hp z2nr@#R$TPrwZFKfzrTca0-EOAckg16oIfvJil!2F^7QE^pFf{AukqP^CI)sretl9~ z+rKu&H0<%?Edf)*BO{gZXCyK+G6Xhl+7I7#vU>csuy)49M>8`stzW-xv9`7@c_{qa zpWac*RpKYAC?_Xp>*%OyXc!I}W2CNVWtD(eR$!TL-@Y}4yfZvFI5}{YZ9Xk7&p`LA z&2@jShnYTYXCAk~%Qy&N*ToQ7s|xi~F5HJ0$f@Ed645wCv$b8lW=+Ia#XA?kV9^^G z_}Xjp^&pp=+i3{lI9kS78WYp*rg%_cy|hEz9WLp}k^l*FJV<7r#?f3bqTe2vM<f|c zPDF`Y)>injpbJDlefqR&rXtBmsja=e8)l$(rfmc)lGwn&jpiQb&OQ42_3K`@$%|Z) zwg++4-Jhb$wr<>h>H09Fz8*lAD&YCH{(e0ZPIh)rZBvxu!m_ez0Q=V7UJDg@yUu%v zb}^3`r&A|SUik9rn)cJPmu+dZp1uPf9=W!*wkgIX>w{#RPhecg@3-H6i({gb)TUd9 zAhqp!KMUh=2URq`N(>Neqp~mGSPon4fxmyu@USh4hdF5D-F;gwafwjO{80RC?CkdT z_xFnj#Ky+nH%!pUJiogxQD460>(}$Ve(=rw5yR$shYw%3ac~f46oYy9Sj48JR7_4z zJ}mMTI#b6jX7)gZlTN3LVrP4%m#h=~m1>OWs!Bn0a8lc09_0T@nEiREoZ7MR;NX6E zHl{8-M9#hTh=}5yjm{6xEiB)UO2mwULD=xa;>A_|>&zfZ6+5}OR0Zg71uFU7)3a_= zR9JZE@)RZ1w6`g*>4*q)q*Us}@NiP8G^PW5gxeEo9B2zI{ry%5T`OJypxTFM-r5tc z81OD`etFfh+h5u-;iu)x6Vf>lntOtz9mO$h74Pi`?WLu&^Mc3hm}bBzrd#~+3;?IW zEx$eIuHdPg3>A;xe!F7YkTu43i*lBgab<bMRw*fM@O$wyx;Mit8nd-LIe4m|A(8aJ zDF)0Slw#lwfIY{I>2zSER}9G;ulrmB2{~Im_xtth+WJr{RT&l@HVITROT^m+QMa(9 z#QW{}XI8?(!pP$)xD!0*c`19ZMLwj>nSvyuZmTOOT-MRm?b(%#ut#~RDs0MjPJaCO zaYe4XtJbcyB_Z<qU^P`UIPhf~_YEp=@Vqdv8*Ena_Vta$o@LQG?ufZ#9%qM;P&EvZ zu=nqew6(Sh2nf9QXF;eE*s-H9GU3S+0TN@tErrF!DUuwZbEHt)O^;imeC*h^?GRc) z7M+#L8~yn4Z+KaEeP(ir$8*PlDnw13CAwg3N)3oX{`&Qy1<#IbJjUVZ=olchnF=9S zC*BsbbByrj=CgRGX?I1mKNOgaTwL2`?KfPov3a(bgF6n#(N~|T&Lw4cZ?xL769o5} zb7O2)mNugz+nE-$&*nGeFk`fO;c+oDEBV-+@y09^!%o-$sm=+p*F!>@UN2pAiLZP$ zt+dI=#KauChw_%jTrfC*T1a}SWn*>?pF1A#sInJ5fBt+0j3FY~$hda5w~t?BqlA7S z?e67E1%B1Qr_@9({jVV*b)-rP2pmyRP$2CA`L}w-ibu`O&Fp3a7`9M%bITsfNl#B7 zn~BWGI1N>eR8rUyDXtS%4H-5E?%cT}g$%b>jXn45G<Vy*%lyk1NR?j#H?YSLJcdN; z+@2xqWAJ9u(O?y&&^lC|U0q{gKcH7VjWz)0Bxnc>*;IA~eKj1hk3%i*%!3BmoLppj z1Jw=Ce6*7c|G*IEqZ>DFxb#5E&w2IwHR)FYqwqX5=H}*Hpb5p9Ldj3vzH0sY7zShI zHhz9q_3IzX_h)Bk?{ykH3oG}ib@tHVoU%vNZ-LBy`spX6RWR-!+CXh;?0A2Jxw;j4 z2RV71T_wU_syh85jizp45qI9eK+Jv0u7`OH#T17OwOz1sI#C?hr4spN@o{{hiS^>$ zg3r<m?V4BMDUai-s;X6a{1`a?1RLhTy?bYAM$fN%dU`^kTD@vj7(9s!fI`FL;~G%n z&Aq)&?wmar5fKq>-IRM&MWu3QV)hn97?-?<4t5<7(wp};{Y{vaV~ZjGUx2Q=@^i<G zj3P=y<b|kJ`W5CHR*TQduQw~zO0)3$_~}!s=_kHLKm5=M=Zj&O<DHTZbe{qpOiI(x z17_z*JfoPMEi=~NaOT>zYY2d++r94Je}p@ULH}4U?~&c#l$Qw5b^<V9nT!a9J0H(l z0c=Fvxbfp|fU`%B9+_3e`GF!!&rY`7{h5O!(Xk1`@iBmJ!0GA-5B_cpTLn@?59AGX zl$-Isx+p+bQy{0Ys`~)}FkE&VI(#?=^l%FDertO>AF3{#*X_c>$Ej4R_OtWXpdm-2 zDic_Y4_q&N{&zgu^a#035&-SH_L+fansrm-uLl;%aNnKY?DT@WGWA%1aDH!n=BKQ) zMDKyc?)&WYuKb3Eh9B731B#241Fk&~5fM?o%PnisJqg4`9>7&TcI=)c_&Qr^-@t%1 znv%swUcW0>7WPbMt#y$Ge&m#KdYXC{i!g?qSv7;*|9m7S)yU9L=JWmq7k*KZ?3kRM zmM`z{fFzp_LscioCACzpEJ(WA$=R9n=G*ftSFTh>NW>rV2c%FreflA3Kq#I-L${pY ztXNon4@YSYdR~G7FY~~_z*}$63u72g7?8Cc3V!~>tTNeK!hC#u`LNlFr{*xdBlhy; z%lv|Zf=MRRFD@>QZuvYk^u@u$hy5{R%8B5?RFA&8j*dl1p6f&cRZUG|!GZ<k^_ZS` z`R46gKRiXKRRlWMh(7xvBQw(vNXRQFsAP+cli9vsM_4RFaAt4<Wq=s;xd>=)a{UPu z7H4e6H?YtCcd%dom1ruIOI*$KPsLIw4wooI;-^rg_h0({r%n6G01V~D^1pF~zux{2 nT!FlV(o0!DSn&VbmtnDKtP$H_o_$Z2+zIuN>cNEH&R+Rn#&aY) literal 0 HcmV?d00001 diff --git a/intervene_test/enhancer_percentages.py b/intervene_test/enhancer_percentages.py new file mode 100644 index 0000000..c4b1048 --- /dev/null +++ b/intervene_test/enhancer_percentages.py @@ -0,0 +1,62 @@ +import pandas as pd +import matplotlib.pyplot as plt +import numpy as np + + +raw_data = {'Cell': ['H3K4me1', 'H3K27ac', 'GRO-seq'], + '1 method': [float((153327)),float((12510)),float((3754))], + '2 method': [float((55923)),float((55820)),float((611))], + '3 method': [float((609)),float((609)),float((609))]} + +df = pd.DataFrame(raw_data, columns = ['Cell', '1 method', '2 method', '3 method']) + +# Create a figure with a single subplot +f, ax = plt.subplots(1, figsize=(10,5)) + +# Set bar width at 1 +bar_width = 0.5 + +# positions of the left bar-boundaries +bar_l = [i for i in range(len(df['1 method']))] + +# positions of the x-axis ticks (center of the bars as bar labels) +tick_pos = [i+(bar_width/2) for i in bar_l] + +# Create the total enhancers +totals = [i+j+k for i,j,k in zip(df['1 method'], df['2 method'], df['3 method'])] + +# Create the percentage of the total unmarked enhancers value for each cell was +meth_1 = [i / float(j) * 100 for i,j in zip(df['1 method'], totals)] + +# Create the percentage of the total H3K4me1 alone enhancers value for each cell +meth_2 = [i / float(j) * 100 for i,j in zip(df['2 method'], totals)] + +# Create the percentage of the total H3K4me1+, H3K27ac+ enhancers value for each cell +meth_3 = [i / float(j) * 100 for i,j in zip(df['3 method'], totals)] + + +N = 3 +ind = np.arange(N) # the x locations for the groups +width = 0.5 # the width of the bars: can also be len(x) sequence + +p1 = plt.bar(ind, meth_1, width, color='#6DA2DB') +p2 = plt.bar(ind, meth_2, width, + bottom=meth_1,color='#D2D5D4') +p3 = plt.bar(ind, meth_3, width, + bottom=[i+j for i,j in zip(meth_1, meth_2)],color='#D52114') + +plt.ylabel('% total enhancers') +plt.title('Stage') +plt.xticks(ind, ('H3K4me1', 'H3K27ac', 'GRo-seq')) +plt.yticks(np.arange(0, 110, 10)) + +plt.savefig('Enhancer_percentages.png') +plt.clf() + +H3K27ac_enhancers=12510, +H3K27ac_enhancers&H3K4me1_enhancers=55566, +H3K4me1_enhancers=153327, +GRO-seq_enhancers&H3K27ac_enhancers&H3K4me1_enhancers=609, +GRO-seq_enhancers&H3K27ac_enhancers=254, +GRO-seq_enhancers=3754, +GRO-seq_enhancers&H3K4me1_enhancers=357 -- GitLab