R/plot_Freq_samps.R

In jafernandez01/cndcR: Critical Nitrogen Dilution Curves in R

#' @name fig3b
#' @title
#' Plot figure 3b from Fernandez et al. (2022)
#' @description
#' This function returns figure 3b in Fernandez et al.(2022): Frequency distribution of number of
#' samplings used for fitting CNDC in literature.

fig3b <- function() {

  figFreq_samps <- cndcR::biblioCNDC %>%
  dplyr::mutate(min = as.numeric(sub("to.*","",.data$nSamp)),
         max = as.numeric(sub(".*to","",.data$nSamp)),
         range = max - min) %>%
  dplyr::mutate(typeStudy = dplyr::case_when(!is.na(range) & range == 0 ~ "Single sampling-design",
                                             is.na(range) ~ NA_character_,
                                             TRUE ~ "Multiple sampling-designs"  )) %>%
  tidyr::drop_na(.data$typeStudy) %>%
  dplyr::mutate(nSamp = dplyr::if_else(.data$typeStudy == "Single sampling-design", min,
                                       min + round((range/2),0 ) )) %>%
  ggplot2::ggplot(ggplot2::aes(x = .data$nSamp)) +
  ggplot2::geom_histogram(ggplot2::aes(fill = .data$typeStudy),binwidth = 2, boundary = 0,
                          color = "black") +

  ggplot2::stat_bin(binwidth=2, geom="text", colour="black", size=3.5, fontface = 3,
           ggplot2::aes(label=.data$..count.., y=(.data$..count..)+2 ), boundary = 0) +
  ggplot2::ylab("Number of CNDCs") +
  ggplot2::ggtitle("Number of sampling times") + ggplot2::xlab(NULL) +
  ggplot2::scale_fill_manual(values = c("#3d3c3c","#c4c2c2")) +
  ggplot2::scale_x_continuous(breaks = seq(0,20,2)) +
  ggplot2::scale_y_continuous(limits = c(0,50)) +

  ggplot2::theme(panel.grid.minor = ggplot2::element_blank(),
        panel.grid.major = ggplot2::element_line(color = "#e8e9eb"),
        panel.background = ggplot2::element_rect(fill = "#f5f5f5"),
        panel.border = ggplot2::element_rect(colour = "black", fill = NA),

        text = ggplot2::element_text(size = 14),
        legend.title = ggplot2::element_blank(),
        legend.position = c(0.7,0.8),
        legend.background = ggplot2::element_blank(),

        axis.ticks.length=ggplot2::unit(-0.15, "cm"),
        axis.text.x = ggplot2::element_text(margin=ggplot2::unit(c(0.5,0.5,0.5,0.5), "cm")),
        axis.text.y = ggplot2::element_text(margin=ggplot2::unit(c(0.5,0.5,0.5,0.5), "cm"),
                                            size = 13))

  return(figFreq_samps)
}