从分组数据中选择两个随机且连续的行

在下面的数据中（包括dput），我对三个人（IndIDII）进行了重复观察（纬度和经度）。请注意，每个人都有不同数量的位置，并且它们是按以下顺序排列的：IndYear.

  IndIDII      IndYear  WintLat  WintLong
1 BHS_265 BHS_265-2015 47.61025 -112.7210
2 BHS_265 BHS_265-2016 47.59884 -112.7089
3 BHS_770 BHS_770-2016 42.97379 -109.0400
4 BHS_770 BHS_770-2017 42.97129 -109.0367
5 BHS_770 BHS_770-2018 42.97244 -109.0509
6 BHS_377 BHS_377-2015 43.34744 -109.4821
7 BHS_377 BHS_377-2016 43.35559 -109.4445
8 BHS_377 BHS_377-2017 43.35195 -109.4566
9 BHS_377 BHS_377-2018 43.34765 -109.4892

我想filter并制作一个新的df每个都有两个连续的行IndIDII。在我的较大数据集中，所有个体至少有 2 个观察值（即行），每个个体有 2 到 4 个观察值。显然，对于只有 2 行的个人，代码将返回仅有的 2 行可用。随着更多数据，第 1 行和第 2 行，or2和3，or3 和 4 将被随机选择。行的顺序并不重要，只要它们是连续的（即可以返回 3 和 4or4 和 3)。

一如既往，非常感谢！

Dat <- structure(list(IndIDII = c("BHS_265", "BHS_265", "BHS_770", "BHS_770", 
"BHS_770", "BHS_377", "BHS_377", "BHS_377", "BHS_377"), IndYear = c("BHS_265-2015", 
"BHS_265-2016", "BHS_770-2016", "BHS_770-2017", "BHS_770-2018", 
"BHS_377-2015", "BHS_377-2016", "BHS_377-2017", "BHS_377-2018"
), WintLat = c(47.6102519805014, 47.5988417247191, 42.9737859090909, 
42.9712914772727, 42.9724390816327, 43.3474354347826, 43.3555934579439, 
43.3519543396226, 43.3476466990291), WintLong = c(-112.720994832869, 
-112.708887595506, -109.039964727273, -109.036693522727, -109.050923061224, 
-109.482114456522, -109.444522149533, -109.45659254717, -109.489241553398
)), class = "data.frame", row.names = c(NA, -9L))

您可以使用ave。在每个组中，创建一个行索引 (i <- seq_along(x)）。要获取要保留的行的第一个索引，请从除最后一个行索引之外的所有行中采样一行（sample(head(i, -1), 1)。还包括下一行（+ 0:1）。检查采样行中有哪些行索引（i %in% ...）。将结果强制返回逻辑子集数据。

Dat[as.logical(ave(Dat$IndIDII, Dat$IndIDII, FUN = function(x){
  i <- seq_along(x)
  i %in% (sample(head(i, -1), 1) + 0:1)
})), ]

#   IndIDII      IndYear  WintLat  WintLong
# 1 BHS_265 BHS_265-2015 47.61025 -112.7210
# 2 BHS_265 BHS_265-2016 47.59884 -112.7089
# 4 BHS_770 BHS_770-2017 42.97129 -109.0367
# 5 BHS_770 BHS_770-2018 42.97244 -109.0509
# 7 BHS_377 BHS_377-2016 43.35559 -109.4445
# 8 BHS_377 BHS_377-2017 43.35195 -109.4566

以类似的方式，但更简洁，data.table及其内置行索引 (.I）和每组的行数（.N)

library(data.table)
setDT(Dat)
Dat[Dat[ , (sample(.I[-.N], 1)) + 0:1, by = IndIDII]$V1]