library(data.table) library(ggplot2) library(knmipluim)
We want to compute the annual cycle of return levels of maximum (and minimum) daily temperature, using the following approach.
windowSize <- 20L
r windowSize
days prior and after.For the maxima at De Bilt this looks like:
maxDat <- txDat[name == "tx260_1981_2010"] #maxDat <- txDat[name == "tx260_WH_2050"] #maxDat[, name := NULL] maxRetLevel <- ReturnLevelClimatology(maxDat, "tx", "max", windowSize = windowSize, kLoc = 15, kScale = 8, kShape = 4)
PlotEstimates(maxRetLevel)
PlotReturnLevels(maxRetLevel)
And for the minima it looks like:
minDat <- tnDat[name == "tn260_1981_2010"] #minDat <- tnDat[name == "tn260_WH_2050"] #minDat[, name := NULL] minRetLevel <- ReturnLevelClimatology(minDat, "tn", "min", windowSize = windowSize, kLoc = 15, kScale = 8, kShape = 4)
PlotEstimates(minRetLevel)
PlotReturnLevels(minRetLevel)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.