Search
Question: Classes and methods for working with segmented copy number data
1
gravatar for Sean Davis
22 months 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 14 months ago by Marcel Ramos ♦♦ 140 • written 22 months ago by Sean Davis21k
2
gravatar for Marcel Ramos
22 months ago by
Marcel Ramos ♦♦ 140
United States
Marcel Ramos ♦♦ 140 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 22 months ago by Marcel Ramos ♦♦ 140
2
gravatar for Marcel Ramos
14 months ago by
Marcel Ramos ♦♦ 140
United States
Marcel Ramos ♦♦ 140 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 14 months ago by Marcel Ramos ♦♦ 140

Update: RaggedExperiment is now in Bioconductor devel!

ADD REPLYlink modified 14 months ago • written 14 months ago by Marcel Ramos ♦♦ 140
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: 177 users visited in the last hour