3

为模拟样本分配出生日期；以下作品，但忽略闰年。

想知道是否有更精确（和优雅）的 r 方法？

# Simulate 10 persons with age evenly distributed 0 to 21
age <- runif(10, 0, 21)
# calc age in seconds
agesecs <- age*365*24*60*60
# subtract from right now to establish 'birthdate'
bday <- as.Date(Sys.time() - agesecs)
bday
[1] "2008-03-28" "1998-06-12" "2010-05-02" "2007-01-11" "2007-06-07"
[6] "1999-05-22" "2004-01-29" "2013-03-29" "1998-06-01" "2006-10-14"
4

1 回答 1

0

在difftime object日期和时间算术方面做得很好。问题是它承认的单位不包括年份。因此，我使用了周，而不是随机数年，一个单位 difftime 接受。

    # Simulate 10 persons with age evenly distributed 0 to 21 (use weeks)
    agewks <- runif(10, 0, (21*52))
    # convert to difftime object
    agedt <- as.difftime(agewks, units="weeks")
    # above could be combined into single step
    # agewks <- as.difftime(runif(10, 0, (21*52)), units='weeks')
    # subtract from right now to establish a 'birthdate' for our simulated persons
    bday <- as.POSIXct(Sys.time() - agedt)
    bday
    [1] "1997-05-26 13:23:07 EDT" "2003-02-24 13:07:48 EST"
    [3] "2006-12-20 12:38:04 EST" "2002-01-02 15:17:14 EST"
    [5] "1993-10-07 15:49:19 EDT" "2001-05-04 04:05:29 EDT"
    [7] "2003-09-28 09:35:30 EDT" "1996-05-17 20:58:15 EDT"
    [9] "2008-08-09 14:17:24 EDT" "2011-05-09 23:26:04 EDT"
    # to create a date object use
    bday <- as.Date(Sys.time() - agedt)

感谢卡尔和其他人引导我完成上述工作。欢迎替代和更好的方法。

于 2013-08-29T00:55:52.887 回答