Question: How to convert segmented aCGH probe-based to gene-based?
gravatar for sjorsvanheuveln
5.0 years ago by
sjorsvanheuveln50 wrote:

I'm looking for the easiest method in R to convert my segmented aCGH probes into log ratios per gene (they're not called). Furthermore I only have chromosome number and start position available. What is the best way to achieve this goal?


R conversion probe acgh gene • 1.1k views
ADD COMMENTlink modified 4.0 years ago by Biostar ♦♦ 20 • written 5.0 years ago by sjorsvanheuveln50
gravatar for arno.guille
5.0 years ago by
arno.guille400 wrote:


Here is my solution:

I use bedtools intersect to find overlaps between segdata output and genes coordinates (refseq from ucsc).

bedtools intersect -a segData.bed -b genes.bed -loj > overlap.bed

Then with R

res = read.table("overlap.bed",sep="\t",header=F)
abs_max = tapply(res$V4,as.character(res$V8),function(x){ x[which.max(abs(x))] })

As you can see, I take the maximum in absolute value

res$V4 correspond to the log2ratio value
res$V8 correspond to the Symbol

If you need more help, let me know.

ADD COMMENTlink modified 10 months ago by RamRS30k • written 5.0 years ago by arno.guille400

Thanks a million Arno, this worked for me! 

ADD REPLYlink written 5.0 years ago by sjorsvanheuveln50
Please log in to add an answer.


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