fixed similar timezone issues that RMevents was having. Fixed in a different way - uses lubridate to find the timezone of the input data, then maintain that tz later on. This avoids R using the timezone of the computer when tz gets dropped through the loop.

Author: limnoliver

NAMESPACE

@@ -10,4 +10,5 @@ export(RMintensity)
 export(RMprep)
 export(discharge_events)
 import(dplyr)
+importFrom(lubridate,tz)
 importFrom(rlang,sym)

R/RMevents_sample.R

@@ -15,6 +15,7 @@
 #' @param bdate character column name in dfsamples for the beginning of the sampling period
 #' @param edate character column name in dfsamples for the ending of the sampling period
 #' @return list of storms and storms2
+#' @importFrom lubridate tz
 #' @export
 #' @examples
 #' RDB <- CedarRRain
@@ -48,6 +49,8 @@ RMevents_sample <- function(df,
   startRainDates <- numeric()
   endRainDates <- numeric()
   tipsbystorm <- data.frame()
+  
+  rain_timezone <- lubridate::tz(df[,time])
 
   for (i in 1:nrow(dfsamples)){
     beginRow <- max(which(df[,time]<dfsamples[i,bdate])+1)
@@ -109,8 +112,8 @@ RMevents_sample <- function(df,
       }
       tipsbystorm <- sub_tips
     } else {
-      startRainDates <- as.POSIXct(c(startRainDates,BD), origin = "1970-01-01")
-      endRainDates <- as.POSIXct(c(endRainDates,ED), origin = "1970-01-01")
+      startRainDates <- as.POSIXct(c(startRainDates,BD), origin = "1970-01-01", tz = rain_timezone)
+      endRainDates <- as.POSIXct(c(endRainDates,ED), origin = "1970-01-01", tz = rain_timezone)
       if (nrow(sub_tips) > 0) {
         event <- event + 1
       } else {