R——投入产出表，直接消耗系数矩阵，完全消耗系数矩阵

library(aTSA)
Q=read.csv("D:Q2017.csv")
X=read.csv("D:X2017.csv")
#直接消耗系数A
Q1<-as.matrix(Q)
X_1<-diag(1/X$total)
A<-Q1%*%X_1

#完全消耗系数B
I<-diag(length(X$total))
B<-solve(I-A)-I
write.csv(B,file="D:B2017.csv")

数据原始

 

Q2017中间投入表

 

X2017总投入