Demonstrate the goodness of fit of the selected non-extreme marginal distribution

Plots demonstrating the goodness of fit of a selected (not truncated) non-extreme marginal distribution to a dataset.

Usage

Diag_Non_Con_Sel(
  Data,
  Omit = NA,
  x_lab = "Data",
  y_lim_min = 0,
  y_lim_max = 1,
  Selected
)

Arguments

Data

Numeric vector containing realizations of the variable of interest.

Omit

Character vector specifying any distributions that are not to be tested. Default "NA", all distributions are fit.

x_lab

Numeric vector of length one specifying. Label on the x-axis of histogram and cumulative distribution plot.

y_lim_min

Numeric vector of length one specifying the lower y-axis limit of the histogram.

y_lim_max

Numeric vector of length one specifying the upper y-axis limit of the histogram.

Selected

Character vector of length one specifying the chosen distribution. Options are: Gaussian "Gaus", Gumbel "Gum", Laplace "Lapl", logistic "Logis", and reverse Gumbel "RGum".

Value

Panel consisting of three plots. Upper plot: Plots depicting the AIC of the two fitted distributions. Middle plot: Probability Density Functions (PDFs) of the selected distributions superimposed on a histogram of the data. Lower plot: Cumulative distribution function (CDFs) of the selected distribution overlaid on a plot of the empirical CDF.

See also

Diag_Non_Con

Examples

S20.Rainfall<-Con_Sampling_2D(Data_Detrend=S20.Detrend.df[,-c(1,4)],
                              Data_Declust=S20.Detrend.Declustered.df[,-c(1,4)],
                              Con_Variable="Rainfall",Thres=0.97)
Diag_Non_Con(Data=S20.Rainfall$Data$OsWL,x_lab="O-sWL (ft NGVD 29)",
             y_lim_min=0,y_lim_max=1.5)
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced

#> $AIC
#>   Distribution      AIC
#> 1         Gaus 201.1158
#> 2          Gum 396.8524
#> 3         Lapl 209.5831
#> 4        Logis 188.5013
#> 5         RGum 167.1707
#> 
#> $Best_fit
#> [1] "RGum"
#> 
Diag_Non_Con_Sel(Data=S20.Rainfall$Data$OsWL,x_lab="O-sWL (ft NGVD 29)",
                 y_lim_min=0,y_lim_max=1.5,Selected="Logis")
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced