library(openxlsx) data1<-read.xlsx("http://kanggc.iptime.org/book/data/ss-fdi-e.xlsx") k2012<-data1$k2012 k2017<-data1$k2017 jj2012<-data1$jj2012 jj2017<-data1$jj2017 A<-matrix(data=NA, nrow=10, ncol=7, byrow=T) for(i in 1:10) { for(j in 1:4) { A[i,1]<-(k2017[i]-k2012[i]) A[i,2]<-(A[i,1]/k2012[i]) A[i,3]<-(jj2017[i]-jj2012[i]) A[i,4]<-(A[i,3]/jj2012[i]) } } for(i in 1:9) { for(j in 5:7) { A[i+1,5]<-(jj2012[i+1])*A[1,2] A[i+1,6]<-jj2012[i+1]*(A[i+1,2]-A[1,2]) A[i+1,7]<-jj2012[i+1]*(A[i+1,4]-A[i+1,2]) } } CS<-colSums(A[,5:7], na.rm=T) A[1,5]<-CS[1] A[1,6]<-CS[2] A[1,7]<-CS[3] AA<-round(A,digits=3) AA