R/key.sector.diff.R

In ignaciomsarmiento/ioanalysis: Input-Output Analysis

key.sector.diff<-function(mip, X.j, X.i, key=TRUE, cutoff=1, write.xlsx=TRUE, name="Key_sector.xlsx"){
  if(dim(mip)[2]!=dim(X.j)[1])stop("check dimensions")
  if(dim(mip)[1]!=dim(X.i)[2])stop("check dimensions")
  nosect<-length(X.j)
  for (i in 1:nosect) if (X.j[i]==0) X.j[i]=0.1
  aij=mip/X.j
  bij_leon<-data.frame(solve(diag(nosect)-aij))
  V<-sum(bij_leon)/nosect
  b.j<-rowSums(bij_leon)
  for (i in 1:nosect) if (b.j[i]==1)b.j[i]=0
  blink<-b.j/V
  aij_gosh=mip/X.i
  inv<-solve(diag(nosect)-aij_gosh)
  bij_gosh<-data.frame(t(inv))
  b.i<-rowSums(bij_gosh)
  for (i in 1:nosect) if (b.i[i]==1)b.i[i]=0
  flink<-b.i/V
  sector.numb<-seq(1:nosect)
  b<-data.frame(sector.numb,blink,flink)
  rownames(b)<-names(mip)
  colnames(b)<- c("Sector Number", "Backward Linkage", "Forward Linkage")
  if(write.xlsx==TRUE) write.xlsx2(b, file=name, sheetName="Sheet1", showNA=TRUE)
  if(key==TRUE){
    b$key<-ifelse(b$Forward>=cutoff&b$Backward>=cutoff,"Key",ifelse(b$Forward>=cutoff&b$Backward<cutoff,"II",ifelse(b$Forward<cutoff&b$Backward<cutoff,"III","IV")))
    b<-b[order(-b$Forward,-b$Backward),] 
  }
  b
}