Update both hake eco summary figures; as submitted.

pbs-assess · Jan 29, 2025 · 9fc1894 · 9fc1894
1 parent b02381f
commit 9fc1894
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 12 deletions.
diff --git a/R-temp/ecosystem-summary-hake-glorys.R b/R-temp/ecosystem-summary-hake-glorys.R
@@ -133,10 +133,8 @@ ecosystem_summary_hake_glorys <- function(max_year = 2024,   # though does to 20
   plot(temp_spawn_index, lwd = lwd_index,
        xlim = x_lim,
        xlab = "",
-       ylab = "",
-       ylim = rev(range(temp_spawn_index$anomaly)),
-       y_axis_reverse = TRUE)
-  mtext("Mean temperature during spawning - fish less likely to spawn when higher",
+       ylab = "")
+  mtext("Mean temperature during spawning - when higher fish are less likely spawn but larvae grow quicker*",
         side = 3, adj = 0, cex = 0.7, line = 0.3)
 
   plot(ast_eggs_index, lwd = lwd_index,
@@ -157,7 +155,6 @@ ecosystem_summary_hake_glorys <- function(max_year = 2024,   # though does to 20
   mtext("Mean mixed layer depth during yolk stage (Jan-Apr) - shallower reduces recruitment",
         side = 3, adj = 0, cex = 0.7, line = 0.3)
 
-  # TODO not sure which way
   plot(ssh_jac_index, lwd = lwd_index,
        xlim = x_lim,
        xlab = "",
@@ -168,9 +165,7 @@ ecosystem_summary_hake_glorys <- function(max_year = 2024,   # though does to 20
   plot(mld_late_larv_index, lwd = lwd_index,
        xlim = x_lim,
        xlab = "",
-       ylab = "",
-       ylim = rev(range(mld_late_larv_index$anomaly)),
-       y_axis_reverse = TRUE)
+       ylab = "")
   mtext("Mean mixed layer depth during late larval stage (Mar-Jun) - shallower reduces recruitment",
         side = 3, adj = 0, cex = 0.7, line = 0.3)
 

diff --git a/R-temp/ecosystem-summary-hake.R b/R-temp/ecosystem-summary-hake.R
@@ -115,10 +115,10 @@ ecosystem_summary_hake <- function(max_year = 2024,
   class(npgo_pre_index) <- class(oni)  # so a pacea_index for plotting
 
 
-  # shift years by 1 like herring, take log (as in Vestfals)  and then restandardise
+  # no shifting, take log (as in Vestfals)  and then restandardise
   hake_log_age1_index <- dplyr::filter(hake_total_biomass_age_1,
-                              year >= min_year - 1) %>%
-    dplyr::mutate(year = year + 1,
+                              year >= min_year) %>%
+    dplyr::mutate(year = year,
                   log_median = log(median),
                   anomaly = standardise(log_median))
   class(hake_log_age1_index) <- class(oni)
@@ -181,7 +181,7 @@ ecosystem_summary_hake <- function(max_year = 2024,
        ylim = rev(range(hake_log_age1_index$anomaly)),
        y_axis_reverse = TRUE)
 
-  mtext("Hake total log biomass of age-1 fish - predation on age-0 fish",
+  mtext("Hake total log biomass of age-1 fish - more predation on age-0 fish lowers recruitment that year",
         side = 3, adj = 0, cex = 0.7, line = 0.3)
 
   mtext("Year of hake recruitment", side = 1, line = 2)