заполнение фрейма данных в список с помощью цикла for - PullRequest
Новая
       38

заполнение фрейма данных в список с помощью цикла for

0 голосов
/ 18 октября 2018

Я выполняю кросс-корреляционный анализ и пытаюсь провести цикл анализа по парам столбцов в кадре данных.Я создал список для заполнения полученного результата, но он не заполняется правильно.Хотя функция печати работает.

Я провел некоторые исследования, но до сих пор не нашел решения этой проблемы.Пожалуйста помоги.

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)
...