#!/usr/bin/env python # -*- coding: utf-8 -*- import pymisca.header as pyhead execfile(pyhead.base__file('headers/header__import.py')) figs = pyutil.collections.OrderedDict() keyDF = pyutil.readData('/home/feng/meta/key_ath.csv') HOST = pyhead.base__check('HOST') gconf = pyutil.readBaseFile( 'results/0219__prepareGCONF__Ath/gconf.npy').tolist() #### Load Tracks chipTks = pyutil.readBaseFile( 'results/0224__chipTracks__CAMICE/tracks.npy').tolist() # url = 'http://172.26.114.34:81/static/figures/0122__prepareRNA__CAMICE/tracks.npy' url = 'results/0224__prepareRNA__CAMICE/tracks.npy' tks = rnaTks = rnaTracks = pyutil.readBaseFile(url).tolist() for v in vars(rnaTks).values(): v.height = 4. panel_kw = dict(figsize=[20, 20], show_axa=1) vdf = pd.concat([ tks.cold_wt, tks.cold_mut, tks.chx_wt, tks.chx_mut,
#!/usr/bin/env python # -*- coding: utf-8 -*- import pymisca.header as pyheader pyheader.base__check() execfile(pyheader.base__file('headers/header__import.py')) figs = pyutil.collections.OrderedDict() chipTks = pyutil.readBaseFile( 'results/0224__chipTracks__CAMICE/tracks.npy').tolist() url = 'results/0224__prepareRNA__CAMICE/tracks.npy' tks = rnaTks = rnaTracks = pyutil.readBaseFile(url).tolist() bkdIndex = rnaTks.rnaDF.index # chipTks['198C-CAMTA1-camta123-CHX-R2'] def plotter(xlab, ylab, ax=None): # ax = plt.gca() # if ax is None # ax = plt.gca() # xlab = 'CAMTA-CHX194' # ylab = 'ICE1-CHX' left = chipTks[xlab] right = chipTks[ylab] indDF, ax = pyvis.qc_index(left.index, right.index, silent=0, ax=ax, xlab=xlab, ylab=ylab) ALIAS = '%s_%s' % (xlab, ylab) figs[ALIAS] = ax.figure
#!/usr/bin/env python2 # -*- coding: utf-8 -*- import pymisca.header as pyhead pyhead.base__check() execfile(pyhead.base__file('headers/header__import.py')) figs = pyutil.collections.OrderedDict() rnaDF = pyext.readBaseFile('results/0427-dump-rna/pk/read-topped-log2.pk') rnaDF = pyext.readBaseFile('results/0427-dump-rna/pk/read-topped-log2.pk') rnaCurr = pyext.readBaseFile('results/0427-dump-rna/rnaCurr.csv') clu = pyext.readBaseFile('results/0408-freezerCluster/0427-dump-rnapkread-topped-log2-meannormpk/baseDist-vmfDistribution_seed-0/clu.csv') chipTarg = pyext.readBaseFile('results/0427-heatmap-brachy/withScore__chipTarg__withCHIP__allTracks.csv') rnaDF.columns = rnaDF.columns.str.split('_').str.get(-1) rnaCurr = rnaCurr.reindex(rnaDF.columns) rnaCurr.columns = rnaCurr.columns.str.upper() qc__subsets = pyvis.qc__subsets mcurr = rnaCurr.copy() _dict = pyext.collections.OrderedDict # COLS = None # COLS = ['AGE','GTYPE','LIGHT'] rnaDFC =rnaDF.copy() res = pyvis.qc_index(chipTarg.index, rnaDFC.index, xlab='ELF3-differentially-bound', ylab='significantly-expressed',
#!/usr/bin/env python # -*- coding: utf-8 -*- # execfile("/home/feng/meta/header__script2figure.py") NCORE = 1 import pymisca.header as pyheader pyheader.base__check() execfile(pyheader.base__file('headers/header__import.py')) HOST = pyheader.base__check('HOST') # fname = 'metaio/simpleModel/Tag/0322-chipTarg-198C.json' # res = pyutil.readBaseFile(fname,ext='json', baseFile=HOST) # jobConfs = res = pyutil.read__buffer(res['fields']['text'],ext='json') # jobConfs = pyutil.readBaseFile('CAMICE/0322-chipTarg-198C.json') jobConfs = pyutil.readBaseFile('CAMICE/0411-chipTarg-CBF2-198C.json') bwMeta = pyutil.readBaseFile('meta/meta_chip.tsv') figs = pyutil.collections.OrderedDict() worker = funcy.partial(sjob.job__chipTargPaired, bwMeta=bwMeta) # assert 0 # for conf in jobConfs: def _worker(conf): try: for k in conf.keys(): if k not in worker.func.func_code.co_varnames: print(k, conf.pop(k))