- 博客(6)
- 收藏
- 关注
RSS订阅原创 pysam筛选reads写bam
使用pysam从bam中筛选readsbamfile = 'Final.bam'bam = pysam.AlignmentFile(bamfile, "r")outfile = pysam.AlignmentFile('out.bam', "wb", template=bam)ALL = set()for reads in bam: status = '''some filters''' if status: outfile.write(reads) print
2021-07-10 16:37:13
853
原创 线性代数--01 方程组的几何解释
麻省理工公开课《线性代数导论》学习笔记以及个人总结文章目录1. 二元一次方程组与矩阵1. 二元一次方程组与矩阵{2x−y=0−x+2y=3⟹[2−1−12]∗[xy]=[03](1) \begin{cases}2x -y = 0 \\-x + 2y = 3\end{cases} \Longrightarrow \left[\begin{matrix} 2 & -1 \\ -1 & 2 \\ \end{matrix}\right]* \left[\b.
2021-01-23 21:01:36
330
原创 什么是p-value
什么是p-value人们经常认为p-value与probability是一个概念，这是不对的。但是他们之间确实有一定的联系。probability表示的是某个事件发生的概率，而p-value表示的该事件与等价事件以及更rarer事件发生的概率之和，以表示当前事件是否是一个罕见的事件（数值越小表示越罕见）。在投硬币实验中，一个硬币连续扔5次，那么4H1T（Head/Tail）的概率为4H与1T的全排列的个数除以252^525，即525=0.15625\frac{5}{2^5}=0.15625255=0.
2021-01-14 16:32:25
9423
原创 使用pandas同时画两个箱线图
pandas 画箱线图测试文件：https://codechina.csdn.net/u014467816/ktools/-/blob/master/%E4%B8%A4%E4%B8%AA%E7%AE%B1%E7%BA%BF%E5%9B%BE.csvimport pandas as pdfrom matplotlib import pyplot as plt# 读取文件信息allInfo = pd.read_csv('test.csv',sep='\t')# 根据type类型分组，分别画出两个
2021-01-06 21:48:05
2059
原创 python处理时间
文章目录使用python处理时间信息1. 获取当前时间使用python处理时间信息1. 获取当前时间import time# 获取当前时间戳print('1:', time.time())# 将时间戳转换成时间数组d = time.time()print('2:', time.localtime(d))# 将时间数组，转成需要格式d = time.localtime(...
2019-09-27 13:00:19
120
原创 用python查看安装的模块
如何用python查看安装的python模块呢？开门见山…import sysimport osfor path in sys.path: print('模块存放路径： %s' % path) ''' 仅仅展示 site-packages文件夹中的模块 如果需要展示其他文件夹中的内容，可以把if 判断改为os.path.isdir(path) ''...
2019-09-25 10:09:54
654
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人