R/utils_internal.R

Defines functions recombine

### Function to colapse into a single data.frame data stored in multiple lists
### Writen by Gustavo Paterno (2017)
recombine <- function(sensi.list, slot1, slot2 = NULL) {
  ### Nulling variables:
  is <- NULL
  
  ### One level indexing list[[slot1]]:
  if (is.null(slot2)) {
    x <- lapply(sensi.list, function(x)
      x[[slot1]])
    
    if (is(x[[1]], "list"))
      stop("Please also provide a value for slot2")
    if (is.na(match(class(x[[1]]), c("data.frame", "matrix"))))
      stop("Check if slot indexing returns a valid data.frame or numeric value withn your list")
    
    res <-
      data.frame(iteration = rep(names(x), sapply(x, nrow)), do.call(rbind, x))
    res$info <- unlist(lapply(x, rownames))
    return(res)
  }
  
  ### Two levels indixing list[[slot1]][[slot2]]
  if (is.null(slot2) == FALSE) {
    x <- lapply(sensi.list, function(x)
      x[[slot1]][[slot2]])
    
    if (inherits(x[[1]], "data.frame") |
        inherits(x[[1]], "matrix") |
        inherits(x[[1]], "numeric")) {
      n.rows <- nrow(x[[1]])
      
      ### If there is only a single value (e.g. AIC, optpar)
      if (is.null(n.rows)) {
        nam <- as.numeric(names(sensi.list))
        n.rows <- length(x[[1]])
        est.name <- names(sensi.list[[1]][[slot1]])[[slot2]]
        N <- rep(nam, each = n.rows)
        res <- data.frame(iteration = N, do.call("rbind", x))
        colnames(res)[2] <- est.name
        return(res)
      }
      
      else
        res <-
          data.frame(iteration = rep(names(x), sapply(x, nrow)), do.call(rbind, x))
      res$info <- unlist(lapply(x, rownames))
      return(res)
    }
    else
      stop("Check if slot indexing returns a valid data.frame or numeric value withn your list")
  }
}

Try the sensiPhy package in your browser

Any scripts or data that you put into this service are public.

sensiPhy documentation built on April 14, 2020, 7:15 p.m.