R/fdvac.R

In DarekScigala/wiotrs: Global input-output data analysis

#'Final demand decomposition for selected country
#'
#' @param country A number 1:41 for version A, 1:44 for version B, see WIOT$countries.
#' @param calc A boolean, for calc = TRUE result as a share of total final demand, for
#'    calc = FALSE result in mln USD.
#' @return A list contains title, country, year and the final demand decomposition.
#' @examples
#' get_io(year = 1998, version = "A")
#' construct()
#' fdvac(10)
#'
#' fdvac(22, calc = FALSE)
#'
#' @export
#'


fdvac <- function(country=10, calc=TRUE){
  # Check if WIOT list exists

  if(exists("WIOT") == FALSE){
    message("WIOT data is missing? :-( Please use get_io() function.")
    stop()
  }

  # Extracting the variables out of the WIOT list

  for(i in 1:length(WIOT)){
    ##extracting the object value
    tempobj=WIOT[[i]]
    ##creating new variable with the original name of the list item
    eval(parse(text=paste(names(WIOT)[[i]],"= tempobj")))
  }

  # Verification of function input

  if(!( country %in% c(1:length(countries)))){
    message("Wrong country ID. Please see WIOT$countries")
    stop()
  }

  # Defining additional variable

  SN <- S*N # SNxSN is the size of intercountry input-output matrix
  CN <- C*N # CNxSN is the size of the intercountry final demand matrix

  # Checking if construct() function was initiated

  if(length(FD) == 1){
    message("Please run the construct() function first.")
    stop()
  }

  FD1 <- matrix(0,  nrow = SN, ncol = CN)   #creating empty matrix SNxCN
  l1 <- (country-1)*S+1                     #counters for selected country
  l2 <- country * S
  FD1[l1:l2,] <- FD[l1:l2,]                 #selection of the final demand submatrix for specified country

   # Calculation of value added amounts of all countries in the sectoral final demand of selected country

  FD1 <- as.matrix(rowSums(FD1))
  q1 <- VA %*% LeoInv %*% FD1



  # Summation and ordering

  kol <- seq_along(q1)
  q2 <- tapply(q1,rep(kol,each =S )[kol],FUN = sum)

  if(calc == TRUE){
    q2 <- 100*round(q2 / sum(q2), 7)
  }


  rownames(q2) <- countries

  q2 <- q2[order(q2, decreasing=TRUE)]


  # Preparing description of output

  q3 <- list("Decomposition of Final Demand in respect to value added from source countries", countries[country],
             year, q2)

  names(q3) <- c("Title", "Country","Year","FDVAC")


  return(q3)

}