R/graphics.R

In Anaquin: Statistical analysis of sequins

#
#  Copyright (C) 2017 - Garvan Institute of Medical Research
#

.m2str <- function(m)
{
    eq <- substitute(italic(y) == a + b * italic(x)*','~~italic(r)^2~'='~r2, 
                     list(a  = format(coef(m)[1], digits = 2), 
                          b  = format(coef(m)[2], digits = 2), 
                          r2 = format(summary(m)$r.squared, digits = 3)))
    as.character(as.expression(eq));
}

.lm2str <- function(data)
{
    return (.m2str(lm(y~x, data)))
}

.transformPlot <- function(p)
{
    p <- p + theme(axis.title.x=element_text(face='bold', size=12))
    p <- p + theme(axis.title.y=element_text(face='bold', size=12))
    p <- p + theme(panel.border=element_rect(colour="black", fill=NA, size=1))
    p <- p + theme(plot.title=element_text(size=15, face="bold"))
    p <- p + theme(legend.title=element_text(face="bold"))
    p <- p + theme(legend.key=element_blank())

    #
    # Calculating the aspect ratio to make the plot a square.
    #
    
    build <- ggplot_build(p)
    
    minX <- build$layout$panel_ranges[[1]]$x.range[[1]]
    maxX <- build$layout$panel_ranges[[1]]$x.range[[2]]
    minY <- build$layout$panel_ranges[[1]]$y.range[[1]]
    maxY <- build$layout$panel_ranges[[1]]$y.range[[2]]
    
    if (is.null(minX) && is.null(maxX) && is.null(minY) && is.null(maxY))
    {
        minX <- build$layout$panel_params[[1]]$x.range[[1]]
        maxX <- build$layout$panel_params[[1]]$x.range[[2]]
        minY <- build$layout$panel_params[[1]]$y.range[[1]]
        maxY <- build$layout$panel_params[[1]]$y.range[[2]]
    }
    
    stopifnot(!is.null(minX))
    stopifnot(!is.null(maxX))
    stopifnot(!is.null(minY))
    stopifnot(!is.null(maxY))    
    
    xrange <- maxX - minX
    yrange <- maxY - minY

    p <- p + coord_fixed(ratio=xrange/yrange)
    
    return (p)    
}