I_1 {PosRatioDist}R Documentation

Lemma

Description

Technical Lemmas for calculating quotient of random variables conditioned to the positive quadrant.For more detailed information please read the first reference paper section 2.2.

Usage

I_1(a, b)

I_2(a, b)

I_3(a, b)

J_1(a, b, c, alpha)

J_2(a, b, c, alpha)

J_3(a, b, c, alpha)

Arguments

a

parameter

b

parameter

c

parameter

alpha

parameter

Details

InI_n Type I Integration

In(a,b)=0ynexp(ay2by)dyI_n (a, b) = \int_0^\infty y^n \exp \left( -a y^2 - b y \right) dy

For <a<,<b<-\infty < a < \infty,-\infty < b < \infty,where n is positive integer.

In particular,for a>0a > 0,we have expressions below

I1(a,b)=πb4a3/2exp(b24a) erfc(b2a)+12aI_1 (a, b) = -\frac {\sqrt{\pi} b}{4 a^{3 / 2}}\exp \left( \frac {b^2}{4 a} \right)\ {\rm erfc} \left( \frac {b}{2 \sqrt{a}} \right) + \frac {1}{2 a}

I2(a,b)=π4a3/2exp(b24a) erfc(b2a)+πb28a5/2exp(b24a) erfc(b2a)b4a2I_2 (a, b) = \frac {\sqrt{\pi}}{4 a^{3 / 2}}\exp \left( \frac {b^2}{4 a} \right)\ {\rm erfc} \left( \frac {b}{2 \sqrt{a}} \right) +\frac {\sqrt{\pi} b^2}{8 a^{5 / 2}} \exp \left( \frac {b^2}{4 a} \right)\ {\rm erfc} \left( \frac {b}{2 \sqrt{a}} \right) - \frac {b}{4 a^2}

I3(a,b)=3πb8a5/2exp(b24a) erfc(b2a)πb316a7/2exp(b24a) erfc(b2a)+12a2+b28a3I_3 (a, b) = -\frac {3 \sqrt{\pi} b}{8 a^{5 / 2}}\exp \left( \frac {b^2}{4 a} \right)\ {\rm erfc} \left( \frac {b}{2 \sqrt{a}} \right) -\frac {\sqrt{\pi} b^3}{16 a^{7 / 2}}\exp \left( \frac {b^2}{4 a} \right)\ {\rm erfc} \left( \frac {b}{2 \sqrt{a}} \right) + \frac {1}{2 a^2} + \frac {b^2}{8 a^3}

JnJ_n Type J Integration

Jn(a,b,c,α)=0yn(ay2+by+c)αdyJ_n (a, b, c, \alpha) = \int_0^\infty y^n \left( a y^2 + b y + c \right)^{-\alpha} dy

In particular,for a>0,b2<4ac,1<n<2α1a > 0,b^2 < 4ac, -1 < n < 2\alpha - 1,we have expressions below

J1(a,b,c,α)=a1c1αB(2,2α2) 2F1(1,α1;α+12;1b24ac)J_1 (a, b, c, \alpha) = a^{-1} c^{1 - \alpha} B \left( 2, 2 \alpha - 2 \right) \ {}_2F_1 \left( 1, \alpha - 1; \alpha + \frac {1}{2}; 1 - \frac {b^2}{4 a c} \right)

J2(a,b,c,α)=a32c32αB(3,2α3) 2F1(32,α32;α+12;1b24ac)J_2 (a, b, c, \alpha) = a^{-\frac {3}{2}} c^{\frac {3}{2} - \alpha} B \left( 3, 2 \alpha - 3 \right) \ {}_2F_1 \left( \frac {3}{2}, \alpha - \frac {3}{2}; \alpha + \frac {1}{2}; 1 - \frac {b^2}{4 a c} \right)

J3(a,b,c,α)=a2c2αB(4,2α4) 2F1(2,α2;α+12;1b24ac)J_3 (a, b, c, \alpha) = a^{-2} c^{2 - \alpha} B \left( 4, 2 \alpha - 4 \right) \ {}_2F_1 \left( 2, \alpha - 2; \alpha + \frac {1}{2}; 1 - \frac {b^2}{4 a c} \right)

Value

I_1 gives value of Type I integration with n=1n = 1

I_2 gives value of Type I integration with n=2n = 2

I_3 gives value of Type I integration with n=3n = 3

J_1 gives value of Type J integration with n=1n = 1

J_2 gives value of Type J integration with n=2n = 2

J_3 gives value of Type J integration with n=3n = 3

Invalid arguments will return an error message.

Author(s)

Saralees Nadarajah & Yuancheng Si siyuanchengman@gmail.com

References

Yuancheng Si and Saralees Nadarajah and Xiaodong Song, (2020). On the distribution of quotient of random variables conditioned to the positive quadrant. Communications in Statistics - Theory and Methods, 49, pp2514-2528.

Balakrishna, N. and Shiji, K. (2014). On a class of bivariate exponential distributions.Statistics and Probability Letters, 85, pp153-160.

Arnold, B. C. and Strauss, D. (1988).Pseudolikelihood estimation.Sankhya B , 53, pp233-243.

Caginalp, C. and Caginalp, G. (2018).The quotient of normal random variables and application to asset price fat tails.Physica A—Statistical Mechanics and Its Applications, 499, pp457-471.

Louzada, F., Ara, A. and Fernandes, G. (2017).The bivariate alpha-skew-normal distribution.Communications in Statistics - Theory and Methods, 46, pp7147-7156.

Nadarajah, S. (2009).A bivariate Pareto model for drought.Stochastic Environmental Research and Risk Assessment, 23, pp811-822.

Nadarajah, S. and Kotz, S. (2006).Reliability models based on bivariate exponential distributions.Probabilistic Engineering Mechanics, 21, pp338-351.

Nadarajah, S. and Kotz, S. (2007).Financial Pareto ratios.Quantitative Finance, 7, pp257-260.

Examples

I_1(1,2)
I_2(1,2)
I_3(1,2)
J_1(1,2,3,3)
J_2(1,2,3,3)
J_3(1,2,3,3)










[Package PosRatioDist version 1.2.1 Index]