Code for analysing SST climatic stability
This folder can be used to recreate the analysis and plots in 'Faster ocean warming threatens richest areas of marine biodiversity'.
To use, simply unzip the folder to your desired location, and then double click on the OceanStability.Rproj file to open the RStudio project. There are two datasets that need to be sourced before you can recreate the analysis.
StableClim: Data available from https://doi.org/10.25909/5ea59831121bc
and
Marine Biogeographic realms: Shapefile available from https://doi.org/10.17608/k6.auckland.5596840.v1
You can then open the Rmarkdown files in RStudio. The files are in Code/, e.g. Code/01_StableClimExtract.Rmd will recreate the first part of the analysis which extracts the climate data from StableClim.
You will need to edit some paths present in the files to get them to work properly.
For example, in 01_StableClimExtract.Rmd, changes need to be made on the following lines:
Line 56: scPath <- "/mnt/c/StableClim", change this to the location of the StableClim data
In 02_sppRichness.Rmd, the location of the marine biogeographic realms needs to be changed:
Line 62: regions <- read_sf("/mnt/c/Shapes/MarineRealmsShapeFile/MarineRealms.shp"), change this to the location where you have extracted the realms.
You can then either knit the documents or run through them interactively