Как записать вывод R с неравной длиной в Excel - PullRequest
Новая
       45

Как записать вывод R с неравной длиной в Excel

1 голос
/ 27 января 2020

У меня есть этот R-код для имитации данных временного ряда, в серии есть 5 колон и 35 слагаемых 

N <- c(15L, 20L, 30L, 50L, 100L)
SD = c(1, 2, 3, 4, 5) ^ 2
theta = c(0.2, 0.4, 0.6, 0.8, 0.9, 0.95, 0.99)

res <- vector('list', length(N))
names(res) <- paste('N', N, sep = '_')

set.seed(123L)
for (i in seq_along(N)){
  res[[i]] <- vector('list', length(SD))
  names(res[[i]]) <- paste('SD', SD, sep = '_')

  ma <- matrix(NA_real_, nrow = N[i], ncol = length(theta)) 

  for (j in seq_along(SD)){
    wn <- rnorm(N[i], mean = 0, sd = SD[j])
    ma[1:2, ] <- wn[1:2]

    for (k in 3:N[i]){
      ma[k, ] <- wn[k - 1L] * theta + wn[k]
    }
    colnames(ma) <- paste('ma_theta', theta, sep = '_')
    res[[i]][[j]] <- ma
  }
}

res

** Я ОЖИДАЮ ЭТО **

v1  v2   v3  v4...v15  
1   4    3    1...1
2   6    3    2...2
3   8    3    12...5
4   4    4    8...5
5   6    4    4...5
6   8    4    0...5
7   4    5    2...5
8   6    5    1...2
9   8    5    2...2
10  11   12   13...2
11  12   13   14...4
8   6    5    1...4
9   8    5    2...4
10  11   12   13...4
11  12   13   14...4

Я хочу запишите res в любой файл Excel. Это может быть файл .csv. Я хочу 1 в одной ячейке, 2 в одной ячейке и т. Д. c, а не 12345 вместе взятых в ячейке.

После этого я записываю его в файл Excel

1 Ответ

1 голос
/ 27 января 2020

res является вложенным list из matrices, мы можем rbind его после зацикливания на list

res1 <- do.call(rbind, lapply(res, function(x) do.call(rbind, x))) 
str(res1)
# num [1:1075, 1:7] -0.56 -0.23 1.513 0.382 0.143 ...
# - attr(*, "dimnames")=List of 2
#  ..$ : NULL
#  ..$ : chr [1:7] "ma_theta_0.2" "ma_theta_0.4" "ma_theta_0.6" "ma_theta_0.8" ...

Затем write.table работает как есть один matrix

Если нам нужно cbind внутренние list элементы, так как количество строк различается, опция cbind.fill из rowr

library(rowr)
out2 <- do.call(rbind, lapply(res, function(x) 
      do.call(cbind.fill, c(as.data.frame(x), fill = NA))))
colnames(out2) <-  make.unique(rep(colnames(res[[1]][[1]]), length.out = ncol(out2)))

Обновление

Если нам нужно несколько выходных данных data.frame, cbind элементы list 

res1 <- lapply(res, function(dat) do.call(cbind,  dat))
sapply(names(res1), function(nm) write.csv(res1[[nm]], 
      file = paste0(nm, ".csv"), row.names = FALSE, quote = FALSE))
...