Bivariate Sinh-Normal Distribution and A Related Model
|
|
- Martina McKinney
- 6 years ago
- Views:
Transcription
1 Bivariate Sinh-Normal Distribution and A Related Model Debasis Kundu 1 Abstract Sinh-normal distribution is a symmetric distribution with three parameters. In this paper we introduce bivariate sinh-normal distribution, which has seven parameters. Due to presence of seven parameters it is a very flexible distribution. Different properties of this new distribution has been established. The model can be obtained as a bivariate Gaussian copula also. Therefore, using the Gaussian copula property, several properties of this proposed distribution can be obtained. Maximum likelihood estimators cannot be obtained in closed forms. We propose to use two step estimators based on Copula, which can be obtained in a more convenient manner. One data analysis has been performed to see how the proposed model can be used in practice. Finally, we consider a bivariate model which can be obtained by transforming the sinh-normal distribution and it is a generalization of the bivariate Birnbaum-Saunders distribution. Several properties of the bivariate Birnbaum-Saunders distribution can be obtained as special cases of the proposed generalized bivariate Birnbaum-Saunders distribution. AMS Subject Classification: Primary 6E15, Secondary 6H10. Key Words and Phrases: Birnbaum-Saunders distribution; bivariate Birnbaum-Saunders distribution; log-birnbaum-saunders distribution; maximum likelihood estimators; copula; Two stage estimators; Total positivity of order two. 1 Department of Mathematics and Statistics, Indian Institute of Technology Kanpur, Pin 08016, India. kundu@iitk.ac.in, Phone no , Fax no
2 1 Introduction Rieck and Nedelman [13], introduced sinh-normal (SHN) distribution, as a three parameter symmetric model with the following cumulative distribution function (CDF) F Y (y;α,σ,µ) = Φ(a(y;α,σ,µ)); y R. (1) Here Φ( ) is the CDF of a standard normal random variable, α > 0, σ > 0, < µ <, sinh(x) is the hyperbolic sine function of x where sinh(x) = (e x e x )/, and a(y;α,σ,µ) = ( ) y µ α sinh. σ Note that a(y;α,σ,µ) is an increasing function of y, and as y, a(y;α,σ,µ), and as y, a(y;α,σ,µ), therefore, F Y ( ) as defined in (1) is a proper distribution function. From now on a three-parameter sinh-normal random variable with CDF (1) will be denoted by SHN(α,σ,µ). It has several important properties, and a brief review of SHN distribution will be provided in Section. The main aim of this paper is to introduce bivariate sinh-normal (BSHN) distribution, which is a natural extension of SHN distribution of one dimension to two dimensions. Recently, Kundu et al. [7] introduced a bivariate Birnbaum-Saunders distribution. It is observed that the bivariate log-birnbaum-saunders distribution can be obtained as a special case of the BSHN distribution. The proposed BSHN distribution has seven parameters, and due to presence of two location, two shapes, two scales and one correlation parameters, it is a very flexible model. After proper normalization, as the shape parameters converge to zero, it approaches to a standard bivariate normal distribution. The probability density function (PDF) of BSHN, can be both unimodal or bimodal depending on the shape parameters. Due to presence of the two shape parameters, it is more flexible than the bivariate normal distribution. The marginals of BSHN are SHN
3 3 distributions, and the conditional distribution also can be obtained in closed form. The generation from a BSHN distribution is quite straight forward, hence simulation experiments can be performed in a routine matter. The proposed BSHN model can be obtained as a bivariate Gaussian copula, hence different copula properties can be easily incorporated for BSHN distribution. The model has the total positivity of order two (TP ) or reverse rule of order two (RR ) properties depending on the correlation parameter. The maximum likelihood estimators of the unknown parameters cannot be obtained in closed form, as expected. Nonlinear optimization method needs to be used to compute the MLEs. We propose to use two stage estimators as proposed by Joe [6]. One data analysis has been performed for illustrative purposes. Finally we introduce a new bivariate distribution, which can be obtained by transforming the BSHN random variable. The new distribution is a generalization of the bivariate Birnbaum-Saunders distribution introduced by Kundu et al. [7], and we call it as the bivariate generalized Birnbaum-Saunders distribution. We establish different properties of the generalized Birnbaum-Saunders distribution, and it is observed that several properties of the bivariate Birnbaum-Saunders distribution can be obtained as special cases of the proposed distribution. Finally we indicate multivariate generalization of the proposed models. Rest of the paper is organized as follows. In Section, we provide some preliminaries. The proposed BSHN distribution is defined and its properties are discussed in Section 3. The statistical inference of the unknown parameters are discussed in Section 4. The analysis of a data set is provided in Section 5. In Section 6, we proposes generalized Birnbaum-Saunders distribution, and finally we conclude the paper in Section 6.
4 4 Preliminaries.1 Sinh-Normal Distribution Suppose Y SHN(α,σ,µ), then the CDF of Y is provided in (1), and the corresponding probability density function (PDF) becomes f Y (y;α,σ,µ) = φ(a(y;α,σ,µ))a(y;α,σ,µ); y R, here φ( ) is the PDF of a standard normal random variable, and A(y;α,σ,µ) = d dy a(y;α,σ,µ) = ( ) y µ ασ cosh = e(y µ)/σ +e (y µ)/σ. σ ασ By simple transformation, it follows that if Y SHN(α,σ,µ), then Z = ( ) y µ α sinh N(0,1). () σ Here N(0,1) denotes a standard normal random variable. From (), it follows that if Z N(0,1), then Y = σ arcsinh ( ) αz +µ SHN(α,σ,µ), (3) where arcsinh (x) = ln(x + x +1). Using (3), random samples from SHN distribution can be easily generated. Moreover, from (3), it follows that if Y SHN(α,σ,µ), then cy SHN(α, c σ, cµ). Rieck [1] has established the following properties of a SHN distribution. The PDF of SHN is symmetric about the location parameter µ. The distribution is strongly unimodal for α, and it is bimodal if α >. It can be easily seen using L Hospital s rule, that if Y SHN(α,σ,µ), then (Y µ) ασ N(0,1) as α 0. (4) Another important property of the SHN distribution is that if Y SHN(α,e µ,), then Y has a log-birnbaum-saunders distribution. The moment generating function of a SHN distribution cannot be obtained in closed form.
5 5. Birnbaum-Saunders Distribution Birnbaum and Saunders [1] introduced a two-parameter failure time distribution for fatigue failure due to cyclic loading. The cumulative distribution function(cdf) of a two-parameter Birnbaum-Saunders random variable T for α > 0 and β > 0, can be written as [ { (t ) 1/ ( ) }] 1/ 1 β F T (t;α,β) = Φ. (5) α β t It is an absolute continuous distribution with the PDF [ (β ) 1/ ( ) ] 3/ [ 1 β f T (t;α,β) = + exp 1 ( t παβ t t α β + β )] t, t > 0. Here α is the shape parameter and β is the scale parameter. From now on a two-parameter Birnbaum-Saunders random variable with CDF (5) will be denoted by BS(α, β). If the following transformation is used, i.e., { (T Z = 1 ) 1/ α β ( ) } 1/ β, (6) T then Z follows standard normal distribution. Several properties and several measures of a Birnbaum-Saunders random variable can be obtained using the transformation (6)..3 Copula It is well known that the dependence among two or more random variables, say X 1,...,X p is completely determined by its joint distribution function F X1,...,X p (x 1,...,x p ). The idea of separating F X1,...,X P (x 1,...,x p ) in two parts - the one which describes the dependence structure, and the other one which describes only the marginal behavior, leads to the concept of copula. A p-variate copula defined on [0,1] p is a multivariate distribution with univariate marginals on [0,1]. Let X 1,...,X p be p random variables with continuous distribution functions F X1 ( ),...,F Xp ( ), respectively, then according to Sklar s theorem, F X1,...,X p (x 1,...,x p )
6 6 has a unique copula representation; Moreover, F X1,...,X p (x 1,...,x p ) = C(F X1 (x 1 ),...,F Xp (x p )). C(u 1,...,u p ) = F X1,...,X p (F 1 X 1 (u 1 ),...,F 1 X p (u p )). It is well known that many dependence properties of a multivariate distribution depend only on the corresponding copula. Therefore, many dependence properties of a multivariate distribution can be obtained by studying the corresponding copula. The bivariate Gaussian copula is defined as follows; C G (u,v) = Φ 1 (u) Φ 1 (v) φ (x,y;ρ)dxdy = Φ (Φ 1 (u),φ 1 (v);ρ). Here φ( ), Φ( ) are same as defined before, and φ (u,v;ρ) denotes the standard bivariate normal density function; { } 1 φ (u,v;ρ) = π 1 1 ρ exp (1 ρ ) (u +v ρuv). The bivariate Gaussian copula density can be obtained as c G (u,v;ρ) = = u v C G(u,v;ρ) = φ (Φ 1 (u),φ 1 (v);ρ) φ(φ 1 (u))φ(φ 1 (v)) ( ) 1 ρφ 1 exp (u)φ 1 (v) ρ (Φ 1 (u) +Φ 1 (v) ). 1 ρ (1 ρ ) We recall that a non-negative function g defined in R is total positivity of order two, abbreviated by TP if for all x 1 < x, y 1 < y, with x,y R g(x 1,y 1 )g(x,y ) g(x,y 1 )g(x 1,y ). (7) If the equality (7) is reversed, it is called reverse rule of order two (RR ). The following result will be useful for future development, and it may have some independent interest also.
7 7 Result : The Gaussian copula density is (a) TP for 0 < ρ < 1, (b) RR if 1 < ρ < 0. Proof: To prove (a), we need to show that for all a 1,a,b 1,b R, if 0 < a 1 < a < 1, and 0 < b 1 < b < 1, then for 0 < ρ < 1, c G (a 1,b 1 ;ρ)c G (a,b ;ρ) c G (a,b 1 ;ρ)c G (a 1,b ;ρ). (8) Proving (8) is equivalent to prove that for all < x 1 < x <, < y 1 < y < φ (a 1,b 1 ;ρ)φ (a,b ;ρ) φ (a,b 1 ;ρ)φ (a 1,b ;ρ). (9) Now considering all possible six cases, like (i) x 1 < x < y 1 < y, (ii) x 1 < y 1 < x < y, etc., (9) can be easily verified for 0 < ρ < 1. Similarly, (b) can be easily obtained along the same line. 3 Bivariate Sinh-Normal Distribution In this section, we introduce bivariate sinh-normal distribution, and study its different properties. The bivariate random vector (Y 1,Y ) is said to have bivariate sinh-normal distribution with parameters α 1,σ 1,µ 1,α,σ,µ,ρ, if the CDF of Y 1 and Y is P(Y 1 y 1,Y y ) = Φ (a(y 1 ;α 1,σ 1,µ 1 ),a(y ;α,σ,µ );ρ); y 1,y R. (10) Here α 1 > 0, α > 0, σ 1 > 0, σ > 0, < µ 1 <, < µ <, 1 < ρ < 1, and Φ (u,v,ρ) is the CDF of a standard bivariate normal vector (Z 1,Z ), with correlation coefficient ρ. The corresponding PDF of Y 1 and Y is given by f Y1,Y (y 1,y ) = φ (a(y 1 ;α 1,σ 1,µ 1 ),a(y ;α,σ,µ );ρ) A(y i ;α i,σ i,µ i ); y 1,y R. Here a( ),A( ) and φ ( ) are same as defined before. From now on a BSHN random variable with CDF (10) will be denoted by BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ).
8 8 It is clear that µ 1, µ are the two location parameters, and σ 1, σ are the two scale parameters. If µ 1 = µ = 0 and σ 1 = σ = 1, it is called a standard BSHN distribution. A standard BSHN random variable has the PDF which is centered at (0,0), and it is symmetric around (0,0). It is observed that the surface of the PDF can be unimodal, bimodal or multimodal depending on the parameter values. It is immediate that if (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), and Z i = a(y i ;α i,σ i,µ i ); i = 1,, then (Z 1,Z ) has a standard bivariate normal distribution with correlation coefficient ρ. Moreover, if (Z 1,Z ) has a standard bivariate normal distribution with correlation coefficient ρ, and if Y i = σ i arcsinh ( ) αi Z i +µ i ; i = 1,, (11) then (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ). The representation (11) can be used to generate samples from a BSHN distribution. We present the following simple algorithm to generate samples from BSHN distribution. Algorithm Step 1: Generate independent U 1 and U from N(0,1). Step : Compute Z 1 = Z = 1+ρ+ 1 ρ 1+ρ 1 ρ U 1 + U 1+ρ 1 ρ 1+ρ+ 1 ρ U 1 + U. Step 3: Y i = σ i arcsinh ( ) αi Z i +µ i ; i = 1,.
9 { ( ) ( )} Y1 µ 1 Y µ ItisfurtherfollowsusingL Hospital srulethatasα 1 0,α 0,, α 1 σ 1 α σ converges to a standard bivariate normal distribution with correlation coefficient ρ. 9 The following theorem provides the marginal and conditional distributions of the BSHN distribution. Theorem 3.1: If (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), then (a) Y i SHN(α i,σ i,µ i ); i = 1,. (b) The conditional CDF of Y 1, given Y = y, is given by { } a(y : α 1,σ 1,µ 1 ) ρa(y ;α,σ,µ ) P(Y 1 y 1 Y = y ) = Φ 1 ρ (c) The conditional PDF of Y 1, given Y = y is { } a(y;α 1,σ 1,µ 1 ) ρa(y ;α,σ,µ ) 1 f Y1 Y =y (y 1 ) = φ A(y 1;α 1,σ 1,µ 1 ). 1 ρ 1 ρ (d) The conditional hazard of Y 1 at y 1, given Y = y is an increasing function of y 1, for y 1 > 0, for all values of y and ρ. Proof: The proof of (a) can be obtained from the definition. The proof of (b) can be obtained by taking the transformation u = a(y 1 ;α 1,σ 1,µ 1 ) and v = u ρa(y ;α,σ,µ ). 1 ρ Part (c) can be obtained from (b) by differentiating with respect to y 1. To prove (d), note that the conditional hazard of Y 1 at y, given Y = y can be written as h Y1 Y =y (y) = φ(b(y)) 1 Φ(b(y)) 1 1 ρ A(y;α 1,σ 1,µ 1 ), where b(y) = a(y;α 1,σ 1,µ 1 ) ρa(y ;α,σ,µ ). 1 ρ
10 10 Since b(y) is an increasing function of y, A(y;α 1,σ 1,µ 1 ) is an increasing function of y for y > 0 and the hazard function of a normal distribution is an increasing function, the result follows. then Consider the following random variable U(y ) = {Y 1 Y = y }. If c = ρa(y ;α,σ,µ ), 1 ρ P(U(y ) y) = Φ(a(y;α 1 1 ρ,σ 1,µ 1 )+c). Therefore, U(y ) has the following representation ( ) α 1 1 ρ V U(y ) = σ 1 arcsinh +µ 1, (1) where V N(c,1). Using (1), different properties of the conditional distribution can be obtained. The following theorem provides the symmetric properties of the BSHN random variable. Theorem 3.: If (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), and 1 if x > 0 sgn(x) = 0 if x = 0 1 if x < 0 then (a) (c 1 Y 1,c Y ) BSHN(α 1, c 1 σ 1,c 1 µ 1,α, c σ,c µ,sgn(c 1 c )ρ). (b) (Y 1, Y ) BSHN(α 1,σ 1,µ 1,α,σ, µ, ρ). (c) ( Y 1,Y ) BSHN(α 1,σ 1, µ 1,α,σ,µ, ρ). (d) ( Y 1, Y ) BSHN(α 1,σ 1, µ 1,α,σ, µ,ρ). Proof: The proof of (a) follows from the joint PDF of (Y 1,Y ) and using the appropriate transformation. The results of (b), (c) and (d) follow from (a).
11 11 Since F Yi (y i ) = Φ(a(y i ;α i σ i µ i )), for i = 1,, therefore, F Y1,Y (y 1,y ) = Φ (a(y 1 ;α 1,σ 1,µ 1 ),a(y ;α,σ,µ );ρ) = C G (Φ 1 (Φ(a(y 1 ;α 1,σ 1,µ 1 ))),Φ 1 (Φ(a(y ;α,σ,µ )));ρ) = C G (F Y1 (y 1 ),F Y (y );ρ). Therefore, BSHN distribution can be obtained as a bivariate Gaussian copula, with marginals as SHN distributions. We immediately have the following results. Theorem 3.3: If (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), then for ρ > 0, it is TP and for ρ < 0, it is RR. Proof: Since both TP and RR are copula properties, the results immediately follow from Result. Theorem 3.4: If (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), then for ρ > 0 (a) T 1 is stochastically increasing in T (b) T is stochastically increasing in T 1, for all values of α 1, β 1, α and β. Proof: We will prove (a), (b) follows along the same line. Note that the results can be established if we can show that C G (u,v;ρ) is a concave function in u for fixed v when ρ > 0, see Nelsen ([10], page 197). It is equivalent to prove that u C G(u,v;ρ) is a decreasing function in u for fixed v. Using Meyer [9], we have ( ) u C(u,v;ρ) = Φ Φ 1 (v) ρφ 1 (u). 1 ρ Clearly, for ρ > 0, the right hand side is a decreasing function in u for fixed v and the result follows. The bivariate hazard rate of Y 1 and Y is defined as h(y 1,y ) = (, ) lnp(y 1 > y 1,Y > y ) = (h 1 (y 1,y ),h (y 1,y )) y 1 y
12 1 see Marshall [8]. Moreover this bivariate vector valued function h(y 1,y ) uniquely determines the probability distribution. We have the following result. Theorem 3.4: If (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), then h 1 (y 1,y ) is an increasing function of y 1, for y 1 > 0, and for all values of y. Similarly, h (y 1,y ) is an increasing function of y for y > 0 for all values of y 1. Proof: Using similar technique as in Gupta and Gupta [3], it can be easily shown that { ( )} a(y 1 Φ ;α,σ,µ ) ρa(y 1 :α 1,σ 1,µ 1 ) φ(a(y 1 ;α 1,σ 1,µ 1 ) 1 ρ h 1 (y 1,y ) = Φ (a(y 1 ;α 1,σ 1,µ 1 ),a(y ;α,σ,µ );ρ) A(y 1;α 1,σ 1,µ 1 ). Since a(y 1 ;α 1,σ 1,µ 1 ) is an increasing function of y 1, and the hazard gradient of a bivariate normal distribution is an increasing function, see Gupta and Gupta [3], it follows that the function within [ ] is an increasing function of y 1. The first part of the result follows by observing the fact that the function A(y 1 ;α 1,σ 1,µ 1 ) is an increasing function of y 1 for y 1 > 0. Similarly, the second part also can be established. From the Gaussian copula, it follows that if (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), then for all values of α 1, α, σ 1, σ, µ 1 and µ, (a) Kendall s tau and (b) Spearman s rho becomes τ = π sin 1 (ρ), and ρ S = 6 π sin 1 ( ρ ), respectively. Now we provide the medial correlation of the BSHN distribution using the copula property. If M Y1 and M Y denote the medians of Y 1 and Y, respectively, then M Y1,Y, the medial correlation coefficient of Y 1 and Y is M Y1,Y = P[(Y 1 M Y1 )(Y M Y ) > 0] P[(Y 1 M Y1 )(Y M Y ) < 0]. It has been shown by Nelsen [10], that the medial correlation coefficient is a copula property
13 13 and M Y1,Y = 4C M Y1,Y = 4C G ( 1, 1 ( 1, 1 ) - 1. Therefore, of (Y 1,Y ) follows BSHN, then ) ( ) 1 1 = 4Φ (0,0,ρ) 1 = sin 1 ρ π 1 = sin 1 ρ. π The concept of bivariate tail dependence relates to the amount of dependence in the upper-quadrant (or lower quadrant) tail of bivariate distribution, see Joe ([5], page 33). In terms of original random variables Y 1 and Y, the upper tail dependence is defines as χ = lim z 1 P(Y F 1 Y (z) Y 1 F 1 Y 1 (z)). Intuitively, the upper-tail dependence exists when there is a positive probability that some positive outliers may occur jointly. Since the Gaussian copula is upper tail and lower tail independent, it follows that if (Y 1,Y ) BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ), then Y 1 and Y are upper tail and lower tail independent. 4 Inference 4.1 Maximum Likelihood Estimators In this section we discuss the estimation of the unknown parameters based on a sample of size n, {(y 11,y 1 ),...,(y n1,y n )} from BSHN(α 1,σ 1,µ 1,α,σ,µ,ρ). Let us denote θ = (α 1,σ 1,µ 1,α,σ,µ,ρ). The log-likelihood function becomes l(θ) = lnφ (a(y i1 ;α 1,σ 1,µ 1 ),a(y i ;α,σ,µ );ρ)+ j=1 lna(y ij ;α j,σ j,µ j ). (13) The MLEs of the unknown parameters can be obtained by maximizing (13) with respect to the unknown parameters. It is a seven-dimensional optimization problem. Clearly, they cannot be obtained in explicit forms. To avoid solving the seven dimensional optimization problem, we propose the profile log-likelihood function approach. Since ) )} ( )} {(e Y 1 µ 1 σ 1 e Y 1 µ 1 σ 1, (e Y µ σ e Y µ σ α N {(0,0), 1 α 1 α ρ α 1 α ρ α,
14 14 for fixed σ 1,µ 1,σ,µ, the MLEs of α 1, α and ρ are ( 1/ 1 α 1 (σ 1,µ 1 ) = g (y i1 ;σ 1,µ 1 )), α (σ,µ ) = n and where ρ(σ 1,µ 1,σ,µ ) = ( 1 n g (y i ;σ,µ ) n g(y i1;σ 1,µ 1 )g(y i ;σ,µ ) n g (y i1 ;σ 1,µ 1 ) n g (y i ;σ,µ ), ) 1/ g(y;σ,µ) = e y µ σ e y µ σ. (14) Observe that α 1 (σ 1,µ 1 ) is a function of σ 1 and µ 1, and α (σ.µ ) is a function of σ and µ only. Finally, the MLEs of σ 1, µ 1, σ and µ can be obtained by maximizing the profile log-likelihood function l profile (σ 1,µ,σ,µ ) = l( α 1 (σ 1,µ 1 ),σ 1,µ 1, α (σ.µ ),σ,µ, ρ(σ 1,µ 1,σ,µ )), (15) with respect to σ 1, µ 1, σ and µ. Note that the maximization of (15) cannot be performed explicitly. Newton-Raphson or some iterative procedure is needed to solve this problem. If the MLEs of σ 1, µ 1, σ and µ are denoted by σ 1, µ 1, σ and µ, respectively, then the MLEs of α 1, α and ρ become α 1 = α 1 ( σ 1, µ 1 ), α = α ( σ, µ ), ρ = ρ( σ 1, µ 1, σ, µ ). The exact distribution of the MLEs cannot be obtained, since BSHN is a regular family, it satisfies all the required conditions for the MLEs to be consistent and asymptotically normally distributed. We have the following results. Theorem 4.1: If θ is the MLE of θ, then n( θ θ) d N 7 (0,I 1 ), here d means convergence in distribution, N 7 (0,I 1 ) denotes the 7-variate normal distribution with mean vector 0 and the covariance matrix I 1. Here the matrix I is the Fisher information matrix.
15 15 It may be mentioned that the elements of the Fisher information matrix can be obtained in a routine matter, although they are not in explicit forms. All the elements are in the double integration form and they are not presented here. 4. Two-Stage Estimators: In the previous subsection we have seen that the MLEs of unknown parameters can be obtained by maximizing the profile log-likelihood function, and it needs solving a four dimensional optimization problem. To avoid that we propose to use the estimation of the unknown parameters using the copula structure based on two stage estimation method. Using the copula structure, the log-likelihood function can be written as l(θ) = lnc G (F Y1 (y i1 ;α 1,σ 1,µ 1 ),F Y (y i ;α,σ,µ );ρ) + lnf Y1 (y i1 ;α 1,σ 1,µ 1 )+ lnf Y (y i ;α,σ,µ ). Two step estimation procedure can be incorporated as follows. First compute the estimates of α 1, σ 1 and µ 1, say α 1, σ 1 and µ 1, respectively, by maximizing h 1 (α 1,σ 1,µ 1 ) = lnf Y1 (y i1 ;α 1,σ 1,µ 1 ) with respect to α 1, σ 1 and µ 1. Similarly, obtain the estimates of α, σ and µ, say α, σ and µ, respectively, by maximizing h (α,σ,µ ) = lnf Y (y i ;α,σ,µ ) with respect to α, σ and µ. Finally, obtain the estimate of ρ, say ρ, by maximizing h(ρ) with respect to ρ, where and h(ρ) = lnc G (u i1,u i ;ρ), u i1 = F Y1 (y i1 ; α 1, σ 1, µ 1 ) and u i = F Y (y i ; α, σ, µ ).
16 16 Now we indicate how to compute α 1, σ 1 and µ 1. Note that h 1 ( ) can be written as h 1 (α 1,σ 1,µ 1 ) = lnφ(a(y i1 ;α 1,σ 1,µ 1 ))+ ln(a(y i1 ;α 1,σ 1,µ 1 )). (16) Hence, for fixed σ 1 and µ 1, the maximization of (16) with respect to α 1 can be obtained at α 1 (σ 1,µ 1 ) = ( 1 n 1/ g (y i1 ;σ 1,µ 1 )), where g( ) is same as defined (14). Further σ 1 and µ 1 can be obtained by maximizing h 1 ( α 1 (σ 1,µ 1 ),σ 1,µ 1 ) = lnφ(a(y i1 ; α 1 (σ 1,µ 1 ),σ 1,µ 1 ))+ ln(a(y i1 ; α 1 (σ 1,µ 1 ),σ 1,µ 1 )), with respect to σ 1 and µ 1. Similarly, we can obtain α, σ and µ. Finally, obtain the ρ by maximizing h(ρ) with respect to ρ, which is equivalent to maximize lnφ (z i1,z i ;ρ), where for i = 1,...,n, z i1 = Φ 1 (u i1 ) = a(y i1 ; α 1, σ 1, µ 1 ) and z i = Φ 1 (u i ) = a(y i ; α, σ, µ ). Therefore, ρ = n z i1z i n n z i1 z i. A two-stage estimation process involves solving two independent two-dimensional optimization problems rather than one four dimensional optimization problem. Clearly it saves a significant amount of computational time. The estimators obtained by the two-stage procedure are consistent estimators of the unknown parameters. The asymptotic distribution of the two-stage estimators can be obtained in a routine matter as it has been obtained in Joe [6], and it is not pursued here.
17 17 5 Data Analysis This data set represents two different measurements of stiffness, Shock and Vibration of each of 30 boards. The first measurement (Shock) involves sending a shock wave down the board and the second measurement(vibration) is determined while vibrating the board. The data set was originally from William Galligan, and it is reported in Johnson and Wichern ([4], page 16). We divide all the data points by 1000, just for computational purposes, it is not going to make any difference in the statistical inference. We fit the SHN distributions to both the marginals and the estimates of the unknown parameters, namely α 1, σ 1, µ 1, α, σ and µ, are as follows: α 1 = , σ 1 = , µ 1 = , α = , σ = , µ 1 = We want to see how good the SHN distribution fits the marginals. The Kolmogorov-Smirnov distance between the empirical distribution function and the fitted distribution function, and the associated p value (reported in brackets) for Y 1 and Y are (0.8965) and (0.9535), respectively. Therefore, based on the Kolmogorov-Smirnov distances, it is clear that the SHN distribution fits the marginals quite well. Finally, we obtain ρ = The 95% confidence intervals of α 1, σ 1, µ 1, α, σ, µ and ρ become, (0.0017,0.0048), ( , ),(0.7781,.1539),( , ),( , ),(0.9874,.4504), (0.8449,1.0000), respectively. For comparison purposes, we have fitted the bivariate normal distribution to the data set, and the estimates are as follows: µ 1 = , σ 1 = 0.987, µ = 1.749, σ = , ρ = The Kolmogorov-Smirnov distance between the empirical distribution function and the fitted distribution function, and the associated p value (reported in brackets) for Y 1 and Y
18 18 are (0.7641) and (0.933), respectively. Therefore, for the marginals SHN distribution fits better than the normal distribution. We also perform the bivariate goodness of fit test. In case of bivariate normal distribution the chi-square value is , and the associated p value is between 0.95 and For BSHN distribution, we transform it to bivariate normal and perform the bivariate normal goodness of fit test. The chi-square value in this case is , and the associated p value is also between 0.95 and Since for BSHN, the chi-square value is lower, it indicates that BSHN provides marginally a better fit than the bivariate normal distribution. 6 Related Model In this section we consider bivariate generalized Birnbaum-Saunders distribution, which can be obtained from BSHN by transformation. Owen [11], see also Diaz-Garcia and Dominguez- Molina [], proposed a generalization of the Birnbaum-Saunders distribution, which can be described as follows. A random variable T, is said to have a generalized Birnbaum-Saunders distribution with parameters α > 0, β > 0 and λ > 0, if the CDF of T is F T (t) = P(T t) = Φ(b(t;α,β,λ)); t > 0. where The corresponding PDF becomes; { (t b(t;α,β,λ) = 1 ) λ α β ( ) } λ β. t f T (t) = φ(b(t;α,β,λ)b(t;α,β,λ); t > 0, where [ (t B(t;α,β,λ) = d dt b(t;α,β,λ) = λ ) λ + αt β ( ) ] λ β. t
19 19 We denote this random variable as GBS(α,β,λ). It can be easily seen that if T GBS(α,β,λ) lnt SHN(α,λ 1,lnβ). Now we are in a position to introduce the bivariate generalized Birnbaum-Saunders distribution as follows. A bivariate random vector (T 1,T ) is said to have bivariate generalized Birnbaum-Saunders distribution with parameters α 1, β 1, λ 1, α, β, λ and ρ, if (T 1,T ) = (e Y 1,e Y ), where (Y 1,Y ) BSHN(α 1,λ 1 1,lnβ 1,α,λ 1,lnβ,ρ). The joint CDF of (T 1,T ) becomes P(T 1 t 1,T t ) = Φ (b(t 1 ;α 1,β 1,λ 1 ),b(t ;α,β,λ );ρ); t 1 > 0,t > 0. (17) A random vector (T 1,T ) with CDF (17), will be denoted by BGBS(α 1,β 1,λ 1,α,β,λ,ρ). The joint PDF of (T 1,T ) becomes f T1,T (t 1,t ) = φ (b(t 1 ;α 1,β 1,λ 1 ),b(t ;α,β,λ );ρ) B(t i ;α i,β i,λ i ); t 1 > 0,t > 0. (18) Notethatwhenλ 1 = λ = 1/,itmatcheswiththebivariateBirnbaum-Saundersdistribution introduced by Kundu et al. [7]. Therefore, most of the results obtained in Kundu et al. [7] can be obtained as special cases of the results of this paper. We present the results without proofs, as they can be obtained along the same line as the proofs provided in Section 3. Theorem 6.1: If (T 1,T ) BGBS(α 1,β 1,λ 1,α,β,λ,ρ), then (a) T i GBS(α i,β i,λ i ); i = 1,. (b) The conditional CDF of T 1 given T = t, is given by [ ] b(t 1 ;α 1,β 1,λ 1 ) ρb(t ;α,β,λ ) P(T 1 t 1 T = t ) = Φ 1 ρ (c) The conditional PDF of T 1, given T = t is { } b(t 1 ;α 1,β 1,λ 1 ) ρb(t ;α,β,λ ) 1 f T1 T =t (t 1 ) = φ B(t 1;α 1,β 1,λ 1 ). 1 ρ 1 ρ
20 0 Theorem 6.: If (T 1,T ) BGBS(α 1,β 1,λ 1,α,β,λ,ρ), then (a) For c 1 > 0,c > 0, (c 1 T 1,c T ) BGBS(α 1,c 1 β 1,λ 1,α,c β,λ,ρ). (b) (T 1,T 1 ) BGBS(α 1,β 1,λ 1,α,β 1,λ, ρ). (c) (T 1 1,T ) BGBS(α 1,β 1 1,λ 1,α,β,λ, ρ). (d) (T 1 1,T 1 ) BGBS(α 1,β 1 1,λ 1,α,β 1,λ,ρ). Theorem 6.3: If (T 1,T ) BGBS(α 1,β 1,λ 1,α,β,λ,ρ), then for ρ > 0, it is TP and for ρ < 0, it is RR. Note that the results of Theorems 3.1, 3., and 3.3 of Kundu et al. [7] can be obtained as special cases of the results of Theorems 6.1, 6. and 6.3, respectively. Due to copula property, Blomqvist s beta, Kendall s tau and Spearman s rho will be same as the bivariate sinh-normal distribution. The maximum likelihood estimators of the unknown parameters cannot be obtained in closed form. In this case also directly the two stage estimators can be used as before, alternatively, after making the data transformation, the two stage procedure as described in Section 4, can be used quite easily. 7 Conclusion In this paper we consider a new bivariate distribution which is a natural extension of the univariate sinh-normal distribution. The proposed model has seven unknown parameters, andweobtaindifferentpropertiesofthemodel. Itisobservedthatthemodelcanbeobtained as a Gaussian copula, hence several properties of the Gaussian copula can be easily obtained for this model. The maximum likelihood estimators of the unknown parameters cannot be obtained in closed form, and they can be obtained by using profile log-likelihood method. It involves solving a four dimensional optimization problem. To avoid that we propose to use
21 1 two step estimation procedure using the copula structure. The estimators obtained by two stage procedure can be obtained by solving two independent two-dimensional optimization problems. They are consistent, and asymptotically normally distributed. Computationally they can be obtained in a much easier manner than the MLEs. Further, we have considered bivariate generalized Birnbaum-Saunders distribution, and studied its different properties. It is observed that many properties of the bivariate Birnbaum-Saunders distribution can be obtained as special cases of the proposed model. It may be mentioned that although in this paper we have mainly concentrated for the bivariate model, the same concept can be extended for any p dimensional model. Several properties can be obtained along the same line as the bivariate model. Clearly, it is more flexible than the multivariate normal model. Different other characteristics are under investigation, it will be reported later. Acknowledgements The author would like to thank the reviewers for their constructive suggestions. References [1] Birnbaum, Z.W., Saunders, S.C. (1969), A new family of life distributions, Journal of Applied Probability, vol. 6, [] Diaz-Garcia, J.A. and Dominguez-Molina, J.R. (006), Some generalizations of Birnbaum-Saunders and sinh-normal distributions, International Mathematical Forum 1, no , [3] Gupta, P.L. and Gupta, R.C. (1997), On the multivariate normal hazard, Journal of Multivariate Analysis, vol. 6,
22 [4] Johnson, R.A. and Wichern, D.W. (199), Applied multivariate statistical analysis, 3-rd edition, Prentice Hall, New Jersey. [5] Joe, H. (1997), Multivariate model and dependence concepts, Chapman and Hall, London. [6] Joe, H. (005), Asymptotic efficiency of the two-stage estimation method for copulabased models, Journal of Multivariate Analysis, vol. 94, [7] Kundu, D., Balakrishnan, N. and Jamalizadeh, A. (010), Bivariate Birnbaum- Saunders distribution and its associated inference, Journal of Multivariate Analysis, vol. 101, [8] Marshall, A. W. (1975), Some comments on the hazard gradient, Stochastic Processes and its Applications, vol. 3, [9] Meyer, C. (009), The bivariate normal copula, arxiv: [math.pr]. [10] Nelsen, R.B. (006), An introduction to copulas, Springer, New York. [11] Owen, W.J. (004), Another look at the Birnbaum-Saunders distribution, / [1] Rieck, J.R. (1989), Statistical analysis for the Birnbaum-Saunders fatigue life distribution, Ph.D. thesis, Clemson University, Department of Mathematical Sciences, Canada. [13] Rieck, J.R. and Nedelman, J.R. (1991), A log-linear model for the Birnbaum-Saunders distribution, Technometrics, vol. 33,
On Bivariate Birnbaum-Saunders Distribution
On Bivariate Birnbaum-Saunders Distribution Debasis Kundu 1 & Ramesh C. Gupta Abstract Univariate Birnbaum-Saunders distribution has been used quite effectively to analyze positively skewed lifetime data.
More informationMultivariate Birnbaum-Saunders Distribution Based on Multivariate Skew Normal Distribution
Multivariate Birnbaum-Saunders Distribution Based on Multivariate Skew Normal Distribution Ahad Jamalizadeh & Debasis Kundu Abstract Birnbaum-Saunders distribution has received some attention in the statistical
More informationOn Weighted Exponential Distribution and its Length Biased Version
On Weighted Exponential Distribution and its Length Biased Version Suchismita Das 1 and Debasis Kundu 2 Abstract In this paper we consider the weighted exponential distribution proposed by Gupta and Kundu
More informationGeneralized Exponential Geometric Extreme Distribution
Generalized Exponential Geometric Extreme Distribution Miroslav M. Ristić & Debasis Kundu Abstract Recently Louzada et al. ( The exponentiated exponential-geometric distribution: a distribution with decreasing,
More informationAn Extension of the Freund s Bivariate Distribution to Model Load Sharing Systems
An Extension of the Freund s Bivariate Distribution to Model Load Sharing Systems G Asha Department of Statistics Cochin University of Science and Technology Cochin, Kerala, e-mail:asha@cusat.ac.in Jagathnath
More informationWeighted Marshall-Olkin Bivariate Exponential Distribution
Weighted Marshall-Olkin Bivariate Exponential Distribution Ahad Jamalizadeh & Debasis Kundu Abstract Recently Gupta and Kundu [9] introduced a new class of weighted exponential distributions, and it can
More informationBivariate Geometric (Maximum) Generalized Exponential Distribution
Bivariate Geometric (Maximum) Generalized Exponential Distribution Debasis Kundu 1 Abstract In this paper we propose a new five parameter bivariate distribution obtained by taking geometric maximum of
More informationContents 1. Contents
Contents 1 Contents 6 Distributions of Functions of Random Variables 2 6.1 Transformation of Discrete r.v.s............. 3 6.2 Method of Distribution Functions............. 6 6.3 Method of Transformations................
More informationA New Two Sample Type-II Progressive Censoring Scheme
A New Two Sample Type-II Progressive Censoring Scheme arxiv:609.05805v [stat.me] 9 Sep 206 Shuvashree Mondal, Debasis Kundu Abstract Progressive censoring scheme has received considerable attention in
More informationOn the Comparison of Fisher Information of the Weibull and GE Distributions
On the Comparison of Fisher Information of the Weibull and GE Distributions Rameshwar D. Gupta Debasis Kundu Abstract In this paper we consider the Fisher information matrices of the generalized exponential
More informationOn Sarhan-Balakrishnan Bivariate Distribution
J. Stat. Appl. Pro. 1, No. 3, 163-17 (212) 163 Journal of Statistics Applications & Probability An International Journal c 212 NSP On Sarhan-Balakrishnan Bivariate Distribution D. Kundu 1, A. Sarhan 2
More informationBayes Estimation and Prediction of the Two-Parameter Gamma Distribution
Bayes Estimation and Prediction of the Two-Parameter Gamma Distribution Biswabrata Pradhan & Debasis Kundu Abstract In this article the Bayes estimates of two-parameter gamma distribution is considered.
More informationGeometric Skew-Normal Distribution
Debasis Kundu Arun Kumar Chair Professor Department of Mathematics & Statistics Indian Institute of Technology Kanpur Part of this work is going to appear in Sankhya, Ser. B. April 11, 2014 Outline 1 Motivation
More informationCopulas. MOU Lili. December, 2014
Copulas MOU Lili December, 2014 Outline Preliminary Introduction Formal Definition Copula Functions Estimating the Parameters Example Conclusion and Discussion Preliminary MOU Lili SEKE Team 3/30 Probability
More informationDependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.
Practitioner Course: Portfolio Optimization September 10, 2008 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y ) (x,
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationA MULTIVARIATE BIRNBAUM-SAUNDERS DISTRIBUTION BASED ON THE MULTIVARIATE SKEW NORMAL DISTRIBUTION
J. Japan Statist. Soc. Vol. 45 No. 1 2015 1 20 A MULTIVARIATE BIRNBAUM-SAUNDERS DISTRIBUTION BASED ON THE MULTIVARIATE SKEW NORMAL DISTRIBUTION Ahad Jamalizadeh* and Debasis Kundu** The Birnbaum-Saunders
More informationOn Bivariate Discrete Weibull Distribution
On Bivariate Discrete Weibull Distribution Debasis Kundu & Vahid Nekoukhou Abstract Recently, Lee and Cha (2015, On two generalized classes of discrete bivariate distributions, American Statistician, 221-230)
More informationarxiv: v4 [stat.co] 18 Mar 2018
An EM algorithm for absolute continuous bivariate Pareto distribution Arabin Kumar Dey, Biplab Paul and Debasis Kundu arxiv:1608.02199v4 [stat.co] 18 Mar 2018 Abstract: Recently [3] used EM algorithm to
More informationDiscriminating Between the Bivariate Generalized Exponential and Bivariate Weibull Distributions
Discriminating Between the Bivariate Generalized Exponential and Bivariate Weibull Distributions Arabin Kumar Dey & Debasis Kundu Abstract Recently Kundu and Gupta ( Bivariate generalized exponential distribution,
More informationLecture 25: Review. Statistics 104. April 23, Colin Rundel
Lecture 25: Review Statistics 104 Colin Rundel April 23, 2012 Joint CDF F (x, y) = P [X x, Y y] = P [(X, Y ) lies south-west of the point (x, y)] Y (x,y) X Statistics 104 (Colin Rundel) Lecture 25 April
More informationAn Extension of the Generalized Exponential Distribution
An Extension of the Generalized Exponential Distribution Debasis Kundu and Rameshwar D. Gupta Abstract The two-parameter generalized exponential distribution has been used recently quite extensively to
More informationHybrid Censoring; An Introduction 2
Hybrid Censoring; An Introduction 2 Debasis Kundu Department of Mathematics & Statistics Indian Institute of Technology Kanpur 23-rd November, 2010 2 This is a joint work with N. Balakrishnan Debasis Kundu
More informationOn the Hazard Function of Birnbaum Saunders Distribution and Associated Inference
On the Hazard Function of Birnbaum Saunders Distribution and Associated Inference Debasis Kundu, Nandini Kannan & N. Balakrishnan Abstract In this paper, we discuss the shape of the hazard function of
More informationFinancial Econometrics and Volatility Models Copulas
Financial Econometrics and Volatility Models Copulas Eric Zivot Updated: May 10, 2010 Reading MFTS, chapter 19 FMUND, chapters 6 and 7 Introduction Capturing co-movement between financial asset returns
More informationThe Delta Method and Applications
Chapter 5 The Delta Method and Applications 5.1 Local linear approximations Suppose that a particular random sequence converges in distribution to a particular constant. The idea of using a first-order
More informationMultivariate Non-Normally Distributed Random Variables
Multivariate Non-Normally Distributed Random Variables An Introduction to the Copula Approach Workgroup seminar on climate dynamics Meteorological Institute at the University of Bonn 18 January 2008, Bonn
More informationApproximation of multivariate distribution functions MARGUS PIHLAK. June Tartu University. Institute of Mathematical Statistics
Approximation of multivariate distribution functions MARGUS PIHLAK June 29. 2007 Tartu University Institute of Mathematical Statistics Formulation of the problem Let Y be a random variable with unknown
More informationA measure of radial asymmetry for bivariate copulas based on Sobolev norm
A measure of radial asymmetry for bivariate copulas based on Sobolev norm Ahmad Alikhani-Vafa Ali Dolati Abstract The modified Sobolev norm is used to construct an index for measuring the degree of radial
More informationLecture 1: August 28
36-705: Intermediate Statistics Fall 2017 Lecturer: Siva Balakrishnan Lecture 1: August 28 Our broad goal for the first few lectures is to try to understand the behaviour of sums of independent random
More informationMultiple Random Variables
Multiple Random Variables This Version: July 30, 2015 Multiple Random Variables 2 Now we consider models with more than one r.v. These are called multivariate models For instance: height and weight An
More informationMultivariate Statistics
Multivariate Statistics Chapter 2: Multivariate distributions and inference Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2016/2017 Master in Mathematical
More informationModified slash Birnbaum-Saunders distribution
Modified slash Birnbaum-Saunders distribution Jimmy Reyes, Filidor Vilca, Diego I. Gallardo and Héctor W. Gómez Abstract In this paper, we introduce an extension for the Birnbaum-Saunders (BS) distribution
More informationA Goodness-of-fit Test for Copulas
A Goodness-of-fit Test for Copulas Artem Prokhorov August 2008 Abstract A new goodness-of-fit test for copulas is proposed. It is based on restrictions on certain elements of the information matrix and
More informationThe Instability of Correlations: Measurement and the Implications for Market Risk
The Instability of Correlations: Measurement and the Implications for Market Risk Prof. Massimo Guidolin 20254 Advanced Quantitative Methods for Asset Pricing and Structuring Winter/Spring 2018 Threshold
More informationModelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich
Modelling Dependence with Copulas and Applications to Risk Management Filip Lindskog, RiskLab, ETH Zürich 02-07-2000 Home page: http://www.math.ethz.ch/ lindskog E-mail: lindskog@math.ethz.ch RiskLab:
More informationPh.D. Qualifying Exam Friday Saturday, January 6 7, 2017
Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017 Put your solution to each problem on a separate sheet of paper. Problem 1. (5106) Let X 1, X 2,, X n be a sequence of i.i.d. observations from a
More information2 Functions of random variables
2 Functions of random variables A basic statistical model for sample data is a collection of random variables X 1,..., X n. The data are summarised in terms of certain sample statistics, calculated as
More informationMultivariate Distribution Models
Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is
More informationEstimation of the Bivariate Generalized. Lomax Distribution Parameters. Based on Censored Samples
Int. J. Contemp. Math. Sciences, Vol. 9, 2014, no. 6, 257-267 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2014.4329 Estimation of the Bivariate Generalized Lomax Distribution Parameters
More information1 Introduction. On grade transformation and its implications for copulas
Brazilian Journal of Probability and Statistics (2005), 19, pp. 125 137. c Associação Brasileira de Estatística On grade transformation and its implications for copulas Magdalena Niewiadomska-Bugaj 1 and
More informationModelling Dependent Credit Risks
Modelling Dependent Credit Risks Filip Lindskog, RiskLab, ETH Zürich 30 November 2000 Home page:http://www.math.ethz.ch/ lindskog E-mail:lindskog@math.ethz.ch RiskLab:http://www.risklab.ch Modelling Dependent
More informationDependence. MFM Practitioner Module: Risk & Asset Allocation. John Dodson. September 11, Dependence. John Dodson. Outline.
MFM Practitioner Module: Risk & Asset Allocation September 11, 2013 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y
More informationp-birnbaum SAUNDERS DISTRIBUTION: APPLICATIONS TO RELIABILITY AND ELECTRONIC BANKING HABITS
p-birnbaum SAUNDERS DISTRIBUTION: APPLICATIONS TO RELIABILITY AND ELECTRONIC BANKING 1 V.M.Chacko, Mariya Jeeja P V and 3 Deepa Paul 1, Department of Statistics St.Thomas College, Thrissur Kerala-681 e-mail:chackovm@gmail.com
More informationA Convenient Way of Generating Gamma Random Variables Using Generalized Exponential Distribution
A Convenient Way of Generating Gamma Random Variables Using Generalized Exponential Distribution Debasis Kundu & Rameshwar D. Gupta 2 Abstract In this paper we propose a very convenient way to generate
More informationCopula modeling for discrete data
Copula modeling for discrete data Christian Genest & Johanna G. Nešlehová in collaboration with Bruno Rémillard McGill University and HEC Montréal ROBUST, September 11, 2016 Main question Suppose (X 1,
More informationCopula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011
Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011 Outline Ordinary Least Squares (OLS) Regression Generalized Linear Models
More informationBivariate Degradation Modeling Based on Gamma Process
Bivariate Degradation Modeling Based on Gamma Process Jinglun Zhou Zhengqiang Pan Member IAENG and Quan Sun Abstract Many highly reliable products have two or more performance characteristics (PCs). The
More information4. Distributions of Functions of Random Variables
4. Distributions of Functions of Random Variables Setup: Consider as given the joint distribution of X 1,..., X n (i.e. consider as given f X1,...,X n and F X1,...,X n ) Consider k functions g 1 : R n
More informationMULTIVARIATE PROBABILITY DISTRIBUTIONS
MULTIVARIATE PROBABILITY DISTRIBUTIONS. PRELIMINARIES.. Example. Consider an experiment that consists of tossing a die and a coin at the same time. We can consider a number of random variables defined
More informationAn Extended Weighted Exponential Distribution
Journal of Modern Applied Statistical Methods Volume 16 Issue 1 Article 17 5-1-017 An Extended Weighted Exponential Distribution Abbas Mahdavi Department of Statistics, Faculty of Mathematical Sciences,
More informationarxiv: v3 [stat.co] 22 Jun 2017
arxiv:608.0299v3 [stat.co] 22 Jun 207 An EM algorithm for absolutely continuous Marshall-Olkin bivariate Pareto distribution with location and scale Arabin Kumar Dey and Debasis Kundu Department of Mathematics,
More informationTail dependence in bivariate skew-normal and skew-t distributions
Tail dependence in bivariate skew-normal and skew-t distributions Paola Bortot Department of Statistical Sciences - University of Bologna paola.bortot@unibo.it Abstract: Quantifying dependence between
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationUnivariate and Bivariate Geometric Discrete Generalized Exponential Distributions
Univariate and Bivariate Geometric Discrete Generalized Exponential Distributions Debasis Kundu 1 & Vahid Nekoukhou 2 Abstract Marshall and Olkin (1997, Biometrika, 84, 641-652) introduced a very powerful
More informationIntroduction of Shape/Skewness Parameter(s) in a Probability Distribution
Journal of Probability and Statistical Science 7(2), 153-171, Aug. 2009 Introduction of Shape/Skewness Parameter(s) in a Probability Distribution Rameshwar D. Gupta University of New Brunswick Debasis
More informationA New Generalized Gumbel Copula for Multivariate Distributions
A New Generalized Gumbel Copula for Multivariate Distributions Chandra R. Bhat* The University of Texas at Austin Department of Civil, Architectural & Environmental Engineering University Station, C76,
More informationAnalysis of Type-II Progressively Hybrid Censored Data
Analysis of Type-II Progressively Hybrid Censored Data Debasis Kundu & Avijit Joarder Abstract The mixture of Type-I and Type-II censoring schemes, called the hybrid censoring scheme is quite common in
More informationBivariate Weibull-power series class of distributions
Bivariate Weibull-power series class of distributions Saralees Nadarajah and Rasool Roozegar EM algorithm, Maximum likelihood estimation, Power series distri- Keywords: bution. Abstract We point out that
More informationINFERENCE FOR BIRNBAUM-SAUNDERS, LAPLACE AND SOME RELATED DISTRIBUTIONS UNDER CENSORED DATA
INFERENCE FOR BIRNBAUM-SAUNDERS, LAPLACE AND SOME RELATED DISTRIBUTIONS UNDER CENSORED DATA INFERENCE FOR BIRNBAUM-SAUNDERS, LAPLACE AND SOME RELATED DISTRIBUTIONS UNDER CENSORED DATA By Xiaojun Zhu A
More informationSome generalisations of Birnbaum-Saunders and sinh-normal distributions
International Mathematical Forum, 1, 006, no. 35, 1709-177 Some generalisations of Birnbaum-Saunders and sinh-normal distributions José A.Díaz-García Universidad Autónoma Agraria Antonio Narro Department
More informationContinuous Random Variables
1 / 24 Continuous Random Variables Saravanan Vijayakumaran sarva@ee.iitb.ac.in Department of Electrical Engineering Indian Institute of Technology Bombay February 27, 2013 2 / 24 Continuous Random Variables
More informationRobust Parameter Estimation in the Weibull and the Birnbaum-Saunders Distribution
Clemson University TigerPrints All Theses Theses 8-2012 Robust Parameter Estimation in the Weibull and the Birnbaum-Saunders Distribution Jing Zhao Clemson University, jzhao2@clemson.edu Follow this and
More informationExact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring
Exact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring A. Ganguly, S. Mitra, D. Samanta, D. Kundu,2 Abstract Epstein [9] introduced the Type-I hybrid censoring scheme
More informationON THE NORMAL MOMENT DISTRIBUTIONS
ON THE NORMAL MOMENT DISTRIBUTIONS Akin Olosunde Department of Statistics, P.M.B. 40, University of Agriculture, Abeokuta. 00, NIGERIA. Abstract Normal moment distribution is a particular case of well
More informationSTAT 512 sp 2018 Summary Sheet
STAT 5 sp 08 Summary Sheet Karl B. Gregory Spring 08. Transformations of a random variable Let X be a rv with support X and let g be a function mapping X to Y with inverse mapping g (A = {x X : g(x A}
More informationIntroduction to Normal Distribution
Introduction to Normal Distribution Nathaniel E. Helwig Assistant Professor of Psychology and Statistics University of Minnesota (Twin Cities) Updated 17-Jan-2017 Nathaniel E. Helwig (U of Minnesota) Introduction
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2008 Prof. Gesine Reinert 1 Data x = x 1, x 2,..., x n, realisations of random variables X 1, X 2,..., X n with distribution (model)
More informationA Brief Introduction to Copulas
A Brief Introduction to Copulas Speaker: Hua, Lei February 24, 2009 Department of Statistics University of British Columbia Outline Introduction Definition Properties Archimedean Copulas Constructing Copulas
More informationOn Parameter-Mixing of Dependence Parameters
On Parameter-Mixing of Dependence Parameters by Murray D Smith and Xiangyuan Tommy Chen 2 Econometrics and Business Statistics The University of Sydney Incomplete Preliminary Draft May 9, 2006 (NOT FOR
More informationSongklanakarin Journal of Science and Technology SJST R1 Sukparungsee
Songklanakarin Journal of Science and Technology SJST-0-0.R Sukparungsee Bivariate copulas on the exponentially weighted moving average control chart Journal: Songklanakarin Journal of Science and Technology
More informationEstimation theory. Parametric estimation. Properties of estimators. Minimum variance estimator. Cramer-Rao bound. Maximum likelihood estimators
Estimation theory Parametric estimation Properties of estimators Minimum variance estimator Cramer-Rao bound Maximum likelihood estimators Confidence intervals Bayesian estimation 1 Random Variables Let
More informationBivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories
Entropy 2012, 14, 1784-1812; doi:10.3390/e14091784 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories Lan Zhang
More informationMA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems
MA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems Review of Basic Probability The fundamentals, random variables, probability distributions Probability mass/density functions
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Kundu, Debasis] On: 4 November 2009 Access details: Access Details: [subscription number 9655482] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationNew properties of the orthant convex-type stochastic orders
Noname manuscript No. (will be inserted by the editor) New properties of the orthant convex-type stochastic orders Fernández-Ponce, JM and Rodríguez-Griñolo, MR the date of receipt and acceptance should
More informationAnalysis of Middle Censored Data with Exponential Lifetime Distributions
Analysis of Middle Censored Data with Exponential Lifetime Distributions Srikanth K. Iyer S. Rao Jammalamadaka Debasis Kundu Abstract Recently Jammalamadaka and Mangalam (2003) introduced a general censoring
More informationIGD-TP Exchange Forum n 5 WG1 Safety Case: Handling of uncertainties October th 2014, Kalmar, Sweden
IGD-TP Exchange Forum n 5 WG1 Safety Case: Handling of uncertainties October 28-30 th 2014, Kalmar, Sweden Comparison of probabilistic and alternative evidence theoretical methods for the handling of parameter
More informationBivariate Paired Numerical Data
Bivariate Paired Numerical Data Pearson s correlation, Spearman s ρ and Kendall s τ, tests of independence University of California, San Diego Instructor: Ery Arias-Castro http://math.ucsd.edu/~eariasca/teaching.html
More informationLatin Hypercube Sampling with Multidimensional Uniformity
Latin Hypercube Sampling with Multidimensional Uniformity Jared L. Deutsch and Clayton V. Deutsch Complex geostatistical models can only be realized a limited number of times due to large computational
More informationProbability Distributions and Estimation of Ali-Mikhail-Haq Copula
Applied Mathematical Sciences, Vol. 4, 2010, no. 14, 657-666 Probability Distributions and Estimation of Ali-Mikhail-Haq Copula Pranesh Kumar Mathematics Department University of Northern British Columbia
More information1. Density and properties Brief outline 2. Sampling from multivariate normal and MLE 3. Sampling distribution and large sample behavior of X and S 4.
Multivariate normal distribution Reading: AMSA: pages 149-200 Multivariate Analysis, Spring 2016 Institute of Statistics, National Chiao Tung University March 1, 2016 1. Density and properties Brief outline
More informationMultivariate Random Variable
Multivariate Random Variable Author: Author: Andrés Hincapié and Linyi Cao This Version: August 7, 2016 Multivariate Random Variable 3 Now we consider models with more than one r.v. These are called multivariate
More information0, otherwise. U = Y 1 Y 2 Hint: Use either the method of distribution functions or a bivariate transformation. (b) Find E(U).
1. Suppose Y U(0, 2) so that the probability density function (pdf) of Y is 1 2, 0 < y < 2 (a) Find the pdf of U = Y 4 + 1. Make sure to note the support. (c) Suppose Y 1, Y 2,..., Y n is an iid sample
More informationRegression and Statistical Inference
Regression and Statistical Inference Walid Mnif wmnif@uwo.ca Department of Applied Mathematics The University of Western Ontario, London, Canada 1 Elements of Probability 2 Elements of Probability CDF&PDF
More informationDiscussion of Maximization by Parts in Likelihood Inference
Discussion of Maximization by Parts in Likelihood Inference David Ruppert School of Operations Research & Industrial Engineering, 225 Rhodes Hall, Cornell University, Ithaca, NY 4853 email: dr24@cornell.edu
More informationBurr Type X Distribution: Revisited
Burr Type X Distribution: Revisited Mohammad Z. Raqab 1 Debasis Kundu Abstract In this paper, we consider the two-parameter Burr-Type X distribution. We observe several interesting properties of this distribution.
More informationA Probability Review
A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in
More informationEstimation of Copula Models with Discrete Margins (via Bayesian Data Augmentation) Michael S. Smith
Estimation of Copula Models with Discrete Margins (via Bayesian Data Augmentation) Michael S. Smith Melbourne Business School, University of Melbourne (Joint with Mohamad Khaled, University of Queensland)
More informationConstruction and estimation of high dimensional copulas
Construction and estimation of high dimensional copulas Gildas Mazo PhD work supervised by S. Girard and F. Forbes Mistis, Inria and laboratoire Jean Kuntzmann, Grenoble, France Séminaire Statistiques,
More informationBirnbaum-Saunders Distribution: A Review of Models, Analysis and Applications
Birnbaum-Saunders Distribution: A Review of Models, Analysis and Applications N. Balakrishnan and Debasis Kundu Abstract Birnbaum and Saunders [43, 44] introduced a two-parameter lifetime distribution
More informationHybrid Censoring Scheme: An Introduction
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline 1 2 3 4 5 Outline 1 2 3 4 5 What is? Lifetime data analysis is used to analyze data in which the time
More informationCopulas. Mathematisches Seminar (Prof. Dr. D. Filipovic) Di Uhr in E
Copulas Mathematisches Seminar (Prof. Dr. D. Filipovic) Di. 14-16 Uhr in E41 A Short Introduction 1 0.8 0.6 0.4 0.2 0 0 0.2 0.4 0.6 0.8 1 The above picture shows a scatterplot (500 points) from a pair
More informationA New Method for Generating Distributions with an Application to Exponential Distribution
A New Method for Generating Distributions with an Application to Exponential Distribution Abbas Mahdavi & Debasis Kundu Abstract Anewmethodhasbeenproposedtointroduceanextraparametertoafamilyofdistributions
More informationOn a simple construction of bivariate probability functions with fixed marginals 1
On a simple construction of bivariate probability functions with fixed marginals 1 Djilali AIT AOUDIA a, Éric MARCHANDb,2 a Université du Québec à Montréal, Département de mathématiques, 201, Ave Président-Kennedy
More informationCopulas and Measures of Dependence
1 Copulas and Measures of Dependence Uttara Naik-Nimbalkar December 28, 2014 Measures for determining the relationship between two variables: the Pearson s correlation coefficient, Kendalls tau and Spearmans
More informationinferences on stress-strength reliability from lindley distributions
inferences on stress-strength reliability from lindley distributions D.K. Al-Mutairi, M.E. Ghitany & Debasis Kundu Abstract This paper deals with the estimation of the stress-strength parameter R = P (Y
More informationOutline of GLMs. Definitions
Outline of GLMs Definitions This is a short outline of GLM details, adapted from the book Nonparametric Regression and Generalized Linear Models, by Green and Silverman. The responses Y i have density
More informationLecture Quantitative Finance Spring Term 2015
on bivariate Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 07: April 2, 2015 1 / 54 Outline on bivariate 1 2 bivariate 3 Distribution 4 5 6 7 8 Comments and conclusions
More informationParameter Estimation of Power Lomax Distribution Based on Type-II Progressively Hybrid Censoring Scheme
Applied Mathematical Sciences, Vol. 12, 2018, no. 18, 879-891 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.8691 Parameter Estimation of Power Lomax Distribution Based on Type-II Progressively
More information