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
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.