首页 > 其他分享 >PyDESeq2使用

PyDESeq2使用

时间:2022-12-23 20:34:27浏览次数:58  
标签:... 01 df seconds 使用 PyDESeq2 condition

PyDESeq2

Python版的DESeq2已上线,以后就可以使用Python来做差异分析了。目前文章还在bioRxiv。我来简单尝尝鲜。

安装

使用mamba或者conda来新建一个虚拟环境,然后使用pip安装。

mamba create -n pydeseq2 python
mamba activate pydeseq2
pip install pydeseq2

用法

作为Python版的DESeq2, 用法和R里差不多。

数据读取

count_file = "test_counts.csv"
condition_file = "test_clinical.csv"

counts_df = pd.read_csv(count_file, index_col=0).T
condition_df = pd.read_csv(condition_file, index_col=0)

数据集是一个100个样本,每个样本10个基因的小测试集。而其中50个样本属于条件A,另50个样本属于条件B。

>>> counts_df.shape
(100, 10)
>>> counts_df.head()
         gene1  gene2  gene3  gene4  gene5  gene6  gene7  gene8  gene9  gene10
sample1     12     22      2    187     15      2     13     57     56       6
sample2     10      6     20     99     55      0     35     96     43       1
sample3      0     28      3     96     38      2      9     54     27      14
sample4      7     28     10    170     16     10     17     38     18      16
sample5      2     31      5    126     23      2     19     53     31      18

>>> condition_df
          condition
sample1           A
sample2           A
sample3           A
sample4           A
sample5           A
...             ...
sample96          B
sample97          B
sample98          B
sample99          B
sample100         B

[100 rows x 1 columns]

构建DeseqDataSet 对象

和DESeq2类似

# 构建DeseqDataSet 对象
dds = DeseqDataSet(counts_df, condition_df, design_factor="condition")
# 离散度和log fold-change评估.
dds.deseq2()
# Fitting size factors...
# ... done in 0.00 seconds.
# Fitting dispersions...
# ... done in 0.64 seconds.
# Fitting dispersion trend curve...
# ... done in 0.03 seconds.
# Fitting MAP dispersions...
# ... done in 0.63 seconds.
# Fitting LFCs...
# ... done in 0.66 seconds.
# Refitting 0 outliers.

统计分析

差异表达统计检验分析

res = DeseqStats(dds)
# 执行统计分析并返回结果
res_df = res.summary()

结果如下

>>> res_df
         baseMean  log2FoldChange     lfcSE       stat        pvalue          padj
gene1   10.306788        1.007045  0.225231   4.471161  7.779603e-06  2.593201e-05
gene2   24.718815       -0.059670  0.165606  -0.360311  7.186146e-01  7.186146e-01
gene3    4.348135       -0.166592  0.325445  -0.511891  6.087275e-01  6.763639e-01
gene4   98.572300       -2.529204  0.136752 -18.494817  2.273125e-76  2.273125e-75
gene5   38.008562        1.236663  0.151824   8.145377  3.781028e-16  1.890514e-15
gene6    4.734285        0.212656  0.304487   0.698408  4.849222e-01  6.061527e-01
gene7   30.011855       -0.445855  0.150575  -2.961017  3.066249e-03  5.110415e-03
gene8   59.330642        0.372080  0.118911   3.129070  1.753603e-03  3.507207e-03
gene9   46.779546        0.547280  0.124922   4.380966  1.181541e-05  2.953853e-05
gene10  11.963156        0.494775  0.229494   2.155940  3.108836e-02  4.441194e-02

最后

PyDESeq2目前主要的核心功能差异分析已完成,虽然相比DESeq2而言还在起步阶段, 不过很明显可以发现Python的组学分析生态正逐渐完善~

参考

https://github.com/owkin/PyDESeq2

标签:...,01,df,seconds,使用,PyDESeq2,condition
From: https://www.cnblogs.com/huanping/p/17001567.html

相关文章

  • Java中资源文件的使用(properties)
    properties文件是一种配置文件,主要用于表达配置信息,文件类型为*.properties,格式为文本文件,文件的内容是格式是"键=值"的格式,在properties文件中,可以用"#"来作注释。 一......
  • 使用感受
    三丰云是一个免费的云服务器平台!云主机质量非常好,推荐大家体验。很适合学生和初学者。现在互联网上有很多免费的主机。我用过很多,对各种空间类型和性能的容忍度都是最高的。......
  • 【Flutter】packages思维以及使用Java添加Android平台特定的实现在Flutter框架里的体
    packages思维即使软件包未在Pub。对于未用于公共发布的特殊插件或尚未准备好发布的软件包,可以使用其他依赖选项:dependencies:flutter:sdk:flutter路径依赖性:Fl......
  • 一文了解 Go fmt 标准库的常用占位符及其简单使用
    耐心和持久胜过激烈和狂热。哈喽大家好,我是陈明勇,今天分享的内容是Gofmt标准库的常用占位符及其简单使用。如果本文对你有帮助,不妨点个赞,如果你是Go语言初学者,不妨点个......
  • ag-grid-vue在排班组件中的使用
    在很多医务系统中,会有医生护士排班的业务,比如下面这种:   可是实现排班数据的上移,下移,并且能够通过右键菜单实现排班等操作;直接上源码:首先安装依赖:npminstall a......
  • 使用 Fail2Ban 保护 Gitea
    使用Fail2Ban保护GiteaFail2Ban是一个入侵检测系统框架,其工作原理是检查客户端登录日志,自动标记多次登录尝试失败的客户端IP,并且在一段时间内阻止其访问。对与互联网......
  • Libevent使用例子,从简单到复杂
    本文从简单到复杂,展示如何使用libevent。网上的许多例子都是只有服务器端的,本文里面客户端和服务器端都有,以飨读者。     关于libevent编程时的一些疑问可以阅读......
  • 使用 REST 访问 Neo4j 数据
    本指南将引导您完成创建应用程序的过程,该应用程序通过基于超媒体休息的前端。您将构建什么您将构建一个Spring应用程序,该应用程序允许您创建和检索存储在​​Person​​Ne......
  • 使用 REST 访问 MongoDB 数据
    本指南将引导您完成创建应用程序的过程,该应用程序通过基于超媒体休息的前端。您将构建什么您将构建一个Spring应用程序,该应用程序允许您创建和检索存储在​​Person​​蒙......
  • 使用 REST 访问 Pivotal GemFire 中的数据
    本指南将引导您完成创建访问存储在中的数据的应用程序的过程ApacheGeode​通过基于超媒体 休息前端.您将构建的内容您将构建一个 Spring Web应用程序,该应用程序允许您......