R/zoom_highlight.R

In bschilder/echoplot: echoverse module: Locus plot creation for fine-mapping and colocalization studies

zoom_highlight <- function(gg,
                           dat,
                           zoom="5x",
                           alpha=.15,
                           verbose=TRUE){
    requireNamespace("ggplot2")
    x <- y <- NULL;
    
    messager("+ Highlighting zoom origin...",v=verbose)
    genomic_units <- guess_genomic_units(gg = gg)
    xlims_zoom <- echoannot::get_window_limits(dat = dat,
                                    zoom = zoom,
                                    genomic_units = genomic_units,
                                    verbose = FALSE)
    rect_dat <- data.frame(x=c(xlims_zoom, rev(xlims_zoom)), 
                           y=c(-Inf,-Inf,Inf,Inf))
    gg2 <- gg +
        ggplot2::geom_polygon(data = rect_dat,
                              ggplot2::aes(x=x, y=y),
                              color="transparent",
                              alpha=alpha,
                              fill="blue")
    return(gg2)
}