跳到正文
cccskills
免费在 GitHub 公开

bioconductor-dnacopy

Implements the circular binary segmentation (CBS) algorithm to segment DNA copy number data and identify genomic regions with abnormal copy number.

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

查看安装方法

包含的文件(1)

  • SKILL.md5.5 KB

SKILL.md(原文)

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

DNAcopy

Dependencies & Environment

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

  • Version: 1.86.0 · Bioconductor: 3.23 · R: ≥ 4.6
  • Depends: Imports
  • System requirements: URL
  • Install: BiocManager::install("DNAcopy")

When to Use

  • Copy Number Segmentation: Finding change-points in array DNA copy number data (e.g., array CGH) using the segment() function.
  • Identifying Gains and Losses: Determining regions of gained and lost copy number by plotting segment means with plot().
  • Outlier Smoothing: Smoothing single point outliers in log-ratio data prior to segmentation using smooth.CNA().

When NOT to Use

  • For allele-specific copy number estimation or tumor purity/ploidy correction, use PureCN because DNAcopy only segments total copy number.
  • For single-cell DNA copy number profiling where specialized noise models are required, use AneuFinder because DNAcopy assumes bulk-level signal distributions.
  • For multi-sample joint segmentation, use copynumber because DNAcopy processes samples independently.

Data Requirements

  • Input Object: Data must be converted into a CNA object using the CNA() function.
  • Required Columns: A matrix or vector of log-ratio data, a vector of Chromosome identifiers, and a vector of Position (map locations).
  • Data Type: Typically specified as data.type="logratio".

Key Parameters

  • undo.splits ("sdundo"): Method to get rid of unnecessary change-points due to local trends.
  • undo.SD (3): The number of standard deviations used when undo.splits="sdundo" to remove splits.
  • plot.type ("w", "s", "c", "p"): Determines the layout of the plot() function (e.g., whole genome, by chromosome, across studies, or ordered by plateau).
  • p.method ("perm"): Method for p-value computation. The default is a faster hybrid approach, but full permutations can be forced.
  • nperm (10000): Number of permutations used in the hybrid method.
  • alpha: Significance level for the test to accept change-points.
  • verbose (1): Controls the printing of progress during the segment() run.

Best Practices

  • Smooth Data First: Always run smooth.CNA() on your CNA object to mitigate the effect of single point outliers before running segment().
  • Undo Local Trends: Use undo.splits="sdundo" with undo.SD=3 in segment() to remove false-positive change-points caused by local data trends.
  • Subset Large Data: Use subset.CNA() to subset by chromosome and sample so that segmentation does not have to be run on a whole data set at once.
  • Determine Thresholds: Use plot() with plot.type="p" to order segments by their chromosome means, helping to visually determine the plateaus for calling gains (e.g., 0.4) and losses (e.g., -0.6).

Common Pitfalls

  • Over-segmentation from Local Trends: The algorithm may find change-points due to local trends rather than true copy number shifts. Fix: Apply the undo.splits="sdundo" parameter in segment().
  • Slow Computation on Large Datasets: Running the full permutation algorithm takes O(N^2) computations. Fix: Rely on the default hybrid approach (Gaussian approximation + permutation) rather than setting p.method='perm', or reduce nperm.
  • Plotting Clutter: Plotting multiple chromosomes at once can be hard to read. Fix: Use plot.type="s" to plot by chromosome within a study.

Alternatives

  • copynumber: For multi-sample segmentation and joint segmentation of copy number data.
  • PureCN: For joint estimation of tumor purity, ploidy, and absolute copy number.
  • AneuFinder: For copy number analysis of single-cell sequencing data.

Citations

  • Olshen, A. B., Venkatraman, E. S., Lucito, R., and Wigler, M. (2004). Circular binary segmentation for the analysis of array-based dna copy number data. Biostatistics, 5:557–72.
  • Venkatraman, E. S. and Olshen, A. B. (2007). A faster circular binary segmentation algorithm for the analysis of array cgh data. Bioinformatics, 23:657–63.

References

<!-- biomate-cta -->

Run this on BioMate

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

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

报告此技能的问题