**Theoretical Economics Letters**

Vol.07 No.02(2017), Article ID:74457,12 pages

10.4236/tel.2017.72021

Statistical Inference for a Novel Health Inequality Index

Cuizhen Niu^{1}, Shaoxin Hong^{2}, Xuejun Jiang^{3*}^{ }

^{1}School of Statistics, Beijing Normal University, Beijing, China

^{2}Department of Mathematics and Statistics, University of North Carolina, Charlotte, USA

^{3}Department of Mathematics, Southern University of Science and Technology of China, Shenzhen, China

Copyright © 2017 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: January 16, 2017; Accepted: February 25, 2017; Published: February 28, 2017

ABSTRACT

In this paper, we develop statistical inference for an important health inequality index proposed by Lv, Wang and Xu [1] for ordinal data. Asymptotic distributions of the indices are established. This allows us to make inference for the indices. Generalizations of the indices to multiple population setting are also studied. We demonstrate the effectiveness of our procedure using the health inequality data of several areas in Switzerland, and our results classify these areas into three classes based on their health inequalities.

**Keywords:**

Health Inequality Index, Multiple Comparison, Self-Reported Health Status, Statistical Inference

1. Introduction

The qualitative nature of SRHS data prevents the straightforward use of conventionally developed indices for measuring income inequality. A reasonable index for SRHS data should be invariant to rescalings of variables which preserve the order of categories.

Assessment on health inequality for ordered data has received attention in the last ten years, [2] and [3] developed median-based concept of inequality. [4] pro- posed polarization measures, which are also median based. These methods are invariant to cardinal scaling on the categories. [5] proposed a method using income-health matrix to measure socioeconomic inequality in health. [6] introduced a family of sub-group decomposable indices and investigated the decomposability of the indices. [7] conducted an empirical study of the health inequality index for ordinal data from China. Reference [8] considered the tools and choices to be made when measuring socioeconomic inequalities with rank-de- pendent inequality indices. [9] made an empirical comparison with several ordinal and cardinal measures of health inequality. [10] proposed a new measure for ordinal health data to monitor income-related health differences between regions in Great Britain. [11] defined a new ratio-scale health status variable and developed positional stochastic dominance conditions that could be implemented in a context of multi-dimensionality categorical variables. [12] examined the measurement of social polarization with categorical and ordinal data. [13] introduced two approaches to measure social polarization in the case where the distance between groups is based on an ordinal variable, such as self-assessed health status. More examples on ordinal inequality measurements can be seen in [14] , [15] and so on. For statistical inference of these recent developed health inequality indices, some authors (e.g. [4] , [16] ) have derived standard errors for the inequality indices they have introduced. [17] presented a unified methodology for the estimation of inequality indices of the cumulative distribution function.

Recently, [1] proposed a class of measures of health inequality, which are easy to compute and have some desirable properties, such as additivity, invariance of parallel shifts, normalization and simple aversion to median-preserving spreads. However, it is designed only for one population and has not developed statistical inference for the index. This motivates us to work along this topic. In this paper, we establish asymptotic distributions of the indices introduced by [1] and extend the indices to multiple population settings. Our procedures allow dependence between the considered populations and different sample sizes. In particular, we answer several important questions, for example, whether the health inequality of one population is the same as others and is there a linear relationship among the health inequalities of different populations?

The reminder of the paper is organized as follows. In Section 2, we review the indices developed by [1] and derive asymptotic distribution of the indices. In Section 3, we develop the indices for multiple populations. Empirical results are reported in Section 4. Section 5 concludes the paper.

2. Inference for the Health Inequality Index

2.1. Review of the Indices

According to [1] , denote $V={\left({V}_{1},\cdots ,{V}_{i},\cdots ,{V}_{n}\right)}^{\text{T}}$ as the health statuses of $n$ individuals. Let $h={\left({h}_{1},\cdots ,{h}_{m}\right)}^{\text{T}}$ be a finite given set of health categories with $m\ge 2$ . Assume that health categories represent various health statuses and satisfy ${h}_{1}<\cdots <{h}_{m}$ . The values of ${h}_{1},\cdots ,{h}_{m}$ are ordinally significant so that, if ${h}_{j}<{h}_{k}$ , then ${h}_{j}$ represents a lower health status than ${h}_{k}$ .

Let $d={\left({d}_{1},\cdots ,{d}_{m}\right)}^{\text{T}}$ and $f={\left({f}_{1},\cdots ,{f}_{m}\right)}^{\text{T}}$ be the empirical and population frequencies of the health categories, respectively, while ${d}_{i}={n}^{-1}{\displaystyle {\sum}_{l=1}^{n}I\left({V}_{l}={h}_{i}\right)}$ represents the relative frequency of individuals with health statuses equal to ${h}_{i}$ . Further let

${\Omega}_{0}\left(f\right)=\text{diag}\left(f\right)-f{f}^{\text{T}},$ (1)

