Asymptotic generalized bivariate extreme with random index
M.A. Abd Elgawad, A.M.Elsawah, Hong Qin, Ting Yan

TL;DR
This paper characterizes the limit distributions of bivariate extreme generalized order statistics with random sample sizes, providing conditions for their convergence and illustrating the results with examples.
Contribution
It fully characterizes the limit distributions of bivariate extreme GOS with random sample sizes and establishes conditions for their weak convergence.
Findings
Characterization of limit distribution functions of bivariate extreme GOS with random sizes
Necessary and sufficient conditions for weak convergence when sample size is independent
Sufficient conditions and forms of limit distributions when size and variables are dependent
Abstract
In many biological, agricultural, military activity problems and in some quality control problems, it is almost impossible to have a fixed sample size, because some observations are always lost for various reasons. Therefore, the sample size itself is considered frequently to be a random variable (rv). The class of limit distribution functions (df's) of the random bivariate extreme generalized order statistics (GOS) from independent and identically distributed RV's are fully characterized. When the random sample size is assumed to be independent of the basic variables and its df is assumed to converge weakly to a non-degenerate limit, the necessary and sufficient conditions for the weak convergence of the random bivariate extreme GOS are obtained. Furthermore, when the interrelation of the random size and the basic rv's is not restricted, sufficient conditions for the convergence and…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsStatistical Distribution Estimation and Applications · Probabilistic and Robust Engineering Design · Hydrology and Drought Analysis
Asymptotic generalized bivariate extreme with random index
M. A. Abd Elgawada,b, A. M. Elsawaha,c,d,111Corresponding author. E-mail: [email protected], [email protected], [email protected] , Hong Qina and Ting Yana
a *Faculty of Mathematics and Statistics, Central China Normal University, Wuhan 430079, China
b* *Department of Mathematics, Faculty of Science, Benha University, Benha 13518, Egypt
c* *Department of Mathematics, Faculty of Science, Zagazig University, Zagazig 44519, Egypt
d* Division of Science and Technology, BNU-HKBU United International College, Zhuhai 519085, China
Abstract
In many biological, agricultural, military activity problems and in some quality control problems, it is almost impossible to have a fixed sample size, because some observations are always lost for various reasons. Therefore, the sample size itself is considered frequently to be a random variable (rv). The class of limit distribution functions (df’s) of the random bivariate extreme generalized order statistics (GOS) from independent and identically distributed rv’s are fully characterized. When the random sample size is assumed to be independent of the basic variables and its df is assumed to converge weakly to a non-degenerate limit, the necessary and sufficient conditions for the weak convergence of the random bivariate extreme GOS are obtained. Furthermore, when the interrelation of the random size and the basic rv’s is not restricted, sufficient conditions of the convergence and the forms of the limit df’s are deduced. Illustrative examples are given which lend further support to our theoretical results.
Keywords: Weak convergence; Random sample size; Generalized order statistics; Generalized bivariate extreme.
1 Introduction
The concept of generalized order statistics (GOS) have been introduced by Kamps (1995). It’s enable a unified approach to ascendingly ordered random variables (rv’s) as ordinary order statistics (oos), sequential order statistics (sos), order statistics with non integral sample size, progressively type II censored order statistics (pos), record values, th record values and Pfeifer’s records. Let and Then the rv’s are called GOS based on the distribution function (df) with density function which are defined by their probability density function (pdf)
[TABLE]
where
In this work, we consider a wide subclass of GOS , by assuming This subclass is known as GOS. Clearly many important practical models of GOS are included such as oos, order statistics with non integer sample size and sos. The marginal df’s of the th and th GOS (Nasri- Roudsari, 1996 and Barakat, 2007) are represented by and respectively, where denotes the incomplete beta ratio function, and Moreover, by using the results of Kamps (1995), we can write explicitly the joint df’s of the th and th GOS, as:
[TABLE]
where and is the usual gamma function. Recently Barakat et al. (2014a) studied the limit df’s of joint extreme GOS, for a fixed sample size. Moreover, the asymptotic behavior for bivariate df of the lower-lower (l-l), upper-upper (u-u) and lower-upper (l-u) extreme GOS in Barakat et al. (2014b).
In the last few years much efforts had been devoted to investigate the limit df’s of independent rv’s with random sample size. The appearance of this trend is naturally because many applications require the consideration of such problem. For example, in many biological, agricultural and in some quality control problems, it is almost impossible to have a fixed sample size because some observations always get lost for various reasons. Therefore, the sample size itself is considered frequently to be a rv where is independent of the basic variables (i.e., the original random sample) or in some applications the interrelation of the basic variables and the random sample size is not restricted. Limit theorems for extremes with random sample size indexes have been thoroughly studied in the above mentioned two particular cases :
The basic variables and sample size index are independents (see, Barakat, 1997). 2. 2.
The interrelation of the basic variables and the random sample size is not restricted (see, Barakat and El Shandidy, 1990, Barakat, 1997 and Barakat et al., 2015a).
Our aim in this paper is to characterize the asymptotic behavior of the bivariate df’s of the (l-l), (u-u) and (l-u) extreme GOS with random sample size. When the random sample size is assumed to be independent of the basic variables and its df is assumed to converge weakly to a non-degenerate limit, the necessary and sufficient conditions for the weak convergence of the random bivariate extreme GOS are obtained. Furthermore, when the interrelation of the random size and the basic rv’s is not restricted, sufficient conditions of the convergence and the forms of the limit df’s are deduced. An illustrative examples are given which lend further support to our theoretical results. Throughout this paper the convergence in probability and the weak convergence, as respectively, denoted as "{\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle p}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}" and "{\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle w}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}".
2 Asymptotic random bivariate extreme under GOS
2.1 Random sample size and basic rv’s are independents
In this subsection we deal with the weak convergence of bivariate df’s of the (u-u), (l-l) and (l-u) extreme GOS are fully characterized in Theorems 2.1, 2.2 and 2.3, respectively. When the sample size itself is a rv which is assumed to be independent of the basic variables
Theorem 2.1. Consider the following three conditions :
[TABLE]
[TABLE]
[TABLE]
[TABLE]
Then any two of the above conditions imply the remaining one, where are suitable normalizing constants, is a non-degenerate df, is a df with
[TABLE]
denotes the incomplete gamma ratio function, and
Remark 2.1. The continuity of the limit df in implies the continuity of the limit Hence the convergence in is uniform with respect to and
Remark 2.2. It is natural to look for the limitations on under which we get the relation In view of Theorem 2.1, the last equation is satisfied if and only if the df is degenerate at one, which means the asymptotically almost randomlessness of We assume, due to Remark 2.2, that is a non-degenerate df and i.e., continuous at zero.
Proof of the implication : First, we note that can be written in the form (see, Theorem 2.3 in Barakat et al., 2014b)
[TABLE]
where Now by using the total probability rule we get,
[TABLE]
Assume that and where denotes the greatest integer part of Thus, the relation (2.1) show that the sum term in (2.2) is a Riemann sum of the integral
[TABLE]
where, for sufficiently large we have
[TABLE]
where Appealing to the condition Theorem 2.3 in Barakat et al. (2014b) and Remark 2.2, we get
[TABLE]
where the convergence is uniform with respect to and over any finite interval of
Now, let be a continuity point of such that ( is an aribtary small value). Then, we have
[TABLE]
In view of the condition we get, for sufficiently large that
[TABLE]
On the other hand, by the triangle inequality, we get
[TABLE]
[TABLE]
[TABLE]
where the convergence in (2.4) is uniform over the finite interval Therefore, for arbitrary and for sufficiently large we have
[TABLE]
In order to estimate the second difference on the right hand side of (2.7), we construct Riemann sums which are close to the integral there. Let be a fixed number and be continuity points of Furthermore, let and be such that
[TABLE]
and
[TABLE]
Since, by the assumption ~{}H_{n}(n\xi_{i}){\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle w}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}H(\xi_{i}),~{}0\leq i\leq T, the two Riemann sums are closer to each other than for all sufficiently large. Thus, once again by the triangle inequality, the absolute value of the difference of the integrals is smaller than Combining this fact with (2.8), the left hand side of (2.7) becomes smaller than for all large Therefore, in view of (2.5), (2.6) and (2.4), we have
[TABLE]
[TABLE]
This completes the proof of the first part of the theorem.
Proof of the implication : Starting with the relation (2.3), we select a subsequence of for which converges weakly to an extended df (i.e., and such a subsequence exists by the compactness of df’s). Then, by repeating the first part of the theorem for the subsequence with the exception that we choose so that we get Since the function is a df, we get which implies that is a df. Now, if did not converge weakly, then we can select two subsequences and such that H_{n^{\prime}}(n^{\prime}z){\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle w}}{{\longrightarrow}}\\ \scriptstyle{n^{\prime}}\end{array}}H^{\prime}(z) and H_{n^{\prime\prime}}(n^{\prime\prime}z){\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle w}}{{\longrightarrow}}\\ \scriptstyle{n^{\prime\prime}}\end{array}}H^{\prime\prime}(z), where and are df’s. In this case, we get
[TABLE]
Thus, let (), we get
[TABLE]
Appealing to equation and by using the same argument which is applied in the proof of the second part of Theorem 2.1 in Barakat, 1997, we can easily prove This complete the proof of the second part.
Proof of the implication : For proving this part, we need first present the following lemma.
Lemma 2.1. For all we have
[TABLE]
and
[TABLE]
where 0<\rho_{i,N}~{},\sigma_{i,N}~{}{\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle}}{{\longrightarrow}}\\ \scriptstyle{N}\end{array}}~{}0 (or equivalently, as and
Proof. Since the proof of the lemma will immediately follow from the result of Smirnov (1952) (Theorem 3, p. 133, or Lemma 2.1 in Barakat, 1997).
We now turn to the proof of the last part of Theorem 2.1. In view of Remark 2.1, we can assume, without any loss of generality, that the df is continuous. Therefore, the condition will be satisfied for all univariate marginals of i.e., we have
[TABLE]
where is the marginals df’s of We shall now prove
[TABLE]
In view of Lemma 2.1, we first show that the sequence is stochastically bounded (see, Feller, 1979). If we assume the contrary, we would find such that at least one of the two following relations
(a) \qquad{\begin{array}[]{c}{}\hfil\\ {}\hfil\\ {}\hfil\\ {\overline{\lim}}\\ \scriptstyle{n\to\infty}\end{array}}~{}P(Z^{\langle n\rangle}_{\grave{r_{i}}:n}\geq x_{i})\geq\varepsilon_{i,1}>0,~{}~{}\forall~{}~{}x_{i}>0,~{}i=1,2,~{}
(b) \qquad{\begin{array}[]{c}{}\hfil\\ {}\hfil\\ {}\hfil\\ {\overline{\lim}}\\ \scriptstyle{n\to\infty}\end{array}}P(Z^{\langle n\rangle}_{\grave{r_{i}}:n}<x_{i})\geq\varepsilon_{i,2}>0,~{}~{}\forall~{}~{}x_{i}<0,~{}i=1,2~{}
is satisfied. The assertions (a) and (b) mean that the sequence is not stochastically bounded at the left and at the right respectively. Let the assumption (a) be true. Since is non-degenerate df, we find and such that
[TABLE]
Using the following well known inequality, for
[TABLE]
We thus get the following inequalities, for sufficiently large
[TABLE]
[TABLE]
(note that ). Therefore,
[TABLE]
Now, if we find such that {\begin{array}[]{c}{}\hfil\\ {}\hfil\\ {}\hfil\\ {\overline{\lim}}\\ \scriptstyle{n\to\infty}\end{array}}P(Z^{\langle n\rangle}_{\grave{r_{i}}:[n\beta]}\geq x_{i})\geq\varepsilon^{\prime}_{i,1}>0, we get {\begin{array}[]{c}{}\hfil\\ {}\hfil\\ {}\hfil\\ {\overline{\lim}}\\ \scriptstyle{n\to\infty}\end{array}}P(Z^{\langle n\rangle}_{\grave{r_{i}}:\nu_{n}}\geq x_{i})\geq\epsilon_{0}\varepsilon^{\prime}_{i,1}>0, which contradicts the right stochastic boundedness of the sequence and consequently contradicts the relation (2.12). However, if such an does not exist we have {\begin{array}[]{c}{}\hfil\\ {}\hfil\\ {}\hfil\\ {\overline{\lim}}\\ \scriptstyle{n\to\infty}\end{array}}P(Z^{\langle n\rangle}_{\grave{r_{i}}:[n\beta]}\geq x_{i})=0, which in view of Lemma 2.1 (relation (2.10)) leads to the following chain of implications () (since which contradicts the assumption (a). Consider the assumption (b). Since is a df we can find a positive integer and real number such that
[TABLE]
Therefore, in view of (2.16) and the inequality (2.15), we have
[TABLE]
[TABLE]
Hence, we get {\begin{array}[]{c}{}\hfil\\ {}\hfil\\ {}\hfil\\ {\overline{\lim}}\\ \scriptstyle{n\to\infty}\end{array}}P(Z^{\langle n\rangle}_{\grave{r_{i}}:\nu_{n}}<x_{i})\geq\alpha{\begin{array}[]{c}{}\hfil\\ {}\hfil\\ {}\hfil\\ {\overline{\lim}}\\ \scriptstyle{n\to\infty}\end{array}}P(Z^{\langle n\rangle}_{\grave{r_{i}}:\delta n}<x_{i}).~{} By using Lemma 2.1 (relation (2.11)) and applying the same argument as in the case (a), it is easy to show that the last inequality leads to a contradiction (the last inequality, in view of the assumption (b)), which yields that the sequences is not stochastically bounded at the left. This completes the proof that the sequences are stochastically bounded. Now, if did not converge weakly, then we could select two subsequences and such that would converge weakly to and to another limit df In this case we get (by repeating the first part of Theorem 2.1 for the univariate case and for the two subsequences )
[TABLE]
However, Lemma 3.2 in Barakat (1997) shows that the last equalities, cannot hold unless Hence the relation (2.13) is proved. Hence, the proof of Theorem 2.1 is completed.
Let and be the classes of all possible limit df’s in and respectively. The class is fully determined by Barakat et al. (2014b). Furthermore, let and be the necessary and sufficient conditions for the validity of the relations and respectively. The following corollary characterizes the class
Corollary 2.1. For every df in there exists a unique df in such that is uniquely determined by the representation Moreover,
Proof of corollary 2.1. Let us first prove the implication If we assume the contrary, we get while Appealing to the first part of Theorem 2.1 , we get The last equalities, as we have seen before, from Lemma 3.2 in Barakat (1997), cannot hold unless Therefore, Corollary 2.1 is followed as a consequence of Theorem 2.1 and the last implication. This completes the proof of Corollary 2.1.
Theorem 2.2. Consider the following three conditions :
[TABLE]
[TABLE]
[TABLE]
[TABLE]
Then any two of the above conditions imply the remaining one, where are suitable normalizing constants, is a non-degenerate df, is a df with
[TABLE]
and
Theorem 2.3. Consider the following three conditions :
[TABLE]
[TABLE]
[TABLE]
[TABLE]
Then any two of the above conditions imply the remaining one, where are suitable normalizing constants, are non-degenerate df’s, is a df with
[TABLE]
Proof of Theorems 2.2 and 2.3. Without significant modifications, the method of the proof of Theorems 2.2 and 2.3 are the same as that Theorem 2.1, except only the obvious changes. Hence, for brevity the details of the proof are omitted.
2.2 The interrelation of and the basic rv’s is not
restricted
When the interrelation between the random index and the basic variables is not restricted, parallel theorem of Theorem 2.1 may be proved by replacing the condition by a stronger one. Namely, the weak convergence of the df must be replaced by the convergence in probability of the rv to a positive rv However, the key ingredient of the proof of this parallel result is to prove the mixing property, due to Rényi (see, Barakat and Nigm, 1996) of the sequence of order statistics under consideration. In the sense of Rényi a sequence of rv’s is called mixing if for any event of positive probability, the conditional df of under the condition converges weakly to a non-degenerate df, which does not depend on as The following lemma proves the mixing property for the sequence
Lemma 2.2. Under the condition in Theorem 2.1 the sequence is mixing.
Proof. The lemma will be proved if one shows the relation P(Z_{\grave{r},\grave{s}:n}^{(n)}<\mathbf{x}\mid Z_{\grave{r},\grave{s}:l}^{(l)}<\mathbf{x}){\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle w}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}\hat{\Phi}_{{r},{s}}^{(m,k)}({x},y), for all integers The sufficiency of the above relation can easily be proved as a direct multivariate extension of Lemma 6.2.1, of Galambos (1987). However, this relation is equivalent to
[TABLE]
where is the survival function of the limit df i.e.,
[TABLE]
Therefore, our lemma will be established if one proves the relation (2.17). Now, we can write
[TABLE]
[TABLE]
Bearing in mind that all are i.i.d rv’s, the first term in (2.18) can be written in the form
[TABLE]
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
Therefore, in view of (2.18), we have
[TABLE]
where By using the well-known inequalities and for any three events for which we get
[TABLE]
On the other hand, by virtue of the condition in Theorem 2.1, it is easy to prove that
[TABLE]
(note that N\overline{L}_{m}(x_{in}){\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}\kappa_{i}^{m+1}\Rightarrow(N-l)\overline{L}_{m}(x_{in}){\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}\kappa_{i}^{m+1},~{}\forall~{}x_{i}^{\prime}s~{} for which ). By combining the relations (2.19)-(2.21), the proof of the relation (2.17) follows immediately. Hence the required result.
Considering the facts that the normalizing constants, which may be used in the bivariate extreme case are the same as those for the univariate case, and the limit df is continuous, we can easily by using Lemma 2.2, show that the proof of the following theorem follows without any essential modifications as a direct multivariate extension of the proof of Theorem 2.1 in Barakat and El Shandidy (1990), except only the obvious changes.
Theorem 2.4. Consider the condition
[TABLE]
where is a positive rv. Under the conditions of Theorem 2.1 , we have the implication
[TABLE]
3 Illustrative examples
The range and midrange are widely used, particularly in statistical quality control as an estimator of the dispersion tendency and in setting confidence intervals for the population standard deviation as well as in Monte Carlo methods. In fact, the range itself is a very simple measure of dispersion, gives a quick and easy to estimate indication about the spread of data. Let us defines the random generalized ranges and the random generalized midranges and The normalized generalized ranges and the normalized generalized midranges are defined by and respectively, where and are suitable normalizing constants. In this section, some illustrative examples for the most important distribution functions are obtained, which lend further support to our theoretical results. In the following examples we consider an important practical situation when has a geometric distribution with mean In this case we can easily show that P(\frac{\nu_{n}}{n}<z){\begin{array}[]{c}\stackrel{{\scriptstyle\textstyle w}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}H(z)=1-e^{-z}(z\geq 0).
Example 3.1 (standard Cauchy distribution). Let Then In view of Theorems 1.1 and 2.1, Part 1, in Barakat et al. (2015b), we get, after some algebra,
[TABLE]
The random generalized ranges and midranges, for standard Cauchy distribution are given by, if
[TABLE]
and
[TABLE]
respectively, with Moreover, if
[TABLE]
and
[TABLE]
respectively, with Finally, when
[TABLE]
and the df of converge weakly to the same limit, with
Example 3.2 (Pareto distribution). It can be shown that, for the Pareto distribution Therefore, in view of Theorems 1.1 and 2.1, Part 1, in Barakat et al. (2015b), since ~{}~{}\frac{c_{n}}{a_{n}}{\begin{array}[]{c}\stackrel{{\scriptstyle}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}0,\forall~{}m. Then
[TABLE]
and the df of converge weakly to the same limit, where
Example 3.3 (uniform distribution). For the uniform distribution, by using Theorem 2.1, Part 1, in Barakat et al. (2015b), since \frac{a_{n}}{c_{n}}{\begin{array}[]{c}\stackrel{{\scriptstyle}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}1, if Then
[TABLE]
and
[TABLE]
respectively, with
Example 3.4 (Beta distribution). For the beta distribution Therefore, in view of Theorem 2.2, Part 5, in Barakat et al. (2015b), if Then
[TABLE]
and
[TABLE]
respectively, where Clearly, the same result holds for the power distribution
Example 3.5 (standard normal, logistic, Laplace, and log-normal distributions). After some algebra, we get,
[TABLE]
Moreover, for the standard normal, logistic, and Laplace distributions, we get
[TABLE]
and
[TABLE]
for the standard normal distribution Then
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
and
[TABLE]
Finally, for the log-normal distribution,
P({\cal R}^{(n)}_{\nu_{n}}(m,k)\leq r){\begin{array}[]{c}\stackrel{{\scriptstyle w}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}\int_{0}^{\infty}\left[1-\Gamma_{\ell}(ze^{-r(m+1)})\right]e^{-z}dz and the df of converge weakly to the same limit.
Example 3.6 (exponential and Rayleigh distributions). In view of Theorem 2.2, Part 4 in Barakat et al. (2015b). Then
P({\cal R}^{(n)}_{\nu_{n}}(m,k)\leq r){\begin{array}[]{c}\stackrel{{\scriptstyle w}}{{\longrightarrow}}\\ \scriptstyle{n}\end{array}}\int_{0}^{\infty}\left[1-\Gamma_{\ell}(ze^{-r(m+1)})\right]e^{-z}dz and the df of converge weakly to the same limit, for exponential and Rayleigh distributions.
Acknowledgements
Elsawah’s work was partially supported by the UIC GRANT R201409 and the Zhuhai Premier Discipline Grant and Qin’s work was partially supported by the National Natural Science Foundation of China (Nos. 11271147, 11471135, 11471136).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Barakat, H. M. (1997). Asymptotic properties of bivariate random extremes. J. Stat. Plann. Inference, 61, 203-217.
- 2[2] Barakat, H. M. (2007). Limit theory of generalized order statistics. J. Statist. Plann. Inference. Vol. 137, No. 1,1-11.
- 3[3] Barakat, H. M. and El Shandidy, M. A. (1990). On the limit distribution of the extremes of a random number of independent random variables. J. Statist. Plann. Inference, 26, 353-361.
- 4[4] Barakat, H. M. and Nigm, E. M. (1996). The mixing property of order ststistics with some applications. Bull. Malaysian Math. Soc. (Second Series) 19, 39-52.
- 5[5] Barakat, H. M., Nigm, E. M. and Abd Elgawad, M. A. (2014 a). Limit theory for joint generalized order statistics, REVSTAT Statistical Journal 12(3), 199-220.
- 6[6] Barakat, H. M., Nigm, E. M. and Abd Elgawad, M. A. (2014 b). Limit theory for bivariate extreme generalized order statistics and dual generalized order statistics, ALEA, Lat. Am. J. Probab. Math. Stat. 11 (1), 331-340.
- 7[7] Barakat, H. M., Nigm, E. M. and Al-Awady, M. A. (2015 a). Asymptotic properties of multivariate order statistics with random index. Bull. Malays. Math. Sci. Soc. 38(1), 289-301.
- 8[8] Barakat, H. M., Nigm, E. M. and Elsawah, A. M. (2015 b). Asymptotic distributions of the generalized range, midrange, extremal quotient, and extremal product, with a comparison study. Comm. Statist. Theory Methods, 44, 900-913.
