inst/doc/Ungauged.R

In martindurocher/floodStat: FloodNet Regional Frequency Analysis

## ------------------------------------------------------------------------
library(floodnetRfa)
data(canFlood)
attach(canFlood)

xdf <- data.frame(q100   = 0,
              area   = scale(log(BASIN_AREA)),
              wb     = scale(log(WB_AREA_PC)),
              stream = scale(log(STREAM_DEN)),
              map    = scale(log(PPTAVG_BAS)),
              lon    = ELON,
              lat    = ELAT)

## ------------------------------------------------------------------------
library(lmomco)
lmm <- split(as.matrix(canFlood[,12:14]), 1:nrow(canFlood))
lmm <- lapply(lmm, vec2lmom, lscale = FALSE)
para <- lapply(lmm, lmom2par, 'gev')
xdf$q100  <- log(sapply(para, function(z) par2qua(.99,z)))