快捷導(dǎo)航

python數(shù)據(jù)可視化繪制火山圖示例

更新時(shí)間：2022年06月25日 17:30:52 作者：Z_Y_H

這篇文章主要為大家介紹了python數(shù)據(jù)可視化繪制火山圖示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步，早日升職加薪

導(dǎo)入模塊

import numpy as np
import pandas as pd

1.讀取測(cè)試數(shù)據(jù)

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

2.查看數(shù)據(jù)

data.head()

3.篩選差異基因

# 3.嘗試寫(xiě)循環(huán)篩選上下調(diào)基因分類賦值給 "up" 和 "down" 和 "nosig" 加入pvalue條件
###loc函數(shù)：通過(guò)行索引 "Index" 中的具體值來(lái)取行數(shù)據(jù)（如取"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.查看數(shù)據(jù),發(fā)現(xiàn)多了type這一列

data.head()

5.統(tǒng)計(jì)個(gè)數(shù)

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
# 對(duì)padj取個(gè)-log10對(duì)數(shù)
data['-logpadj']=-data.padj.apply(math.log10)
# 查看
data[['log2FoldChange','padj','type','-logpadj']].head()

# 先設(shè)置一下自己的顏色
colors = ["#01c5c4","#ff414d", "#686d76"]
sns.set_palette(sns.color_palette(colors))
# 繪圖
ax=sns.scatterplot(x='log2FoldChange', y='-logpadj',data=data,
                hue='type',#顏色映射
                edgecolor = None,#點(diǎn)邊界顏色
                s=8,#點(diǎn)大小
                )
# 標(biāo)簽
ax.set_title("vocalno")
ax.set_xlabel("log2FC")
ax.set_ylabel("-log10(padj)")
#移動(dòng)圖例位置
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數(shù)據(jù)可視化繪制火山圖示例的詳細(xì)內(nèi)容，更多關(guān)于python數(shù)據(jù)可視化火山圖的資料請(qǐng)關(guān)注腳本之家其它相關(guān)文章！

您可能感興趣的文章:

相關(guān)文章

Python之多線程退出與停止的一種實(shí)現(xiàn)思路
這篇文章主要介紹了Python之多線程退出與停止的一種實(shí)現(xiàn)思路,具有很好的參考價(jià)值,希望對(duì)大家有所幫助,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教
2024-02-02
對(duì)python抓取需要登錄網(wǎng)站數(shù)據(jù)的方法詳解
今天小編就為大家分享一篇對(duì)python抓取需要登錄網(wǎng)站數(shù)據(jù)的方法詳解，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧
2018-05-05
解析Tensorflow之MNIST的使用
這篇文章主要介紹了解析Tensorflow之MNIST的使用，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧
2020-06-06
python編寫(xiě)五子棋游戲
這篇文章主要為大家詳細(xì)介紹了python編寫(xiě)五子棋游戲，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下
2021-05-05
Python3中str、bytes、bytearray轉(zhuǎn)化
本文主要介紹了Python3中str、bytes、bytearray轉(zhuǎn)化，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧
2023-04-04
python中安裝Scrapy模塊依賴包匯總
Scrapy的安裝有一些依賴包,沒(méi)有這些包是會(huì)安裝失敗的，下面我們就來(lái)詳細(xì)探討下
2017-07-07
詳解如何優(yōu)雅的用PyQt訪問(wèn)http
這篇文章主要我打開(kāi)詳細(xì)介紹了如何優(yōu)雅的用PyQt實(shí)現(xiàn)訪問(wèn)http,文中的示例代碼講解詳細(xì),具有一定的借鑒價(jià)值,感興趣的小伙伴可以了解下
2024-11-11
Python的IDEL增加清屏功能實(shí)例
下面小編就為大家?guī)?lái)一篇Python的IDEL增加清屏功能實(shí)例。小編覺(jué)得挺不錯(cuò)的，現(xiàn)在就分享給大家，也給大家做個(gè)參考。一起跟隨小編過(guò)來(lái)看看吧
2017-06-06
PyTorch 中的傅里葉卷積實(shí)現(xiàn)示例
這篇文章主要介紹了PyTorch 中的傅里葉卷積實(shí)現(xiàn)示例，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧
2020-12-12
使用Python3?Boto3包刪除AWS?CloudFormation的棧(Stacks)
這篇文章主要介紹了如何使用Python3?Boto3刪除AWS?CloudFormation的棧（Stacks）,本文給大家介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值，需要的朋友可以參考下
2024-01-01