On the Relationship between Pearson Correlation Coefficient and Kendall’s Tau under Bivariate Homogeneous Shock Model

This paper studies the relationship between Kendall's tau and Pearson correlation coefficient under the so-called bivariate homogeneous shock (BHS) model. We find Caperaa-Genest-type inequality may not hold for general BHS model. Computational simulations suggest that the Denials' inequality is likely to be true.


Introduction
Pearson's correlation coefficient r is the mostly used nonparametric measure of association for two random variables. Besides it, the Spearman's ρ and Kendall's τ are two very useful measures of association. As we all know, the Spearman's ρ is the ordinary Pearson's correlation coefficient.
The relationship between ρ and τ has received considerable attention recently. For instance, Hutchinson and Lai 1 conjecture that −1 √ 1 3τ ≤ ρ ≤ min{3τ/2, 2τ − τ 2 } for stochastically increasing random variables. Hürlimann 2 has shown that the entire Hutchinson and Lai conjecture holds for bivariate extreme value distributions. Munroe et al. 3 show that the Hutchinson and Lai conjecture, 1 3τ ≤ 1 ρ 2 , does not hold when X and Y are stochastically increasing. Fredricks and Nelsen 4 show that, under mild regularity conditions, the limit of the ratio ρ/τ is 3/2 as the joint distribution of the random variables approaches to independence. Capéraà and Genest 5 have shown that ρ ≥ τ ≥ 0 when one variable is simultaneously left-tail decreasing and the other right-tail increasing. Genest and Nešlehová 6 give a short analytical proof for classical Daniels' inequality |3τ −2ρ| ≤ 1, found by Daniels 7 .

ISRN Probability and Statistics
However, relatively little attention has been paid on the relationships between r and τ. One example on this track is Edwardes 8 ,in which the author has shown that r τ under a bivariate exponential BVE model introduced by Marshall and Olkin 9 .
Based on a two-component series system subjected to some fatal shocks, Wang and Li 10 introduced a more general bivariate model, which is referred to as bivariate homogeneous shock BHS model. In this paper, we study the relationships between r and τ under BHS model. We find that the most revealed relationships between ρ and τ will no longer hold for the relationships between r and τ. However, computational simulations suggest that the Daniels-type inequality, |3τ − 2r| ≤ 1, is likely to be true.

Main Results
For any pair of random variables X, Y , let F x, y be its bivariate cumulative distribution function. The classical Pearson correlation coefficient r of X, Y is defined as follows: and Kendall's τ is defined as follows: A prominent feature of BHS model is its singularity. More specifically, even though U, V , and W all are continuous random variables, the joint distribution of X and Y is usually discontinuous.
Denote the survival functions of U, V , and W as u x for some baseline survival function u 0 x . Clearly, the BVE model is a proportional hazard model with u 0 x e −x , x ≥ 0. We refer the constants λ, α, and β as to shape constants since they are determined by the structural representation of system. If an association measure depends only on the structural representation of system, we then say it is fully structure determined. Under proportional hazard model, by 2.3 , we can easily obtain τ β λ α β .

2.6
Thus, the Kendall's τ is fully structure determined. The correlation coefficient r, as we can verify, is not fully structure determined. The exact expression of r is not so easy to obtain in general, except for two special cases when u 0 x e −x , x ≥ 0, or u 0 x 1 − x, 0 ≤ x ≤ 1, that is, when the baseline distribution is exponential or uniform. For our convenience, we refer to the first model as model A and the second one as model B. Clearly, model A is just the BVE model. Now we list the main results in the following theorem. i r ≥ τ.
ii The limit of the ratio of r/τ can be any number larger than 1 when the model approaches to independence situation.

2.11
ISRN Probability and Statistics 5 Clearly, we have, max{u, v} ≤ w ≤ u v. With a little bit notational confusion, we relabel 2/u, 2/v and 2/w as u, v, and w, respectively. Then, we have, w ≤ min{u, v}. Without loss of generality, we assume u ≥ v, and then,

2.12
As we can see, the equality holds only when u v w, that is, α λ β 0. Hence, when the three parameters are not all zero, Γ > 1, and thus r > τ. Consider In a similar way, we can show that Φ > 1. We can show that Φ can be any number that is larger than 1. Let α kλ, then, As k → 0, Φ → √ 1 λ, which can be any number that is larger than 1. Denote Ψ 2r − 3τ, then, Since Ψ is symmetric about α and λ, the minimum or maximum of Ψ will be attained on α λ. So we just need to show that the minimum or maximum of Ψ will be between −1 and 1. where N − β 3 β 2 3α 2 β 2α 2 − 2α , D β 3 β 2 5α 2 β 8α 2 6α 4α 3 4α.

2.17
We have,

Computational Simulations
In order to investigate the relationships between r and τ under general BHS model, we conduct some computational simulations. For the sample data { x i , y i , i 1, 2, · · · , n}, the sample's r is computed as follows:

3.1
While the sample's Kendall's τ is computed as follows: where a ij 1 if x i < x j , and −1 if x i > x j , and b ij 1 if y i < y j , and −1 if y i > y j . We compute the sample values of r, τ, r/τ, and 2r − 3τ under several cases. The sample size is set as n 100, and for each computation, the iteration is 100. Table 1 gives the results.
In case 1, we set U, V , and W as uniform variables on 0, 1 . By 2.6 and 2.9 , we obtain, r 0.4 and τ 1/3, and the ratio r/τ 1.2. In case 2, the variables U, V , and W all follow exponential distributions with means 1, 2, and 3, respectively. In this case, r τ 1/3 / 1 1/2 1/3 2/11 0.1818, so the ratio is exactly 1. Based on the computation results for these two cases, we can see that the numerical computations are quite precise. In case 3, we set all variables to follow standard normal distribution. In case 4, we set U Z 2 , V Z 0.2 2 , and W Z 1.0 2 , where Z is the variable of standard normal. In case 5, we set U E 2 1 , V E 2 2 , and W E 2 3 , where E 1 ∼ exp 1 , E 2 ∼ exp 2 , and E 3 ∼ exp 3 . Surprisingly, in this case, the ratio r/τ is less than 1. In all the cases, we find the Daniels' inequality |2r − 3τ| ≤ 1 holds.

Concluding Remarks
The relationship between Spearman's ρ, which is the ordinary Pearson correlation coefficient, and Kendall's τ has been of interest for a long time. However, little attention has been paid on the relationship of Pearson correlation coefficient and Kendall's τ. In this paper, we investigate their relationship under the so-called BHS model. We find that even though for some typical BHS models, the Capéraà-Genest type inequality, r ≥ τ ≥ 0, holds, but for general BHS model, the inequality may not hold. Our simulation studies suggest that the Daniels-type inequality, −1 ≤ 2r − 3τ ≤ 1, will hold under BHS model. We thus conjecture that the Daniels-type inequality will be valid in general. However, theoretical confirmation for such conjecture merits further study.