pydeseq2

Perform differential expression analysis on bulk RNA-seq count data.

52|6|Updated Nov 24, 2025
One-click install
npx skills add https://github.com/ovachiever/droid-tings --skill pydeseq2
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: pydeseq2
Source: https://github.com/ovachiever/droid-tings/tree/main/skills/pydeseq2
Command: npx skills add https://github.com/ovachiever/droid-tings --skill pydeseq2

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires pydeseq2, pandas, numpy, scipy, and includes scripts (resource) and references (resource) components.

What problem does it solve?

PyDESeq2 provides a Python implementation of DESeq2 for differential expression analysis on bulk RNA-seq data, enabling Wald tests, multiple testing correction, and LFC shrinkage in Python.

Core Features & Use Cases

  • DESeq2-style normalization, dispersion estimation, and LFC fitting
  • Wald tests with BH FDR control
  • Optional apeGLM LFC shrinkage for visualization
  • Full workflow from counts + metadata to results

Quick Start

Load counts and metadata, run DESeq2 pipeline, and extract significant genes.

Frequently Asked Questions about pydeseq2

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I perform differential expression analysis on RNA-seq count data?

Differential expression analysis identifies genes with significant expression changes between conditions. PyDESeq2 applies DESeq2-style normalization, dispersion estimation, Wald tests with BH-FDR correction, and log2 fold-change calculation to bulk RNA-seq counts, producing a results dataframe with p-values and adjusted significance for each gene.

Does PyDESeq2 support multi-factor experimental designs?

Yes, PyDESeq2 handles both single-factor and multi-factor designs. It applies Wald tests and BH-FDR correction across complex experimental layouts, integrating with pandas and AnnData to manage metadata and design matrices for end-to-end workflows.

What's the best way to identify and handle outliers in RNA-seq analysis?

PyDESeq2 uses Cook's distance to detect outlier samples and observations. This helps flag potentially artifactual counts before downstream interpretation, improving the reliability of differential expression results.

Can I use PyDESeq2 with pandas and AnnData for my RNA-seq workflow?

Yes, PyDESeq2 integrates directly with pandas DataFrames for count matrices and metadata, and with AnnData for annotated data objects. This enables seamless end-to-end workflows from raw counts through results export.

What does log fold-change shrinkage do, and when should I use it?

Log fold-change shrinkage via apeGLM reduces noise in LFC estimates, particularly for lowly-expressed genes. It improves visualization and interpretation of results by producing more stable effect-size estimates for downstream analysis.

How do I export and interpret differential expression results from PyDESeq2?

PyDESeq2 generates a results dataframe (results_df) containing p-values, adjusted p-values, log fold-changes, and per-gene statistics. Export this for filtering by FDR threshold, fold-change cutoff, or downstream functional annotation.