copBasic (version 2.0.1)

statTn: The Tn Statistic of a Fitted Copula to the Empirical Copula

Description

Compute the $T_n(p)$ statistic of Genest et al. (2011) that is defined as $$T_n(p) = \sum_{i=1}^n \big|\mathbf{C}_n(u_i, v_i) - \mathbf{C}_{\Theta_n}(u_i, v_i)\big|^p\mbox{,}$$ where $\mathbf{C}_n(u,v)$ is the empirical copula, $\mathbf{C}_{\Theta_n}(u,v)$ is the copula with estimated parameters $\Theta_n$ from the sample of size $n$. The $T_n$ for $p = 2$ is reported by those authors to be of general purpose and overall performance in large scale simulation studies. The extension here for arbitary exponent $p$ is made for flexibility. Alternatively the definition could be associated with the statistic $T_n(p)^{1/p}$ in terms of a root $1/p$ of the summation as shown above.

The $T_n$ statistic is obviously a form of deviation between the empirical (nonparametric) and parametric fitted copula. The distribution of this statistic through Monte Carlo simulation could be used for inference that a chosen parametric model is suitably close to the empirical copula. This statistic has an advantage of being relatively straightforward to understand and explain to stakeholders and decision makers, is attractive for being suitable in a wide variety of circumstances, but intuitively might have limited statistical power in some situations for it looks at whole copula structure and note say on tail dependency.

Usage

statTn(uv, cop=NULL, para=NULL, p=2, proot=FALSE, ...)

Arguments

uv

An Rtwo column matrix or data.frame of a sample of nonexceedance probabilities $u$ and $v$;

cop

A copula function;

para

Vector of parameters or other data structure, if needed, to pass to the copula;

p

The value for $p$, and the default follows that of Genest et al. (2011);

proot

A logical controling whether the $T_n$ returned be rooted by $1/p$, and the default follows that of Genest et al. (2011); and

...

Additional arguments to pass to the copula function and (or) the empirical copula.

Value

• The value for $T_n$ is returned dependent on the specification of $p$ and whether rooting of the result is desired.

encoding

utf8

References

Genest, C., Kojadinovic, I., Nešlehová{Neslehova}, J., and Yan, J., 2011, A goodness-of-fit test for bivariate extreme-value copulas: Bernoulli, v. 17, no. 1, pp. 253--275.

See Also

EMPIRcopdf, vuongCOP, kullCOP

Examples

# Example here is just for Tn. For the example below, the PSP copula is radically
# different from the Gumbel-Hougaard copula and thus, the hatTn would be expected
# to be quite different from those of the Gumbel-Hougaard and certainly not too
# near to zero. The Bernstein extension to the empirical copula is CPU intensive.
samUV  <- simCOP(n=40, cop=PSP, graphics=FALSE)   # random sample
hatTau <- cor(samUV$U, samUV$V, method="kendall") # Kendall Tau
hatTn  <- statTn(sampleUV, cop=GHcop, para=GHcop(tau=hatTau)$para,
                 ctype="bernstein") # approximate range 0.0141 -- 0.030
# hatTn in this case is cernby itself is somewhat uninformative and requires
# Monte Carlo to put an individual value into context.