and let $G$ be a $m\times m$ matrix with the $\left(i,j\right)$ th entry being $2g\left(\left|i-j\right|\right)$ , where $g(\cdot )$ is a function of nonnegative integers, such that $0=g\left(0\right)<g\left(1\right)<g\left(2\right)<\cdots <g\left(m-1\right)$ . Lv, Wang and Xu [1] proposed the following classes of health inequalities:

$I\left(f\right)={\displaystyle {\sum}_{i=1}^{m}{\displaystyle {\sum}_{j\ne i}g\left(\left|i-j\right|\right){f}_{i}{f}_{j}}}.$ (2)

Two typical choices for $g(\cdot )$ include the following:

$g\left(i\right)=\frac{2i}{m-1},\text{}i=1,\cdots ,m-1;\text{}g\left(i\right)=2{\alpha}^{m-1-i},\text{}i=1,\cdots ,m-1,\text{}0<\alpha <1.$ (3)

Intuitively, the index is estimated by $I\left(d\right)$ , an empirical plug-in estimator in statistics.

2.2. Asymptotic Results

Base on the above indices, we establish the following asymptotic distribution.

Theorem 1. Using the delta method, we can establish that

$\sqrt{n}\left(I\left(d\right)-I\left(f\right)\right)\Rightarrow N\left(0,{\sigma}^{2}\right),$

where ${\sigma}^{2}={f}^{\text{T}}{G}^{\text{T}}{\Omega}_{0}\left(f\right)Gf$ .

In practice, ${\sigma}^{2}$ is unknown and must be estimated. Given that $d$ is a consistent estimator of $f$ , the asymptotic variance can be estimated by ${\widehat{\sigma}}^{2}={d}^{\text{T}}{G}^{\text{T}}{\Omega}_{0}\left(d\right)Gd$ . Based on the asymptotic result, the two-sided symmetric $100\left(1-\alpha \right)\%$ asymptotic confidence interval for the health inequality index $I\left(f\right)$ can be constructed as

$\left(I\left(d\right)-{z}_{1-\alpha /2}\sqrt{{\widehat{\sigma}}^{2}/n},I\left(d\right)+{z}_{1-\alpha /2}\sqrt{{\widehat{\sigma}}^{2}/n}\right),$

where ${z}_{1-\alpha /2}$ is the $1-\alpha /2$ quantile of the standard normal distribution.

3. Extension to Multiple Populations

3.1. Testing for Equivalence

We first consider two populations with ${f}_{i}={\left({f}_{i1},\cdots ,{f}_{im}\right)}^{\text{T}}$ , ${d}_{i}={\left({d}_{i1},\cdots ,{d}_{im}\right)}^{\text{T}}$ and ${V}_{i}={\left({V}_{i1},\cdots ,{V}_{i{n}_{i}}\right)}^{\text{T}},\text{}i=1,2$ . Our analysis considers the cases of mutually de-

pendent samples and independent samples, with the former being relevant in examining the evolution of health inequalities in a single group (e.g., changes in health inequality over time), while the latter being relevant in comparing health inequality between two groups (e.g., cross-national). The sampling is performed independently within each group.

Lemma 1. Using the delta method, we have

$\begin{array}{l}\sqrt{n}\left(I\left(d\right)-I\left(f\right)\right)=\sqrt{n}{f}^{\text{T}}{G}^{\text{T}}\left(d-f\right)+{o}_{p}\left(1\right)\\ =\frac{1}{\sqrt{n}}{\displaystyle {\sum}_{l=1}^{n}{f}^{\text{T}}{G}^{\text{T}}\left(I\left({V}_{l}=h\right)-f\right)+{o}_{p}\left(1\right)}\end{array}$

Theorem 2. Let ${\sigma}_{ij}$ be the $\left(i,j\right)$ th entry of two populations’ covariance matrix. Denote $\tilde{n}=\mathrm{max}\left\{{n}_{1},{n}_{2}\right\}$ , ${\lambda}_{i}={\mathrm{lim}}_{\tilde{n}\to \infty}\left({n}_{i}/\tilde{n}\right),\text{}0<{\lambda}_{i}\le 1,\text{}i=1,2$ . The asymptotic distribution of $I\left({d}_{1}\right)-I\left({d}_{2}\right)$ is

$\sqrt{\tilde{n}}\left[\left(I\left({d}_{1}\right)-I\left({d}_{2}\right)\right)-\left(I\left({f}_{1}\right)-I\left({f}_{2}\right)\right)\right]\Rightarrow N\left(0,{\sigma}_{11}^{2}/{\lambda}_{1}+{\sigma}_{22}^{2}/{\lambda}_{2}-2{\sigma}_{12}^{2}/\sqrt{{\lambda}_{1}{\lambda}_{2}}\right).$

Now we consider hypothesis testing problem,

${H}_{0,1}:I\left({f}_{1}\right)=I\left({f}_{2}\right)\text{}v.s.\text{}{H}_{a,1}:I\left({f}_{1}\right)\ne I\left({f}_{2}\right).$

We introduce the following Wald statistic:

