diff --git a/R-temp/ecosystem-summary-hake-glorys.R b/R-temp/ecosystem-summary-hake-glorys.R index b0e9082c..6f1f5c26 100644 --- 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 index 4fec9c39..d7869a53 100644 --- 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)