当前位置:天才代写 > tutorial > R语言教程 > R语言教程之时间序列问题汇总实现代码教学

R语言教程之时间序列问题汇总实现代码教学

2018-05-24 08:00 星期四 所属: R语言教程 浏览:707

把表明时间序列的字符串转成时刻类型数据,在R言语里面有两个根本的函数:as.POSIXlt() 和 as.POSIXct()。两者都是S3泛型函数,依据参数的数据类型挑选不同的变换办法,除字符串外还能够变换数字、因子等数据类型,适应性很强:

as.POSIXlt("2010/01/01")
## [1] "2010-01-01 CST"
as.POSIXlt("2010/01/01 10:30:30")
## [1] "2010-01-01 10:30:30 CST"
as.POSIXlt("10:30:30 2010/01/01", format = "%H:%M:%S %Y/%m/%d")
## [1] "2010-01-01 10:30:30 CST"
tt <- Sys.time() as.POSIXlt(3600 * 24 * 1001, origin = tt)
## [1] "2019-07-20 07:15:10 CST"
as.POSIXlt(3600 * 24 * 1001, origin = tt, tz = "UTC")
## [1] "2019-07-19 23:15:10 UTC"
as.POSIXct("2010/01/01")
## [1] "2010-01-01 CST"
as.POSIXct("2010/01/01 10:30:30")
## [1] "2010-01-01 10:30:30 CST"
as.POSIXct("10:30:30 2010/01/01", format = "%H:%M:%S %Y/%m/%d")
## [1] "2010-01-01 10:30:30 CST"
as.POSIXct(3600 * 24 * 1001, origin = tt)
## [1] "2019-07-20 07:15:10 CST"
as.POSIXct(3600 * 24 * 1001, origin = tt, tz = "UTC")
## [1] "2019-07-19 23:15:10 UTC"

非标准次序日期时刻字符串需求运用format参数。

那么,什么是POSIXct和POSIXlt?它们间有什么不同?咱们需求了解它们的界说:POSIXct用某个时刻到UNIX元年(1970-01-01 00:00:00)所经过的秒数来记载时刻,即经过计数的方法表明时刻(count time);而POSIXlt用列表方法表明时刻(list time),时刻的每一部分都是列表的一个元素:

(ct <- as.POSIXct("1980-09-30 10:11:12"))
## [1] "1980-09-30 10:11:12 CST"
(lt <- as.POSIXlt("1980-09-30 10:11:12"))
## [1] "1980-09-30 10:11:12 CST"
unlist(ct)
## [1] "1980-09-30 10:11:12 CST"
unlist(lt)
##    sec    min   hour   mday    mon   year   wday   yday  isdst   zone 
##   "12"   "11"   "10"   "30"    "8"   "80"    "2"  "273"    "0"  "CST" 
## gmtoff 
##     NA

lt时刻的年份是从1900年开端计数的。在没有任何辅助工具函数的条件下,从一个lt时刻中很容易取得年月日时分秒和星期等数据,而ct时刻不行:

c(1900 + lt$year, lt$mon, lt$mday, lt$hour, lt$min, lt$sec)
## [1] 1980    8   30   10   11   12

但另一方面,ct时刻的存储只需求用一个双精度数,lt则需求一个长列表,相比之下当然是ct节约资源。

尽管终究可能表明相同的时刻,但ct时刻和lt时刻类型数据不是随意能够兼并的:

c(ct, lt)
##                                                 sec 
## "1980-09-30 10:11:12 CST" "1970-01-01 08:00:12 CST" 
##                       min                      hour 
## "1970-01-01 08:00:11 CST" "1970-01-01 08:00:10 CST" 
##                      mday                       mon 
## "1970-01-01 08:00:30 CST" "1970-01-01 08:00:08 CST" 
##                      year                      wday 
## "1970-01-01 08:01:20 CST" "1970-01-01 08:00:02 CST" 
##                      yday                     isdst 
## "1970-01-01 08:04:33 CST" "1970-01-01 08:00:00 CST" 
##                      zone                    gmtoff 
##                        NA                        NA
c(lt, ct)
## [1] "1980-09-30 10:11:12 CST" "1980-09-30 10:11:12 CST"
c(NULL, ct)
## [1] 339127872
c(NULL, lt)
## $sec
## [1] 12
## 
## $min
## [1] 11
## 
## $hour
## [1] 10
## 
## $mday
## [1] 30
## 
## $mon
## [1] 8
## 
## $year
## [1] 80
## 
## $wday
## [1] 2
## 
## $yday
## [1] 273
## 
## $isdst
## [1] 0
## 
## $zone
## [1] "CST"
## 
## $gmtoff
## [1] NA
c(ct, NULL)
## [1] "1980-09-30 10:11:12 CST"
c(lt, NULL)
## Error in as.POSIXct.default(X[[i]], ...): do not know how to convert 'X[[i]]' to class "POSIXct"


2 时刻的算术运算

ct和lt能够混合运算,但只能求差值,得到的是天数。时刻和数字的加减运算是把数字当成秒对待,在次序上也有考究,详细看例子:

tt1 <- as.POSIXct("2010-11-12 10:20:30") tt2 <- as.POSIXlt("2010-11-21 10:30:40") tt1 - tt2
## Time difference of -9.00706 days
tt1 + tt2
## Error in `+.POSIXt`(tt1, tt2): binary '+' is not defined for "POSIXt" objects
tt1 - 10
## [1] "2010-11-12 10:20:20 CST"
tt1 + 10
## [1] "2010-11-12 10:20:40 CST"
tt1 - 10
## [1] "2010-11-12 10:20:20 CST"
10 + tt1
## [1] "2010-11-12 10:20:40 CST"
10 - tt1
## Error in `-.POSIXt`(10, tt1): can only subtract from "POSIXt" objects