${T}_{n}=\frac{\tilde{n}{\left(I\left({d}_{1}\right)-I\left({d}_{2}\right)\right)}^{2}}{{\sigma}_{11}^{2}/{\lambda}_{1}+{\sigma}_{22}^{2}/{\lambda}_{2}-2{\sigma}_{12}^{2}/\sqrt{{\lambda}_{1}{\lambda}_{2}}}.$ (4)

Then under the null hypothesis ${H}_{0}$ , ${T}_{n}\Rightarrow {\mathcal{X}}_{1}^{2}$ as $n\to \infty $ . The corresponding $p$ -value can be computed by the following formula:

$p=Pr\left(T\ge {t}_{obs}\right)=1-{F}_{{\mathcal{X}}_{1}^{2}}\left({t}_{obs}\right),$ (5)

where ${F}_{{\mathcal{X}}_{1}^{2}}(\cdot )$ represents the cumulative distribution function of the chi-squared variable with one degree of freedom.

These results are general, an assumption of independent populations is not required, this implies that our test work with the unbalanced designs case. If these two populations are treated as independent, then $\text{Cov}\left\{I\left({V}_{1l}=h\right),{I}^{\tau}\left({V}_{2l}=h\right)\right\}=0$ and thus ${\sigma}_{12}^{2}=0$ . For a particular circumstance, when the sample sizes of these two populations are equal, ${n}_{1}={n}_{2}=n$ , we can have ${\lambda}_{1}={\lambda}_{2}=1$ and then the asymptotic distribution in Theorem 2 reduces to

$\sqrt{n}\left[\left(I\left({d}_{1}\right)-I\left({d}_{2}\right)\right)-\left(I\left({f}_{1}\right)-I\left({f}_{2}\right)\right)\right]\Rightarrow N\left(0,{\sigma}_{11}^{2}+{\sigma}_{22}^{2}-2{\sigma}_{12}^{2}\right).$ (6)

We propose statistical inference procedures to test the equality between samples in terms of their health inequality indices. This equality issue often emerges when checking for the similarity of the health inequalities in the whole country or in a specified region. For example, China, a country consists of many administrative regions, such as Eastern China, North China, and Central Region, with each region having several provinces. Those provinces in the same region have similar economic and/or social behaviors. Therefore, those provinces in the same region are assumed to have the same health inequalities. We also examine whether the health inequality index of a province is the same as the average index of the entire region. The above two testing problems lead to another application. If the preceding analysis reveals that the provinces within each region have equal indices, then we can check whether the common means in two regions are also the same. Accordingly, we cluster the regions based on the test results. In other words, if several regions have the same health inequality, then we can view these regions as one cluster.

3.2. Global Test

Suppose there are $r\left(r\ge 3\right)$ populations with ${f}_{i}={\left({f}_{i1},\cdots ,{f}_{im}\right)}^{\text{T}}$ ,

${d}_{i}={\left({d}_{i1},\cdots ,{d}_{im}\right)}^{\text{T}}$ and ${V}_{i}={\left({V}_{i1},\cdots ,{V}_{i{n}_{i}}\right)}^{\text{T}},\text{}i=1,\cdots ,r$ . For the dependent sam-

ples, we can obtain the similar results as those presented in Section 2. However, the covariance structure becomes too complex to be practical when more samples are used. We only consider independent samples for simplification. A global test can be constructed as:

${H}_{0,2}:I\left({f}_{1}\right)=\cdots =I\left({f}_{r}\right)\text{}v.s.\text{}{H}_{a,2}:I\left({f}_{i}\right)\ne I\left({f}_{j}\right)$ (7)

for some $i\ne j\left(i,j=1,2,\cdots ,r\right)$ .

Define the matrix $R=\left[{I}_{\left(r-1\right)\times \left(r-1\right)},-{L}_{\left(r-1\right)\times 1}\right]$ , where ${I}_{\left(r-1\right)\times \left(r-1\right)}$ is an identity matrix with $r-1$ dimension, and ${L}_{\left(r-1\right)\times 1}$ is a $r-1$ dimensional vector with all the elements being 1. Then, Hypothesis in (7) can be rewritten as follows:

${H}_{0,3}:RI\left(f\right)=0,\text{}v.s.\text{}{H}_{a,3}:RI\left(f\right)\ne 0,$

where $I\left(f\right)={\left(I\left({f}_{1}\right),\cdots ,I\left({f}_{r}\right)\right)}^{\text{T}}$ .

Define $\tilde{n}=\mathrm{max}\left\{{n}_{i},i=1,\cdots ,r\right\}$ , ${\lambda}_{i}={\mathrm{lim}}_{\tilde{n}\to \infty}\left({n}_{i}/\tilde{n}\right),\text{}0<{\lambda}_{i}\le 1,\text{}i=1,\cdots ,r$ . Given the independence of the $r$ groups of samples, we can obtain

$\sqrt{\tilde{n}}\left(I\left(d\right)-I\left(f\right)\right)\Rightarrow N\left(0,\Sigma \right),$ (8)

where $\Sigma =\text{diag}\left({\sigma}_{11}^{2}/{\lambda}_{1},\cdots ,{\sigma}_{rr}^{2}/{\lambda}_{r}\right)$ . Therefore,

