【问题标题】:Determine birthdate from age in r根据 r 中的年龄确定出生日期
【发布时间】:2013-08-27 19:04:31
【问题描述】:

为模拟样本分配出生日期;以下工作,但忽略闰年。

想知道是否有更精确(和优雅)的 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"

【问题讨论】:

  • 用365.25会更准确
  • 看看这些链接是否有帮助,您可能需要找到一种方法来“逆向工程”它们的工作。 link1link2
  • 网络上充斥着那些没有正确计算日期的人的恐怖故事(更不用说一种叫做“Y2K”的东西了)。对于您的示例案例,只需将 Sys.time 转换为年-月-日-小时-分钟-秒格式,然后将您的随机“年”转换为相同的格式,您就可以使用内置的时间函数来获取生日。程序员的规则编号[某事]:永远不要重新发明以前肯定由其他人完成的功能。
  • DWin 的解决方案对于我的目的来说无疑是快速且“足够好”的。回答 Fernando 的问题的目的是创建一个模拟记录集,这些记录集与实际记录集相似,但绝不是直接派生的,也不是实际记录掩码的结果。卡尔的建议让我找到了解决方案。
  • 它必须是base R吗?使用 lubridate 软件包可以解决问题

标签: r date-arithmetic


【解决方案1】:

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)

感谢 Carl 和其他人引导我完成上述工作。欢迎使用其他更好的方法。

【讨论】:

    猜你喜欢
    • 2013-10-31
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2014-04-04
    • 2013-08-16
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多