xiAnova {flood}

Homogeneity test for the shape

Description

A test for assumption H_0: "shape parameter is equal for all d GEV margins" with test statistic based on (trimmed) L-moments.

Usage

xiAnova(x, leftrim = 0, rightrim = 1)

Arguments

x

matrix of observations (rows: observations, d columns: stations).

leftrim

integer indicating lower trimming parameter (\ge 0).

rightrim

integer indicating upper trimming parameter (\ge 0).

Value

p-value of the test.

Examples

library("evd")
# sample observations of 100 years at 5 stations:
set.seed(1053)
x19 <- matrix(rgev(400, 2, 1, 0.1), ncol=4) # 4 stations with the same shape
x10 <- rgev(100, 2, 1, 0.4) # one station with a different shape
x <- cbind(x19, x10)
xiAnova(x=x, leftrim=0, rightrim=1)