S3 `alogLik`

method to perform loglikelihood adjustment for fitted
extreme value model objects returned from the functions
`gev`

, `gpd`

and `pot`

in the evir package.
If `x`

was returned from `pot`

then the model will
need to be re-fitted using `pot_refit`

.

- x
A fitted model object with certain associated S3 methods. See

**Details**.- cluster
A vector or factor indicating from which cluster the respective log-likelihood contributions from

`loglik`

originate. The length of`cluster`

must be consistent with the`estfun`

method to be used in the estimation of the 'meat'`V`

of the sandwich estimator of the covariance matrix of the parameters to be passed to`adjust_loglik`

. In most cases,`cluster`

must have length equal to the number of observations in data. The exception is the GP (only) model (`binom = FALSE`

), where the`cluster`

may either contain a value for each observation in the raw data, or for each threshold exceedance in the data.If

`cluster`

is not supplied (is`NULL`

) then it is assumed that each observation forms its own cluster. See**Details**for further details.- use_vcov
A logical scalar. Should we use the

`vcov`

S3 method for`x`

(if this exists) to estimate the Hessian of the independence loglikelihood to be passed as the argument`H`

to`adjust_loglik`

? Otherwise,`H`

is estimated inside`adjust_loglik`

using`optimHess`

.- ...
Further arguments to be passed to the functions in the sandwich package

`meat`

(if`cluster = NULL`

), or`meatCL`

(if`cluster`

is not`NULL`

).

An object inheriting from class `"chandwich"`

. See

`class(x)`

is a vector of length 5. The first 3 components are

`c("lax", "chandwich", "evir")`

.
The remaining 2 components depend on the model that was fitted.
If `gev`

was used then these components are

`c("gev", "stat")`

.
If `gpd`

was used then these components are

`c("gpd", "stat")`

.
If `pot_refit`

was used then these components are

`c("potd", "stat")`

.

See `alogLik`

for details.

If `pot`

was used then `x`

does not contain the
raw data that `alogLik`

needs. The model will need to be
re-fitted using `pot_refit`

and the user will be prompted to
do this by an error message produced by `alogLik`

.

Chandler, R. E. and Bate, S. (2007). Inference for clustered
data using the independence loglikelihood. *Biometrika*,
**94**(1), 167-183. doi:10.1093/biomet/asm015

Suveges, M. and Davison, A. C. (2010) Model
misspecification in peaks over threshold analysis, *The Annals of
Applied Statistics*, **4**(1), 203-221.
doi:10.1214/09-AOAS292

Zeileis (2006) Object-Oriented Computation and Sandwich
Estimators. *Journal of Statistical Software*, **16**, 1-16.
doi:10.18637/jss.v016.i09

`alogLik`

: loglikelihood adjustment for model fits.

```
# We need the evir package
got_evir <- requireNamespace("evir", quietly = TRUE)
if (got_evir) {
library(evir)
# An example from the evir::gev documentation
data(bmw)
out <- gev(bmw, "month")
adj_out <- alogLik(out)
summary(adj_out)
# An example from the evir::gpd documentation
data(danish)
out <- gpd(danish, 10)
adj_out <- alogLik(out)
summary(adj_out)
# An example from the evir::pot documentation
# We use lax::pot_refit() to return the input data
out <- pot_refit(danish, 10)
adj_out <- alogLik(out)
summary(adj_out)
}
#>
#> Attaching package: 'evir'
#> The following objects are masked from 'package:eva':
#>
#> dgpd, pgev, pgpd, qgev, qgpd, rgpd
#> The following object is masked from 'package:extRemes':
#>
#> decluster
#> The following objects are masked from 'package:evd':
#>
#> dgev, dgpd, pgev, pgpd, qgev, qgpd, rgev, rgpd
#> MLE SE adj. SE
#> xi 0.4972 0.1359 0.1435
#> sigma 1.1610 0.6888 0.7086
#> mu -1.6960 3.4380 3.4270
```