标题: 为R代码提速的2条小技巧 [打印本页] 作者: 晚风 时间: 2011-7-27 20:14 标题: 为R代码提速的2条小技巧 一、如果需要保存数据到向量或矩阵，首选方法是：新建一个向量或矩阵，然后为新建向量或矩阵赋值。备选方法是：通过连接当前向量或矩阵，以生成新的向量或矩阵。在数据量较少的情况下，上述两种方法没有明显区别，在遇到大量数据时，前者可以提高60倍效率。
小试验：
n <- 10^5
a <- NULL
#连接
system.time(for(i in 1:n){a <- c(a,1)})
#新建并赋值
b <- rep(0,n)
system.time(for(i in 1:n){b[i] <- 1})
前者耗时31.78秒，后者耗时0.44秒。
二、能避免使用圆括号的，尽量不要使用。
小试验：
f <- function(n){ for(i in 1:n){ x <- 1/(1+x) }}
g <- function(n){ for(i in 1:n){ x <- (1/(1+x)) }}
h <- function(n){ for(i in 1:n){ x <- (1+x)^(-1) }}
j <- function(n){ for(i in 1:n){ x <- {1/{1+x}} }}
k <- function(n){ for(i in 1:n){ x <- 1/{1+x} }}
x <- 1
system.time(f(10^6))
system.time(g(10^6))
system.time(h(10^6))
system.time(j(10^6))
system.time(k(10^6))
分别耗时：1.36、1.51、1.90、1.44、1.25。