bulk-rna-seq-batch-correction-with-combat

Harmonize batch effects across bulk RNA-seq matrices with ComBat and export raw and corrected matrices.

32|5|Updated Mar 12, 2026
One-click install
npx skills add https://github.com/omicverse/omicclaw --skill bulk-rna-seq-batch-correction-with-combat-omicverse
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: bulk-rna-seq-batch-correction-with-combat
Source: https://github.com/omicverse/omicclaw/tree/main/src/omicverse_skills/skills/bulk-combat-correction
Command: npx skills add https://github.com/omicverse/omicclaw --skill bulk-rna-seq-batch-correction-with-combat-omicverse

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill includes references (resource) components.

What problem does it solve?

Bulk RNA-seq studies often generate batch effects across cohorts; this skill provides a structured workflow to remove batch effects using ComBat, enabling integrated analyses and fair cross-cohort comparisons.

Core Features & Use Cases

  • Import libraries and prepare data: load expression matrices, convert to AnnData objects, and annotate batches.
  • Concatenate and harmonize: merge cohorts on shared genes to preserve common features across batches.
  • Batch correction and export: run ComBat to produce a corrected data layer while preserving raw counts, and export both raw and corrected matrices plus a harmonised AnnData file for downstream analyses.
  • Benchmark and visualize: generate before/after visualizations (e.g., PCA embeddings) to assess mixing and batch removal.

Quick Start

Run the ComBat batch correction workflow on your merged bulk RNA-seq datasets to obtain raw and corrected matrices and a harmonised AnnData file.

Frequently Asked Questions about bulk-rna-seq-batch-correction-with-combat

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

FAQPage Schema
How do I remove batch effects in bulk RNA-seq data before PCA and clustering?

Remove batch effects in bulk RNA-seq data by merging expression matrices on shared genes, annotating batches, and running ComBat to produce a corrected data layer while preserving raw counts for downstream PCA and clustering.

What is the best way to integrate bulk RNA-seq cohorts from diverse experiments?

Integrate bulk RNA-seq cohorts by concatenating matrices on shared genes to preserve common features, then applying ComBat harmonization to ensure balanced integration and fair cross-cohort comparisons.

Can I use anndata to store both raw and ComBat-corrected matrices?

Use anndata to store both raw and ComBat-corrected matrices by running the correction to create a harmonized AnnData file that preserves original counts while exporting the adjusted data for downstream analyses.

How do I visualize batch correction results for bulk RNA-seq?

Visualize batch correction results for bulk RNA-seq by generating before and after PCA embeddings using matplotlib to assess mixing and verify batch removal across merged cohorts.

Do I need pyComBat to harmonize bulk RNA-seq expression matrices?

You need pyComBat to harmonize bulk RNA-seq expression matrices, along with anndata and pandas for data handling, to effectively correct batch effects and export harmonized matrices.