本文へ移動
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(原文)

インストールする前に、エージェントに与えられる指示の中身を確認できます。

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 →

レビュー

まだレビューはありません。使ってみた感想をお寄せください。

同じリポジトリのスキル

概要と使いどころ

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

bioMate-AI のスキルをすべて見る

このスキルの問題を報告する