[問題] list 元素加總,如何不用loop做到

看板R_Language作者 (Benjimine)時間10年前 (2015/01/15 16:17), 10年前編輯推噓8(8015)
留言23則, 9人參與, 最新討論串1/2 (看更多)
[問題類型]: 程式諮詢(我想用R 做某件事情,但是我不知道要怎麼用R 寫出來) [軟體熟悉度]: 入門(寫過其他程式,只是對語法不熟悉) [問題敘述] 想使用apply類型的函數將list內的元素加總,而不使用for loop [程式範例]: mat=matrix(1:4,2,2) mat.list=list(mat,mat,mat) 想加總成單一的矩陣 結構為:matrix(c(3,6,9,12),2,2) 有人知道不用loop的情況可以怎麼做嗎? 謝謝 -- -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.116.152.120 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1421309862.A.0CF.html

01/15 16:18, , 1F
試試看 ?sapply
01/15 16:18, 1F
感覺好像有難度,lapply系列的,都是針對list的每一個物件,物件跟物件間的 元素好像沒辦法@@

01/15 17:38, , 2F
do.call
01/15 17:38, 2F
我試試!謝謝

01/15 18:14, , 3F
apply(mat. list, 1:2, sum)
01/15 18:14, 3F
這個好像無法度

01/15 18:15, , 4F
Reduce也可以,只是do.call最有效率
01/15 18:15, 4F
我試試看,謝謝

01/15 18:44, , 5F
Reduce('+' , mat.list)
01/15 18:44, 5F
感謝!!!解決了!!!! 抱歉,有人能示範在do.call要怎麼做到這件事嗎?我目前沒試出來

01/15 21:49, , 6F
do.call好像不能用`+`,還是先用Reduce吧Orz
01/15 21:49, 6F

01/15 21:54, , 7F
要指定`+`.__C__matrix,但是這個我試不出來
01/15 21:54, 7F
感謝你,以前沒用過do.call,現在看到就能開始想用在哪了~

01/15 22:07, , 8F
如果list中的matrix大小一樣 建議你這樣做
01/15 22:07, 8F

01/15 22:08, , 9F
matrix(rowSums(matrix(unlist(mat.list), nrow=4)), nro
01/15 22:08, 9F

01/15 22:08, , 10F
w=2)
01/15 22:08, 10F
謝謝!!學習了!!!

01/16 02:46, , 11F
謝謝!!!學習了
01/16 02:46, 11F
回報一下,我以Reduce的方法與m大以拆開再拼的做法比較,以長度為12582912的2by2 矩陣相加,Reduce是12.38秒,unlist是1.84秒。謝謝!!

01/16 05:48, , 12F
do.call(.Primitive("+"), mat.list)
01/16 05:48, 12F
感謝!這樣寫在length(mat.list)==2時確實可以,長度>2的還要再試試

01/16 09:26, , 13F
apply和 for的效能應該沒差多少
01/16 09:26, 13F
你好,參考一下,因為最近學校給了份模擬跌代的作業,寫出來的結果,for要約3.7天, 通篇apply2.5天,再配合sapply進行優化有人到1.5天以內(Reduec),接著可能會試M大 的方法。 ※ 編輯: a78998042a (140.116.152.120), 01/16/2015 14:26:33 ※ 編輯: a78998042a (140.116.152.120), 01/16/2015 14:28:09

01/16 15:14, , 14F
for 要用 compiler,不過我太曉得目前jit的情況
01/16 15:14, 14F

01/16 15:14, , 15F
如果要超過小時的code,還是推RCpp吧
01/16 15:14, 15F
謝謝!!又學習了,看來我對R的了解仍是一知半解阿 ※ 編輯: a78998042a (140.116.1.136), 01/16/2015 17:35:47

01/16 18:47, , 16F
推Rcpp, RcppArmadillo.還有若是你每次迴圈或是模擬
01/16 18:47, 16F

01/16 18:48, , 17F
都是獨立的,平行運算也很快.Rcpp+平行運算直接逆天!!
01/16 18:48, 17F

01/16 18:49, , 18F
簡單說,模擬K次,K次間都獨立,可以分散到N個CPU平行~~
01/16 18:49, 18F

01/17 15:31, , 19F
為什麼不用array?
01/17 15:31, 19F

01/17 15:31, , 20F
mat.array = array(mat, 1:3) # 可能是吧, 隨意寫.
01/17 15:31, 20F

01/17 15:32, , 21F
sum(mat.array[,,1:3])
01/17 15:32, 21F

01/17 15:33, , 22F
你的每個list維度都一樣, 用array比較簡單吧.
01/17 15:33, 22F

01/17 19:44, , 23F
我回錯了...apply那個只能用到3維陣列...
01/17 19:44, 23F
文章代碼(AID): #1KjtUc3F (R_Language)
文章代碼(AID): #1KjtUc3F (R_Language)