UIBE在价值链指标体系中有一个基于3x4(3国4部门)的生产分解案例,附有代码,但是第四个脚本中有一部分代码运行出现错误,而且没看明白是什么意思,具体代码如下(红色部分即为出错的代码):
##***************************************************##
## Created by: Fei Wang Date: 05/08/2016 ##
## E-Mail:wangfei64@hotmail.com ##
## School of International Trade and Economics ##
## University of International Business and Economics, China ##
#####################################
## Part 0 = GVC participation index based on forward industrial linkage======
path
pathdata
pathinc
pathout
flnm
load( flnm )
GN
N
flnm = paste0( path,"/index1_PROD/SVAbyFDbyYear.RData" )
load( flnm )
#View( SVAbyFDbyYear[ ,,1] )
GVC_Pat_f
for ( j in 1:length( year ) ){
z
GVC_Pat_f[ ,1,j ]
GVC_Pat_f[ ,2,j ]
GVC_Pat_f[ ,3,j ]
}
colnam
dimnames( GVC_Pat_f )
##== VAX_Fs, DVA_Fs,RCA_DVA_f
VAX_Fs
DVA_Fs
##== RCA_DVA_f ==============================
DVA_G
for ( i in 1:length( year ) ){
#i=1
flnm
load( flnm )
DVA_G[ ,i ]
rm( T01_DVA_FIN); gc( )
flnm
load( flnm )
DVA_G[ ,i ]
rm( T02_DVA_INT); gc( )
flnm
load( flnm )
DVA_G[ ,i ]
rm( T03_DVA_INTrex1); gc( )
flnm
load( flnm )
DVA_G[ ,i ]
rm( T04_DVA_INTrex2 );gc( )
flnm
load( flnm )
DVA_G[ ,i ]
rm( T05_DVA_INTrex3 );gc( )
flnm
load( flnm )
DVA_G[ ,i ]
rm( T06_RDV_FIN1 );gc( )
flnm
load( flnm )
DVA_G[ ,i ]
rm( T07_RDV_FIN2 );gc( )
flnm
load( flnm )
DVA_G[ ,i ]
rm( T08_RDV_INT); gc( )
print( paste( year[ i ], " is done !" ) )
}
VAEXbySVA
rm( DVA_G);gc( )
RCAVAbySVAbyStateYear
WIbyYear
XIbyYear
for ( i in 1:length(year) ){
x
dim( x )
WIbyYear[ ,i ]
for ( s in 1:G ){
x
XIbyYear[(1+(s-1)*N):(N+(s-1)*N),i ]
RCAVAbySVAbyStateYear[(1+(s-1)*N):(N+(s-1)*N),i ]
XIbyYear[(1+(s-1)*N):(N+(s-1)*N),i ]/WIbyYear[ ,i ]
}
}
dimnames( RCAVAbySVAbyStateYear )
RCA_DVA_f
require( reshape )
# == -------------==============================
require( plyr )
nrownam
nrownam2
nrownam2
nrownam2
nrownam
z
rownames( z )
SVAbyFDbyYear
z
rownames( z )
GVC_Pat_f
dim( VAX_Fs )
dimnames( VAX_Fs )
dim( DVA_Fs )
dimnames( DVA_Fs )
dim( RCA_DVA_f )
dimnames( RCA_DVA_f )
flnm
write.csv( cbind(SVAbyFDbyYear,GVC_Pat_f,
VAX_Fs,DVA_Fs,RCA_DVA_f),file=flnm )
运行时显示:
Error in readChar(con, 5L, useBytes = TRUE) : 无法打开链结
此外: Warning message:
In readChar(con, 5L, useBytes = TRUE) :
无法打开压缩文件'D:/GVCinR/inc/2000/T01_DVA_FIN.RData',可能是因为'No such file or directory'
在之前的代码中并没有看到输出T01_DVA_FIN.RData,所以这一步不知道是从哪里来的,麻烦大神指导,十分感谢