R/LQ.R

In jjpwade/ioanalysis: Input-Output Analysis

LQ <- function(io){
  # Preliminaries
  if(class(io) != "InputOutput") stop("io must be of class InputOutput. See ?as.inputoutput")
  A <- io$A
  X <- io$X
  RS_label <- io$RS_label
  # Updating
  Xn <- sum(X)
  regions <- unique(RS_label[, 1])
  location <- 0
  LQ <- rep(1/100, length(X))
  for(r in 1:length(regions)){
    i <- which(RS_label[, 1] %in% regions[r])
    Xr <- sum(X[i])
    sectors <- unique(RS_label[i, 2])
    for(s in 1:length(sectors)){
      j <- which( RS_label[i, 2] %in% sectors[s])
      Xri <- sum(X[i[j]])
      k <- which(RS_label[, 2] %in% sectors[s])
      Xni <- sum(X[k])
      top <- Xri/Xr
      bot <- Xni/Xn
      location <- location + 1
      LQ[location] <- top/bot
    }
  }
  SLQ <- matrix(LQ, ncol = length(X), nrow = length(X))
  i <- which(SLQ >= 1)
  SLQ[i] <- 1
  upd <- A * SLQ
  upd
}