$\sqrt{\tilde{n}}\left(RI\left(d\right)-RI\left(f\right)\right)\Rightarrow N\left(0,R\Sigma {R}^{\text{T}}\right)$

and

$\tilde{n}{\left(RI\left(d\right)-RI\left(f\right)\right)}^{\text{T}}{\left(R\Sigma {R}^{\text{T}}\right)}^{-1}\left(RI\left(d\right)-RI\left(f\right)\right)\Rightarrow {\mathcal{X}}_{r-1}^{2}.$ (9)

Note that under the null hypothesis, $RI\left(f\right)=0$ in (9). Consequently, a Wald type of test statistic can be defined as

${T}_{r}=\tilde{n}{I}^{\text{T}}\left(d\right){R}^{\text{T}}{\left(R\widehat{\Sigma}{R}^{\text{T}}\right)}^{-1}RI\left(d\right),$ (10)

where $\widehat{\Sigma}=\text{diag}\left(\tilde{n}{\widehat{\sigma}}_{11}^{2}/{n}_{1},\cdots ,\tilde{n}{\widehat{\sigma}}_{rr}^{2}/{n}_{r}\right)$ is an estimator of $\Sigma $ . Given the central role of the test statistic ${T}_{r}$ , we state the asymptotic behavior of ${T}_{r}$ under the null hypothesis in the following theorem.

Theorem 3. Let $\tilde{n}=\mathrm{max}\left\{{n}_{i},i=1,\cdots ,r\right\}$ , ${\lambda}_{i}={\mathrm{lim}}_{\tilde{n}\to \infty}\left({n}_{i}/\tilde{n}\right)$ , $0<{\lambda}_{i}\le 1$ , $i=1,\cdots ,r$ , then under the null hypothesis ${H}_{0,1}:I\left({f}_{1}\right)=\cdots =I\left({f}_{r}\right)$ in (7), we have ${T}_{r}\Rightarrow {\mathcal{X}}_{r-1}^{2}.$

The corresponding p-value can be computed by:

$p=Pr\left({T}_{r}\ge {t}_{r,obs}\right)=1-{F}_{{\mathcal{X}}_{r-1}^{2}}\left({t}_{r,obs}\right),$ (11)

where ${F}_{{\mathcal{X}}_{r-1}^{2}}(\cdot )$ represents the cumulative distribution function of the chi-squared variable with $r-1$ degrees of freedom. The equality hypothesis (7) can be regarded as a generalization of the two-sample comparison case. The availability of this hypothesis can be seen clearly in our empirical application.

3.3. Hypothesis Testing within a Cluster

Another interesting problem in the multiple sample case is whether the health inequality of a specified population is the same as the average health inequality of entire population. For instance, one may interest to investigate the health inequality level in Hebei province is higher or lower than the average level of all provinces in the North China region. Accordingly, we propose the following testing hypothesis:

${H}_{0,4}:I\left({f}_{j}\right)=I\left({f}_{0}\right)\text{}v.s.\text{}{H}_{a,4}:I\left({f}_{j}\right)\ne I\left({f}_{0}\right)$ (12)

for some $j\left(j=1,2,\cdots ,r\right)$ . If the null hypothesis ${H}_{0,1}$ in (7) holds, then null hypothesis ${H}_{0,4}$ holds naturally. In other words, hypothesis ${H}_{0,4}$ only becomes meaningful when hypothesis ${H}_{0,1}$ is not true.

Define ${a}_{j}={\left(-1/r,\cdots ,-1/r,1-1/r,-1/r,\cdots ,-1/r\right)}^{\text{T}}$ , that is, ${a}_{j}$ is a $r\times 1$ vector with its $j$ -th element being $1-1/r$ and other elements all being $-1/r$ . Hypothesis (12) can be rewritten as follows

${H}_{0,4}:{a}_{j}^{\text{T}}I\left(f\right)=0,\text{}v.s.\text{}{H}_{a,4}:{a}_{j}^{\text{T}}I\left(f\right)\ne 0.$

Recall that $\sqrt{\tilde{n}}\left(I\left(d\right)-I\left(f\right)\right)\Rightarrow N\left(0,\Sigma \right)$ in (8) holds, we can obtain

$\sqrt{\tilde{n}}\left({a}_{j}^{\text{T}}I\left(d\right)-{a}_{j}^{\text{T}}I\left(f\right)\right)\Rightarrow N\left(0,{a}_{j}^{\text{T}}\Sigma {a}_{j}\right).$

Similar to the derivation of ${T}_{r}$ , we can construct the following test statistic

${T}_{r{a}_{j}}=\tilde{n}{I}^{\text{T}}\left(d\right){a}_{j}{\left({a}_{j}^{\text{T}}\Sigma {a}_{j}\right)}^{-1}{a}_{j}^{\text{T}}I\left(d\right).$ (13)

Under the null hypothesis in (12), ${T}_{r{a}_{j}}\Rightarrow {\mathcal{X}}^{2}\left(1\right)$ . Then the p-value can be determined similarly as that for ${T}_{r}$ .

