跳到正文
cccskills
免费在 GitHub 公开

bioconductor-seqlogo

seqLogo takes the position weight matrix of a DNA sequence motif and plots the corresponding sequence logo as introduced by Schneider and Stephens (1990).

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

查看安装方法

包含的文件(1)

  • SKILL.md4.2 KB

SKILL.md(原文)

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

seqLogo

When to Use

  • Visualizing sequence logos for DNA sequence alignments using seqLogo().
  • Representing position weight matrices (PWM) and their information content profiles using makePWM().
  • Plotting RNA sequence logos by specifying the RNA alphabet and custom fill colors.

When NOT to Use

  • For highly customized ggplot2-based sequence logos, use ggseqlogo instead because seqLogo uses the base grid package for rendering.
  • For protein/amino acid sequence logos, use ggseqlogo instead because seqLogo is currently implemented primarily for DNA/RNA sequences.

Data Requirements

  • A simple matrix or data.frame where column probabilities add up to 1.0.
  • An instance of the pwm class created by makePWM().

Key Parameters

  • pwm: The position weight matrix for which the sequence logo is to be plotted (can be a pwm object, matrix, or data.frame).
  • ic.scale: Logical indicating whether the height of each column is to be proportional to its information content.
  • fill: A named character vector to change the default colors for the nucleotides (e.g., A, C, G, T or U).
  • alphabet: Character string to specify the alphabet, such as "RNA" for RNA logos.

Best Practices

  • Construct a pwm object from a matrix or data frame using makePWM() to automatically check that column probabilities sum to 1.0.
  • Use ic.scale=FALSE if you want all columns to have the same uniform height rather than being scaled by information content.
  • When plotting RNA logos, specify alphabet="RNA" in makePWM() and provide fill colors for either c("A", "C", "G", "U") or c("A", "C", "G", "T").

Common Pitfalls

  • Unnormalized Matrix: Column probabilities do not sum to 1.0. Fix: Ensure the input matrix or data.frame is properly normalized before passing it to makePWM().
  • Incorrect color mapping for RNA logos: The default colors may not map correctly to Uracil. Fix: Provide a named character vector to the fill argument in seqLogo() with specific colors for A, C, G, and T/U.
  • Attempting to plot amino acid sequences: The package will fail or produce incorrect logos. Fix: Use a different package, as the vignette notes seqLogo is designed for DNA/RNA alignments.

Alternatives

  • ggseqlogo: Uses ggplot2 for rendering, allowing for more modern and flexible plot customization.
  • universalmotif: Provides comprehensive motif manipulation and visualization for various alphabets.
  • motifStack: Allows plotting stacked logos and tree-structured logos for multiple motifs.

Citations

  • Schneider, Thomas D., and R.Michael Stephens. 1990. "Sequence logos: a new way to display consensus sequences." Nucleic Acids Research 18 (20): 6097-6100.
  • Schneider, Thomas D., Gary D. Stormo, Larry Gold, and Andrzej Ehrenfeucht. 1986. "Information Content of Binding Sites on Nucleotide Sequences." Journal of Molecular Biology 188 (3): 415-31.

References

<!-- biomate-cta -->

Run this on BioMate

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

▶ Open seqlogo 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-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 的全部技能

报告此技能的问题