跳到正文
cccskills
免费在 GitHub 公开

bioconductor-sechm

sechm provides a simple interface between SummarizedExperiment objects and the ComplexHeatmap package. It enables plotting annotated heatmaps from SE objects, with easy access to rowData and colData columns, and implements a number of features to make the generation of heatmaps easier and more flexible. These functionalities used to be part of the SEtools package.

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

实测评分 未实测

查看安装方法

包含的文件(1)

  • SKILL.md4.5 KB

SKILL.md(原文)

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

sechm

Workflows

Standard Workflow

sechm provides a simple interface between SummarizedExperiment objects and the ComplexHeatmap package. It enables plotting annotated heatmaps from SE objects, with easy access to rowData and colData columns, and implements a number of features to make the generation of heatmaps easier and more flexible.

library(SummarizedExperiment)
library(sechm)

# Load example data
data("Chen2017", package="sechm")
SE <- Chen2017

# 1. Basic scaled heatmap of top variable genes
g <- c("Egr1", "Nr4a1", "Fos", "Egr2", "Sgk1", "Arc", "Dusp1", "Fosb", "Sik1")
sechm(SE, features=g, do.scale=TRUE)

# 2. Heatmap of Log2 Fold Changes with quantile capping (breaks)
sechm(SE, features=g, assayName="logFC", breaks=0.985)

# 3. Customizing colors using global package options
setSechmOption("hmcols", value=c("white","grey","black"))
sechm(SE, features=g, do.scale=TRUE)
resetAllSechmOptions()

Input: A SummarizedExperiment object and a vector of feature names. Output: A Heatmap object from ComplexHeatmap.

When to Use

  • To plot highly customizable, annotated heatmaps directly from SummarizedExperiment objects using sechm().
  • To display row and column annotations automatically extracted from rowData and colData (e.g., top_annotation, left_annotation).
  • To combine multiple heatmaps from different SummarizedExperiment objects with aligned rows and consistent scales using crossHm().
  • To highlight specific genes of interest on a large heatmap using the mark argument.

When NOT to Use

  • For plotting heatmaps from raw matrices or data frames without a SummarizedExperiment container, use ComplexHeatmap::Heatmap or pheatmap directly.

Data Requirements

  • A SummarizedExperiment object containing at least one assay (e.g., "logcpm", "logFC").
  • Row annotations stored in rowData(SE) and column annotations stored in colData(SE).

Key Parameters

  • features: A vector of row/feature names to plot.
  • assayName: Name of the assay to extract data from (e.g., "logFC").
  • do.scale (FALSE): Logical indicating whether to scale rows.
  • top_annotation: Column names from colData to display as top annotations.
  • left_annotation: Column names from rowData to display as left annotations.
  • breaks (TRUE): Controls symmetric scaling and quantile capping (e.g., 0.985 for 98.5% quantile capping).
  • gaps_at: Column name from colData to introduce gaps between columns.
  • uniqueScale (FALSE): Logical indicating whether to enforce a unique color scale across datasets in crossHm().

Best Practices

  • Store default visualization settings (like default assay and annotations) in the object's metadata under metadata(SE)$default_view to simplify repetitive plotting.
  • Use quantile capping (e.g., breaks=0.985) when plotting fold-changes to prevent extreme outlier values from dominating the color scale.
  • Define custom annotation colors in metadata(SE)$anno_colors to ensure they are automatically applied across all heatmaps generated from that object.

Common Pitfalls

  • Extreme values compressing the color scale: Outliers in fold-change data can make normal variations invisible. Fix: Set breaks to a quantile value less than 1 (e.g., breaks=0.985) to cap the color scale.
  • Duplicate heatmap or annotation names when combining heatmaps with +: Combining raw sechm objects can lead to conflicts. Fix: Use crossHm() to safely plot multiple SummarizedExperiment objects side-by-side.

Alternatives

  • ComplexHeatmap: The underlying engine, but requires manual extraction of assays and annotations from SummarizedExperiment.
  • pheatmap: A classic heatmap package, but lacks native integration with SummarizedExperiment metadata.

Citations

  • Chen et al., 2017 (vignette reference for example data).

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-kb7932026年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-kb7932026年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-kb7932026年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-kb7932026年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-kb7932026年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-kb7932026年6月21日 更新

查看 bioMate-AI 的全部技能

报告此技能的问题