# Select the appropiate columns from the example dataset
# Dates_BT and rename column names to make the file compatible
# with the function
library(magrittr)
library(dplyr)
library(lubridate)
Bloom <- Dates_BT %>%
select(Year, sbloom) %>%
rename(Fday=sbloom) %>%
filter(Year==2003)
# Obtain estimated hourly RH from the example dataset Tudela_DW
Weather <- Tudela_DW %>%
filter (Tudela_DW$Year==2003)
RH_h <- hourly_RH(Weather, 42.13132)
# Estimate the number of russet-inducing days for a RH>55\%
# between 30 to 34 days after full bloom for each season
Russet_Risk <-russet(RH_h,Bloom,55,30,34)
Run the code above in your browser using DataLab