#!/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))
