R/stitchData.R

In castels/interpTools: Tools for the Systematic Testing of Interpolation Algorithms

#' Stitch Data 
#' 
#' Function to stitch together data
#' 
#' @param directory The directory location of the .rda files to stitch. Must be devoid of other files, and each .rda must be identical in structure and follow the same naming convention. 


stitchData <- function(directory){
  
  if(!all(grepl(".rda", list.files(directory)))) stop(c("Directory must only contain .rda files that are identical in structure. Please relocate or delete the following files: ",
                                                             
                                                             paste(list.files(directory)[which(!grepl(".rda",list.files(directory)))], 
                                                                   collapse = ", ")))
  
  files <- list.files(directory, pattern=".rda", full.names = TRUE, recursive = TRUE)
  split <- length(files)
  
  fileList <- list()
  for(i in 1:split){
  fileList[[i]] <- load_obj(files[i]) 
  }
  
  names(fileList) <- letters[1:split]
  
  M <- length(fileList[[1]])
  P <- length(fileList[[1]][[1]])
  G <- length(fileList[[1]][[1]][[1]])
  K <- length(fileList[[1]][[1]][[1]][[1]])*split
    
  stitchList <-  lapply(stitchList <- vector(mode = 'list',M), function(x)
    lapply(stitchList <- vector(mode = 'list', P), function(x) 
        x<-vector(mode='list',G)))

    for(m in 1:M){
      for(p in 1:P){
        for(g in 1:G){
          for(i in 1:split){
          stitchList[[m]][[p]][[g]] <- append(stitchList[[m]][[p]][[g]],fileList[[i]][[m]][[p]][[g]])
          }
          names(stitchList[[m]][[p]]) <- names(fileList[[i]][[m]][[p]])
        }
        names(stitchList[[m]]) <- names(fileList[[i]][[m]])
      }
      names(stitchList) <- names(fileList[[i]])
  }

  return(stitchList)
  }