3.4. Hypothesis Testing between Clusters

Further, we discuss the hypothesis testing between clusters. Assume now that our preliminary analysis reveals that the provinces the corresponding region (cluster), such as Eastern China region, have the same health inequality indices. We may then examine whether the health inequalities between two regions are similar. To this end, we choose two representative provinces in each region and then compare their health inequality indices following the proposed approaches in Section 2. However, this method does not employ all information in these groups. To use all underlying information, we compare the common means of these two regions. We consider the following hypothesis:

${H}_{0,5}:I\left({f}_{01}\right)=I\left({f}_{02}\right)\text{versus}{H}_{a,5}:I\left({f}_{01}\right)\ne I\left({f}_{02}\right)$ (14)

where $I\left({f}_{01}\right)={\displaystyle {\sum}_{i=1}^{{\tau}_{1}}I\left({f}_{i}\right)/{r}_{1}}$ and $I\left({f}_{02}\right)={\displaystyle {\sum}_{i={r}_{1}+1}^{{\tau}_{1}+{r}_{2}}I\left({f}_{i}\right)/{r}_{2}}$ .

Without loss of generality, we assume that the first ${r}_{1}$ populations are clustered in one group with a common health inequality $I\left({f}_{01}\right)$ , while the ${r}_{1}+1$ to ${r}_{1}+{r}_{2}$ populations are clustered in another group with another common health inequality $I\left({f}_{02}\right)$ . ${H}_{0,5}$ only becomes meaningful when null hypothesis ${H}_{0,1}$ in (7) is not true. Define

$b={\left(1/{r}_{1},\cdots ,1/{r}_{1},-1/{r}_{2},\cdots ,-1/{r}_{2},0,\cdots ,0\right)}^{\text{T}}.$ .

That is, $b$ is a $r\times 1$ vector with its first ${r}_{1}$ elements being $1/{r}_{1}$ , the ${r}_{1}+1$ to ${r}_{1}+{r}_{2}$ elements being $-1/{r}_{2}$ and the other elements being 0. Similar to the derivation of ${T}_{ra}$ , we can construct the test statistic as follows:

${T}_{rb}=\tilde{n}{I}^{\text{T}}\left(d\right)b{\left({b}^{\text{T}}\Sigma b\right)}^{-1}{b}^{\text{T}}I\left(d\right).$ (15)

Under the null hypothesis in (14), ${T}_{rb}\Rightarrow {\mathcal{X}}^{2}\left(1\right)$ , thus p-value can be determined similarly as that for ${T}_{r}$ .

4. Empirical Application

To illustrate our proposed procedures, we present a real application by using the data of the Swiss Health Survey [SHS] in 2002, conducted by Switzerland's Federal Statistical Office. A total of 19,706 observations were collected from seven areas in Switzerland. The survey respondents were asked to rate their health statuses on a five-point scale ranging from very bad to very good. This dataset was also analyzed by [3] and [6] . We do not include the distributions of SHS in the seven regions in this paper, this information can be found in [6] . We use the health inequality indices proposed by [1] to analyze the survey data and yield new observations. Denote the index with $g\left(\left|i-j\right|\right)=2\left|i-j\right|/\left(m-1\right)$ by F1 and the index with $g\left(\left|i-j\right|\right)=2{\alpha}^{m-1-\left|i-j\right|}$ by F2. For checking the robustness of the results obtained, we choose $\alpha =0.9,0.6$ and 0.3, then these related indices are denoted as F2-1, F2-2 and F2-3, respectively.

Table 1 presents the health inequalities of seven areas in Switzerland based on F1 and F2 with different $\alpha $ . The standard errors are enclosed in parentheses, and the health inequalities are ranked based on the proposed measures. From this table, in all four different measures, we can find that Leman is the region with the highest health inequality value, which implies that the health status exists the most significant difference between Leman citizens. The other regions show ambiguous ranking. Specifically, for F1 and F2-2, Zurich is the region with least difference in health status, and Central has the second-to-the-lowest inequality. However, for F2-3, Central is identified as the least imbalanced region in health status, while Zurich has the second-to-the-lowest inequality. East and Ticino show the similar behavior.

Table 1. Health inequality in the seven statistical areas of Switzerland.

Due to the reason of random sampling of the data set, it is natural to ask questions, like, do East and Ticino have different health inequalities in fact? Do Central and Zurich have the same health inequality actually? We use statistical inferences to address these problems. To fully answer these questions, various interesting two-sample comparison tests are carried out, the results are reported in Table 2. We set the significance level to $5\%$ . From Table 2, we can conclude that Leman is significantly more imbalanced than Middle-Land in health status. In contrast to the findings in Table 1, Middle-Land and North-West do not show statistically significant differences in their health inequalities. In other words, these two regions have the same health inequality level base on the data set we have. North-West is significantly more unbalanced in health status than East. Except for F2-3, all p-values for North-West and Ticino are all smaller than $5\%$ . Therefore, the difference of health inequality between North-West and Ticino can be confirmed almost. Central and Zurich have the same inequality level, and the same finding has been observed for East and Ticino.

