R/krugman.spec.R

In REAT: Regional Economic Analysis Toolbox

krugman.spec <- function (e_ij, e_il) {

  industries <- nrow(as.matrix(e_ij))

  s_ij <- vector()
  s_il <- vector()
  
  e_j <- sum (e_ij)
  e_l <- sum (e_il)

  for (i in 1:industries) {
    s_ij[i] <- e_ij[i]/e_j
    s_il[i] <- e_il[i]/e_l
  }
  
  K_jl <- sum(abs(s_ij-s_il))

  return(K_jl)
  
}