Skip to content

Commit

Permalink
updating readmes to reflect new model names
Browse files Browse the repository at this point in the history
  • Loading branch information
melofton committed Aug 12, 2021
1 parent dd92172 commit eb24cba
Show file tree
Hide file tree
Showing 3 changed files with 12 additions and 11 deletions.
4 changes: 2 additions & 2 deletions 4.1_JAGS_models/readme.txt
Original file line number Diff line number Diff line change
Expand Up @@ -6,15 +6,15 @@ Refer to Table 2 and Table S4 of the Lofton et al. 2021 manuscript for descripti

AC = null model with autocorrelation term; corresponds to AC model in manuscript

base_DLM = null dynamic linear model; correspond to BaseDLM model in manuscript
base_DLM = null dynamic linear model; correspond to BaseLM model in manuscript

GDD = GDD model in manuscript

precip = Precip model in manuscript

RW = random walk model with uninformed observation error prior; not presented in manuscript

RW_bias = random walk null model with bias term; correspond to BiasRW in manuscript
RW_bias = random walk null model with bias term; correspond to OffsetRW in manuscript

RW_obs = random walk model with informed observation error prior; corresponds to RW model in manuscript

Expand Down
4 changes: 2 additions & 2 deletions 5_Model_output/readme.txt
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@ There are between 80-120 output files for each Bayesian state-space model, depen
The names of the .csv files for each model follow the naming convention below:
[model name]_[sources of uncertainty included in hindcast]_[year of hindcast]_[week of year].csv
For example, the following file:
AR_det.prediction_2015_1.csv contains the 1-4 week hindcasts generated for weeks 1-4 of the 2015 sampling season for the AR model, with no forms of uncertainty included
AC_det.prediction_2015_1.csv contains the 1-4 week hindcasts generated for weeks 1-4 of the 2015 sampling season for the AC model, with no forms of uncertainty included
A second example:
GDD_hindcast.IC.Pa.D.P._2016_12.csv contains the 1-4 week hindcasts generated for weeks 12-15 of the 2016 sampling season for the GDD model, and includes initial conditions, parameter, driver, and process uncertainty

Expand Down Expand Up @@ -46,7 +46,7 @@ There are between 80-120 output files for each Bayesian state-space model, depen
The names of the .csv files for each model follow the naming convention below:
[model name]_[sources of uncertainty included in hindcast]_[year of hindcast]_[week of year].csv
For example, the following file:
AR_det.prediction_2015_1.csv contains the 1-4 week hindcasts generated for weeks 1-4 of the 2015 sampling season for the AR model, with no forms of uncertainty included
AC_det.prediction_2015_1.csv contains the 1-4 week hindcasts generated for weeks 1-4 of the 2015 sampling season for the AC model, with no forms of uncertainty included
A second example:
GDD_hindcast.IC.Pa.D.P._2016_12.csv contains the 1-4 week hindcasts generated for weeks 12-15 of the 2016 sampling season for the GDD model, and includes initial conditions, parameter, driver, and process uncertainty

15 changes: 8 additions & 7 deletions 9_Data_publication/9A_rename_hindcast_files.R
Original file line number Diff line number Diff line change
Expand Up @@ -3,22 +3,23 @@
#Date: 25MAY20
library(tidyverse)

#get sampling dates
obs <- read_csv("./00_Data_files/Covariate_analysis_data/HC_Gechinulata_long.csv") %>%
filter(year %in% c(2015:2016))
dates <- as.character(obs$date)
new_wks <- c("2015-05-14",dates[1:19],"2016-05-19",dates[21:39])

#reset working directory
setwd("./9_Data_publication/")

startingDir <- "./5_Model_output/5.2_Hindcasting/"

#create wks vector
wks <- c(paste0("2015_",1:20),paste0("2016_",1:20))

#get sampling dates
obs <- read_csv("./00_Data_files/Covariate_analysis_data/HC_Gechinulata_long.csv") %>%
filter(year %in% c(2015:2016))
dates <- as.character(obs$date)
new_wks <- c("2015-05-14",dates[1:19],"2016-05-19",dates[21:39])

#create vector of model names
model_names <- c("RW_obs","RW_bias","AC","base_DLM","wtrtemp_min","wtrtemp_min_lag","wtrtemp_MA7","GDD","wnd_dir_2day_lag","schmidt_max_lag","schmidt_and_wind","temp_and_wind","wind_and_GDD")
new_model_names <- c("RW","BiasRW","AC","BaseDLM","MinWaterTemp","MinWaterTempLag","WaterTempMA","GDD","WindDir","SchmidtLag","SchmidtAndWind","TempAndWind","WindAndGDD")
new_model_names <- c("RW","OffsetRW","AC","BaseLM","MinWaterTemp","MinWaterTempLag","WaterTempMA","GDD","WindDir","SchmidtLag","SchmidtAndWind","TempAndWind","WindAndGDD")

for(n in 13:length(model_names)){

Expand Down

0 comments on commit eb24cba

Please sign in to comment.