Based on the above analysis, we classify North-West and Middle-Land, East and Ticino, and Central and Zurich into three groups. However, can we combine two groups, such as the East and Ticino group with the Central and Zurich group? The question is equivalent to ask whether the average health inequality of the East and Ticino group is the same as that of the other group. The p-values of tests by using the above four measures are 0.5505, 0.1778, 0.7105 and 0.0140, respectively, which are all larger than $5\%$ except for F2-3. Therefore, East, Ticino, Central, and Zurich may be clustered into one group. We also check whether these four regions have the same health inequality levels. The p-values for this global equality hypothesis testing are 0.8805, 0.1824, 0.8946 and 0.0942, respectively, which suggest that these regions have the same inequality levels. We then examine whether this four-member group can be enlarged by including the North-West and Middle-Land group? We propose two hypotheses to investigate this question. First, are the average inequalities of North-West and Middle-Land similar to those of the other groups? Second, do these six regions have the same health inequality levels? For these two hypotheses, all the p-values resulting from tests with the four measures are significantly smaller than $5\%$ , which indicate that the average health inequality of the North-West and Middle-Land group is different from that of the four-member group. We then examine whether the

Table 2. p-values for two-sample comparison problems.

health inequality level of Leman is the same as the average level of the North- West and Middle-Land group. The p-values of all four measures are strongly smaller than $5\%$ , which indicate that the health inequality level of Leman is different from the average level of the North-West and Middle-Land group. In sum, we classify these seven regions into three groups, that is, Leman, North- West and Middle-Land, and the other four regions.

5. Conclusion

In this paper, we propose several statistical inference procedures for the novel health inequality indices introduced in [1] . We consider one-, two-, and multiple-sample cases. Given that health surveys generally cover multiple regions, the health inequalities of multiple sample cases must be tested. The health inequality in various regions of Switzerland validates the availability of our proposed tools. Seven regions covered by SHS can be categorized into three groups after the numerical study; Leman has the highest health inequality followed by the North-West and Middle-Land group. The other four regions (i.e., Central, East, Ticino, and Zurich) have the same health inequality. Our proposed procedures can also be applied to other recently proposed health inequality indices. The subjective well-being is influenced by many factors such as health inequality, education, environment and so on. The statistical inference on multi-dimen- sionality well-being inequality can be investigated ongoing.

Acknowledgements

This research was funded by the Fundamental Research Funds for the Central Universities, China Postdoctoral Science Foundation (2016M600951), National Natural Science Foundation of China (11101432) and Natural Science Foundation of Guangdong Province, China (2016A030313856). The authors thank the editor, the associate editor and the anonymous referee for their constructive comments and suggestions which led to a substantial improvement of an early manuscript. All correspondence should be addressed to Xuejun, Jiang, Department of Mathematics, Southern 6 University of Science and Technology, Shenzhen, China, E-mail: jiangxj@sustc.edu.cn.

Cite this paper

Niu, C.Z., Hong, S.X. and Jiang, X.J. (2017) Statistical Inference for a Novel Health Inequality Index. Theoretical Economics Letters, 7, 251-262. https://doi.org/10.4236/tel.2017.72021

References

