added plot_trans_dates to replace

- plot_sst_transition_dates
- plot_sst_tranition_length

NAMESPACE

@@ -20,8 +20,7 @@ export(plot_seasonal_bt_anomaly_gridded)
 export(plot_seasonal_oisst_anom)
 export(plot_seasonal_sst_anomaly_gridded)
 export(plot_setup)
-export(plot_sst_transition_dates)
-export(plot_sst_transition_dates_length)
+export(plot_trans_dates)
 export(plot_zoo_diversity)
 export(plot_zoo_regime)
 export(theme_facet)

R/plot_habitat_diversity.R

@@ -4,7 +4,7 @@
 #'
 #' @param shadedRegion Numeric vector. Years denoting the shaded region of the plot (most recent 10)
 #' @param report Character string. Which SOE report ("MidAtlantic", "NewEngland")
-#' @param Var Character string. Which Variable to plot ("Diversity","Richness")
+#' @param varName Character string. Which Variable to plot ("Diversity","Richness")
 #'
 #' @return ggplot object
 #'

R/plot_setup.R

@@ -10,7 +10,7 @@
 #'
 #' @export
 #'
-plot_setup <- function(shadedRegion=shadedRegion,
+plot_setup <- function(shadedRegion = shadedRegion,
                        report = report){
 
   #library(sf) # now needed in function for filter to work on sf object? https://github.com/r-spatial/sf/issues/1381
@@ -59,6 +59,10 @@ plot_setup <- function(shadedRegion=shadedRegion,
     x.shade.min <- shadedRegion[1]
     x.shade.max <- shadedRegion[2]
 
+    ## Additional parameters for 4 panel maps eg. plot_seasonal_sst_anomaly_gridded
+    ###################################################
+    ###################################################
+
     #Define constants for figure plot
     series.col <- c("indianred","black")
 

R/plot_trans_dates.R

@@ -0,0 +1,90 @@
+#' plot SST transition dates and timing
+#'
+#' uses ecodata::trans_dates
+#'
+#' @param shadedRegion Numeric vector. Years denoting the shaded region of the plot (most recent 10)
+#' @param report Character string. Which SOE report ("MidAtlantic", "NewEngland")
+#' @param varName Character string. Which Variable to plot ("timing","length")
+#'
+#' @return ggplot object
+#'
+#' @export
+
+plot_trans_dates <- function(shadedRegion = shadedRegion,
+                                      report="MidAtlantic",
+                                      varName = "timing") {
+
+
+  setup <- ecodata::plot_setup(shadedRegion = shadedRegion,
+                               report=report)
+
+  if (report == "MidAtlantic") {
+    filterEPUs <- c("MAB")
+  } else {
+    filterEPUs <- c("GB", "GOM")
+  }
+
+  if (varName == "timing") {
+
+    fix <- ecodata::trans_dates |>
+      dplyr::filter(EPU %in% filterEPUs,
+                    Var %in% c("falltrans", "sprtrans","maxday"),
+                    !Value == "NA",
+                    !Var == "NA",
+                    !Time == "NA") |>
+      dplyr::mutate(Var = dplyr::recode(Var, "falltrans"="Fall",
+                                 "sprtrans" = "Spring",
+                                 "maxday" = "Max"))
+
+    p <-  fix |>
+
+      ggplot2::ggplot(ggplot2::aes(x= Time, y = Value, color = Var)) +
+      ggplot2::annotate("rect", fill = setup$shade.fill, alpha = setup$shade.alpha,
+                        xmin = setup$x.shade.min , xmax = setup$x.shade.max,
+                        ymin = -Inf, ymax = Inf) +
+      ggplot2::geom_point() +
+      ggplot2::geom_line() +
+      ecodata::geom_gls() +
+      # ggplot2::theme(strip.text=ggplot2::element_text(hjust=0),
+      #                plot.title = ggplot2::element_text(size = 12)) +
+      ecodata::theme_title("") +
+      ggplot2::ylab("Day of Year") +
+      ggplot2::facet_wrap(.~EPU) +
+      ecodata::theme_facet() +
+      ggplot2::xlab(ggplot2::element_blank()) +
+    ggplot2::ggtitle(paste0(report,": SST transition dates"))
+
+  } else if (varName == "length") {
+    fix <- ecodata::trans_dates |>
+      dplyr::filter(EPU %in% filterEPUs,
+                    Var == "sumlen",
+                    !Value == "NA")
+    p <- fix |>
+      #dplyr::filter(Var %in% season) %>%
+      ggplot2::ggplot(ggplot2::aes(x= Time, y = Value))+
+      ggplot2::annotate("rect", fill = setup$shade.fill, alpha = setup$shade.alpha,
+                        xmin = setup$x.shade.min , xmax = setup$x.shade.max,
+                        ymin = -Inf, ymax = Inf) +
+      ggplot2::geom_point()+
+      ggplot2::geom_line()+
+      ecodata::geom_gls() +
+      # ggplot2::theme(strip.text=ggplot2::element_text(hjust=0),
+      #                plot.title = ggplot2::element_text(size = 12))+
+      ecodata::theme_title()+
+      ggplot2::ylab("Number of Days")+
+      ggplot2::ggtitle(paste0(report,": Number of days between spring and fall transition dates")) +
+      #ggplot2::xlab(ggplot2::element_blank())+
+      ecodata::theme_ts()+
+      ggplot2::facet_wrap(.~EPU)+
+      ecodata::theme_facet()
+
+
+
+
+
+  } else {
+    stop("Transition date `timing` or `length` allowed")
+  }
+
+  return(p)
+}