中的R

Question

滑動窗口我有以下所示的數據幀DF，有兩列A和B：

A     B     
1     0    
3     0    
4     0     
2     1      
6     0      
4     1      
7     1     
8     1      
1     0 

的滑動窗口方法如下所示進行。對於列3，在滑動窗口大小爲3的情況下，使用以下公式計算平均值B：B（寬度= 3，寬度= 1）。每個窗口的平均值顯示在左側。現在

A:         1    3    4    2    6    4    7    8    1             
    B:        0    0    0    1    0    1    1    1    0         
       [0    0    0]                                             0 
        [0    0    1]             0.33 
          [0    1    0]                        0.33 
           [1    0    1]                         0.66 
             [0    1    1]             0.66 
              [1    1    1]           1 
               [1    1    0]           0.66 
output:       0   0.33 0.33 0.66  0.66 1  1 1 0.66 

，對於每一行/列座標，含有該座標被認爲是和應當保留的最高平均值其給出如圖列「輸出」的結果的所有窗口。

我需要獲得如上所示的輸出。輸出應該是這樣的：

A     B     Output 
1     0      0 
3     0      0.33 
4     0      0.33 
2     1      0.66 
6     0      0.66 
4     1      1 
7     1      1 
8     1      1 
1     0     0.66 

R的任何幫助？

Answer 1

試試這個：

# form input data 
library(zoo) 
B <- c(0, 0, 0, 1, 0, 1, 1, 1, 0) 

# calculate 
k <- 3 
rollapply(B, 2*k-1, function(x) max(rollmean(x, k)), partial = TRUE) 

最後一行返回：

[1] 0.0000000 0.3333333 0.3333333 0.6666667 0.6666667 1.0000000 1.0000000 
[8] 1.0000000 0.6666667 

如果有你可能想試試這個NA值：

k <- 3 
B <- c(1, 0, 1, 0, NA, 1) 
rollapply(B, 2*k-1, function(x) max(rollapply(x, k, mean, na.rm = TRUE)), partial = TRUE) 

其中最後行給出了這個：

[1] 0.6666667 0.6666667 0.6666667 0.5000000 0.5000000 0.5000000 

擴大它這些構成爲：

c(mean(B[1:3], na.rm = TRUE), ## 
max(mean(B[1:3], na.rm = TRUE), mean(B[2:4], na.rm = TRUE)), ## 
max(mean(B[1:3], na.rm = TRUE), mean(B[2:4], na.rm = TRUE), mean(B[3:5], na.rm = TRUE)), 
max(mean(B[2:4], na.rm = TRUE), mean(B[3:5], na.rm = TRUE), mean(B[4:6], na.rm = TRUE)), 
max(mean(B[3:5], na.rm = TRUE), mean(B[4:6], na.rm = TRUE)), ## 
mean(B[4:6], na.rm = TRUE)) ## 

如果你不想在每端（標有上述##）下降partial = TRUE的k-1組件。

Answer 2

將R庫TTR具有用於計算平均值在滑動窗口

SMA =簡單移動平均

data$sma <- SMA(data$B, 3) 

更多的文檔許多功能是這裏http://cran.r-project.org/web/packages/TTR/TTR.pdf