R/plotScore.R

In ababaian/RRNAframework: R interface for RNAframework

# plotScore
#' Scatter plot of Score vs. Ratio for each nt in mod.df
#'
#' @param mod.df modification data.frame object
#' @param log.scale Plot log(score, base = 10) instead of linear. Default = T
#' @return geom_hex() object
#' @keywords RNAframework rf-modcall
#' @examples
#' plot <- plotScore(RTstop.wt1)
#' plot <- plot + xlab('log(RT stop score)') + ylab('ratio')
#' plot
#'
#' @seealso \code{\link{geom_hex}}
#' @export
plotScore <- function( mod.df, log.scale = T ){

  score <- as.numeric(unlist(mod.df$score))
  ratio <- as.numeric(unlist(mod.df$ratio))
  log.score <- log( score, base = 10)
  log.ratio <- log( ratio, base = 10)

  plot.df <- data.frame(score, ratio, log.score, log.ratio)

  if (log.scale){
    # Plot Exp Data
    plot <- ggplot( data = plot.df, aes(log.score, ratio) )
    plot <- plot + geom_hex(bins = 100, aes(fill = log(..count.., base = 10)))
    plot <- plot + theme_bw()
    #plot

  } else {
    # Plot Linear Data
    plot <- ggplot( data = plot.df, aes(score, ratio) )
    plot <- plot + geom_hex(bins = 100, aes(fill = log(..count.., base = 10)))
    plot <- plot + theme_bw()
    #plot

    return(plot)
  }
}