add load_tables_vectors function

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: decompr
 Type: Package
 Title: Export Decomposition (Wang-Wei-Zhu and source)
-Version: 1.2.1
+Version: 1.3
 Date: 2014-09-01
 Authors@R: c( person("Bastiaan", "Quast", email = "[email protected] ", role =
     c("aut", "cre") ), person("Fei", "Wang", role = "aut"), person("Victor",

NAMESPACE

@@ -3,4 +3,5 @@
 export(decomp)
 export(kung_fu)
 export(load_tables)
+export(load_tables_vectors)
 export(wwz)

NEWS

@@ -1,3 +1,9 @@
+decompr 1.3
+=======================
+
+* add load_tables_vectors to input in simple form
+
+
 decompr 1.2.1
 =======================
 

R/decompr.R

@@ -4,7 +4,7 @@
 #' @author
 #' Bastiaan Quast \email{bquast@@gmail.com}
 #' Fei Wang
-#' Victor Kümmritz
+#' Victor Kummritz
 #' @references {Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods", WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
 #' 
 #' {Wang, Zhi, Shang-Jin Wei, and Kunfu Zhu. Quantifying international production sharing at the bilateral and sector levels. No. w19677. National Bureau of Economic Research, 2013.}
@@ -22,4 +22,25 @@ NULL
 #' @description the intermediate demand data
 #' @references {Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods",
 #' WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
+NULL
+#' @name region_names
+#' @docType data
+#' @title World Input-Output Database 2011
+#' @description the names of the regions data
+#' @references {Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods",
+#' WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
+NULL
+#' @name industry_names
+#' @docType data
+#' @title World Input-Output Database 2011
+#' @description the names of the industries data
+#' @references {Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods",
+#' WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
+NULL
+#' @name output
+#' @docType data
+#' @title World Input-Output Database 2011
+#' @description final output
+#' @references {Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods",
+#' WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
 NULL

R/load_tables_vectors.R

@@ -0,0 +1,215 @@
+#' Load the Input-Output and Final demand tables
+#' 
+#' This function loads the demand tables
+#' and defines all variables for the decomposition
+#' 
+#' @param x the intermediate demand table, it has dimensions GN x GN (G = no. of country, N = no. of industries),
+#'  excluding the first row and the first column which contains the country names,
+#'  and the second row and second column which contain the industry names for each country.
+#'  In addition, an extra row at the end should contain final demand.
+#' @param y the final demand table it has dimensions GN x MN,
+#'  excluding the first row and the first column which contains the country names,
+#'  the second column which contains the industry names for each country,
+#'  and second row which contains the five decomposed final demands (M).
+#' @param k is a vector of country of region names
+#' @param i is a vector of sector or industry names
+#' @param o is a vecotr of final outputs
+#' @return a decompr class object
+#' @author Bastiaan Quast
+#' @details Adapted from code by Fei Wang.
+#' @export
+#' @examples
+#' # load World Input-Output Database for 2011
+#' data(wiod2)
+#' 
+#' # create intermediate object (class decompr)
+#' decompr_object <- load_tables_vectors(intermediate_demand, final_demand, region_names, industry_names, output)
+#' str(decompr_object)
+
+
+load_tables_vectors <- function(x, y, k, i, o) {
+
+  # Part 1: getting the rownames etc.
+  # GN     <- length(x) - 2
+  # regnam <- unique(x[3:(GN+2),1])
+  # G      <- length(regnam)
+  # N      <- GN / G
+  # secnam <- unique( x[3:(N+2),2] )
+  # rownam <- paste( x[3:(GN+2),1], ".", x[3:(GN+2),2], sep="" )
+
+
+  G      <- length(k)
+  N      <- length(i)
+  GN     <- G * N
+  rownam <- as.vector(t(outer(k, i, paste, sep="."))) 
+
+  # making regions' names
+  z <- rownam
+  dim( z ) <- c( N, G )
+
+  # making the big rownames: bigrownam
+  z1 <-  t( array( rownam,dim=c( GN,G ) ) )
+  tot <- rep( "sub", times=GN )
+  z01 <- rbind( z1,tot )
+  dim(z01 ) <- c( (G+1)*GN,1 )
+
+  z2 <-  c( k, "TOTAL" )
+  z02 <- rep( z2,times=GN )
+
+  bigrownam <- paste( z01, z02, sep="." )
+
+  # define dimensions
+  Ad   <- array( 0,dim=c( GN,GN ) )
+  Am   <- array( 0,dim=c( GN,GN ) )
+  Bd   <- array( 0,dim=c( GN,GN ) )
+  Bm   <- array( 0,dim=c( GN,GN ) )
+  Y    <- array( 0,dim=c( GN,length( k ) ) )
+  Yd   <- array( 0,dim=c( GN,length( k ) ) )
+  Ym   <- array( 0,dim=c( GN,length( k ) ) )
+  ESR  <- array( 0,dim=c( GN,length( k ) ) )
+  Eint <- array( 0,dim=c( GN,length( k ) ) )
+  Efd  <- array( 0,dim=c( GN,length( k ) ) )
+
+  X <- o
+
+  #### this might not be the best way to construct V
+  V <- o - colSums( x )
+
+  A <- t(t(x)/o)
+  A[ is.na( A ) ] <- 0
+  A[ A==Inf ] <- 0
+  II <- diag(GN)
+  B <- solve( II-A )
+  Bm <- B
+  Am <- A
+
+  for (j in 1:length(k) )  {
+    m=1+(j-1)*N
+    n=N+(j-1)*N
+
+    Ad[m:n,m:n]  <- A[m:n,m:n]
+    Bd[m:n,m:n]  <- B[m:n,m:n]
+    Bm[m:n,m:n]  <- 0
+    Am[m:n,m:n]  <- 0
+  }
+
+  L <- solve( II-Ad )
+  Vc <- V/o
+  Vc[ is.na( Vc ) ] <- 0
+  Vc[ Vc==Inf ] <- 0
+
+  Vhat <- diag(GN)
+  diag(Vhat) <- Vc
+
+
+  # contruct final demand components
+  fdc <- dim(y)[2] / G
+
+  # Part 2: computing final demand: Y
+  for ( j in 1:length(k) ){
+    m = 1 + (j-1) * fdc
+    n = fdc + (j-1) * fdc
+    Y[ ,j ] <- rowSums( y[,m:n ] )
+  }
+  Ym <- Y
+
+
+  # Part 3: computing export: E, Esr
+  E <- cbind( x, y )
+  rm( x, y )
+  gc()
+
+  for (j in 1:length(k) )  {
+    m=1+(j-1)*N
+    n=N+(j-1)*N
+    E[m:n,m:n]  <- 0   # intermediate demand for domestic goods
+  }
+
+  for (j in 1:length(k))  {
+    m=1+(j-1)*N
+    n=N+(j-1)*N
+    s=GN+1+(j-1)*fdc
+    r=GN+fdc+(j-1)*fdc
+    E[m:n,s:r]  <- 0  # final demand for domestic goods
+    Yd[ m:n,j ] <- Y[m:n,j]
+    Ym[ m:n,j ] <- 0
+  }
+
+  z <- E
+  E <- rowSums( E )
+  E <- as.matrix( E )
+
+  for (j in 1:length(k))  {
+    m = 1 + (j-1) * N
+    n = N + (j-1) * N
+    s = GN + 1 + (j-1) * fdc
+    r = GN + fdc + (j-1) * fdc
+    ESR[ ,j ] <- rowSums( z[,m:n] ) + rowSums( z[ ,s:r ] )
+    Eint[ ,j ] <- rowSums( z[,m:n] )
+    Efd[ ,j ] <- rowSums( z[ ,s:r ] )
+  }
+
+  Exp <- diag(GN)
+  diag(Exp) <- rowSums(ESR)
+
+
+  # Part 4: naming the rows and columns in variables
+  colnames(A)     <-  rownam
+  rownames( A )   <- rownam
+  dimnames(B)     <- dimnames( A )
+  dimnames(Bm)    <- dimnames( A )
+  dimnames(Bd)    <- dimnames( A )
+  dimnames(Ad)    <- dimnames( A )
+  dimnames(Am)    <- dimnames( A )
+  dimnames(L)     <- dimnames( A )
+  names(Vc)       <- rownam
+  names(o)        <- rownam
+  colnames(Y)     <- k
+  rownames( Y )   <- rownam
+  dimnames(Ym)    <- dimnames( Y )
+  names(E)        <- rownam
+  colnames(ESR)   <- k
+  rownames( ESR ) <- rownam
+  dimnames(Eint)  <- dimnames( ESR )
+  dimnames(Efd)   <- dimnames( ESR )
+
+
+  # Part 5: creating decompr object
+  out <- list( Exp = Exp,
+               Vhat = Vhat,
+               A = A,
+               B = B,
+               Ad = Ad,
+               Am = Am,
+               Bd = Bd,
+               Bm = Bm,
+               L = L,
+               Vc = Vc,
+               X = o,
+               Y = Y,
+               Yd = Yd,
+               Ym = Ym,
+               E = E,
+               ESR = ESR,
+               Eint = Eint,
+               Efd = Efd,
+               G = G,
+               N = N,
+               GN = GN,
+               bigrownam = bigrownam,
+               k <- k,
+               rownam = rownam,
+               tot = tot,
+               z = z,
+               z01 = z01,
+               z02 = z02,
+               z1 = z1,
+               z2 = z2
+               )
+
+  class(out) <- 'decompr'
+
+  # Part 6: returning object
+  return(out)
+
+}

R/wwz.R

@@ -283,7 +283,7 @@ wwz <- function( x ) {
   rm( EEr, z, VrBrs )
   gc()
 
-  dimnames( ALL )  <-  list( x$rownam, x$regnam, decomp19)
+  dimnames( ALL )  <-  list( x$rownam, x$k, decomp19)
 
 
 

decompr.Rproj

@@ -13,6 +13,5 @@ RnwWeave: knitr
 LaTeX: pdfLaTeX
 
 BuildType: Package
-PackageUseDevtools: Yes
 PackageInstallArgs: --no-multiarch --with-keep.source
 PackageRoxygenize: rd,namespace

man/decompr.Rd

@@ -10,7 +10,7 @@ Export Decomposition using the Wang-Wei-Zhu and source decompositions algorithms
 \author{
 Bastiaan Quast \email{bquast@gmail.com}
 Fei Wang
-Victor Kümmritz
+Victor Kummritz
 }
 \references{
 {Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods", WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }

man/industry_names.Rd

@@ -0,0 +1,13 @@
+% Generated by roxygen2 (4.0.2): do not edit by hand
+\docType{data}
+\name{industry_names}
+\alias{industry_names}
+\title{World Input-Output Database 2011}
+\description{
+the names of the industries data
+}
+\references{
+{Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods",
+WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
+}
+

man/load_tables_vectors.Rd

@@ -0,0 +1,46 @@
+% Generated by roxygen2 (4.0.2): do not edit by hand
+\name{load_tables_vectors}
+\alias{load_tables_vectors}
+\title{Load the Input-Output and Final demand tables}
+\usage{
+load_tables_vectors(x, y, k, i, o)
+}
+\arguments{
+\item{x}{the intermediate demand table, it has dimensions GN x GN (G = no. of country, N = no. of industries),
+excluding the first row and the first column which contains the country names,
+and the second row and second column which contain the industry names for each country.
+In addition, an extra row at the end should contain final demand.}
+
+\item{y}{the final demand table it has dimensions GN x MN,
+excluding the first row and the first column which contains the country names,
+the second column which contains the industry names for each country,
+and second row which contains the five decomposed final demands (M).}
+
+\item{k}{is a vector of country of region names}
+
+\item{i}{is a vector of sector or industry names}
+
+\item{o}{is a vecotr of final outputs}
+}
+\value{
+a decompr class object
+}
+\description{
+This function loads the demand tables
+and defines all variables for the decomposition
+}
+\details{
+Adapted from code by Fei Wang.
+}
+\examples{
+# load World Input-Output Database for 2011
+data(wiod2)
+
+# create intermediate object (class decompr)
+decompr_object <- load_tables_vectors(intermediate_demand, final_demand, region_names, industry_names, output)
+str(decompr_object)
+}
+\author{
+Bastiaan Quast
+}
+

man/output.Rd

@@ -0,0 +1,13 @@
+% Generated by roxygen2 (4.0.2): do not edit by hand
+\docType{data}
+\name{output}
+\alias{output}
+\title{World Input-Output Database 2011}
+\description{
+final output
+}
+\references{
+{Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods",
+WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
+}
+

man/region_names.Rd

@@ -0,0 +1,13 @@
+% Generated by roxygen2 (4.0.2): do not edit by hand
+\docType{data}
+\name{region_names}
+\alias{region_names}
+\title{World Input-Output Database 2011}
+\description{
+the names of the regions data
+}
+\references{
+{Timmer, Marcel P. (ed) (2012), "The World Input-Output Database (WIOD): Contents Sources and Methods",
+WIOD Working Paper Number 10, downloadable at http://www.wiod.org/publications/papers/wiod10.pdf }
+}
+