• Standard QC workflow with adjustable thresholds (all cells filtered the same way)
• Batch processing multiple datasets
• Quick exploratory analysis
• User wants the "just works" solution

Requirements: anndata, scanpy, scipy, matplotlib, seaborn, numpy

Parameters:

Customize filtering thresholds and gene patterns using command-line parameters:

• --output-dir - Output directory
• --mad-counts, --mad-genes, --mad-mt - MAD thresholds for counts/genes/MT%
• --mt-threshold - Hard mitochondrial % cutoff
• --min-cells - Gene filtering threshold
• --mt-pattern, --ribo-pattern, --hb-pattern - Gene name patterns for different species

Use --help to see current default values.

Outputs:

All files are saved to <input_basename>_qc_results/ directory by default (or to the directory specified by --output-dir):

• qc_metrics_before_filtering.png - Pre-filtering visualizations
• qc_filtering_thresholds.png - MAD-based threshold overlays
• qc_metrics_after_filtering.png - Post-filtering quality metrics
• <input_basename>_filtered.h5ad - Clean, filtered dataset ready for downstream analysis
• <input_basename>_with_qc.h5ad - Original data with QC annotations preserved

If copying outputs for user access, copy individual files (not the entire directory) so users can preview them directly.

Workflow Steps

The script performs the following steps:

1. Calculate QC metrics - Count depth, gene detection, mitochondrial/ribosomal/hemoglobin content
2. Apply MAD-based filtering - Permissive outlier detection using MAD thresholds for counts/genes/MT%
3. Filter genes - Remove genes detected in few cells
4. Generate visualizations - Comprehensive before/after plots with threshold overlays

Approach 2: Modular Building Blocks (For Custom Workflows)

For custom analysis workflows or non-standard requirements, use the modular utility functions from scripts/qc_core.py and scripts/qc_plotting.py:

# Run from scripts/ directory, or add scripts/ to sys.path if needed
import anndata as ad
from qc_core import calculate_qc_metrics, detect_outliers_mad, filter_cells
from qc_plotting import plot_qc_distributions  # Only if visualization needed

adata = ad.read_h5ad('input.h5ad')
calculate_qc_metrics(adata, inplace=True)
# ... custom analysis logic here

When to use this approach:

• Different workflow needed (skip steps, change order, apply different thresholds to subsets)
• Conditional logic (e.g., filter neurons differently than other cells)
• Partial execution (only metrics/visualization, no filtering)
• Integration with other analysis steps in a larger pipeline
• Custom filtering criteria beyond what command-line params support

Available utility functions:

From qc_core.py (core QC operations):

…

customer-research

围绕客户问题进行多来源调研与溯源，快速整理背景并支持准确回复。

MCP Analyst

帮助用户直接分析本地 CSV 或 Parquet 大数据文件并生成洞察。