dt <- ch_rfa_extractamax(flow~date, CAN01AD002)$date
ch_rfa_seasonstat(dt)
## Illustration of the analysis of multiple sites
F0 <- function(ii) data.frame(site = ii, dt = sample(dt, replace = TRUE))
x <- lapply(1:10, F0)
x <- do.call(rbind, x)
st <- ch_rfa_seasonstat(dt ~ site, x)
ch_rfa_julianplot()
points(y ~ x, st, col = 2, pch = 16)
Run the code above in your browser using DataLab