R/PlotParametricSpace.R

In KamilSJaron/ConjunctionStats: The parser of simulations of hybrid zones performed by simulator Conjunction (http://github.com/KamilSJaron/Conjunction)

#' @title PlotParametricSpace
#'
#' @description
#' \code{PlotParametricSpace} plots range of selections, betas, lambdas and
#' number of chromosomes in one set of simulations
#'
#' @param GradTable
#'
#' @param global_margins set margins around plot
#'
#' @param denisity_haching set density of lines in rectangles
#'
#' @param cex is size of all letters in plot (2 or 3 are reasonable)
#'
#' @author Kamil Jaron \email{kamiljaron at gmail.com}
#'
#' @export

PlotParametricSpace <- function(GradTable, global_margins = 0.1, # in inches
                                denisity_haching = 10, cex = 2){

  par(mfrow=c(4,1), tcl=-0.5, family="serif", las=1, omi=c(global_margins,
                                                           global_margins,
                                                           global_margins,
                                                           global_margins))
  # SELECTION (min, max, by)
  PlotParametricTemplate(c(0,1), 0.05, cex)
  rect(min(GradTable$s), 0, max(GradTable$s), 1, density = denisity_haching)
  mtext("s", side=2, line=0.5, cex = cex)

  # BETA (min, max, by)
  PlotParametricTemplate(c(-4, 4), NA, cex)
  beta_axis <- 2^(seq(-4,4,1))
  axis(1, labels=beta_axis, at=log2(beta_axis), lwd=0, lwd.ticks=1, cex.axis=cex)
  if(min(GradTable$b) == max(GradTable$b)){ext_width = 0.03} else {ext_width = 0}
  rect(log2(min(GradTable$b)) - ext_width, 0, log2(max(GradTable$b)) + ext_width, 1,
       density = denisity_haching)
  mtext(expression(beta), side=2, line=0.5, cex = cex)

  # LAMBDA (min, max, by)
  PlotParametricTemplate(c(0, 2), 0.2, cex)
  if(min(GradTable$r) == max(GradTable$r)){ext_width = 0.02} else {ext_width = 0}
  rect(min(GradTable$r) - ext_width, 0 + ext_width, max(GradTable$r), 1,
       density = denisity_haching)
  mtext(expression(lambda), side=2, line=0.5, cex = cex)

  # CHROMOSOMES (min, max, by)
  PlotParametricTemplate(c(1, 20), NA, cex)
  axis(1, at = c(1,5,10,20), lwd=0, lwd.ticks=1, cex.axis=cex)
  # -0.1, +0.1 is defined for situation of one simulated value only
  if(min(GradTable$C) == max(GradTable$C)){ext_width = 0.1} else {ext_width = 0}
  rect(min(GradTable$C) - ext_width, 0, max(GradTable$C) + ext_width, 1,
       density = denisity_haching)
  mtext("C", side=2, line=0.5, cex = cex)
}