# Daily data with separate min temperature
daily_data <- data.frame(
date = seq(as.Date("2000-01-01"), as.Date("2002-12-31"), by = "day"),
tmin = rnorm(1096, mean = 18, sd = 5)
)
calculate_TR20(
df = daily_data,
frequency = "daily",
time_col = "date",
tmin_col = "tmin"
)
# Hourly data (will be aggregated to daily min temperature)
hourly_data <- data.frame(
datetime = seq(
as.POSIXct("2000-01-01 00:00", tz = "UTC"),
as.POSIXct("2000-01-31 23:00", tz = "UTC"),
by = "hour"
),
temperature = rnorm(31*24, mean = 16, sd = 3)
)
calculate_TR20(
df = hourly_data,
frequency = "hourly",
time_col = "datetime",
temp_col = "temperature",
threshold = 20
)
Run the code above in your browser using DataLab