calc_returnValueDiff_fevd {climextRemes}

Calculates return value difference for two sets of covariates and standard error of difference given return period(s) of interest

Description

Calculates difference in return values (also known as return levels) for two sets of covariates given return period(s) of interest, using model fit from extRemes::fevd. Standard error is obtained via the delta method. The return value is the value for which the expected number of blocks until an event that exceeds that value is equal to the return period. Differences and standard errors are returned for as many contrasts between covariate sets as provided.

Usage

calc_returnValueDiff_fevd(
  fit,
  returnPeriod,
  covariates1,
  covariates2,
  getSE = TRUE
)

Arguments

fit

fitted object from extRemes fevd

returnPeriod

value(s) for which return value difference is desired

covariates1

matrix of covariate values, each row a set of covariates for which the return value difference relative to the corresponding row of covariates2 is desired

covariates2

matrix of covariate values, each row a set of covariates

getSE

logical indicating whether standard error is desired, in addition to the point estimate

Details

This is designed to calculate differences in return values and associated standard errors for different covariate values based on the same model fit. It is not designed for differences based on separate model fits, although it may be possible handle this case by fit two models in a single model specification using dummy variables.