Evaluates the h-function (or its inverse) corresponding to a kdecopula object. H-functions are conditional distribution functions obtained by integrating the copula density w.r.t. to one of its arguments (see also VineCopula::BiCopHfunc().

hkdecop(u, obj, cond.var, inverse = FALSE)

Arguments

u

\(n x 2\) matrix of evaluation points.

obj

kdecopula object.

cond.var

integer; cond.var = 1 conditions on the first variable, cond.var = 2 on the second.

inverse

logical; indicates whether the h-function or its inverse shall be calculated.

Value

A length \(n\) vector of the (inverse) h-function evaluated at u.

Examples

## load data and transform with empirical cdf data(wdbc) udat <- apply(wdbc[, -1], 2, function(x) rank(x) / (length(x) + 1)) ## estimation of copula density of variables 5 and 6 fit <- kdecop(udat[, 5:6]) plot(fit)
## evaluate h-function estimate and its inverse at (u1|u2) = (0.123 | 0.321) hkdecop(c(0.123, 0.321), fit, cond.var = 2)
#> [1] 0.116027
hkdecop(c(0.123, 0.321), fit, cond.var = 2, inverse = TRUE)
#> [1] 0.1283875