feat: get_saved_model_runs for new MSE run save pattern

Author: Ovec8hkin

R/utils.R

@@ -168,8 +168,8 @@ generate_annual_frequency <- function(frequency, len){
 #'
 #' @example
 #'
-get_saved_model_runs <- function(om_order=NULL, hcr_order=NULL){
-    fs <- list.files(file.path(here::here(), "data", "active"), full.names = TRUE)
+get_saved_model_runs_old <- function(om_order=NULL, hcr_order=NULL){
+    fs <- list.files(file.path(here::here(), "data", "mse_runs_good"), full.names = TRUE)
     model_runs <- unlist(lapply(fs, function(x){
         m <- readRDS(x)
         mse <- m$mse_objects
@@ -202,6 +202,58 @@ get_saved_model_runs <- function(om_order=NULL, hcr_order=NULL){
 
 }
 
+#' Load Saved MSE Model Runs from Disk
+#' 
+#' Read all saved RDS files present in data/active and coerce into proper
+#' model_runs list object. Also setup correctly specified extra_columns
+#' object for use with bind_mse_outputs. 
+#'
+#' @param om_order vector of correct order of OMs (used to set OM factor level)
+#' @param hcr_order vector of correct order of HCRs (used to set HCR factor level)
+#'
+#' @export get_saved_model_runs
+#'
+#' @example
+#'
+get_saved_model_runs <- function(om_order=NULL, hcr_order=NULL){
+    fs <- list.files(file.path(here::here(), "data", "active"), full.names = TRUE)
+    model_runs <- lapply(seq_along(fs), function(i){
+        x <- fs[i]
+        m <- readRDS(x)
+        mse <- m$mse_objects
+        mse[[length(mse)]]
+    })
+
+    om_names <- sapply(fs, function(x){
+        m <- readRDS(x)
+        m$om$name
+    })
+    names(om_names) <- NULL
+
+    hcr_names <- sapply(fs, function(x){
+        m <- readRDS(x)
+        m$hcr$name
+    })
+    names(hcr_names) <- NULL
+
+    extra_columns2 <- data.frame(om=om_names, hcr=hcr_names)
+
+    if(!is.null(om_order))
+        extra_columns2$om <- factor(extra_columns2$om, levels=om_order)
+    
+    if(!is.null(hcr_order))
+        extra_columns2$hcr <- factor(extra_columns2$hcr, levels=hcr_order)
+
+    return(listN(model_runs, extra_columns2))
+
+}
+
+
+
+
+
+
+
 #' Get maximum value without considering infinite values
 #' 
 #' Wrapper around max that ignores infinite values