logo
publist
写文章

简介

该用户还未填写简介

擅长的技术栈

可提供的服务

暂无可提供的服务

Python中pairplot函数可视化探索数据特征间的关系之seaborn

import matplotlib.pylab as pltimport seaborn as snsdf = sns.load_dataset('iris')my_dpi=96plt.figure(figsize=(480/my_dpi, 480/my_dpi), dpi=my_dpi)sns.pairplot(df)plt.show()import matplotlib.pylab as pl

#python
TCGA癌症中英文及缩写名称对应表

本文转自:https://blog.csdn.net/Taylent/article/details/89575854主要是为了记录一下TCGA数据库的癌症名称信息。缩写英文名称中文名称ACCAdrenocortical carcinoma肾上腺皮质癌BLCABladder Urothelial Carcinoma膀...

Python画散点图和三维图之seaborn

import matplotlib.pyplot as pltimport numpy as npimport seaborn as snsimport pandas as pddf=pd.DataFrame({'x': range(1,101), 'y': np.random.randn(100)*15+range(1,101) })my_dpi=96plt.figure(figsize=(48

#python
samtools的学习之路(一)

1.本文主要是针对samtools view 命令的用法和相关参数的说明。2.view命令的主要功能是查看bam和sam文件的内容。3.view命令的用法和常用参数。samtools view -hsamtools view: No input provided or missing option argument.Usage: samtools view [options] ...

TCGA中RNA-counts数据中ensemble_ID转换成gene_ID的方法

1.测试数据的准备和相关包的安装。library(stringr)> d1 <- read.table('test.txt', sep = '\t', header = TRUE)> d1tagtcga1 ENSG00000000003.13 2969 4725 1350 16672ENSG00...

Python画核密度估计图之seaborn

1.核密度估计图。import matplotlib.pylab as pltimport seaborn as snsdf = sns.load_dataset('iris')p1=sns.kdeplot(df['sepal_width'])plt.show()import matplotlib.pylab as pltimport seaborn as snsdf = sns.load_dat

#python
C语言学习之求1到100之间奇数和

1.求1到100之间的奇数之和。# include <stdio.h>int main(void){int i;int sum = 0;for (i=1; i<101; ++i){if (i%2 == 1)sum += i; // sum = sum + i;}printf("sum = %d\n", sum)return...

去除fastq文件中不足四行的

1.下机的fastq文件中四行表示一个reads,但是有的reads 是不足四行的,需要我们去过滤这些reads并且也过滤掉第三行碱基数不等的。下面的python脚本专门处理该fastq文件。#!/usr/bin/env python#-*- encoding=UTF-8 -*-import sysfa = open(sys.argv[1], "r")fb = open(sys.a...

#python
python简单画图示例

1.使用python中matplotlib模块简单的画图脚本示例。#!/usr/bin/env python# -*- coding: UTF-8 -*-import matplotlib.pyplot as pltimport pandas as pdfrom pandas import DataFramedata1 = pd.read_csv("/test/king.xls...

#python#matplotlib
根据提供的fastq文件列表,针对单个fastq文件并行计算reads的长度,全部完成后计算N50

#!/usr/bin/env python#-*- encoding=UTF-8 -*-import multiprocessingimport gzipimport sys## this function cal the sequence length per filedef cal_length(f):seq = []reads = 0tt = gzip.open(f...

#python
    共 12 条
  • 1
  • 2
  • 请选择