ShixiangWang / sigminer

🌲 An easy-to-use and scalable toolkit for genomic alteration signature (a.k.a. mutational signature) analysis and visualization in R https://shixiangwang.github.io/sigminer/reference/index.html
https://shixiangwang.github.io/sigminer/
Other
141 stars 18 forks source link

show_cn_freq_circos() y axis #416

Closed rajeha closed 1 year ago

rajeha commented 1 year ago

Hi @ShixiangWang,

I hope you are doing well! Would you please clarify what the y-axis of the show_cn_freq_circos() function output represents? It is not clear in the documentation.

Thank you!

ShixiangWang commented 1 year ago

@rajeha

This function draws copy number amplification and deletion across the genome, you can check the example code. For a specific region, the height shows the proportion with AMP/DEL in all samples. If you provide a group list, each group will be drawn and the ratio will be calculated separately. (The reference, i.e., normal copy number is 2, >2 is defined as AMP, <2 is defined as DEL)

For example code:

load(system.file("extdata", "toy_copynumber.RData",
  package = "sigminer", mustWork = TRUE
))

show_cn_freq_circos(cn)

image

ss <- unique(cn@data$sample)
show_cn_freq_circos(cn, groups = list(a = ss[1:5], b = ss[6:10]), cols = c("red", "green"))

image

rajeha commented 1 year ago

Super helpful - thanks!