差距怎么这么大,这个图是用了楼主给的数据,用前十个行业画的。
不过这种方法,说实在话,没有什么意义,应该是研究部门为了交差的 应付之作。
data = fread("E:/data/data.csv")
dat = data[-c(1,2),2:11,with=FALSE]
setDT(dat)[, names(dat) := lapply(.SD, function(x) diff(log(as.numeric(x))))]
setDT(dat)[, names(dat) := lapply(.SD, function(x) as.numeric(x))]
cp = sapply(1:nrow(dat),function(i) {dat = dat[1:i,]
dat = cor(dat)
dat = dat[upper.tri(dat)]
return(mean(dat)/sd(dat))
})
dat1 = data.frame(x = 1:nrow(dat),y = cp)
plot(dat1,type = 'l')