# extRemes - Analýza extrémních hodnot
# metoda blokovych maxim

# nastaveni praconvniho adresare
getwd()
# setwd("C:/data/R_data/MFGcviceni")
# setwd("\\\\ha-ntc.ics.muni.cz/profiles/XXX/Desktop/cviceni")
getwd()

# install.packages("extRemes")
library(extRemes)

# rocni hodnoty maximalnich dennich uhrnu srazek, Brno, Turany
pX <- read.table('turany_p_max_annual_1958_2018.csv', header = T, dec = '.', sep = ';')
names(pX)
plot(pX$yeart, pX$maxpt, type="l", col = "darkblue", lwd = 1.5, cex.lab = 1.25,
xlab = "rok", ylab = "Suma srazek [mm]", main="Maximalni denni sumy srazek")

# ?fevd()
model1 <- fevd(pX$maxpt)
plot(model1)

ci(model1)
ci(model1, type = "parameter")
return.level(model1)
return.level(model1, do.ci = TRUE)
ci(model1, return.period = c(100, 150, 200))