我有一个像下面这样包含两个循环的代码。该代码读取每月流数据并使其成为多重复制。循环太慢了。我想知道是否有其他方法可以让它更快?
library(xlsx)
library(data.table)
a <- read.xlsx("streamflow.xlsx",sheetName = "Sheet1", header = TRUE)
b=matrix(nrow=129792,ncol=17)
b= data.frame(b)
i=0
for (j in -11:1236)
{
for (k in 1:104)
{
i=i+1
j=j+12
j[j > 1248] <-j-1248
b[i,] <-a[j,]
}
}
由于
答案 0 :(得分:0)
我相信这是你的双for循环到矢量化代码的正确翻译。它应该大大提高速度。此外,无需将b
声明为矩阵并将其转换为data.frame,可以从a
获取值。
j_iter <- -11:1236
k_iter <- 1:104
k <- seq(12, length(k_iter) * 12, 12)
k <- rep(k, times=length(j_iter))
j <- rep(j_iter, each=length(k_iter))
j <- j + k
j[j > 1248] <- j[j > 1248] - 1248
b <- a[j,]