brry / extremestat

Fit, plot and compare several (extreme value) distributions by means of a plot with return periods on a linear scale

intro

Fit, plot and compare several (extreme value) distributions. Can also compute (truncated) distribution quantile estimates and draw a plot with return periods on a linear scale.

See the Vignette for an introduction to the package.

installation

Install / load the package and browse through the examples:

install.packages("extremeStat")
library(extremeStat)
vignette(extremeStat)

# update to the current development version, incl. vignette:
remotes::install_github("brry/extremeStat", build_vignettes=TRUE)

trouble

If normal installation doesn't work, click on Code - Download ZIP (topright), unzip the file to some place, then

setwd("that/path")
dd <- dir("extremeStat-master/R", full=T)
dummy <- sapply(dd, source)

This creates all R functions as objects in your globalenv workspace (and overwrites existing objects of the same name!).

extremestat's Issues

is there a way to make sure the fitted curve go through the max observation data?

Hi,would you please help me with the problem that when fitting GEV distribution to some kind of block maxima rainfall data?The question is that" how to make sure the fitted curve go through the max observation rainfalldata using extremeStat?",many thanks!

log scale on y axis in plotLextreme

I plot the return periods graphs with 17 distribution. The only issue is, i want log scale on y axis. As in code only "log=TRUR" which shows log scale on x axis. But i want log scale on y axis also.
As the peak values of different station are too vary.

NAs still result in a selected distribution

Hello, this package is very useful. I am using the distLfit function on a gridded dataset. The gridded dataset include grid cells with missing/NA values. When I loop this function through each grid cell, however, the cells with NAs still include a selected distribution (and it's usually the first one listed alphabetically (e.g., c("gev", "glo", "nor"). so when plotting the results, grid cells that should have a missing value, it now says 'gev'. why is this?

Statistical Distributions names

Dear Sir,
I am using your extremeStat package for the prediction of different Rivers using the available 17 different statistical distribution, which is mentioned in this package (extremeStat). it is a kind request that please provide the full names of each distribution, as the distributions are in the short form (i.e. wak, gev, nor, lap, etc). I googled most of them for the full name but some were found for more than one distribution.

Best Regards
Muhammad Irfan

bootstrapping sampling technique

Calculating confidence intervalls of return periods using extremeStat::distLexBoot I stumbled across the bootstrapping sampling.
I assumed bootstrapping samples are taken with replacement, eg. (sample(x, size, replace = TRUE)).
Currently sampling is done without replacement, is this on intend?