我正在进行互相关分析,并且试图通过数据帧中的成对列来循环进行分析。我已经创建了一个列表来填充结果输出,但是不能正确填充。打印功能可以工作。
我已经进行了一些研究,但仍未找到解决该问题的方法。请帮忙。
library(tseries)
blahthresholdtest <- matrix(rnorm(100, 50, 1), nrow = 20, ncol = 4)
thislist <- list()
test <- for(i in seq(1, ncol(blahthresholdtest), 2)){
for(k in seq(1, ncol(blahthresholdtest), 1)) {
vector <- numeric(ncol(blahthresholdtest))
ccftime <- ccf(blahthresholdtest[, i], blahthresholdtest[, i+1],
type="correlation", na.action=na.omit, plot=FALSE)
crosscorr <- cbind(ccftime$acf, ccftime$lag)
crosscorr <- as.data.frame(crosscorr)
colnames(crosscorr) <- c("CCF", "lag")
vector[k+1] <- with(crosscorr, lag[CCF == max(CCF)])
thislist[[k]] <- vector
}
print(vector[k])
print(vector[k + 1])
}
thislist
do.call(rbind, thislist)