2013-08-05 39 views
2

我需要一些幫助來理解和解決與as.POSIXct相關的一些奇怪行爲。POSIXct的奇怪轉換問題

我正在從三個不同的米以三分鐘的時間間隔收集三年的光數據。爲了將此數據與其他數據集匹配,我試圖使用動物園以5分鐘的間隔創建插值光測量。

當我將日期和時間轉換爲as.POSIXct時會出現奇怪的問題。該轉換對大部分數據集都能正常工作,但始終將時間設置爲三個特定時間段的一小時:03/13/11 2:00:00 - 2:45:00,03/11/12 2:00 :00 - 2:45:00和03/10/13 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 

回答

2

這是由夏令時造成的。您需要設置的參數tz=as.POSIXct調用匹配的數據被記錄在時區例如,有沒有問題,如果你設置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"