4 字符串变换为时刻:

除了as.POSIXlt() 和 as.POSIXct()函数外,还有其他一些函数。如果载入了R言语的其他包,会有更多用法。

4.1 as.Date

用法很简单:

as.Date("2010-10-30")
## [1] "2010-10-30"
as.Date("2010/10/30")
## [1] "2010-10-30"
as.Date("20101030")
## Error in charToDate(x): character string is not in a standard unambiguous format
as.Date("20101030", format = "%Y%m%d")
## [1] "2010-10-30"

format参数用于识别输入字符串的摆放次序。

4.2 strptime与strftime

这两者称号和参数都很像,但这儿面有乌龙。先看参数:

args(strptime)
## function (x, format, tz = "") 
## NULL
args(strftime)
## function (x, format = "", tz = "", usetz = FALSE, ...) 
## NULL
strptime("2010-10-10")
## Error in strptime("2010-10-10"): 短少参数"format",也没有缺省值
strptime("2010-10-10", format = "%Y-%m-%d")
## [1] "2010-10-10 CST"
strftime("2010-10-10")
## [1] "2010-10-10"
strftime("2010-10-10", usetz = TRUE)
## [1] "2010-10-10 CST"
  • format参数在strptime中为必选,在strftime中为可选

  • strptime强制包括时区,而strftime默许不设置时区。如果strftime设置usetz=TRUE,输出成果就和strptime相同

但是,它们返回值的数据类型是彻底不相同的:

str(strptime("2010-10-10", format = "%Y-%m-%d"))
##  POSIXlt[1:1], format: "2010-10-10"
str(strftime("2010-10-10", usetz = TRUE))
##  chr "2010-10-10 CST"

strptime得到的是时刻类型数据,strftime得到字符串。检查它们的函数源代码会了解得更清楚:

strptime
## function (x, format, tz = "") 
## {
##     y <- .Internal(strptime(as.character(x), format, tz))
##     names(y$year) <- names(x)
##     y
## }
####
strftime
## function (x, format = "", tz = "", usetz = FALSE, ...) 
## format(as.POSIXlt(x, tz = tz), format = format, usetz = usetz, 
##     ...)
####

strptime函数运用的是内部C言语函数变换字符串为时刻;而strftime本质上是字符串格式化函数,先用as.POSIXlt将字符串转为时刻,然后用format函数输出字符串。

再说format参数吧,strftime尽管终究运用了strptime,但format参数并没有传递下去,而是作为format函数的参数,用于格式化输出字符串:

strptime("30-10-2010", format = "%d-%m-%Y")
## [1] "2010-10-30 CST"
strftime("30-10-2010", format = "%d-%m-%Y", tz = "CST", usetz = TRUE)
## [1] "20-10-30 CST"
strptime("2010-10-10", format = "%A")
## [1] NA
strftime("2010-10-10", format = "%A")
## [1] "周日"

strptime得到的是lt时刻。再检查as.POSIXlt对字符串处理的源代码会发现,它其实就是对strptime函数的封装。

<div id="outline-container-sec-5" class="outline-2" font-size:14px;background-color:#ffffff;"="" style="word-wrap: break-word; font-family: "sans serif", tahoma, verdana, helvetica; font-size: 12px; white-space: normal; padding: 0px; margin: 0px; color: rgb(51, 51, 51);">

5 格式化日期/时刻的字符串

经过格式化时刻,能够得到表明时刻的不同部分,也能够把时刻依照需求的方法展现出来。在处理日期时刻类型数据上,as.character、strftime和format的用法和成果是相同的:

as.character(Sys.Date(), format = "%A")
## [1] "星期六"
strftime(Sys.Date(), format = "%A")
## [1] "星期六"
format(Sys.Date(), format = "%A")
## [1] "星期六"
format(Sys.Date(), format = "%Y年%m月%d日 %A")
## [1] "2016年10月22日 星期六"

strftime灵敏一些,能直接处理代表日期时刻的字符串:

strftime("2010-10-10", format = "%A")
## [1] "周日"
as.character(as.POSIXlt("2010-10-10"), format = "%A")
## [1] "周日"
format(as.POSIXlt("2010-10-10"), format = "%A")
## [1] "周日"

不同的格式化字符代码不同的含义,详细可用“?format.POSIXct”检查:

strftime("2010-10-30 01:02:03", format = "%Y")
## [1] "2010"
strftime("2010-10-30 01:02:03", format = "%m")
## [1] "10"
strftime("2010-10-30 01:02:03", format = "%d")
## [1] "30"
strftime("2010-10-30 01:02:03", format = "%H")
## [1] "01"
strftime("2010-10-30 01:02:03", format = "%M")
## [1] "02"
strftime("2010-10-30 01:02:03", format = "%S")
## [1] "03"
strftime("2010-10-30 01:02:03", format = "%A")
## [1] "星期六"
strftime("2010-10-30 01:02:03", format = "%B")
## [1] "十月"
strftime("2010-10-30 01:02:03", format = "%c")
## [1] "2010年10月30日 星期六 01时02分03秒"

字符串的输出与locale有关,可经过LC_TIME设置以取得想要的格式。

xx <- Sys.getlocale("LC_TIME") Sys.setlocale("LC_TIME", "C")
## [1] "C"
strftime("2010-10-30 01:02:03", format = "%A")
## [1] "Saturday"
strftime("2010-10-30 01:02:03", format = "%B")
## [1] "October"
strftime("2010-10-30 01:02:03", format = "%c")
## [1] "Sat Oct 30 01:02:03 2010"
Sys.setlocale("LC_TIME", xx)
## [1] "zh_CN.utf8"

 

    关键字:

天才代写-代写联系方式