R/exit.list.r

In PABalland/EconGeo: Computing Key Indicators of the Spatial Distribution of Economic Activities

#' Generate a data frame of exit events from multiple regions - industries matrices (same matrix composition for the different periods)
#'
#' This function generates a data frame of exit events from multiple regions - industries matrices (different matrix compositions are allowed). In this function, the maximum number of periods is limited to 20.
#' @param mat1 An incidence matrix with regions in rows and industries in columns (period 1 - mandatory)
#' @param mat2 An incidence matrix with regions in rows and industries in columns (period 2 - mandatory)
#' @param mat... An incidence matrix with regions in rows and industries in columns (period ... - optional)
#' @keywords diversification
#' @export
#' @examples
#' ## generate a first region - industry matrix in which cells represent the presence/absence
#' ## of a RCA (period 1)
#' set.seed(31)
#' mat1 <- matrix(sample(0:1,20,replace=T), ncol = 4)
#' rownames(mat1) <- c ("R1", "R2", "R3", "R4", "R5")
#' colnames(mat1) <- c ("I1", "I2", "I3", "I4")
#'
#' ## generate a second region - industry matrix in which cells represent the presence/absence
#' ## of a RCA (period 2)
#' mat2 <- mat1
#' mat2[2,1] <- 0
#'
#' ## run the function
#' exit.list (mat1, mat2)
#'
#' ## generate a third region - industry matrix in which cells represent the presence/absence
#' ## of a RCA (period 3)
#' mat3 <- mat2
#' mat3[5,1] <- 0
#'
#' ## run the function
#' exit.list (mat1, mat2, mat3)
#'
#' ## generate a fourth region - industry matrix in which cells represent the presence/absence
#' ## of a RCA (period 4)
#' mat4 <- mat3
#' mat4[5,3] <- 0
#'
#' ## run the function
#' exit.list (mat1, mat2, mat3, mat4)
#' @author Pierre-Alexandre Balland \email{p.balland@uu.nl} \cr
#' Wolf-Hendrik Uhlbach \email{w.p.uhlbach@students.uu.nl}
#' @seealso \code{\link{entry}}, \code{\link{exit}}, \code{\link{entry.list}}
#' @references Boschma, R., Balland, P.A. and Kogler, D. (2015) Relatedness and Technological Change in Cities: The rise and fall of technological knowledge in U.S. metropolitan areas from 1981 to 2010, \emph{Industrial and Corporate Change} \strong{24} (1): 223-250 \cr
#' \cr
#' Boschma, R., Heimeriks, G. and Balland, P.A. (2014) Scientific Knowledge Dynamics and Relatedness in Bio-Tech Cities, \emph{Research Policy} \strong{43} (1): 107-114

"exit.list" <- function (mat1, mat2, mat3, mat4, mat5, mat6, mat7, mat8,
                                mat9, mat10, mat11, mat12, mat13, mat14, mat15,
                                mat16, mat17, mat18, mat19, mat20)
{
  a <- list(mat1, if (missing(mat2)) NA else mat2,
            if (missing(mat3)) NA else mat3,
            if (missing(mat4)) NA else mat4,
            if (missing(mat5)) NA else mat5,
            if (missing(mat6)) NA else mat6,
            if (missing(mat7)) NA else mat7,
            if (missing(mat8)) NA else mat8,
            if (missing(mat9)) NA else mat9,
            if (missing(mat10)) NA else mat10,
            if (missing(mat11)) NA else mat11,
            if (missing(mat12)) NA else mat12,
            if (missing(mat13)) NA else mat13,
            if (missing(mat14)) NA else mat14,
            if (missing(mat15)) NA else mat15,
            if (missing(mat16)) NA else mat16,
            if (missing(mat17)) NA else mat17,
            if (missing(mat18)) NA else mat18,
            if (missing(mat19)) NA else mat19,
            if (missing(mat20)) NA else mat20)

  cols <- sort(unique(unlist(lapply(a, colnames))))
  rows <- sort(unique(unlist(lapply(a, rownames))))
  unimat <- array(0, dim=c(length(rows), length(cols)), dimnames=list(rows,cols))
  unimat

  {
    n = 2
    nam <- paste("bim", 1, sep = "")
    object = get(paste0("mat", n - 1))

    unimat1 <- unimat
    col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
    row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
    unimat1[row, col] <- object[row, col]
    object <- unimat1

    assign(nam, object)
    nam2 <- paste("bim", 2, sep = "")
    object = get(paste0("mat", n))

    unimat2 <- unimat
    col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
    row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
    unimat2[row, col] <- object[row, col]
    object <- unimat2

    assign(nam2, object)
    bim1[bim1 == 0] <- NA
    exit = bim1 + bim2
    exit[exit == 2] <- 0
    exit = get.list(exit)
    colnames(exit) <- c("region", "industry", "exit")
    exit$period <- n
    if (missing(mat3)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat4)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat5)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat6)) {
      return(exit)

    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat7)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat8)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat9)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat10)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat11)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat12)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat13)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat14)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat15)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat16)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat17)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat18)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat19)) {
      return(exit)
    }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
    }
    if (missing(mat20)) {
      return(exit)

      }
    else {
      n = n + 1
      nam <- paste("bim", 1, sep = "")
      object = get(paste0("mat", n - 1))

      unimat1 <- unimat
      col <- colnames(unimat1)[colnames(unimat1) %in% colnames(object)]
      row <- rownames(unimat1)[rownames(unimat1) %in% rownames(object)]
      unimat1[row, col] <- object[row, col]
      object <- unimat1

      assign(nam, object)
      nam2 <- paste("bim", 2, sep = "")
      object = get(paste0("mat", n))

      unimat2 <- unimat
      col <- colnames(unimat2)[colnames(unimat2) %in% colnames(object)]
      row <- rownames(unimat2)[rownames(unimat2) %in% rownames(object)]
      unimat2[row, col] <- object[row, col]
      object <- unimat2

      assign(nam2, object)
      bim1[bim1 == 0] <- NA
      exit2 = bim1 + bim2
      exit2[exit2 == 2] <- 0
      exit2 = get.list(exit2)
      colnames(exit2) <- c("region", "industry", "exit")
      exit2$period <- n
      exit <- rbind(exit, exit2)
      return(exit)
    }
  }
}