R/create_figure_alignment_qualities.R

#' Create figure 'alignment_qualities'
#' @param alignments data frame with alignments, as returned by 'read_collected_alignments'
#' @param filename the name of the file the figure is saved to
#' @return nothing, it saves a file
#' @author Richel Bilderbeek
#' @export
create_figure_alignment_qualities <- function(
  alignments,
  filename
) {

  status <- NULL; rm(status) # nolint, should fix warning: no visible binding for global variable
  n <- NULL; rm(n) # nolint, should fix warning: no visible binding for global variable

  df <- data.frame(
    status = c("ok", "zero", "na"),
    n = c(
      sum(alignments$n_alignments_ok),
      sum(alignments$n_alignments_zeroes),
      sum(alignments$n_alignments_na)
    )
  )
  df$status <- as.factor(df$status)

  ggplot2::ggplot(
    data = df,
    ggplot2::aes(x = status, y = n, fill = status)
  ) + ggplot2::geom_col() +
    ggplot2::ylab(label = "count") +
    ggplot2::labs(
      fill = "quality",
      title = "Alignment qualities",
      caption  = "Figure 220"
    ) +
    ggplot2::theme(plot.title = ggplot2::element_text(hjust = 0.5))

  ggplot2::ggsave(file = filename, width = 7, height = 7)
}
richelbilderbeek/wiritttea documentation built on May 27, 2019, 8:02 a.m.