如何利用Python数据可视化画火山图，代码是什么-群英

上一篇：Python列表初始化的方法及代码是什么下一篇：用Python画散点图和直方图的操作分别是什么

今天这篇我们来学习和了解“如何利用Python数据可视化画火山图，代码是什么”，下文的讲解详细，步骤过程清晰，对大家进一步学习和理解“如何利用Python数据可视化画火山图，代码是什么”有一定的帮助。有这方面学习需要的朋友就继续往下看吧！

导入模块

import numpy as np
import pandas as pd

1.读取测试数据

data=pd.read_csv(r'E:\ZYH\R.project\rna-seq\lianxi1\exon_level\df.csv')

2.查看数据

data.head()

3.筛选差异基因

# 3.尝试写循环筛选上下调基因分类赋值给 "up" 和 "down" 和 "nosig" 加入pvalue条件
###loc函数：通过行索引 "Index" 中的具体值来取行数据（如取"Index"为"A"的行）
data.loc[(data.log2FoldChange>1)&(data.padj<0.05),'type']='up'
data.loc[(data.log2FoldChange<-1)&(data.padj<0.05),'type']='down'
data.loc[(abs(data.log2FoldChange)<=1)|(data.padj>=0.05),'type']='nosig'

4.查看数据,发现多了type这一列

data.head()

5.统计个数

data.type.value_counts()
up      123
down    103
Name: type, dtype: int64

6.绘火山图

import seaborn as sns
import math
import matplotlib.pyplot as plt
import matplotlib as mpl
%matplotlib inline
# 对padj取个-log10对数
data['-logpadj']=-data.padj.apply(math.log10)
# 查看
data[['log2FoldChange','padj','type','-logpadj']].head()

# 先设置一下自己的颜色
colors = ["#01c5c4","#ff414d", "#686d76"]
sns.set_palette(sns.color_palette(colors))
# 绘图
ax=sns.scatterplot(x='log2FoldChange', y='-logpadj',data=data,
                hue='type',#颜色映射
                edgecolor = None,#点边界颜色
                s=8,#点大小
                )
# 标签
ax.set_title("vocalno")
ax.set_xlabel("log2FC")
ax.set_ylabel("-log10(padj)")
#移动图例位置
ax.legend(loc='center right', bbox_to_anchor=(0.95,0.76), ncol=1)

7.保存图片

fig = ax.get_figure()
fig.savefig('./python_vocalno.pdf')

关于“如何利用Python数据可视化画火山图，代码是什么”就介绍到这了，如果大家觉得不错可以参考了解看看，如果想要了解更多，欢迎关注群英网络，小编每天都会为大家更新不同的知识。

标签： python绘制火山图

免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：mmqy2019@163.com进行举报，并提供相关证据，查实之后，将立刻删除涉嫌侵权内容。

上一篇：Python列表初始化的方法及代码是什么下一篇：用Python画散点图和直方图的操作分别是什么