lubridate包

R语言的基础包中提供了两种类型的时间数据，一类是Date日期数据，它不包括时间和时区信息，另一类是POSIXct/POSIXlt类型数据，其中包括了日期、时间和时区信息。一般来讲，R语言中建立时序数据是通过字符型转化而来，但由于时序数据形式多样，而且R中存贮格式也是五花八门，例如Date/ts/xts/zoo/tis/fts等等。用户很容易被一系列的数据格式所迷惑，所以时序数据的转化和操作并不是非常方便。所幸的是，我们有了lubridate包。lubridate包主要有两类函数，一类是处理时点数据（time instants），另一类是处理时段数据（time spans）

时点类函数，它包括了解析、抽取、修改。

library(lubridate)

##
## Attaching package: 'lubridate'

## The following object is masked from 'package:base':
##
## date

# 从字符型数据解析时间，会自动识别各种分隔符
x <- ymd('2010-04-08')
# 观察x日期是一年中的第几天
yday(x)

## [1] 98

# 修改x日期中的月份为5月
month(x) <- 5

时段类函数，它可以处理三类对象，分别是： interval：最简单的时段对象，它由两个时点数据构成。 duration：去除了时间两端的信息，纯粹以秒为单位计算时段的长度，不考虑闰年和闰秒，它同时也兼容基本包中的difftime类型对象。 period：以较长的时钟周期来计算时段长度，它考虑了闰年和闰秒，适用于长期的时间计算。以2012年为例，duration计算的一年是标准不变的365天，而period计算的一年就会变成366天。 有了时点和时段数据，就可以进行各种计算了。

# 从两个时点生成一个interval时段数据
y <- interval(x,now())
# 从interval格式转为duration格式
z <- as.duration(as.integer(as.difftime(y,units = "secs")))
# 时点+时段生成一个新的时点
now() + as.duration(z)

## [1] "2024-07-16 11:58:30 CST"

# 10天后的时间数据
now() + ddays(10)

## [1] "2017-06-21 21:59:15 CST"

# 下面的两条语句很容易看出duration和period的区别，dyears(1)表示duration对象的一年，它永远是365天。而years(1)表示period对象的一年，它识别出2012是闰年，它有366天，所以得到正确的时点。
ymd('20120101') + dyears(1)

## [1] "2012-12-31"

ymd('20120101') + years(1)

## [1] "2013-01-01"

为了处理时区信息，lubridate包提供了三个函数： tz：提取时间数据的时区 with_tz：将时间数据转换为另一个时区的同一时间 force_tz：将时间数据的时区强制转换为另一个时区

# 输入欧洲杯决赛在乌克兰的开场时间，再转为北京时间
eurotime <- ymd_hms('2012-07-01 21:45:00',tz='EET')
with_tz(eurotime,tzone='asia/shanghai')

## [1] "2012-07-01 18:45:00 asia"

最后来玩玩股票指数作为结束吧，在金融市场中谣传着一种日历效应。即在一周的第一天或者一年的第一个月份，股票会出现不错的上涨。让我们用热图来观察一下。首先是获取上证指数数据，然后根据不同的月份和星期数，将收益率汇集到不同的组中。将该组收益率的中位数映射到图形颜色上去。可以从下图看到，似乎并不存在明显周一效应或是一月效应。

library(quantmod)

## Loading required package: xts

## Loading required package: zoo

##
## Attaching package: 'zoo'

## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric

## Loading required package: TTR

## Version 0.4-0 included new data defaults. See ?getSymbols.

library(ggplot2)
library(lubridate)
# 读取上证指数历史数据
getSymbols('^SSEC',src='yahoo',from = '1997-01-01')

## 'getSymbols' currently uses auto.assign=TRUE by default, but will
## use auto.assign=FALSE in 0.5-0. You will still be able to use
## 'loadSymbols' to automatically load data. getOption("getSymbols.env")
## and getOption("getSymbols.auto.assign") will still be checked for
## alternate defaults.
##
## This message is shown once per session and may be disabled by setting
## options("getSymbols.warning4.0"=FALSE). See ?getSymbols for details.

##
## WARNING: There have been significant changes to Yahoo Finance data.
## Please see the Warning section of '?getSymbols.yahoo' for details.
##
## This message is shown once per session and may be disabled by setting
## options("getSymbols.yahoo.warning"=FALSE).

## Warning: ^SSEC contains missing values. Some functions will not work if
## objects contain missing values in the middle of the series. Consider using
## na.omit(), na.approx(), na.fill(), etc to remove or replace them.

## [1] "SSEC"

time <- ymd(as.character(index(SSEC)))
open <- as.numeric(Op(SSEC))
high <- as.numeric(Hi(SSEC))
low <- as.numeric(Lo(SSEC))
close <- as.numeric(Cl(SSEC))
volume <- as.numeric(Vo(SSEC))
# 根据收盘和开盘计算当日收益率
profit <- (close-open)/open
# 提取时间数据中的周数和月份
wday <- wday(time)-1
mday <-month(time)
data <- data.frame(time,wday,mday,profit)
p <- ggplot(data,aes(factor(mday),factor(wday),z=profit))
# 收益率的热图，图中颜色越浅，表示汇集到这个组中的收益率中位数越高
p +stat_summary_2d(fun=function(x) median(x))+ labs(x='月份',y='星期')

## Warning: Removed 172 rows containing non-finite values (stat_summary2d).