Colocalization plot

[sahirbhatnagar]

Hi,
Do you think it's possible to use the existing topr functions to create a colocalization plot? Specifically, multiple locuszoom plots are superimposed on each other with a shared gene track at the bottom. I found a solution, but it's quite hacky. The locuszooms are using standard ggplot2 functions, the gene track is extracted from topr::locuszoom with extract_plots=TRUE, and I combine everything with ggbio::tracks.

[totajuliusd]

Hi, you cannot do that using the current version of topr. You can plot two datasets together with the ntop argument set to 1, like this:

locuszoom(list(R2_CD_UKBB,R2_CD_UKBB), ntop=1, annotate_with_vline = 1e-6)

But then one of the datasets is displayed upside down, and the legend_label is missing. The locuszoom function in topr has been neglected a bit and I already have some improvements lined-up for the next release. I like the colocalization plot you did and it would be nice to include something like that in the next topr release?