Question: Rawcopy copy number analysis
0
gravatar for moritz.feigl
8 months ago by
moritz.feigl0 wrote:

Hallo,

I used the R rawcopy package to analyze my affymetrix SNP 6.0 data. This produced the segmentation and the log ratios and allelic imbalance values. Now I'm trying to analyze those, but I can't find a suitable way. I tried the packages that are stated in the rawcopy paper for further analysis and failed to apply any of them. This was mainly because all of them have practical no useful documentation. Or at least I didn't find one. Right now I started trying to calculate the gains and losses with a simple threshold method myself, but even this is tricky because I'm not sure if I use the correct outputs from rawcopy... If anyone has experience with the rawcopy output and its further usage, I would be extremely happy to get some advice.

copy number rawcopy • 260 views
ADD COMMENTlink modified 8 weeks ago by hu199301210 • written 8 months ago by moritz.feigl0

Is this problem resolved ?

ADD REPLYlink modified 6 months ago • written 6 months ago by sofie_carolina20

Hi jgarce recently, I have run the rawcopy program on my snp6.0 data, and i don't understand the segmentation output(segments.txt). could it be the input data for GISTIC2.0, but it still missing the n.probes values for input data.
follow your sugggestion, i used copynumber package for additional segmentation (use the segments.txt as the input). finally, the segments was achieved, but the segments is so litter. could be something wrong.

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by hu199301210
0
gravatar for jgarces
9 weeks ago by
jgarces10
Spain
jgarces10 wrote:

Hi Moritz,

I'll try to answer you with my (little) experience with rawcopy (and after some emails with the package's author, he's very kind):

  • The correct output, I guess you mean the segmentation output, is the segments.txt file (the other ones are for generating the different plots for the QC).
  • A possible method for establishing a threshold could be those contained within GISTIC2... it performs a double-step procedure: 1) it calculates a significance according the entire genome as background (this'd be the value you're searching) and 2) it searches for recurrent alterations among all your samples and gives them a score.
  • And finally, after some (and not so difficult) code, you can use copynumber package for plotting (and additional segmentation if you don't like rawcopy's one.

I hope these tips help you a bit. Feel free to ask whatever you want. Bests.

ADD COMMENTlink written 9 weeks ago by jgarces10
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.3.0
Traffic: 1502 users visited in the last hour