scripts/05_extract_temp_data.R

In lorecatta/DENVclimate:

# rescale predictor values as you need degree celsius in R0 formula


devtools::load_all()


# define parameters -----------------------------------------------------------


parameters <- list(resample_grid_size = 20,
                   year.i = 2007,
                   year.f = 2014,
                   ppyear = 64)

out_path <- "output"

out_name <- "foi_data_cov_rescaled"


# define variables

# pre processing --------------------------------------------------------------


my_predictors <- predictor_rank$name


# scale covariate values before using R0 model --------------------------------


for (i in seq_along(my_predictors)){

  my_pred <- my_predictors[i]

  scale <- get_covariate_scaling_factor(my_pred, parameters)

  message(scale)

  foi_data[, my_pred] <- foi_data[, my_pred] / scale

}


# save ------------------------------------------------------------------------


write_out_rds(foi_data, out_path, out_name)