按随机间隔在小标题中选择行

时间:2018-12-07 16:54:18

标签: r dplyr tidyverse lubridate

我正在尝试采用一系列日期-从第一个日期开始-通过从正态分布生成的随机数选择后续日期。目前,我有通过随机数选择行号的代码,但是每次都使用相同的数字。在此示例中,它每12天选择一行:

set.seed(123)

library(tidyverse)
library(lubridate)

start_date <- as.Date('2018-03-01')
end_date <- as.Date('2018-07-01')

seq_dates <- seq(ymd(start_date), ymd(end_date), by='1 days')

seq_dates <- seq_dates %>%
  as.tibble()
seq_dates

seq_dates %>% 
  filter(row_number() %% round(rnorm(n=1, mean=14, sd=3), 0) == 1) 

是否可以使用dplyr做到这一点,但是每次都以随机间隔从开始日期中选择一行?因此,从2018-03-01开始,下一个日期可能是12天后,14天后,19天后,依此类推?

1 个答案:

答案 0 :(得分:2)

library(dplyr)

set.seed(10) 
n <- rnorm(50, 14, 3)
rows <- cumsum(round(n, 0))
diff(rows) # random ~normal increments used when selecting your rows
#  [1] 13 10 12 15 15 10 13  9 13 17 16 13 17 16 14 11 13 17 15 12  7 12  8 10 13 12 11 14 13  8 14 17
# [33] 15 10 10 15  9 13 12 17 12 12 17 11 14 15 13 12 16

seq_dates %>% 
  slice(rows[rows <= n()])
# # A tibble: 9 x 1
#   value     
#   <date>    
# 1 2018-03-14
# 2 2018-03-27
# 3 2018-04-06
# 4 2018-04-18
# 5 2018-05-03
# 6 2018-05-18
# 7 2018-05-28
# 8 2018-06-10
# 9 2018-06-19