PTT數位生活區 / R_Language

[問題] 一維陣列相臨值取平均的實作

看板R_Language作者 (討厭有好心推文後刪文者)時間9年前 (2015/07/24 00:28), 9年前編輯推噓7(7推 0噓 23→)
留言30則, 4人參與, 最新討論串1/2 (看更多)
[問題類型]: 程式諮詢 [軟體熟悉度]: 使用者 [問題敘述]: 想寫一個 function 自動求出相鄰元素之平均（或其它統計量） 我想做出一個 funciton，可以做相鄰值的統計量或套用特定的 function。 目前想到的參數有 1. x: 來源 numeric vector 2. n: 取幾個相鄰元素 3. FUN: 想套用的統計量或 function 我的第一個困難是，在頭和在尾的元素在取相鄰元素會有例外。 例如 1:3 的第一個元素是 1，但它沒有上一個元素， 所以就只能往之後的元素納入。 如果是用 for loop，裡頭做例外處理， 我還辦得到，但不知道有沒有更好的寫法。 我的第二個困難是，我想寫成類似 R 中 *apply 系列的 FUN 的風格， 但我不甚了解怎麼撰寫這類風格的 function。 我猜是建出一個 list 再用 lapply 來延伸，不知道好不好？ 舉一個例好了 x <- 1:5 newFun(x, n, FUN) <- function{...} newFun(x, 0, sum) #回傳 1, 2, 3, 4, 5 newFun(x, 1, sum) #回傳 3, 6, 9, 12, 9 # =1+2 =1+2+3 =2+3+4 =3+4+5 =4+5 newFun(x, 2, sum) #回傳 6, 10, 15, 14, 12 newFun(x, 0, function(a){a+1}) #回傳 2,3,4,5,6 如果有什麼想法，歡迎請提供線索給我即可，不用全刻出來沒關係。 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 122.117.37.172 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1437668910.A.2A0.html ※ 編輯: andrew43 (122.117.37.172), 07/24/2015 00:29:42 ※ 編輯: andrew43 (122.117.37.172), 07/24/2015 00:30:59 ※ 編輯: andrew43 (122.117.37.172), 07/24/2015 00:31:54
→
07/24 01:50, , 1F
我給一個想法 明天再寫 頭尾先獨立， 所以剩下只是
07/24 01:50, 1F
→
07/24 01:50, , 2F
組合問題，像是五個元素，n=1的話，中間三個考慮的
07/24 01:50, 2F
→
07/24 01:50, , 3F
是(111), (0), (0)這三個的組合，最後頭尾只要用hea
07/24 01:50, 3F
→
07/24 01:50, , 4F
d(, n)還有tail(, n)做新向量插入就好
07/24 01:50, 4F
→
07/24 01:51, , 5F
n=2的話 就要頭尾兩個獨立，然後算中間那個
07/24 01:51, 5F
→
07/24 01:52, , 6F
所以要做組合的元素個數是length(v)-n*2
07/24 01:52, 6F
→
07/24 01:52, , 7F
頭尾個數都是n
07/24 01:52, 7F
→
07/24 01:53, , 8F
頭尾都要做到2*n-1個
07/24 01:53, 8F
→
07/24 01:54, , 9F
更正頭尾是從 n到2*n-1
07/24 01:54, 9F
→
07/24 01:56, , 10F
你幾乎把過程都寫出來了，謝謝。
07/24 01:56, 10F
推
07/24 01:59, , 11F
我會考慮location matrix(logical) 列數行數都是輸
07/24 01:59, 11F
推
07/24 01:59, , 12F
入向量的長度，預設是FALSE，用迴圈做head, tail做n
07/24 01:59, 12F
推
07/24 01:59, , 13F
個從n到n-1，中間用combinations把，只是要注意comb
07/24 01:59, 13F
推
07/24 01:59, , 14F
inations的順序，要運算的值改為TRUE，再用apply做
07/24 01:59, 14F
推
07/24 01:59, , 15F
就好
07/24 01:59, 15F
推
07/24 02:00, , 16F
更正應該是從n到2*n-1
07/24 02:00, 16F
→
07/24 02:12, , 17F
c兄，多謝。延著你的提示，我突然想到另x<-c(NA,x,NA)
07/24 02:12, 17F
→
07/24 02:12, , 18F
之類的辦法就可以解決不少問題。
07/24 02:12, 18F
→
07/24 02:13, , 19F
不過這還是有很多限制，例如FUN本身要能處理NA。
07/24 02:13, 19F
→
07/24 02:14, , 20F
但對於FUN=sum|mean之類的，這樣就已經足夠了。
07/24 02:14, 20F
→
07/24 02:14, , 21F
例如FUN=sum的情況，頭和尾直接填0就可，連NA處理都免了
07/24 02:14, 21F
→
07/24 02:18, , 22F
這是一個方法只是其他計算可能沒用
07/24 02:18, 22F
→
07/24 02:23, , 23F
是的。要一般化的話不能這麼做。
07/24 02:23, 23F
推
07/24 08:20, , 24F
第一個想到的是moving average, 例如hydroTSM::ma
07/24 08:20, 24F
→
07/24 08:22, , 25F
不過若還要換成別的Function. 應該改成去尾補NA即可.
07/24 08:22, 25F
→
07/24 08:23, , 26F
原來已經討論過了. 0跟NA不同. 真的建議不要補0
07/24 08:23, 26F
→
07/24 19:16, , 27F
celestialgod提到，package zoo裡有 roll*() 一系列的
07/24 19:16, 27F
→
07/24 19:16, , 28F
function，完全是我需要的沒錯。參考看看。
07/24 19:16, 28F
→
07/25 00:54, , 29F
有考慮使用filter這個內建函數嗎?
07/25 00:54, 29F
→
07/25 12:14, , 30F
w兄，我沒搞懂你的意思。
07/25 12:14, 30F
‣ 更多 andrew43 的文章
文章代碼(AID): #1LiHOkAW (R_Language)
R_Language 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 andrew43 的文章
文章代碼(AID): #1LiHOkAW (R_Language)