Build data and AI skills | 50% off

Last chance! 50% off unlimited learning

Sale ends in


copBasic (version 2.2.7)

RAYcop: The Rayleigh Copula

Description

The Rayleigh copula (Boškoskia and others, 2018) is CΘ(u,v)=RAY(u,v;Θ)=1+AB, A=eΘa2a2(ea10Θa2esI0(2a1s)ds1)ds, B=ea10a2esI0(2Θa1s)ds, where a1=log(1u)/(1Θ), a2=log(1v)/(1Θ), Iν(x) is the modified Bessel function of the first kind of order ν (see base::besselI()), and Θ(0,1]. The copula, as Θ0+ limits, to the independence coupla (Π(u,v); P) and as Θ1 limits to the comonotonicity copula (M(u,v); M). Finally, there are formulations of the Rayleigh copula using the Marcum-Q function, but the copBasic developer has not been able to make such work. If the Marcum-Q function could be used, then only one integration and not the two involving the modified Bessel function are possible. Infinite integrations begin occurring in the upper right corner for about Θ>0.995 at which point the M(u,v) copula is called in the source code.

Usage

RAYcop(u, v, para=NULL, rho=NULL, method=c("default"),
             rel.tol=.Machine$double.eps^0.5, ...)

Value

Value(s) for the copula are returned.

Arguments

u

Nonexceedance probability u in the X direction;

v

Nonexceedance probability v in the Y direction;

para

A vector (single element) of parameters---the Θ parameter of the copula;

rho

Value for Spearman Rho from which parameter Θ is computed by polynomial approximation and returned. The estimation appears sufficient for most pratical applications (see Examples);

method

The computational method of integrals associated with the definition of the copula; this is designed for the ability to switch eventually in sources to Marcum-Q function implementation. The definition in January 2023 and default is to call the two Bessel function integrals shown for the definition in this documentation;

rel.tol

Argument of the same name for integrate() call; and

...

Additional arguments to pass.

Author

W.H. Asquith

References

Boškoskia, P., Debenjaka, A., Boshkoskab, B.M., 2018, Rayleigh copula for describing impedance data with application to condition monitoring of proton exchange membrane fuel cells: European Journal of Operational Research, v. 266, pp. 269--277, tools:::Rd_expr_doi("10.1016/j.ejor.2017.08.058").

Zeng, X., Ren, J., Wang, Z., Marshall, S., and Durrani, T., [undated], Copulas for statistical signal processing (Part I)---Extensions and generalization, accessed January 14, 2024, at https://pure.strath.ac.uk/ws/portalfiles/portal/34078849/Copulas_Part1_v2_6.pdf.

Zeng, X., Ren, J., Sun, M., Marshall, S., and Durrani, T., [undated], Copulas for statistical signal processing (Part II)---Simulation, optimal selection and practical applications, accessed January 14, 2024, at https://strathprints.strath.ac.uk/48371/1/Copulas_Part2s_v2_5_2.pdf

See Also

M, P