- 1. Lv, G. M., Wang, Y. and Xu, Y. S. (2015) On a New Class of Measures for Health Inequality Based on Ordinal Data. The Journal of Economic Inequality, 13, 465-477. https://doi.org/10.1007/s10888-014-9289-4
- 2. Allison, R.A. and Foster, J.E. (2004) Measuring Health Inequalities Using Qualitative Data. Journal of Health Economics, 23, 505-524. https://doi.org/10.1016/j.jhealeco.2003.10.006
- 3. Abul Naga, R.H. and Yalcin, T. (2008) Inequality Measurement for Ordered Response Health Data. Journal of Health Economics, 27, 1614-1625. https://doi.org/10.1016/j.jhealeco.2008.07.015
- 4. Apouey, B. (2007) Measuring Health Polarization with Self-Assessed Health Data. Health Economics, 16, 875-894. https://doi.org/10.1002/hec.1284
- 5. Zheng, B. (2011) A New Approach to Measure Socioeconomic Inequality in Health. The Journal of Economic Inequality, 9, 555-577. https://doi.org/10.1007/s10888-010-9148-x
- 6. Kobus, M. and Milos, P. (2012) Inequality Decomposition by Population Subgroups for Ordinal Data. Journal of Health Economics, 31, 15-21. https://doi.org/10.1016/j.jhealeco.2011.11.005
- 7. Wang, H.L. and Yu, Y.W. (2016) Increasing Health Inequality in China: An Empirical Study with Ordinal Data. The Journal of Economic Inequality, 14, 41-61. https://doi.org/10.1007/s10888-015-9315-1
- 8. Erreygers, G. and Van Ourti, T. (2011) Measuring Socioeconomic Inequality in Health, Health Care and Health Financing by Means of Rank-Dependent Indices: A Recipe for Good Practice. Journal of Health Economics, 30, 685-694. https://doi.org/10.1016/j.jhealeco.2011.04.004
- 9. Madden, D. (2010). Ordinal and cardinal measures of health inequality: an empirical comparison. Health Economics, 19, 243-250. https://doi.org/10.1002/hec.1472
- 10. Allanson, P. (2017) Monitoring Income-Related Health Differences between Regions in Great Britain: A New Measure for Ordinal Health Data. Social Science and Medicine, 175, 72-80. https://doi.org/10.1016/j.socscimed.2016.12.033
- 11. Makdissi, P. and Yazbeck, M. (2014) MEASURING Socioeconomic Health Inequalities in Presence of Multiple Categorical information. Journal of Health Economics, 34, 84-95. https://doi.org/10.1016/j.jhealeco.2013.11.008
- 12. Permanyer, I. and D’ambrosio, C. (2015) Measuring Social Polarization with Ordinal and Categorical Data. Journal of Public Economic Theory, 17, 311-327. https://doi.org/10.1111/jpet.12093
- 13. Fusco, A. and Silber, J. (2014) On Social Polarization and Ordinal Variables: The Case of Self-Assessed Health. The European Journal of Health Economics, 15, 841-851. https://doi.org/10.1007/s10198-013-0529-5
- 14. Erreygers, G. (2009) Correcting the Concentration Index. Journal of Health Economics, 28, 504-515. https://doi.org/10.1016/j.jhealeco.2008.02.003
- 15. Lazar, A. and Silber, J. (2013) On the Cardinal Measurement of Health Inequality When Only Ordinal Information Is Available on Individual Health Status. Health Economics, 22, 106-113. https://doi.org/10.1002/hec.1821
- 16. Cowell, F. and Flachaire, E. (2014) Inequality with Ordinal data. Unpublished.
- 17. Abul Naga, R.H. and Stapenhurst, C. (2015) Estimation of Inequality Indices of the Cumulative Distribution Function. Economics Letters, 130, 109-112. https://doi.org/10.1016/j.econlet.2015.03.004

Appendix. Proofs of Theorems

Proof of Theorem 1.

Given that $d$ is a consistent estimator of $f$ , for the empirical frequency of the health categories, we can easily obtain the following:

$\sqrt{n}\left(d-f\right)\Rightarrow N\left(0,{\Omega}_{0}\left(f\right)\right),$ (A.1)

where

${\Omega}_{0}\left(f\right)=\text{diag}\left(f\right)-f{f}^{\text{T}}.$

Note that $g\left(\left|i-j\right|\right){f}_{i}{f}_{j}=g\left(\left|j-i\right|\right){f}_{j}{f}_{i}$ . Therefore,

$I\left(f\right)=2{\displaystyle \sum _{i=1}^{m-1}{\displaystyle \sum _{j=i+1}^{m}g\left(\left|i-j\right|\right){f}_{i}{f}_{j}}}.$

Define $J={\left(\frac{\partial I\left(f\right)}{\partial {f}_{1}},\cdots ,\frac{\partial I\left(f\right)}{\partial {f}_{m}}\right)}^{\text{T}}$ It can be easily shown that

$\frac{\partial I\left(f\right)}{\partial {f}_{i}}=2{\displaystyle \sum _{l\ne i}g\left(\left|l-i\right|\right){f}_{l}}.$ (A.2)

Alternatively, we can have $J=Gf$ . Keeping only the first two terms of the Taylor expansion, we can estimate $I\left(f\right)$ as

$I\left(d\right)\approx I\left(f\right)+{J}^{\text{T}}\cdot \left(d-f\right).$

Then the variance of $I\left(d\right)-I\left(f\right)$ is approximated by

$\begin{array}{c}\text{Var}\left({J}^{\text{T}}\cdot \left(d-f\right)\right)={J}^{\text{T}}\cdot \text{Cov}\left(d-f,d-f\right)\cdot J\\ ={f}^{\text{T}}{G}^{\text{T}}\left({\Omega}_{0}/n\right)Gf.\end{array}$

Also $E\left(I\left(d\right)\right)=E\left(I\left(f\right)\right)$ since $d$ is a consistent estimator of $f$ , it follows that

$\sqrt{n}\left(I\left(d\right)-I\left(f\right)\right)\Rightarrow N\left(0,{\sigma}^{2}\right),$

where ${\sigma}^{2}={f}^{\text{T}}{G}^{\text{T}}{\Omega}_{0}\left(f\right)Gf$ . $$

Proof of Theorem 2.

Define $\tilde{n}=\mathrm{max}\left\{{n}_{1},{n}_{2}\right\}$ , ${\lambda}_{i}={\mathrm{lim}}_{\tilde{n}\to \infty}\left({n}_{i}/\tilde{n}\right),\text{}0<{\lambda}_{i}\le 1,\text{}i=1,2$ . From Lemma 1, we have

