作者andrew43 (讨厌有好心推文後删文者)
看板R_Language
标题[问题] 一维阵列相临值取平均的实作
时间Fri Jul 24 00:28:26 2015
[问题类型]: 程式谘询
[软体熟悉度]: 使用者
[问题叙述]: 想写一个 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://webptt.com/cn.aspx?n=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
1F：→ celestialgod: 我给一个想法 明天再写 头尾先独立， 所以剩下只是 07/24 01:50
2F：→ celestialgod: 组合问题，像是五个元素，n=1的话，中间三个考虑的 07/24 01:50
3F：→ celestialgod: 是(111), (0), (0)这三个的组合，最後头尾只要用hea 07/24 01:50
4F：→ celestialgod: d(, n)还有tail(, n)做新向量插入就好 07/24 01:50
5F：→ celestialgod: n=2的话 就要头尾两个独立，然後算中间那个 07/24 01:51
6F：→ celestialgod: 所以要做组合的元素个数是length(v)-n*2 07/24 01:52
7F：→ celestialgod: 头尾个数都是n 07/24 01:52
8F：→ celestialgod: 头尾都要做到2*n-1个 07/24 01:53
9F：→ celestialgod: 更正头尾是从 n到2*n-1 07/24 01:54
10F：→ andrew43: 你几乎把过程都写出来了，谢谢。 07/24 01:56
11F：推 celestialgod: 我会考虑location matrix(logical) 列数行数都是输 07/24 01:59
12F：推 celestialgod: 入向量的长度，预设是FALSE，用回圈做head, tail做n 07/24 01:59
13F：推 celestialgod: 个从n到n-1，中间用combinations把，只是要注意comb 07/24 01:59
14F：推 celestialgod: inations的顺序，要运算的值改为TRUE，再用apply做 07/24 01:59
15F：推 celestialgod: 就好 07/24 01:59
16F：推 celestialgod: 更正应该是从n到2*n-1 07/24 02:00
17F：→ andrew43: c兄，多谢。延着你的提示，我突然想到另x<-c(NA,x,NA) 07/24 02:12
18F：→ andrew43: 之类的办法就可以解决不少问题。 07/24 02:12
19F：→ andrew43: 不过这还是有很多限制，例如FUN本身要能处理NA。 07/24 02:13
20F：→ andrew43: 但对於FUN=sum|mean之类的，这样就已经足够了。 07/24 02:14
21F：→ andrew43: 例如FUN=sum的情况，头和尾直接填0就可，连NA处理都免了 07/24 02:14
22F：→ celestialgod: 这是一个方法只是其他计算可能没用 07/24 02:18
23F：→ andrew43: 是的。要一般化的话不能这麽做。 07/24 02:23
24F：推 Edster: 第一个想到的是moving average, 例如hydroTSM::ma 07/24 08:20
25F：→ Edster: 不过若还要换成别的Function. 应该改成去尾补NA即可. 07/24 08:22
26F：→ Edster: 原来已经讨论过了. 0跟NA不同. 真的建议不要补0 07/24 08:23
27F：→ andrew43: celestialgod提到，package zoo里有 roll*() 一系列的 07/24 19:16
28F：→ andrew43: function，完全是我需要的没错。参考看看。 07/24 19:16
29F：→ Wush978: 有考虑使用filter这个内建函数吗? 07/25 00:54
30F：→ andrew43: w兄，我没搞懂你的意思。 07/25 12:14