跳到正文
cccskills
免费在 GitHub 公开

bioconductor-scarray

Provides large-scale single-cell omics data manipulation using Genomic Data Structure (GDS) files. It combines dense and sparse matrices stored in GDS files and the Bioconductor infrastructure framework (SingleCellExperiment and DelayedArray) to provide out-of-memory data storage and large-scale manipulation using the R programming language.

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

查看安装方法

包含的文件(1)

  • SKILL.md3.8 KB

SKILL.md(原文)

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

scarray

Workflows

Standard Workflow

Convert single-cell data to GDS format, load it as a GDS-backed SingleCellExperiment, and perform memory-efficient downstream analysis.

library(SCArray)
library(SingleCellExperiment)
library(scuttle)
library(scater)

# Load example SingleCellExperiment
fn <- system.file("extdata", "example.rds", package="SCArray")
sce <- readRDS(fn)

# Convert to GDS file
scConvGDS(sce, "test.gds")

# Load SingleCellExperiment from GDS
sce_gds <- scExperiment("test.gds")

# Normalize counts
sce_gds <- logNormCounts(sce_gds)

# Row and column summarization
col_mean <- colMeans(assays(sce_gds)$counts)
mvar <- scRowMeanVar(assays(sce_gds)$counts)

# Memory-efficient PCA and UMAP
sce_gds <- scRunPCA(sce_gds)
sce_gds <- runUMAP(sce_gds)

# Plot PCA
plotReducedDim(sce_gds, dimred="PCA")

Input: A SingleCellExperiment object or matrix; Output: A GDS-backed SingleCellExperiment object with normalized counts and reduced dimensions.

When to Use

  • To convert large-scale single-cell datasets (matrices or SingleCellExperiment objects) to GDS format using scConvGDS.
  • To load and manipulate single-cell data out-of-memory as a GDS-backed SingleCellExperiment using scExperiment.
  • To perform memory-efficient row/column summarization (e.g., colMeans, rowMeans, scRowMeanVar) on GDS-backed DelayedMatrix objects.
  • To run memory-efficient PCA on large-scale single-cell datasets using scRunPCA.

When NOT to Use

  • When the dataset is small enough to fit comfortably in memory and standard in-memory matrices are preferred.
  • For non-array-oriented genomic data (e.g., raw FASTQ or BAM files).

Data Requirements

  • Input single-cell datasets as a SingleCellExperiment object or a dense/sparse matrix (e.g., dgCMatrix).
  • Output GDS files (.gds) containing compressed assay data.

Key Parameters

  • SCArray.verbose (TRUE/FALSE): Global option to enable displaying debug information during SCArray operations.
  • dimred: The name of the reduced dimension slot to plot in plotReducedDim (e.g., "PCA", "UMAP").

Best Practices

  • Use scRunPCA instead of scater::runPCA for large-scale datasets to avoid high memory usage from realizing the full matrix in memory.
  • Enable verbose debugging with options(SCArray.verbose=TRUE) to track GDS-backed matrix operations.
  • Keep the GDS file connection closed when not in use by calling scClose if manually opened with scOpen.

Common Pitfalls

  • Attempting to run standard memory-intensive functions (like scater::runPCA) on very large GDS-backed objects, which triggers in-memory realization and defeats the out-of-memory benefit. Use scRunPCA instead.
  • Forgetting that GDS-backed matrices are DelayedMatrix objects, meaning operations are delayed until explicitly realized.

Alternatives

  • HDF5Array: For HDF5-backed out-of-memory data storage and manipulation.
  • TileDBArray: For TileDB-backed multi-dimensional array storage.
  • LoomExperiment: For Loom-backed single-cell data storage.

Citations

  • Zheng et al. (GDS format and gdsfmt package).

References

评价

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

评论

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

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

提交后会立即公开。请勿填写邮箱、电话等个人信息。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 的全部技能

报告此技能的问题