$\begin{array}{l}\sqrt{n}\left(I\left(d\right)-I\left(f\right)\right)=\sqrt{n}{f}^{\text{T}}{G}^{\text{T}}\left(d-f\right)+{o}_{p}\left(1\right)\\ =\frac{1}{\sqrt{n}}{\displaystyle {\sum}_{l=1}^{n}{f}^{\text{T}}}{G}^{\text{T}}\left(I\left({V}_{l}=h\right)-f\right)+{o}_{p}\left(1\right).\end{array}$

Let ${\sigma}_{h}^{2}=\text{Cov}\left(I\left({V}_{1l}=h\right){I}^{\text{T}}\left({V}_{2l}=h\right)\right)$ . After that, similarly as the proof of Theorem 1, it can be derived that

$\begin{array}{c}\sqrt{\tilde{n}}\left(\begin{array}{c}I\left({d}_{1}\right)-I\left({f}_{1}\right)\\ I\left({d}_{2}\right)-I\left({f}_{2}\right)\end{array}\right)=\left(\begin{array}{c}\sqrt{\frac{\tilde{n}}{{n}_{1}}}\times \sqrt{{n}_{1}}\left(I\left({d}_{1}\right)-I\left({f}_{1}\right)\right)\\ \sqrt{\frac{\tilde{n}}{{n}_{2}}}\times \sqrt{{n}_{2}}\left(I\left({d}_{2}\right)-I\left({f}_{2}\right)\right)\end{array}\right)\end{array}$ $=\left(\begin{array}{c}\sqrt{\frac{\tilde{n}}{{n}_{1}}}\times \frac{1}{\sqrt{{n}_{1}}}{\displaystyle {\sum}_{l=1}^{{n}_{1}}{f}_{1}^{\text{T}}}{G}^{\text{T}}\left(I\left({V}_{1l}=h\right)-{f}_{1}\right)\\ \sqrt{\frac{\tilde{n}}{{n}_{2}}}\times \frac{1}{\sqrt{{n}_{2}}}{\displaystyle {\sum}_{k=1}^{{n}_{2}}{f}_{2}^{\text{T}}}{G}^{\text{T}}\left(I\left({V}_{2k}=h\right)-{f}_{2}\right)\end{array}\right)$ $+{o}_{p}\left(1\right)\Rightarrow N\left(\left(\begin{array}{c}0\\ 0\end{array}\right),\left(\begin{array}{cc}{\sigma}_{11}^{2}/{\lambda}_{1}& {\sigma}_{12}^{2}/\sqrt{{\lambda}_{1}{\lambda}_{2}}\\ {\sigma}_{21}^{2}/\sqrt{{\lambda}_{1}{\lambda}_{2}}& {\sigma}_{22}^{2}/{\lambda}_{2}\end{array}\right)\right).$ (A.3)

Here,

$\begin{array}{l}{\sigma}_{11}^{2}={f}_{1}^{\text{T}}{G}^{\text{T}}{\Omega}_{0}\left({f}_{1}\right)G{f}_{1},\\ {\sigma}_{12}^{2}={\sigma}_{21}^{2}={f}_{1}^{\text{T}}{G}^{\text{T}}\text{Cov}\left(I\left({V}_{1l}=h\right){I}^{\text{T}}\left({V}_{2l}=h\right)\right)G{f}_{2}={f}_{1}^{\text{T}}{G}^{\text{T}}{\sigma}_{h}^{2}G{f}_{2},\\ {\sigma}_{22}^{2}={f}_{2}^{\text{T}}{G}_{0}^{\text{T}}\left({f}_{2}\right)G{f}_{2}\end{array}$

Plug in the consistent estimators of ${f}_{i}$ and ${\sigma}_{h}^{2}$ by ${d}_{i}$ and

${\widehat{\sigma}}_{h}^{2}={\left\{\mathrm{min}\left({n}_{1},{n}_{2}\right)\right\}}^{-1}{\displaystyle \sum _{l=1}^{\mathrm{min}\left\{{n}_{1},{n}_{2}\right\}}I\left({V}_{1l}=h\right){I}^{\text{T}}\left({V}_{2l}=h\right)-{d}_{1}{d}_{2}^{\text{T}}},$

respectively, thus we can easily estimate ${\sigma}_{i,j}^{2},i,j=1,2$ consistently.

Combining the above result in (A.3), we can obtain

$\begin{array}{l}\sqrt{\tilde{n}}\left[\left(I\left({d}_{1}\right)-I\left({d}_{2}\right)\right)-\left(I\left({f}_{1}\right)-I\left({f}_{2}\right)\right)\right]\\ =\sqrt{\tilde{n}}\left[\left(I\left({d}_{1}\right)-I\left({f}_{1}\right)\right)-\left(I\left({d}_{2}\right)-I\left({f}_{2}\right)\right)\right]\\ \Rightarrow N\left(0,{\sigma}_{11}^{2}/{\lambda}_{1}+{\sigma}_{22}^{2}/{\lambda}_{2}-2{\sigma}_{12}^{2}/\sqrt{{\lambda}_{1}{\lambda}_{2}}\right).\end{array}$ $$