如何在R中解析年份+周数?

时间:2012-02-21 15:42:01

标签: r datetime

有没有一种方法可以将一年+周数转换为R中的日期?我尝试过以下方法:

> as.POSIXct("2008 41", format="%Y %U")
[1] "2008-02-21 EST"
> as.POSIXct("2008 42", format="%Y %U")
[1] "2008-02-21 EST"

根据?strftime

  

%Y年与世纪。请注意,虽然没有零   原始公历,ISO 8601:2004将其定义为有效   (解释为1BC):见http://en.wikipedia.org/wiki/0_(year)。注意   该标准还说明了1582年之前的日历   只应在有关各方同意的情况下使用。

     

%U使用星期日作为十进制数字(00-53)   一周的第一天(通常是第一天的第一天)   年作为第1周的第1天)。美国公约。

3 个答案:

答案 0 :(得分:22)

这有点像另一个问题you may have seen before。 :)

关键问题是:一周的数字指定哪一天?这是一周的第一天吗?最后?这是模棱两可的。我不知道第一周是一年的第一天还是一年的第七天,或者可能是一年的第一个星期日或星期一(这是一个经常的解释)。 (而且比这更糟糕:这些通常看起来是0索引,而不是1索引。)因此,需要指定一周中列举的日期。

例如,试试这个:

as.POSIXlt("2008 42 1", format = "%Y %U %u")

%u指标指定星期几。

附加说明:有关格式转换的各种选项,请参阅?strptime。重要的是要注意几周的枚举,因为这些可以在年末分开,而1天是不明确的:是根据星期日或星期一,还是从第一天开始指定年?这应该在R代码运行的不同系统上进行指定和测试。我不确定Windows和POSIX系统在其中一些转换中是否会唱出相同的曲调,因此我会再次进行测试和测试。

答案 1 :(得分:5)

POSIXlt DateTimesClasses系统中的星期几==零是星期日。不完全符合圣经,也不符合从“1”惯例开始的R索引,但这就是它的本质。第0周是一年中的第一个(部分)周。第一周(但是第零周的一天)从第一个星期日开始。 POSIXlt中的所有其他序列类型都以0为起点。看看强制POSIXlt对象的列表元素做什么很有意思。实际更改POSIXlt日期的唯一方法是更改​​$ year,$ mon或$ mday元素。其他似乎是附带现象。

  today <- as.POSIXlt(Sys.Date())
  today  # Tuesday
#[1] "2012-02-21 UTC"
     today$wday <- 0  # attempt to make it Sunday
     today
# [1] "2012-02-21 UTC"   The attempt fails
 today$mday <- 19
 today
#[1] "2012-02-19 UTC"   Success

答案 2 :(得分:1)

我自己没有提出这个问题(它来自Forester的blog post),但我认为我会将其添加到答案列表中,因为它是ISO 8601周数公约的第一个实现我在R中看过。

毫无疑问,周数是一个非常含糊不清的主题,但我更喜欢通过format(..., "%U")而不是当前实施周数的ISO标准,因为这似乎是大多数人同意的,至少在德国(日历等。)。

我已将实际功能def置于底部,以便首先关注输出。另外,我偶然发现了ISOweek包,也许值得一试。

方法比较

x.days  <- c("Mon", "Tue", "Wed", "Thu", "Fri", "Sat", "Sun")
x.names <- sapply(1:length(posix), function(x) {
    x.day <- as.POSIXlt(posix[x], tz="Europe/Berlin")$wday
    if (x.day == 0) {
        x.day <- 7
    }
    out <- x.days[x.day]
})

data.frame(
    posix, 
    name=x.names,
    week.r=weeknum, 
    week.iso=ISOweek(as.character(posix), tzone="Europe/Berlin")$weeknum
)

# Result

        posix name week.r week.iso
