copula (version 0.99-0)

psiDabsMC: Absolute Value of Generator Derivatives via Monte Carlo

Description

Computes the absolute values of the $d$th generator derivative $\psi^{(d)}$ via Monte Carlo simulation.

Usage

psiDabsMC(t, family, theta, degree = 1, n.MC,
          method = c("log", "direct", "pois.direct", "pois"),
          log = FALSE)

Arguments

t

numeric vector of evaluation points.

family

Archimedean family (name or object).

theta

parameter value.

degree

order $d$ of the derivative.

n.MC

Monte Carlo sample size.

method

different methods: [object Object],[object Object],[object Object],[object Object]

log

if TRUE the logarithm of psiDabs is returned.

Value

• numeric vector of the same length as t containing the absolute values of the generator derivatives.

Details

The absolute value of the $d$th derivative of the Laplace-Stieltjes transform $\psi=\mathcal{LS}[F]$ can be approximated via $$(-1)^d\psi^{(d)}(t)=\int_0^\infty x^d\exp(-tx)\,dF(x)\approx\frac{1}{N}\sum_{k=1}^NV_k^d\exp(-V_kt),\ t > 0,$$ where $V_k\sim F,\ k\in{1,\dots,N}$. This approximation is used where $d=$degree and $N=$n.MC. Note that this is comparably fast even if t contains many evaluation points, since the random variates $V_k\sim F,\ k\in{1,\dots,N}$ only have to be generated once, not depending on t.

References

Hofert, M., Mächler{Maechler}, M., and McNeil, A. J. (2011a), Estimators for Archimedean copulas in high dimensions: A comparison, to be submitted.

See Also

acopula-families.

Examples

t <- c(0:100,Inf)
set.seed(1)
psiDabsMC(t, family="Gumbel", theta=2, degree=10, n.MC=10000, log=TRUE)
## Note: The absolute value of the derivative at 0 should be Inf for
## Gumbel, however, it is always finite for the Monte Carlo approximation