※ 引述《a78998042a (Benjimine)》之銘言：
: [問題類型]:
: 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
: [軟體熟悉度]:
: 使用者(已經有用R 做過不少作品)
: [問題敘述]:
: 我需要將程式向量化
: 使用迴圈、apply系列的函數太慢了
: 但不用apply又做不出來
: 想詢問大家的經驗
: [程式範例]:
: ###vactor的問題
: #我想將一個向量分別加上一個數字
: x = 1:3
: y = rep(0, 3)
: #我想得到類似下面結果，但不用apply, 不用迴圈的做法
: sapply(x, function(x) y + x, simplify = FALSE)
: [[1]]
: [1] 1 1 1
: [[2]]
: [1] 2 2 2
: [[3]]
: [1] 3 3 3
我猶豫了一下要怎麼回你XD
不過我覺得sapply, lapply夠快了，而且for也很好用
我寫了一個簡單的測試，我不認為這會是你程式速度的瓶頸：
PS: 相信你的y不是全部0，不然用mapply + rep最快
y <- 1:3e4
x <- 1:1e4
st <- proc.time()
z1 <- sweep(array(y, c(length(y), length(x))), 2, x, '+')
proc.time() - st
# user system elapsed
# 8.72 0.45 9.19
st <- proc.time()
z2 <- sapply(x, function(xi) y + xi)
proc.time() - st
# user system elapsed
# 3.32 0.48 3.83
all.equal(z1, z2) # TRUE
st <- proc.time()
z12 <- sweep(matrix(rep(y, length(x)), length(y)), 2, x, '+')
proc.time() - st
# user system elapsed
# 8.93 1.10 10.05
all.equal(z1, z12)
st <- proc.time()
z22 <- sapply(x, function(xi) y + xi, simplify = FALSE)
proc.time() - st
# user system elapsed
# 0.58 0.45 1.05
st <- proc.time()
z3 <- lapply(x, function(xi) y + xi)
proc.time() - st
# user system elapsed
# 0.53 0.17 0.70
all.equal(z22, z3) # TRUE
st <- proc.time()
z4 <- vector('list', length(x))
for (i in seq_along(x))
z4[[i]] <- y + x[i]
proc.time() - st
# user system elapsed
# 0.61 0.04 0.65
all.equal(z4, z3) # TRUE
: ###array的問題
: #我想將array不同維度的物件，扣除不同數字
: temp.ar = array(matrix(0, 2, 3), c(2, 3, 2))
: #我想得到類似下面的結果
: sapply(1:2, function(x)temp.ar[,,x] + x,simplify = FALSE)
: [[1]]
: [,1] [,2] [,3]
: [1,] 1 1 1
: [2,] 1 1 1
: [[2]]
: [,1] [,2] [,3]
: [1,] 2 2 2
: [2,] 2 2 2
: #另一個問題是，sapply會將物件轉成list，增加後續處理上的難度
: #例如我無法直接對list取round
: #使用sapply的simplify又會將破壞結構，之後可能需要再重新架構
: #但如果結構是array就可以直接做一些函數運算
: #不知道大家有這方面的經驗可以分享嗎
: #謝謝!
arr <- array(0, c(2, 3, 1e4))
st <- proc.time()
z <- sweep(arr, 3, 1:1e4, '+')
proc.time() - st
# user system elapsed
# 0.00 0.02 0.02