1  2012-01-01  Sun      1  4480458
2  2012-01-02  Mon      1        1
3  2012-01-03  Tue      1        1
4  2012-01-04  Wed      1        1
5  2012-01-05  Thu      1        1
6  2012-01-06  Fri      1        1
7  2012-01-07  Sat      1        1
8  2012-01-08  Sun      2        1
9  2012-01-09  Mon      2        2
10 2012-01-10  Tue      2        2
11 2012-01-11  Wed      2        2
12 2012-01-12  Thu      2        2
13 2012-01-13  Fri      2        2
14 2012-01-14  Sat      2        2
15 2012-01-15  Sun      3        2
16 2012-01-16  Mon      3        3
17 2012-01-17  Tue      3        3
18 2012-01-18  Wed      3        3
19 2012-01-19  Thu      3        3
20 2012-01-20  Fri      3        3
21 2012-01-21  Sat      3        3
22 2012-01-22  Sun      4        3
23 2012-01-23  Mon      4        4
24 2012-01-24  Tue      4        4
25 2012-01-25  Wed      4        4
26 2012-01-26  Thu      4        4
27 2012-01-27  Fri      4        4
28 2012-01-28  Sat      4        4
29 2012-01-29  Sun      5        4
30 2012-01-30  Mon      5        5
31 2012-01-31  Tue      5        5

功能定义

直接来自blog post,我刚刚改变了一些小事。该功能仍然是粗略的(例如,第一个日期的周数很远),但我发现这是一个不错的开始!

ISOweek <- function(
    date, 
    format="%Y-%m-%d", 
    tzone="UTC", 
    return.val="weekofyear"
){
  ##converts dates into "dayofyear" or "weekofyear", the latter providing the ISO-8601 week
  ##date should be a vector of class Date or a vector of formatted character strings
  ##format refers to the date form used if a vector of
  ##  character strings  is supplied

  ##convert date to POSIXt format 
  if(class(date)[1]%in%c("Date","character")){
    date=as.POSIXlt(date,format=format, tz=tzone)
  }

#  if(class(date)[1]!="POSIXt"){
  if (!inherits(date, "POSIXt")) {
    print("Date is of wrong format.")
    break
  }else if(class(date)[2]=="POSIXct"){
    date=as.POSIXlt(date, tz=tzone)
  }
print(date)

  if(return.val=="dayofyear"){
    ##add 1 because POSIXt is base zero
    return(date$yday+1)
  }else if(return.val=="weekofyear"){
    ##Based on the ISO8601 weekdate system,
    ## Monday is the first day of the week
    ## W01 is the week with 4 Jan in it.
    year=1900+date$year
    jan4=strptime(paste(year,1,4,sep="-"),format="%Y-%m-%d")
    wday=jan4$wday

    wday[wday==0]=7  ##convert to base 1, where Monday == 1, Sunday==7

    ##calculate the date of the first week of the year
    weekstart=jan4-(wday-1)*86400  
    weeknum=ceiling(as.numeric((difftime(date,weekstart,units="days")+0.1)/7))

    #########################################################################
    ##calculate week for days of the year occuring in the next year's week 1.
    #########################################################################
    mday=date$mday
    wday=date$wday
    wday[wday==0]=7
    year=ifelse(weeknum==53 & mday-wday>=28,year+1,year)
    weeknum=ifelse(weeknum==53 & mday-wday>=28,1,weeknum)

    ################################################################
    ##calculate week for days of the year occuring prior to week 1.
    ################################################################

    ##first calculate the numbe of weeks in the previous year
    year.shift=year-1
    jan4.shift=strptime(paste(year.shift,1,4,sep="-"),format="%Y-%m-%d")
    wday=jan4.shift$wday
    wday[wday==0]=7  ##convert to base 1, where Monday == 1, Sunday==7
    weekstart=jan4.shift-(wday-1)*86400
    weeknum.shift=ceiling(as.numeric((difftime(date,weekstart)+0.1)/7))

    ##update year and week
    year=ifelse(weeknum==0,year.shift,year)
    weeknum=ifelse(weeknum==0,weeknum.shift,weeknum)

    return(list("year"=year,"weeknum"=weeknum))
  }else{
    print("Unknown return.val")
    break
  }
}