R，在 data.frame 中找到连续时间步长的平均长度

Question

我有以下 data.frame，其中 time 列按升序排序:

colA=c(1,2,5,6,7,10,13,16,19,20,25,40,43,44,50,51,52,53,68,69,77,79,81,82)
colB=rnorm(24)

df=data.frame(time=colA, x=colB)

如何计算并取 time 列中观察到的连续时间步长的平均值？

具体来说，我需要根据连续观察对 time 列中的行进行分组，例如1,2 和 5,6,7 和 19,20 和 43,44 等等...然后取每组长度的平均值。

Answer 1

您可以像这样对连续观察的集群进行分组:

df$group <- c(0, cumsum(diff(df$time) != 1)) + 1

给出:

df
#>    time          x group
#> 1     1  0.7443742     1
#> 2     2  0.1289818     1
#> 3     5  1.4882743     2
#> 4     6 -0.6626820     2
#> 5     7 -1.1606550     2
#> 6    10  0.3587742     3
#> 7    13 -0.1948464     4
#> 8    16 -0.2952820     5
#> 9    19  0.4966404     6
#> 10   20  0.4849128     6
#> 11   25  0.0187845     7
#> 12   40  0.6347746     8
#> 13   43  0.7544441     9
#> 14   44  0.8335890     9
#> 15   50  0.9657613    10
#> 16   51  1.2938800    10
#> 17   52 -0.1365510    10
#> 18   53 -0.4401387    10
#> 19   68 -1.2272839    11
#> 20   69 -0.2376531    11
#> 21   77 -0.9268582    12
#> 22   79  0.4112354    13
#> 23   81 -0.1988646    14
#> 24   82 -0.5574496    14

您可以通过以下方式获得这些组的长度:

rle(df$group)$lengths
#> [1] 2 3 1 1 1 2 1 1 2 4 2 1 1 2

连续组的平均长度为:

mean(rle(df$group)$lengths)
#> [1] 1.714286

和每组内x的平均值使用

tapply(df$x, df$group, mean)
#>          1          2          3          4          5          6          7 
#>  0.4366780 -0.1116876  0.3587742 -0.1948464 -0.2952820  0.4907766  0.0187845 
#>          8          9         10         11         12         13         14 
#>  0.6347746  0.7940166  0.4207379 -0.7324685 -0.9268582  0.4112354 -0.3781571