跳到正文
cccskills
免费在 GitHub 公开

bioconductor-glmgampoi

Fit linear models to overdispersed count data. The package can estimate the overdispersion and fit repeated models for matrix input. It is designed to handle large input datasets as they typically occur in single cell RNA-seq experiments.

中文简介正在准备中。当前显示原文说明。

查看安装方法

包含的文件(1)

  • SKILL.md5.3 KB

SKILL.md(原文)

安装前可以查看交给 AI 代理的指令内容。

glmGamPoi

Dependencies & Environment

Package-intrinsic requirements from the Bioconductor landing page — reproduce in any R environment.

  • Version: 1.24.0 · Bioconductor: 3.23 · R: ≥ 4.6
  • Imports: Rcpp, beachmat, DelayedMatrixStats, matrixStats, MatrixGenerics, SparseArray, S4Vectors, DelayedArray, HDF5Array, Matrix, SummarizedExperiment, SingleCellExperiment, BiocGenerics, rlang, vctrs
  • System requirements: C++17
  • Install: BiocManager::install("glmGamPoi")

When to Use

  • Fitting Gamma-Poisson (negative binomial) generalized linear models on high-dimensional count data using glm_gp().
  • Performing differential expression analysis on single-cell RNA-seq data using quasi-likelihood ratio testing via test_de().
  • Aggregating single-cell counts into pseudobulk samples for multi-condition comparisons using pseudobulk().
  • Processing massive datasets (e.g., millions of cells) efficiently by setting on_disk = TRUE to avoid memory limits.

When NOT to Use

  • For standard bulk RNA-seq with small sample sizes, use DESeq2 or edgeR instead because they are classical methods heavily optimized for low-replicate bulk data.
  • For normalized or log-transformed data, use standard linear models instead because glmGamPoi requires raw count data to fit Gamma-Poisson models.

Data Requirements

  • A raw count matrix (dense matrix, sparse matrix, or on-disk HDF5Array).
  • A SingleCellExperiment or SummarizedExperiment object containing unnormalized integer counts.

Key Parameters

  • design: Formula or model matrix specifying the experimental design (e.g., ~ 1 or ~ condition*cell + ind).
  • on_disk: Boolean indicating whether to force the calculation to happen on-disk (e.g., TRUE or FALSE).
  • size_factor: Method to calculate size factors (e.g., "ratio").
  • reference_level: String specifying the baseline level for factors in the design formula.
  • contrast: The specific comparison to test in test_de() (e.g., using the cond() function or coefficient arithmetic).
  • group_by: Used in pseudobulk() with vars() to quote the grouping columns from colData.
  • sort_by: Column to sort the differential expression results by in test_de() (e.g., pval).

Best Practices

  • Filter out genes with very low counts (e.g., rowSums(counts) > 5) before fitting the model to improve performance and reliability.
  • Use pseudobulk() to aggregate single-cell data by individual/replicate before running glm_gp() to ensure reliable statistical tests and avoid underestimating variance.
  • Use the cond() function within test_de() to cleanly specify complex contrasts (e.g., cond(cell = "B cells", condition = "stim") - cond(...)).

Common Pitfalls

  • Overly optimistic p-values: Applying test_de() directly to single-cell data without pseudobulking treats cells from the same sample as independent replicates; fix this by aggregating data with pseudobulk() first.
  • Empty factor levels: glm_gp() will complain if factor levels have no observations; fix this by calling droplevels() on the metadata columns before fitting.
  • Deprecated pseudobulk workflow: Using the pseudobulk_by argument inside test_de() is deprecated and wastes computation time; fix this by calling pseudobulk() before glm_gp().

Alternatives

  • DESeq2: A classical method for bulk RNA-seq differential expression, though significantly slower on large single-cell datasets.
  • edgeR: Another classical method for RNA-seq that uses negative binomial models, but glmGamPoi is optimized to be faster for single-cell data.

Citations

  • Kang et al. (2018). Multiplexed droplet single-cell RNA-sequencing using natural genetic variation. Nature Biotechnology, 36(1), 89-94.

References

<!-- biomate-cta -->

Run this on BioMate

This skill is the knowledge layer — when, why, and how to use glmgampoi. To run this analysis on your own data with managed compute, automated QC, and reproducible outputs, use BioMate — free to start.

▶ Open glmgampoi on BioMate →

评价

还没有评价。欢迎分享你的使用感受。

评论

无需登录即可发表。欢迎分享使用感受、提问,或说明在哪些环境下无法运行。

还没有评论。来写第一条吧。

提交后会立即公开。请勿填写邮箱、电话等个人信息。URL 最多 2 个。

0 / 1000 字

此表单受 reCAPTCHA 保护,适用 Google 的隐私权政策和服务条款。

同一仓库的技能

The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7972026年6月21日 更新

This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7972026年6月21日 更新

This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7972026年6月21日 更新

apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7972026年6月21日 更新

ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7972026年6月21日 更新

AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7972026年6月21日 更新

查看 bioMate-AI 的全部技能

报告此技能的问题