The evreg
package provides functions that model non-stationary extreme values using univariate extreme value regression modelling, while also performing variable selection to objectively determine the best model for a given dataset. This is achieved by using generalized linear modelling for each parameter, and by fitting models with maximum likelihood estimates.
One of the main functions in the evreg package is gevreg
, which performs the model fitting for generalized extreme value regression modelling. Then, forward_gevreg
can be used to perform forward selection to objectively determine the best model. The following code fits a GEV regression model using the gevreg
function for fremantle
dataset and then uses forward_gevreg
perform forward selection.
library(evreg)
data <- fremantle[ , -which(names(fremantle) %in% c("Year"))]
fit0 <- gevreg(SeaLevel, data)
forward_gevreg(fit0)
#>
#> Call:
#> gevreg(y = SeaLevel, data = data, mu = mu, mustart = mustart,
#> sigmastart = sigmastart, xistart = xistart)
#>
#> Convergence: TRUE
#>
#> Coefficients:
#> mu: (Intercept) mu: Year01 mu: SOI
#> 1.38433 0.19449 0.05452
#> sigma: (Intercept) xi: (Intercept)
#> -2.11418 -0.14999
#>
#> Log-likelihood: 53.9 DF: 5 AIC: -97.8
To get the current released version from CRAN:
install.packages("evreg")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.