copula (version 0.999-19)

evTestK: Bivariate Test of Extreme-Value Dependence Based on Kendall's Distribution

Description

Test of extreme-value dependence based on the bivariate probability integral transformation. The test statistic is defined in Ben Ghorbal, G. Ne<U+0161>lehov<U+00E1>, and Genest (2009).

Usage

evTestK(x, method = c("fsample","asymptotic","jackknife"), ties = NA, N = 100)

Arguments

x

a data matrix.

method

specifies the variance estimation method; can be either "fsample" (finite-sample, the default), "asymptotic" or "jackknife".

ties

logical; if TRUE, the original test is adapted to take the presence of ties in the coordinate samples of x into account; the default value of NA indicates that the presence/absence of ties will be checked for automatically.

N

number of samples to be used to estimate a bias term if ties = TRUE.

Value

An object of class htest which is a list, some of the components of which are

statistic

value of the test statistic.

p.value

corresponding p-value.

Details

The code for this test was generously provided by Johanna G. Ne<U+0161>lehov<U+00E1>. More details are available in Appendix B of Ben Ghorbal, G. Ne<U+0161>lehov<U+00E1> and Genest (2009).

References

Ghorbal, M. B., Genest, C., and G. Ne<U+0161>lehov<U+00E1>, J. (2009) On the test of Ghoudi, Khoudraji, and Rivest for extreme-value dependence. The Canadian Journal of Statistics 37, 1--9.

Kojadinovic, I. (2017). Some copula inference procedures adapted to the presence of ties. Computational Statistics and Data Analysis 112, 24--41, http://arxiv.org/abs/1609.05519.

See Also

evTestC, evTestA, evCopula, gofEVCopula, An.

Examples

# NOT RUN {
set.seed(321)
## Do the data come from an extreme-value copula?
evTestK(Ug <- rCopula(200, gumbelCopula(3))) # not significant => yes, EV
    dim(Uc <- rCopula(200, claytonCopula(3)))
## Clayton:                       tests are highly significant => no, not EV
(K1 <- evTestK(Uc))
(K2 <- evTestK(Uc, method = "asymptotic"))
# }
# NOT RUN {
system.time(print(K3 <- evTestK(Uc, method = "jackknife")))
## slower: 1.06 sec (2015 intel core i7)
# }
# NOT RUN {
<!-- %\dont show -->
# }