代码拉取完成,页面将自动刷新
同步操作将从 Vincent.Gao/甲基化生信分析pipeline 强制同步,此操作会覆盖自 Fork 仓库以来所做的任何修改,且无法恢复!!!
确定后同步将在后台操作,完成时将刷新页面,请耐心等待。
import shelve
import matplotlib.pyplot as plt
# with shelve.open("/home/DUAN/methylation/data_files/ucsc") as db:
# dic = db["ucsc"]
# # new_dic = {i:dic[i].values() for i in dic}
# # print(new_dic)
# new_list = list(dic.keys())[:9]+list(dic.keys())[10:25]
# print(new_list)
# # print(list(dic["chr1"].values()))
# new_dic = {i:list(dic[i].values()) for i in new_list}
# with shelve.open("/home/DUAN/methylation/data_files/dic") as dbs:
# dbs["dic"] = new_dic
# with shelve.open("/home/DUAN/methylation/data_files/dic") as dbs:
# new_dic = dbs["dic"]
# print(new_dic.keys())
# # def func(listTemp, n):
# # for i in range(0, len(listTemp), n):
# # yield listTemp[i:i + n]
# # a=len(range(247754845, 247760556))
# # b = 247754845+620
# # print(a)
# # print(list(range(247754845, 247760556))[0])
# # print((b-list(range(247754845, 247760556))[0])/a)
# help_dict = {"-2kb":0,"10%":0,"20%":0,"30%":0,"40%":0,"50%":0,"60%":0,"70%":0,"80%":0,"90%":0,"100%":0,"+2kb":0}
# plt.plot(x,y)
# plt.show()
# plt.savefig("/home/DUAN/plt.png")
dic = {'-2kb': 0, '10%': 5934, '20%': 5749, '30%': 5679, '40%': 5613, '50%': 5578, '60%': 5666, '70%': 5972, '80%': 5621, '90%': 5796, '100%': 6058, '+2kb': 0}
y = list(dic.values())[1:11]
x = list(dic.keys())[1:11]
plt.plot(x,y)
plt.xlabel('gene position sliced to ten equal parts')
plt.ylabel('methylation count in certain area')
plt.show()
plt.savefig("/home/DUAN/plt.png")
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。