H_inv: Inverse survival function if Monte-Carlo approximation is set...

Description Usage Arguments Value Examples

View source: R/liouville.R

Description

The function is used internally for optimization.

Usage

1
H_inv(u, alphavec, family, theta, MC = 1e+05, TRUNC = FALSE)

Arguments

u

data at which to compute the survival inverse

alphavec

vector of Dirichlet allocations (must be a vector of integers)

family

family of the Liouville copula. Either "clayton", "gumbel", "frank", "AMH" or "joe"

theta

parameter of the corresponding Archimedean copula

MC

number of Monte-Carlo points for evaluation

trunc

whether to truncate at low quantile. Will be based on numerical root finding for the lower 0.025 fraction of the data

Value

Inverse survival function values

Examples

1
2
3
4
5
6
7
8
## Not run: 
u <- rliouv(n = 100, family = "frank", alphavec <- c(2,3), theta = 1)
H_inv(u=u, family="frank", alphavec=c(2,3), theta=2)
#Difference between true value and approximation (can be large depending on family)
sum(abs(H_inv(u=u, family="frank", alphavec=c(2,3), theta=2)-
isliouv_m(u=u, family="frank", alphavec=c(2,3), theta=2)))

## End(Not run)

lcopula documentation built on May 30, 2017, 4:36 a.m.

Search within the lcopula package
Search all R packages, documentation and source code