Calculates log return probability given a return value of interest, using model fit from `extRemes::fevd`

. Standard error is obtained via the delta method. The return probability is the probability of exceeding the return value in a single block. For non-stationary models (those that include covariates for the location, scale, and/or shape parameters, log probabilities and standard errors are returned for as many sets of covariates as provided.

1 | ```
calc_logReturnProb_fevd(fit, returnValue, covariates = NULL)
``` |

`fit` |
fitted object from extRemes |

`returnValue` |
value for which log return probability is desired |

`covariates` |
matrix of covariate values, each row a set of covariates for which the return probability is desired |

Results are calculated (and returned) on log scale as delta-method based standard errors are more accurate for the log probability. Confidence intervals on the probability scale should be calculated by calculating a confidence interval for the log probability and exponentiating the endpoints of the interval.

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.