我需要一些帮助来理解和解决与 as.POSIXct
相关的一些奇怪行为。
我正在处理三年来的光数据,这些数据以 15 分钟的间隔从三个不同的仪表收集。为了将此数据与另一个数据集匹配,我尝试使用 zoo 以 5 分钟为间隔创建插值光测量。
当我使用 as.POSIXct
转换日期和时间时,出现了奇怪的问题。该转换适用于大多数数据集,但始终将三个特定时间段的时间设置为后一小时:2011 年 3 月 13 日 2:00:00 - 2:45:00、2012 年 11 月 3 日 2:00 :00 - 2:45:00 和 2013 年 10 月 3 日 2:00:00 - 2:45:00。我已经在下面发布了 R 控制台的结果。谁能解释这种行为并提出解决方案?
> Ldata <- read.csv("9813306_LI_AIRTEM_HENGILL2_SUM.csv")
> head(Ldata)
X. date time Temp_C Light_lux
1 1 10/29/10 15:30:00 2.195 3444.5
2 2 10/29/10 15:45:00 1.330 3100.0
3 3 10/29/10 16:00:00 1.330 3100.0
4 4 10/29/10 16:15:00 1.221 2927.8
5 5 10/29/10 16:30:00 1.221 2152.8
6 6 10/29/10 16:45:00 1.112 1463.9
> str(Ldata)
'data.frame': 88417 obs. of 5 variables:
$ X. : int 1 2 3 4 5 6 7 8 9 10 ...
$ date : Factor w/ 923 levels "1/1/11","1/1/13",..: 105 105 105 105 105 105 105 105 105 105 ...
$ time : Factor w/ 98 levels "0:00:00","0:15:00",..: 33 34 35 36 37 38 39 40 41 42 ...
$ Temp_C : num 2.19 1.33 1.33 1.22 1.22 ...
$ Light_lux: num 3444 3100 3100 2928 2153 ...
这部分数据给我带来了麻烦
> Ldata[c(12908:12911, 44194:44197, 79136:79139),]
X. date time Temp_C Light_lux
12908 12908 3/13/11 2:00:00 -5.496 0
12909 12909 3/13/11 2:15:00 -5.249 0
12910 12910 3/13/11 2:30:00 -5.995 0
12911 12911 3/13/11 2:45:00 -6.246 0
44194 44197 3/11/12 2:00:00 0.674 0
44195 44198 3/11/12 2:15:00 0.563 0
44196 44199 3/11/12 2:30:00 0.453 0
44197 44200 3/11/12 2:45:00 0.343 0
79136 79139 3/10/13 2:00:00 -2.494 0
79137 79140 3/10/13 2:15:00 -2.610 0
79138 79141 3/10/13 2:30:00 -2.377 0
79139 79142 3/10/13 2:45:00 -2.610 0
> P_datetime <- as.POSIXct(paste(Ldata$date, Ldata$time), format = "%m/%d/%y %H:%M:%S")
> Ldata_m <- cbind(Ldata, P_datetime)
请注意 P_datetime
比时间列早 1 小时。
> Ldata_m[duplicated(Ldata_m$P_datetime),]
X. date time Temp_C Light_lux P_datetime
12908 12908 3/13/11 2:00:00 -5.496 0 2011-03-13 01:00:00
12909 12909 3/13/11 2:15:00 -5.249 0 2011-03-13 01:15:00
12910 12910 3/13/11 2:30:00 -5.995 0 2011-03-13 01:30:00
12911 12911 3/13/11 2:45:00 -6.246 0 2011-03-13 01:45:00
44194 44197 3/11/12 2:00:00 0.674 0 2012-03-11 01:00:00
44195 44198 3/11/12 2:15:00 0.563 0 2012-03-11 01:15:00
44196 44199 3/11/12 2:30:00 0.453 0 2012-03-11 01:30:00
44197 44200 3/11/12 2:45:00 0.343 0 2012-03-11 01:45:00
79136 79139 3/10/13 2:00:00 -2.494 0 2013-03-10 01:00:00
79137 79140 3/10/13 2:15:00 -2.610 0 2013-03-10 01:15:00
79138 79141 3/10/13 2:30:00 -2.377 0 2013-03-10 01:30:00
79139 79142 3/10/13 2:45:00 -2.610 0 2013-03-10 01:45:00
最佳答案
这是由夏令时引起的。您需要在 as.POSIXct
调用中设置 tz=
参数,以匹配记录数据的时区。例如,如果您设置 tz="UTC"
:
> as.POSIXct(paste(Ldata$date, Ldata$time), format="%m/%d/%y %H:%M:%S", tz="UTC")
# [1] "2011-03-13 02:00:00 UTC" "2011-03-13 02:15:00 UTC"
# [3] "2011-03-13 02:30:00 UTC" "2011-03-13 02:45:00 UTC"
# [5] "2012-03-11 02:00:00 UTC" "2012-03-11 02:15:00 UTC"
# [7] "2012-03-11 02:30:00 UTC" "2012-03-11 02:45:00 UTC"
# [9] "2013-03-10 02:00:00 UTC" "2013-03-10 02:15:00 UTC"
# [11] "2013-03-10 02:30:00 UTC" "2013-03-10 02:45:00 UTC"
关于r - POSIXct 的奇怪转换问题,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/18066786/