Search
Question: Classes and methods for working with segmented copy number data
1
gravatar for Sean Davis
2.1 years ago by
Sean Davis21k
United States
Sean Davis21k wrote:

There are dozens of packages for working with copy number data available in Bioconductor.  However, I don’t get the sense that there is a class that unifies the somewhat standardized “processed” output that consists of:

  • sample_id
  • chromosome
  • start_pos
  • end_pos
  • copy_number

I am interested in applying something like gistic over my data, summarizing to gene, etc. Any suggestions on classes (and associated methods) would be much appreciated.

ADD COMMENTlink modified 17 months ago by Marcel Ramos ♦♦ 290 • written 2.1 years ago by Sean Davis21k
2
gravatar for Marcel Ramos
2.1 years ago by
Marcel Ramos ♦♦ 290
United States
Marcel Ramos ♦♦ 290 wrote:

Hi Sean!

We have the `MultiAssayExperiment::RangedRaggedAssay` class that can handle all of the above.

It would require you to have a `GRangesList` object with a "score" metadata column denoting copy_number. The `GRangesList` can be created from a `DataFrame` using the `GenomicRanges::makeGRangesListFromDataFrame` and subsequently converted to RangedRaggedAssay from `MultiAssayExperiment`. 

I hope this helps and feel free to reach out for more details. We have a run-able example in the package via `example("RangedRaggedAssay")`.

 

Marcel 

ADD COMMENTlink written 2.1 years ago by Marcel Ramos ♦♦ 290
2
gravatar for Marcel Ramos
17 months ago by
Marcel Ramos ♦♦ 290
United States
Marcel Ramos ♦♦ 290 wrote:

Update: Please see the latest package submission (RaggedExperiment) for this type of data. https://github.com/Bioconductor/Contributions/issues/339 and the GitHub version of the package https://github.com/Bioconductor/RaggedExperiment/ 

 

Sincerely, 

Marcel

ADD COMMENTlink written 17 months ago by Marcel Ramos ♦♦ 290

Update: RaggedExperiment is now in Bioconductor devel!

ADD REPLYlink modified 17 months ago • written 17 months ago by Marcel Ramos ♦♦ 290
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 160 users visited in the last hour