This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Copy number profile using Control-FREEC

Hello, I'm using Control FREEC to generate a copy number profile(Normalized copy number) of tumor samples. Although I've generated a profile, I have a doubt. Please see the profile in which red color suggests copy number gain. But I don't know what does each dots represent? Could anyone please describe that? I have used the following commands using R:

getwd() [1] "~/FREEC-11.5/output" ploidy = 2 > ratio <- read.table("Tumor.recal.bam_ratio.txt", header = TRUE) > png(filename = "Tumor_plot.png",width = 1180, height = 1180,units = "px", pointsize = 20, bg = "white", res=NA) > plot(ratio$Start,ratio$Ratioploidy,ylim = c(0,3ploidy),xlab = paste ("Position\n Tumor"),ylab = "Normalised copy number profile", pch = ".", col = "black", cex = 4)
gains <- which(ratio$CopyNumber > ploidy ) > points(ratio$Start[gains], ratio$Ratio[gains]ploidy, pch = ".", col = "blue", cex = 4) > losses <- which(ratio$CopyNumber < ploidy) >points(ratio$Start[losses],ratio$Ratio[losses]ploidy,pch = ".", col = "green", cex = 4)

Please see the left panel of the following figure:

https://www.ncbi.nlm.nih.gov/core/lw/2.0/html/tileshop_pmc/tileshop_pmc_inline.html?title=Click%20on%20image%20to%20zoom&p=PMC3&id=3268243_btr670f1.jpg

control-freec cnv copy number profilr

Figure link: http://%20https://www.ncbi.nlm.nih.gov/core/lw/2.0/html/tileshop_pmc/tileshop_pmc_inline.html?title=Click%20on%20image%20to%20zoom&p=PMC3&id=3268243_btr670f1.jpg

1 answer

As you have used _ratio.txt, each dot represents the predicted copy number for each window.

The size of the window will depend on your input configuration file.

Kevin

Log in to answer this question.