martakarass / arctools Goto Github PK

devtools::install_github("martakarass/arctools")
#install.packages("arctools")
library(arctools)
library(lubridate)

# generate a time vector 6am to 6pm next day
start.time <- ymd_hms("2024-02-07 06:00:00 UTC")
acc_ts <- start.time + minutes(0:1439)
class(acc_ts)
# Result is "POSIXct" "POSIXt" 

acc <- rnorm(length(acc_ts))

ts_mtm  <- midnight_to_midnight(acc = acc_ts, acc_ts = acc_ts)
# Warning message:
# 1 failed to parse.

# check that the length is the desired output
length(acc_ts) #1440
length(ts_mtm) #4320 
# in addition to the parsing error, the mtm function adds 1 extra day (1440 minutes) 
# at the end

# the error appears to happen when midnight_to_midnight calls lubridate::ymd_hms

acc_ts2 <- lubridate::ymd_hms(acc_ts)

# print original time stamps around the problematic one
acc_ts[(which(is.na(acc_ts2))-2):(which(is.na(acc_ts2))+2)]

################################################################################
# a "hack" solution: adding 1 second to date times
acc_ts <- start.time + minutes(0:1439) + seconds(1)
ts_mtm  <- midnight_to_midnight(acc = acc_ts, acc_ts = acc_ts)
# no warning

# check that the length is the desired output
length(acc_ts) #1440
length(ts_mtm) #2880 

ymd_hms("2024-02-08 00:00:00 UTC")
# "2024-02-08 UTC", works

ymd_hms(ymd_hms("2024-02-08 00:00:00 UTC"))
# fails

ymd_hms("2024-02-08 00:01:00 UTC")
# "2024-02-08 00:01:00 UTC", works

ymd_hms(ymd_hms("2024-02-08 00:01:00 UTC"))
# "2024-02-08 00:01:00 UTC", works