ó ¦s§]c@s ddlZddlZddlmZddlTddlTdefd„ƒYZdefd„ƒYZ defd „ƒYZ ddddd d „Z ddd „Z ddd „Zd„Zd„Zd„Zd„Zd„Zd„Zdddd„Zd„ZdS(iÿÿÿÿN(tTracker(t*tChromParameterscBs;eZdZd„Zd„Zd„Zd„Zd„ZRS(s9Basic information on chromosome, inferred from input filecCsg|dk rt|ƒ|_n|dk r<t|ƒ|_n|dk rZt|ƒ|_n||_dS(N(tNonetinttminPostmaxPostrestname(tselfRRRR((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt__init__ s   cCst|j|j|jƒdS(sNumber of possible locii(RRRR(R ((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt getLengthscCsE||jks%||j|jkr)dSt||j|jƒSdS(s`Converts genomic coordinate into absolute index. Absolute indexing includes empty (zero) points.N(RRRRR(R tgenCoord((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pytgetAbsoluteIndexs%cCs|j|j|S(s/Converts absolute index into genomic coordinate(RR(R t abs_index((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt getGenCoordscCsE|j|}|j||}|j||}t||||jƒS(s*Creates low-res version of this chromosome(RRRRR(R tresRatiotlowRest lowMinPost lowMaxPos((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt reduceRes#s (t__name__t __module__t__doc__R R R RR(((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyR s    t StructurecBs•eZdZd„Zd„Zd„Zd„Zd„Zd„Zd„Z d„Z d „Z d „Z d „Z d „Zd „Zd„Zd„ZRS(sBIntrachromosomal structure of points or substructures in 3-D spacecCsk||_t|ƒdks'|dkr3g|_n |j|ƒ||_|dk rgt|ƒ|_ndS(Ni(tpointstlenRt structurest setstructurestchromRtoffset(R RRRR((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyR ,s     cCs g|jƒD]}|j^q S(N(t getPointstpos(R tpoint((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt getCoords6scCsAx:t||jƒƒD]#\}}||j||j_qWdS(N(tziptnonzero_abs_indicesRRR (R tcoordstcoordR((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt setCoords9s"cCs)tjg|jƒD]}|j^qƒS(s)Absolute indices for all non-zero points.(tnptarrayRtabsolute_index(R R!((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyR$=scCs$|jƒt|jj|jjƒS(sVNonzero bin numbers with indexing relative to chromosome position 0 (not chrom.minPos)(R$RRRR(R ((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pytnonzero_bins_whole_chromAscCs!|jtj|jdkƒdS(sAll non-zero pointsi(RR(twhere(R ((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyREscCsŒ|j||d!}|jj|ƒ|j_|jj|ƒ|j_x.tj|dkƒdD]}|||_q^W||_|jƒdS(sLSet structure's points to only include start_abs_index through end_abs_indexiiN( RRRRRR(R,R*tset_rel_indices(R tstart_abs_indext end_abs_indexRR((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pytsubsamplePointsIs  cCs)g|jƒD]}|jj|ƒ^q S(s)Non-zero genomic coordinates of structure(R$RR(R R((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt getGenCoordsTscCs|jj|ƒ}|dkr"dS||j8}|dkrw|t|jƒkrw|j|}|dkrmdS|jSndSdS(s0Converts genomic coordinate into relative index.iN(RR RRRRtrelative_index(R R RR!((sm/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2019-09-27-figures/multi_mds/miniMDS-master/data_tools.pyt get_rel_indexXs  !   cCs˜||_tjtg|D]}t|jƒƒ^qƒddtjƒ|_xD|jD]9}x0|jD]%}|dkrg||j|js$    !t %#