Merge branch 'gh-pages' of https://github.com/OHI-Science/ohiprep_v2024…

… into gh-pages

globalprep/fis/v2024/STEP5_RAM_CMSY_combine.Rmd

@@ -316,7 +316,7 @@ comparison_plot
 
 For understanding results: 
 
-[When B/BMSY = 1, then biomass equals BMSY. If B/BMSY falls below 1, biomass is too low to provide maximum sustainable yield.](https://catchshareindicators.org/wp-content/uploads/2013/12/NE_Biomass_Feb-2014.pdf) 
+[When B/BMSY = 1, then biomass equals BMSY. If B/BMSY falls below 1, biomass is too low to provide maximum sustainable yield.](https://catchshareindicators.org/wp-content/uploads/2013/12/NE_Biomass_Feb-2014.pdf).  This is also seen in Christopher Costello's 2016 paper, ["Global fishery prospects under contrasting management regimes."](https://www.pnas.org/doi/10.1073/pnas.1520420113) 
 
 ```{r}
 # read in our final data frame saved earlier
@@ -345,24 +345,39 @@ fish_status_summary_2019 <- fish_status %>%
   summarise(count = n()) %>%
   filter(!is.na(fish_status))
 
-# ggplot(fish_status_summary_2019, aes(x = fish_status, y = count, fill = fish_status)) +
-#   geom_bar(position = "dodge", stat = "identity") +
-#   labs(title = "Comparison of Overfished and Sustainable Stocks by Region in 2019",
-#        x = "B/Bmsy Stock Status",
-#        y = "Number of stocks") +
-#   theme_minimal()
-# 
-# fish_status_summary_all_time <- fish_status %>%
-#   group_by(rgn_id, fish_status) %>%
-#   summarise(count = n()) %>%
-#   filter(!is.na(fish_status))
-# 
-# ggplot(fish_status_summary, aes(x = fish_status, y = count, fill = fish_status)) +
-#   geom_bar(stat = "identity") +
-#   labs(title = "Comparison of Overfished and Sustainable Stocks",
-#        x = "B/Bmsy Stock Status",
-#        y = "Number of stocks") +
-#   theme_minimal()
+ggplot(fish_status_summary_2019, aes(x = fish_status, y = count, fill = fish_status)) +
+  geom_bar(position = "dodge", stat = "identity") +
+  labs(title = "Comparison of Overfished and Sustainable Stocks in 2019",
+       x = "B/Bmsy Stock Status",
+       y = "Number of stocks") +
+  theme_minimal()
+
+fish_status_summary_all_time <- fish_status %>%
+  group_by(rgn_id, fish_status) %>%
+summarise(count = n()) %>%
+  filter(!is.na(fish_status))
+
+ggplot(fish_status_summary_all_time, aes(x = fish_status, y = count, fill = fish_status)) +
+  geom_bar(stat = "identity") +
+  labs(title = "Comparison of Overfished and Sustainable Stocks 2012 - 2019",
+       x = "B/Bmsy Stock Status",
+       y = "Number of stocks") +
+  theme_minimal()
+
+# ---- 10 most overfished stocks in 2019 (B/Bmsy < 1) ----
+overfished_2019 <- fish_status %>% 
+  filter(fish_status %in% "overfished") 
+
+# interactive plotly of b/bmsy for each stock
+over_fis_plot <- overfished_2019 %>%
+  plot_ly(x = ~year, y = ~bbmsy, color = ~stock_id, 
+          type = "scatter", mode = "lines") %>%
+  layout(title = "B/Bmsy by Stock ID",
+         xaxis = list(title = "Year"),
+         yaxis = list(title = "B/Bmsy < 1)"))
+over_fis_plot
+
+# ---- all sustainable stocks in 2019 (B/Bmsy > 1) ----
 
 
 # comparison of how many stocks are over/underfished by region

globalprep/ico/v2024/README.md

@@ -0,0 +1,17 @@
+## Ocean Health Index: Iconic Species (Sense of Place)
+
+See full data prep details: ICO [here](http://ohi-science.github.io/ohiprep_v2023/globalprep/ico/v2023/ico_data_prep.html).
+
+More information about the Sense of Place Goal [here](https://oceanhealthindex.org/images/htmls/Supplement.html#69_Sense_of_Place)
+
+## Layers created
+* ico_spp_iucn_status
+* ico_status_calc
+* ico_trend
+
+Our [data managment SOP](https://rawgit.com/OHI-Science/ohiprep_v2021/Reference/SOP_dataOrganization/dataOrganization_SOP.html) describes how we manage OHI global data, including a description of the file structure.
+
+Please see our [citation policy](https://oceanhealthindex.org/global-scores/data-download/) if you use OHI data or methods.
+
+Thank you!
+