R/plot_qpcr_target.R

In J9DM/j9lablib: What the Package Does (One Line, Title Case)

#' QC plot for qPCR data
#'
#' @param df input data set as tibble, normalized to control condition
#' @param facet_vars variables for faceting, quoted
#' @param color_by variable for coloring, quoted
#' @param pch_by variable for plotting character, quoted
#'
#' @importFrom rlang .data
#'
#' @return A plot of 2^(CT_ctrl - CT) for each target and condition
#' @export
#'
# @examples
plot_qpcr_target <- function(df, facet_vars = c("target", "condition"),
                             color_by, pch_by) {

  facet_vars <- rlang::syms(facet_vars)
  color_by <- rlang::sym(color_by)
  pch_by <- rlang::sym(pch_by)
  ctrl_condition <- df$ctrl %>% unique()

  df %>%
    ggplot2::ggplot(
      ggplot2::aes(x = .data$replicate, y = .data$fold_ctrl)) +
    ggplot2::facet_grid(
      ggplot2::vars(!! facet_vars[[1]]),
      ggplot2::vars(!! facet_vars[[2]])) +
    ggplot2::geom_hline(yintercept = 1, color = "red") +
    ggplot2::geom_boxplot() +
    ggplot2::geom_point(
      ggplot2::aes(color = !! color_by, pch = !! pch_by)) +
    ggplot2::ylab(paste0("Fold change over ", ctrl_condition))
}