python 鉴定H3k27ac峰,其与所有BLUEPRINT数据集中的启动子重叠并随后鉴定转录因子
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了python 鉴定H3k27ac峰,其与所有BLUEPRINT数据集中的启动子重叠并随后鉴定转录因子相关的知识,希望对你有一定的参考价值。
import xmlrpclib
import time
url = "http://deepblue.mpi-inf.mpg.de/xmlrpc"
user_key = "anonymous_key"
server = xmlrpclib.Server(url, allow_none=True)
(status, exps) = server.select_regions("", "GRCh38", "H3k27ac", "", "",
"BLUEPRINT Epigenome", "chr1", None,
None, user_key)
(status, exps_peaks) = server.query_experiment_type(exps, "peaks", user_key)
(status, promoters) = server.select_annotations("promoters", "GRCh38","chr1",
None, None, user_key)
(status, exps_promoters) = server.intersection(exps_peaks, promoters, user_key)
(status, tf) = server.select_regions("","hg19", ["SP1"], "",
"", "ENCODE", "chr1", None, None, user_key)
(status, ts_signals) = server.query_experiment_type(tf, "signal", user_key)
(status, final) = server.intersection(tf, exps_promoters, user_key)
(status, request) = server.get_regions(final,
"CHROMOSOME,START,END,@NAME,@EPIGENETIC_MARK,@BIOSOURCE",
user_key)
# Wait the processing
(status, info) = server.info(request, user_key)
while info[0]["state"] != "done" and info[0]["state"] != "error":
time.sleep(5)
print info
(status, info) = server.info(request, user_key)
(status, regions) = server.get_request_data(request, user_key)
print regions
以上是关于python 鉴定H3k27ac峰,其与所有BLUEPRINT数据集中的启动子重叠并随后鉴定转录因子的主要内容,如果未能解决你的问题,请参考以下文章
R语言学习 - 线图绘制
python 通过H-ESC H3K4me3峰聚集肝脏数据实验的甲基化数据信号
python 总结来自人肝细胞的H3K4me3峰区域的肝组织中的DNA甲基化水平
硬件开源电路STC8H3K64S4开发板
R语言学习 - 线图一步法
Journal of Proteomics Research | 构建用于鉴定蓖麻毒素的串联质谱库