輥的意思是基於R中

Question

我想使用在實施例R.數據移動平均函數（rollmean）條件：

x<-seq(1,48,by=1) 
y<-c(rep(11,12), rep(12,12)) 
z<-data.frame(x,y) 

欲計算與3天滯後移動z[,"x"]的平均用於每個z[,"y"] （即11 & 12）並將它們填充到（新）第三列，其中每個z[,"y"]的前兩行由NAs指定。

Answer 1

使用data.table包使用zoo例如：

library(zoo) 
library(data.table) 
DT <- data.table(z) 
DT[, roll_x := rollmeanr(x, 3, fill = NA), y] 

當然，如果有一些缺失值最好是使用rollapplyr：

DT[, roll_x := rollapplyr(x, 3, mean, fill = NA), y] 

那將是更快的另一個包caTools：

library(caTools) 
DT[, roll_x := runmean(x, 3, align = 'right', endrule = 'NA'), by = y] 

在這個數據幀爲單位工作過：

transform(z, roll_x = ave(x, y, FUN = function(x) rollmeanr(x, 3, fill = NA))) 

Answer 2

您可以使用基本功能filter。

z[ , moving_average := filter(x, rep(1/3, 3), sides = 1), by = "y"] 

注意的3個結果在 NA秒的滯後。