Background {#Sec1}
==========

The determinants of health outcomes include well-known factors such as income, medical technology and education. Some recent empirical studies have focused on the important role of information and communication technology (ICT) \[[@CR1]\], income inequality \[[@CR2]\], ethnic diversity \[[@CR1], [@CR3], [@CR4]\] and public healthcare expenditure \[[@CR5]\] to explain the observed health outcomes of countries. However, most previous research (both time-series and panel data estimation) have used a single variable, such as life expectancy or mortality rates (infant, maternal or under-five) to measure health outcomes \[[@CR5]--[@CR7]\]. Yet other studies have used outcome variables such as the prevalence or incidence of diseases or infections, the number of visits and readmissions to hospitals, hospital length of stay, work absences due to illness, perceived health status, quality-adjusted life years and disability-adjusted life years. Each of these measures and some recently developed health outcome indices (including the 'better life index published' by the Organisation for Economic Co-operation and Development (OECD) and 'urban health index' published by the World Health Organisation (WHO) is justifiable to answer individual research questions (often with cross-sectional data). Although reliable, these measures are inadequate for time series and panel data analysis. This is because continuous data on some of the indicators (e.g. self-reported health, disability, healthcare access) used in these indices are unavailable. For example, none of the health outcome indices provide data for an adequate number of time periods to enable statistically robust estimation to be undertaken. Also, studies estimating health outcomes with panel data are often obliged to use a single variable (mostly mortality variables) for estimation. Unfortunately, even for the OECD countries, which have rich sources of data for both mortality and morbidity variables, a comprehensive health outcome index suitable for panel data estimation is non-existent.

The current study constructs a health outcome index that incorporates both mortality (e.g. life expectancy and infant mortality) and morbidity (e.g. perceived health status) variables. The main objective of this study is to re-examine the factors that determined health outcomes for 30 OECD countries over the period 2004 to 2015. Along with some conventional determinants (e.g. income, education and technology), the estimated model also incorporates the influence of income inequality, ethnic diversity, development in ICT, and public health expenditure on health outcomes into the estimation process.

Past empirical findings have showed that societies with unequal incomes often experience poorer health outcomes \[[@CR8]--[@CR11]\]. However, the hypothesis that income inequality contributes to poorer health outcomes is not definitively determined with mixed empirical findings. Some studies have found a significant relationship between health outcomes (mortality, longevity) and income inequality \[[@CR2], [@CR8], [@CR9], [@CR12], [@CR13]\] and others have concluded no significant association \[[@CR11], [@CR14], [@CR15]\]. Furthermore, Rözer and Volker \[[@CR16]\] found that the impact of income inequality on health reduces to zero after the age of 36 and Li and Zhu \[[@CR10]\] reported a U shaped relationship, which means that at lower-levels of inequality, health outcomes improve with rising inequality and vice-versa. Hence, despite the large volume of literature, it is seemingly impossible to definitively determine the causality of income inequality on public health outcomes \[[@CR17]\]. However, these past studies have mostly used cross-sectional data and a single variable to measure health outcomes. The health impacts of inequality will develop over a long period of time as the negative social consequences of that inequality will gradually take account of an individual's health status \[[@CR18]\]. Therefore, a representative health outcome index and panel data of 12 continuous years have been used to generate a better understanding of the issue.

Another key variable of interest in this study is ethnic diversity. Previous researchers have primarily focused on whether there are disparities in healthcare access, use and health status amongst majority and ethnic minority population \[[@CR3], [@CR4], [@CR19], [@CR20]\]. Most of the empirical findings indicated that ethnic minorities experience poorer access to healthcare services and their health status is significantly inferior to the majority of the population \[[@CR3], [@CR21]\]. Nonetheless, it is still unclear as to whether having higher levels of ethnic diversity reduces the overall health status of a country or not. Moreover, no study has yet examined the influence of public health expenditure on the potential adverse effects of higher levels of ethnic diversity on health outcomes.

Better health outcomes are also firmly linked with ICT development. ICT can support and promote efficient disease management and prevention \[[@CR1]\], effective communication and reduction in the knowledge gap between patient and provider \[[@CR22], [@CR23]\], facilitate long-distance medical consultations \[[@CR24]\], underpin the prevention of and supervision of public health risks \[[@CR25], [@CR26]\] and boost the performance and administration of the health system \[[@CR27]\]. On the contrary, some empirical studies reported adverse effects of ICT on health outcomes. For example, Kim et al. \[[@CR28]\] found that higher usage of the internet is associated with irregular food intake; Rosell et al. \[[@CR29]\] indicated an association between mental distress and ICT use; Booth et al. \[[@CR30]\] linked ICT use with less physical activity; and Punamäki et al. \[[@CR31]\] concluded that ICT use is related to lower perceived health status. The contradictory findings may arise because previous studies used variables of choice that most suited the study purpose- to scale health outcome (mortality or morbidity) or ICT development (adoption or diffusion). Therefore, the overall impact of ICT developments on health outcomes is still unconfirmed \[[@CR24]\]. In particular, what is the impact on overall public health when there is increased access to and use of ICT?

In light of the above discussion, the primary objective of this study is to examine the impact of income inequality, ethnic diversity and ICT development on health outcomes in the OECD countries utilising a large panel data. The secondary objective is to construct a health outcome index with indicators representing both mortality and morbidity variables. This study, therefore, investigates the following key questions: (i) does income inequality significantly affect health outcomes? (ii) do countries with greater levels of ethnic diversity have lower health status than countries with lower levels of diversity? (iii) what impact does ICT development (access, use and skill level) have on health outcomes? and (iv) does public health expenditure significantly improve health outcomes? These questions are empirically investigated by developing a health outcome index and using a panel dataset of OECD countries over the period of 2004 to 2015 (Additional file [1](#MOESM1){ref-type="media"}: Appendix A).

This study will make the following contributions to the existing literature. First, instead of using a single variable to measure health outcome, a health outcome index is constructed using both morbidity and mortality variables in a partial least square estimator based on a structural equation model (PLS-SEM). To the best of the authors' knowledge this is the first study that has used perceived health status as an indicator of health outcome in panel data analysis. In the OECD database perceived health status was measured with the question "How is your health in general" and the responses were recorded in a Likert scale (from very bad to very good). Second, given that the adverse effects of income inequality on health occurs over time, this study assembled income inequality data and constructed a health outcome index over 12 continuous years. Third, the relationship between higher ethnic diversity and overall health outcomes was examined. Previous studies have investigated health care access of ethnic minorities but no previous study has examined health outcome variances across countries with different levels of ethnic diversity. Fourth, the relationship between ICT development and health outcome was re-examined using credible ICT development index data from the International Telecommunication Union (ITU). Furthermore, this study also investigated which aspect of ICT (access, use or skills) development is most strongly related to population health outcomes. Fifth, unlike previous studies, this study adopted a heterogeneous panel study method, rather than use homogenous panel data analysis. This enables control of cross-section dependence in the panel data by using OLS estimation with panel-corrected standard errors.

The paper consists of five sections. Section "[Methods](#Sec2){ref-type="sec"}" outlines the model variables and research method. Section "[Results](#Sec10){ref-type="sec"}" shows the results and section "[Discussion and Policy Implications](#Sec11){ref-type="sec"}" consists of the discussion and finally, section "[Conclusions](#Sec12){ref-type="sec"}" draws conclusions of the study.

Methods {#Sec2}
=======

Data and variables {#Sec3}
------------------

This paper utilises a heterogeneous panel data for 30 OECD countries over the period 2004 to 2015. The rationale behind selecting the countries and the time period is the availability of data, especially for the variables of Gini coefficient and perceived health status. Data on infant mortality rate and income inequality came from World Development Indicators database \[[@CR32]\]; ICT index data from World Telecommunication/ICT indicators database \[[@CR33]\], population ethnicity data from the CIA World Factbook \[[@CR34]\] and the remainder was sourced from the OECD database \[[@CR35]\].

For the purpose of empirical analysis, apart from the measured variable HO*index*, other two key variables of interest are: income inequality (Gini) measured by the *Gini coefficient or Gini index* and the ICT*index* that represents the ICT development index which ranks all countries based on access, use and skills index \[[@CR33]\]. To investigate the role of ethnic diversity (Eth*dum*) and public health expenditure (PHE*dum*), two dummy variables were generated. Ethnic diversity was low when the percentage of the largest ethnic group is high and vice versa. For example, in South Korea, 98% of the population are from the one ethnic background (Korean) yet in Switzerland, it is only 65%. Average ethnic diversity was calculated for the 30 countries (81.30%) and the variable was recorded 1 for countries below the average value and 0 otherwise. This being a time-invariant variable, only 2014 data was used. For PHE*dum* the mean value of public health expenditure as a percentage of total health expenditure was calculated. Countries with values above the average (71.91%) were scaled as 1 (generous levels of public health expenditure) and 0 otherwise. The average value in the year 2004 was 70.72 and 72.97% in 2015.

Model specifications {#Sec4}
--------------------

The estimated model also controlled for the log of GDP per capita (lnGDPpc) expressed at USD current PPP prices; the percentage of obese population (males 15--64 years) (ObstM); the percentage of population (between 15 and 64 years) with tertiary education (TerEdu); public expenditure on research and development (RdExp) as a proxy for technological advancement; growth in actual hours of work per person per year (lnHrWrk); and the suicide rate (SudRt) as total suicides per 100,000 population as proxy for the overall mental health status of a country \[[@CR36], [@CR37]\]. Any missing data were replaced by linear interpolation \[[@CR38]\]. The empirical models are specified as follows,$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \mathrm{HO}{index}_{it}=\alpha +\kern0.5em {\beta}_1\mathrm{ICT}{index}_{it}+{\beta}_2{\mathrm{Gini}}_{it}+{\beta}_3{\mathrm{lnGDPpc}}_{it}+{\beta}_4{ObsM}_{it}+{\beta}_5{TerEdu}_{it}+{\beta}_6\ {RdExp}_{it}+{\beta}_7\ {lnHrWrk}_{it}+{\beta}_8{SudRt}_{it}+{e}_i $$\end{document}$$$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \mathrm{HO}{index}_{it}=\alpha +\kern0.5em {\beta}_1\mathrm{ICT}{index}_{it}+{\beta}_2{\mathrm{Gini}}_{it}+{\beta}_3{\mathrm{lnGDPpc}}_{it}+{\beta}_4{ObsM}_{it}+{\beta}_5{TerEdu}_{it}+{\beta}_6\ {RdExp}_{it}+{\beta}_7\ {lnHrWrk}_{it}+{\beta}_8{SudRt}_{it}+{\beta}_8{Dummy}_{it}+{e}_i $$\end{document}$$$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \mathrm{HO}{index}_{it}=\alpha +\kern0.5em {\beta}_1\mathrm{ICT}{index}_{it}+{\beta}_2{\mathrm{Gini}}_{it}+{\beta}_3{\mathrm{lnGDPpc}}_{it}+{\beta}_4{ObsM}_{it}+{\beta}_5{TerEdu}_{it}+{\beta}_6\ {RdExp}_{it}+{\beta}_7\ {lnHrWrk}_{it}+{\beta}_8{SudRt}_{it}+{\beta}_8 Interaction\ {\mathit{\operatorname{var}}}_{it}+{e}_i $$\end{document}$$

Lastly, some interaction variables were generated to better understand the relationship among the variables in the model and to test some additional hypothesis. They are:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \mathrm{EthPhe}=\mathrm{Eth} dum\ \mathrm{multiplied}\ \mathrm{by}\ \mathrm{PHE} dum $$\end{document}$$$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \mathrm{SklPhe}=\mathrm{ICT} skl\ \mathrm{multiplied}\ \mathrm{by}\ \mathrm{PHE} dum $$\end{document}$$

Equations [1](#Equ1){ref-type=""}, [2](#Equ2){ref-type=""} and [3](#Equ3){ref-type=""} are used for statistical estimation by panel corrected standard errors (PCSE) and by pooled OLS Driscoll and Kraay standard errors (DKSE) models. Initial diagnostic tests showed that the panel data displayed cross-section dependence \[[@CR39]\] and heteroskedasticity \[[@CR40]\] but no serial correlation \[[@CR41]\]. Ignoring these issues will seriously compromise the estimated results. The PCSE estimation method proposed by Beck and Katz \[[@CR42]\] can account for these data difficulties and provide accurate estimates from linear models with panel data \[[@CR43], [@CR44]\]. Driscoll and Kraay \[[@CR45]\] recommended a nonparametric covariance matrix estimator which can generate consistent standard errors in the presence of heteroscedasticity when the errors are spatial and temporal-dependent \[[@CR46]\]. These estimation approaches will provide efficient and robust outcomes, suitable for formulating accurate conclusions.

Construction of the health outcome index {#Sec5}
----------------------------------------

The rationale behind constructing the health outcome index was the unavailability of data for the appropriate indicator variables with which the health outcome indices were developed. This study does not propose that the newly created index is better than the existing ones. However, it made panel estimation possible using a continuous time series dataset. The reliability of the index was tested by estimating the Cronbach coefficient alpha and the robustness and sensitivity of the index was examined by including and excluding indicators and altering weights in the aggregation process \[[@CR47]\]. The findings are reported in Additional file [1](#MOESM1){ref-type="media"}: Appendix A.

The construction of the health index has four steps which include data of 36 countries for the year 2014.

### Variable selection {#Sec6}

The initial idea was to construct a health outcome index with three mortality variables (life expectancy at birth, infant mortality and death due to cancer) and three morbidity variables (perceived health status, absence at work due to illness and low birth weight). However, the partial least square based structural equation model (PLS-SEM) showed that only three variables; life expectancy at birth (LEB), infant mortality rate (IMR) and perceived health status as Bad (PHSB) were significant variables in the model.

To ensure that all three variables have a similar causal direction towards the dependent variable (i.e. health outcome data for the life expectancy at birth) LEB was modified into LEB1 (100- LEB) with the hypothesis that 100 years is the highest achievable life expectancy in the time period. Hence, a country with LEB = 82 had a revised LEB1 value of 18 (100--82) and a country with LEB 75 earned a LEB1 value of 25 (100--75). This ensured that all the indicator variables LEB1, IMR and PHSB are negatively related to health outcome, the latent variable. Additional file [1](#MOESM1){ref-type="media"}: Appendix A illustrates the detail process of the selection of the indicators.

### Calculating weights {#Sec7}

The estimation method for constructing the index is founded on the structural equation model (SEM) proposed by Wright \[[@CR48]\]. The approach incorporates the advantages of both multiple regression and principal components analysis \[[@CR49], [@CR50]\]. Using SEM we have identified the structural model with the set of indicators (LEB1, IMR and PHSB) used as the latent variables. One benefit of using the SEM approach is that it can measure the reliability of the relationship between indicators and the latent variable, and is suitable to control for the multicollinearity problem among the indicators in the model. Urbach and Ahlemann \[[@CR51]\] outlined the detailed methodology of the SEM approach.

This study used the PLS-SEM which aims to estimate the association between indicator variables and the latent variables with a system of autonomous equations based on simple and multiple regressions \[[@CR52]\]. This study employed the reflective measurement model which considers the indicator variables as reflections of the latent variables and the path of causality runs from construct to the indicators \[[@CR52], [@CR53]\]. The reflective measure tests the reliability of the estimation while considering the indicator variables as a linear combination of its latent variable. The estimated model is specified as,$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ {x}_{pq}={\lambda}_{pq}{\xi}_q+{\varepsilon}_{pq} $$\end{document}$$where *x*~*pq*~ indicates indicator variable, *ξ*~*q*~ is the latent variable, *λ*~*pq*~ represents the generic loading coefficients of the *p*-th indicator variable in the q block and *ε*~*pq*~ is the residual terms of the generic indicator variables \[[@CR53]\]. See Chin \[[@CR52]\] for a detail discussion of the method.

SmartPLS (version 2) software was used to run the model.

Lastly, the factor loadings for each indicator were divided by the summation of the factor loadings (0.758 + 0.897 + 0.782 = 2.437) of the three indicators to calculate the weights which are 0.311 (0.758/2.437) for IMR, 0.369 (0.897/2.437) for LEB1 and 0.320 (0.782/2.437) for PHSB (Fig. [1](#Fig1){ref-type="fig"}). This process normalised the corresponding loadings of each indicator and eased the comparability of the weights. All the variables are statistically significant at a 5% confidence interval and an R-Sq = 0.425 indicates that 42.5% of the observed variable (per capita GDP in Fig. [1](#Fig1){ref-type="fig"}) is explained by the indicators in the model.Fig. 1*p-values* of the proposed structural equation model

### Standardization of the indicators {#Sec8}

The aim of this step was to accommodate disparities in the metrics and scales of the indicators. The score of each unit was standardised by dividing the difference between the current value and lowest value with the difference between the highest value and lowest value. This method has been used in previous studies, for example in constructing the 'Urban Health Index' \[[@CR54]\].The simple standardisation formula is as follows,$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ {\mathrm{I}}^{\mathrm{SD}}=\frac{\mathrm{CV}-\mathrm{LV}\ }{\mathrm{HV}-\mathrm{LV}} $$\end{document}$$

where *I*^*SD*^ stands for the standardised value of the indicator, *CV, LV* and *HV* are the current value, lowest value and highest value, respectively, for each of the indicators among the 36 countries in a particular year.

### Estimation of index scores {#Sec9}

After standardisation of the values of the indicators, the next step was to integrate the *I*^*SD*^into a single composite index. To do so, standard value of each indicator was multiplied by corresponding normalised weights. Afterwards, each sub-index scores were summed for individual countries for each year to construct the health index scores. The basic health outcome index is described by the following equation,$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ \mathrm{HO} index=\sum \left(\mathrm{ZLEB}1\ \mathrm{x}\ {\mathrm{W}}_{LEB1\kern0.5em }+\mathrm{ZIMR}\ \mathrm{x}\ {\mathrm{W}}_{IMR}+\mathrm{ZPHSB}\ \mathrm{x}\ {\mathrm{W}}_{PHSB}\right) $$\end{document}$$

Here Z indicates the standardised value of the indicators. As mentioned, all of the indicators are negatively related to health outcome or better health status. Therefore, for ease of interpretation of the findings, we deducted the HO*index* scores from the value of 1 (the maximum attainable value). The modified value indicates that higher the index score, the better the health outcomes or health status of the country. Table [1](#Tab1){ref-type="table"} illustrates the HO*index* scores for the OECD countries for the year 2014.Table 1The calculated health outcome index and indicator variables for 2014RankingCountryHO IndexLEB1IMRPHSB1Sweden0.92088217.72.442Switzerland0.90660916.73.53.83Iceland0.90114917.11.66.34Spain0.85146516.73.68.35Norway0.85006617.82.26.96Ireland0.84728517.93.147Canada0.84627517.94.43.18Finland0.83029918.726.49Netherlands0.82462718.23.35.410Luxembourg0.81179517.71.68.311France0.80676217.23.68.412Italy0.77200416.831213Austria0.76405618.438.914Germany0.75964118.83.2815Slovenia0.74435618.82.21116Belgium0.73628818.63.49.217Denmark0.72894719.237.418UK0.72370618.63.7919Greece0.69722118.53.710.820USA0.68218221.25.72.821Israel0.64861517.83.315.522Korea0.64367817.8315.923Czech Republic0.61780421.12.911.624Portugal0.5564418.8318.325Poland0.46505222.34.513.726Estonia0.43680422.82.516.527Slovak Republic0.39952423.16.112.728Hungary0.33585824.15.315.929Turkey0.2858852212.311.830Latvia0.26899125.77.217.1Note: These 30 countries were selected due to availability of continuous time-series data

Results {#Sec10}
=======

Table [2](#Tab2){ref-type="table"} shows the descriptive statistics of the variables. It includes the mean values and the range (highest minus lowest) of the values. During the study period the mean health outcome index improved slightly from 0.70 to 0.73 however, the disparity in health outcomes amongst the 30 countries reduced significantly from 0.77 in 2004 to 0.54 in 2015. This finding indicates that countries with historically poor health status are gradually closing the gap to countries with traditional better health status. Average income inequality reduced slightly (33.01 to 31.80) and the ICT development index increased from 4.55 to 7.70. Noticeably, the number of obese population (male) has risen alarmingly from 17.09 in 2004 to 22.43 in 2015 and the suicide rate showed a minor decrease from 114.12 to 111.67.Table 2Descriptive analysisVariableMeanMax-MinMeanMax-MinMeanMax-MinMeanMax-MinMeanMax-MinMeanMax-MinYear200420052006200720082009 HO*index*0.700.770.700.740.690.740.690.730.690.700.670.77 Gini33.01116.31332.72117.07332.30118.12332.26117.38332.04117.78331.97117.37 ICT*index*44.5533.6455.2944.9055.6744.3266.0533.7466.4333.9566.6044.03 ICT*acs*55.5855.2266.1044.8966.6144.5777.1344.2477.0233.9877.2033.93 ICT*use*11.5922.8422.3233.5533.0544.2633.7944.9744.4655.3355.0133.98 ICT*skl*88.3233.0788.4333.1788.5433.3188.6433.6088.7133.0088.7633.05 RdExp11.6933.4811.7333.5611.7533.6711.7633.9811.8733.8811.9233.68 lnGDPpc110.1711.78110.2211.76110.3111.75110.3711.74110.4211.69110.3811.68 ObstM117.09222.20117.57223.00118.08223.60118.60224.20119.11224.80119.63225.30 TerEdu330.81441.72332.14441.18333.14441.16334.06441.47335.37442.39336.92446.45 lnHrWk77.4700.5277.4700.5177.4600.5077.4600.4877.4600.4577.4500.46 SucdRt114.12226.20114.04227.70113.41227.90112.83224.10112.72226.70112.77227.00Year201020112012201320142015 HO*index*0.680.770.650.740.660.780.670.630.680.630.730.54 Gini331.91117.56331.93117.13332.11116.14332.10115.47331.78115.47331.80115.60 ICT*index*66.7844.2866.9744.1877.4333.6977.5733.5777.6433.5777.7033.14 ICT*acs*77.3933.9477.3633.7177.4933.8288.0933.6388.1133.6388.1333.49 ICT*use*55.5644.0155.6755.8766.0755.6266.3655.4766.6055.4766.8455.06 ICT*skl*88.8133.1088.8133.0788.8322.6388.9322.8288.9222.8288.9222.63 RdExp11.9033.3411.9933.3522.0533.5022.0033.5422.0233.5422.0833.63 lnGDPpc110.4211.60110.4711.55110.4911.49110.5311.45110.5611.45110.6011.45 ObstM220.09225.80220.55226.10221.02226.40221.50226.80221.95226.80222.43227.60 TerEdu337.76447.57338.58444.95339.92444.69441.00444.68441.80444.68442.05443.82 lnHrWk77.4500.4477.4500.4377.4400.4377.4400.4377.4300.4377.4300.42 SucdRt113.06331.40112.52331.80112.68331.20112.62226.50112.02226.50111.67117.40HO*index* Health outcome index, *Gini* Gini index, ICT*index* Information and communication technology, *acs* Access, skl Skil, *RdExp* Public expenditure on research and development, *lnGDPpc* Log of per capita gross domestic product, *ObstM* The percentage of obese population (males 15--64 years), *TerEdu* Percentage of total population (15 to 64) are tertiary educated, *lnHrWrk* Growth in actual hours of work per person per year, *SudRt* Suicide rate as total suicide per 100,000 population

Table [3](#Tab3){ref-type="table"} illustrates the findings of eq. 1 with four different models estimated by PCSE and Pooled OLS with DKSE approaches. For model 1, we used variable ICT*index;* for model 2 ICT*acs*; for model 3 ICT*use;* and for model 4 ICT*skl*. According to the PCSE approach, ICT development, income inequality, male obesity, growth in hours of work and the suicide rate have a negative impact on health outcomes in model 1. However, income inequality and the suicide rate are not significant. A one unit increase in the ICT development index reduces the health outcome index by 0.051. To further investigate the impact of ICT development, the variable was segregated into three parts. Notably, the findings of models 2, 3 and 4 indicated that although ICT*acs* and ICT*use* are inversely related to health outcome, conversely, improvement in ICT*skl* has a positive and significant effect on health. Moreover, expectedly, growth in income, technological development (R&D expenditure) and a rising percentage of tertiary educated population have a significant and positive influence on health outcomes. A 1% increase in per capita GDP increases the health outcome index by 0.328.Table 3OLS estimation results for models 1 to 4TestPooled OLS DKSEOLS PCSEDep Var. (Health outcome)Model 1Model 2Model 3Model 4Model 1Model 2Model 3Model 4ICT*index*−0.051\*\*\* (0.002)−0.027\*\* (0.007)ICT*acs*−0.033\* (0.017)−0.014\* (0.008)ICT*use*−0.039\*\*\* (0.003)−0.024\*\*\* (0.004)ICT*skl*0.063\*\*\* (0.013)0.022\*\* (0.012)Gini−0.012\*\*\* (0.009)−0.010\*\*\* (0.001)− 0.013\*\*\* (0.001)−0.004\*\*\* (0.001)− 0.002 (0.002)−0.001 (0.002)− 0.003\* (0.002)−0.000 (0.001)*RdExp*0.042\*\*\* (0.003)0.037\*\*\* (0.002)0.041\*\*\* (0.006)0.016\*\*\* (0.002)0.024\*\* (0.011)0.019\* (0.011)0.0 7\*\* (0.010)0.014 (0.011)lnGDPpc0.371\*\*\* (0.042)0.365\*\*\* (0.056)0.372\*\*\* (0.034)0.342\*\*\* (0.016)0.328\*\*\* (0.014)0.311\*\*\* (0.043)0.344\*\*\* (0.039)0.282\*\*\* (0.037)ObstM−0.007\*\*\* (0.001)− 0.009\*\*\* (0.002)− 0.005\*\*\* (0.000)− 0.011 (0.002)− 0.007\*\*\* (0.002)−0.009\*\*\* (0.002)− 0.005\*\*\* (0.002)−0.010\*\*\* (0.002)TerEdu0.004\*\*\* (0.001)0.003\*\*\* (0.001)0.005\*\*\* (0.001)0.0002 (0.001)0.002\*\* (0.001)0.001 (0.001)0.002\*\*\* (0.001)0.000 (0.001)lnHrWrk0.065 (0.026)0.033 (0.029)0.075 (0.025)−0.062\* (0.030)−0.167\* (0087)− 0.196\*\* (0.090)−0.108 (0.078)− 0.226\*\* (0.090)SudRt−0.011\*\*\* (0.002)− 0.011\*\*\* (0.001)−0.011\*\*\* (0.002)− 0.010\*\*\* (0.001)−0.002 (0.002)− 0.002 (0.002)−0.004\*\* (0.002)− 0.002 (0.002)*Obs*339339339339339339339339*N-group*3030303030303030*R-Sq*0.7085.68050.74880.69950.85940.84330.85180.8438*VIF* ~*mean*~2.232.402.072.04*CSD* ~*p-value*~0.0000.0000.0000.000*Hetero* ~*p-value*~0.0000.0000.0000.000*AutoCorr* ~*p-value*~0..170.150.190.13*JB test* ~*p-value*~0.0540.1900.1430.285Notes: Significant at 10, 5, and 1% level is defined by \*, \*\*, and \*\*\*, respectively. Standard error appear in the parentheses. Obs indicates the number of observations, VIF shows the variance inflation factor test results for multicollinearity, Hetero is the test of heteroskedasticity, Autcorr is the first-order serial correlation test \[[@CR41]\], CSD is the cross-section dependence test \[[@CR37]\] and JB test is the Jarque-Bera residual normality test \[[@CR55]\]

The results of the Pooled OLS DKSE estimation indicate similar signs for almost all the coefficients, however, there are some notable differences. The variable income inequality is not only negatively related to health outcome but is also significant. Growth in hours of work is positive but not significant and the suicide rate significantly impacts health outcomes. The only variable that has a different sign is growth in hours of work but it is not significant. For model 1, a one unit increase in income inequality and the suicide rate reduces the health outcome index by 0.012 and 0.011, respectively.

The diagnostic test results showed that the panel dataset exhibits heterogeneity and the errors are cross-sectional dependent. Serial correlation and multicollinearity were not detected and the errors are normally distributed \[[@CR55]\]. All the results were significant at 95% confidence interval.

Table [4](#Tab4){ref-type="table"} shows three alternative versions of model 1 (with Eth*dum*, PHE*dum* and EthPhe) and one alternative of model 4 (with SklPhe). As expected, all the coefficient signs and significance level are similar to the previously estimated results in Table [3](#Tab3){ref-type="table"}. However, there were minor changes in the coefficient values. The PCSE estimation results show that countries with higher ethnic diversity have lower health outcomes than countries with lower ethnic diversity and the result is significant at 1% confidence interval. In addition, countries where health expenditure is generously funded by government enjoy significantly better health outcomes. Moreover, countries with higher ethnic diversity and more generous public health expenditure are significant at 1%. Finally, countries with a higher ICT*skl* index and generous public health expenditure illustrates significantly higher health index values. The Pooled OLS DKSE approach also shows identical results.Table 4OLS estimation results with dummy and interaction variablesTestPooled OLS DKSEOLS PCSEDep var. (Health outcome)Ethnicity*dum*Public Health Exp*dum*Interatction term 1Interaction term 2Ethnicity*dum*Public Health Exp*dum*Interatction term 1Interaction term 2ICT*index*−0.057\*\*\* (0.007)−0.054\*\*\* (0.007)−0.052\*\*\* (0.007)−0.031\*\*\* (0.007)−0.030\*\*\* (0.006)− 0.31\*\*\* (0.006)ICT*skl*0.054\*\*\* (0.013)0.012\* (0.012)Gini−0.008\*\*\* (0.001)−0.008\*\*\* (0.001)− 0.007\*\*\* (0.001)−0.002 (0.001)− 0.001 (0.002)−0.0001 (0.002)− 0.0007 (0.002)−0.002 (0.002)RdExp0.040\*\*\* (0.003)0.046\*\*\* (0.004)0.031\*\*\* (0.003)0.019\*\*\* (0.002)0.025\*\* (0.011)0.031\*\*\* (0.011)0.019\*\* (0.010)0.019\*\* (0.011)lnGDPpc0.373\*\*\* (0.038)0.351\*\*\* (0.0445)0.366\*\*\* (0.039)0.330\*\*\* (0.018)0.332\*\*\* (0.039)0.314\*\*\* (0.042)0.334\*\*\* (0.038)0.265\*\*\* (0.038)ObstM−0.005\*\*\* (0.001)−0.004\*\*\* (0.001)−0.006\*\*\* (0.001)− 0.009\*\*\* (0.002)−0.005\*\* (0.002)− 0.004\*\* (0.027)−0.006\*\*\* (0.002)− 0.008\*\*\* (0.002)TerEdu0.005\*\*\* (0.001)0.004\*\*\* (0.001)0.005\*\*\* (0.001)0.000 (0.001)0.003\*\*\* 0.001)0.002\*\* (0.001)0.003\*\*\* (0.001)0.001 (0.001)lnHrWrk0.008 (0.025)0.207\*\*\* (0.028)0.118\*\*\* (0.023)0.063\*\* (0.026)−0.196 (0.085)−0.012 (0.097)− 0.061 (0.081)−0.083 (0.102)SudRt−0.009\*\*\* (0.002)−0.009\*\*\* (0.001)− 0.009\*\*\* (0.001)−0.009\*\*\* (0.001)− 0.002 (0.002)−0.002 (0.002)− 0.002 (0.001)−0.002 (0.002)Eth*dum*−0.076\*\*\* (0.011)- 0.077\*\*\* (0.020)PHE*dum*0.084\*\*\* (0.013)0.085\*\*\* (0.026)EthPhe−0.044\*\*\* (0.013)− 0.054\*\*\* (0.018)SklPhe0.008\*\*\* (0.001)0.009\*\*\* (0.003)*N-group*339339339339339339*Obs*303030303030*R-Sq*0.73750.72960.74230.85660.85590.8599Notes: Significant at 10, 5, and 1% level is defined by \*, \*\*, and \*\*\*, respectively. Standard error appear in the parentheses. Obs indicates the number of observations

Discussion {#Sec11}
==========

All estimated coefficients have the a priori anticipated signs. This confirms the validity of the health outcome index constructed in this study. Thus the method was then used investigate the key determinants of health outcome using longitudinal data while controlling for heteroskedasticity and cross-section dependence.

The results of the Pooled OLS with DKSE and PCSE estimation showed that income inequality negatively impacts health outcome for the OECD countries. However, this outcome is significant only in the Pooled OLS test. One observation is that to establish the relationship beyond doubt, maybe, time series data on a longer time scale is required \[[@CR56]\] as income inequality often affects health through other components such as harmful health behaviour \[[@CR57]\] and stress-related diseases \[[@CR58]\]. Nonetheless, this study supports the income inequality hypothesis that the more unequal the society, the higher the probability of poorer health controlling for per capita GDP. Moreover, past empirical works often produced mixed conclusions, where some researchers supported this hypothesis \[[@CR2], [@CR8], [@CR9], [@CR13], [@CR59]\] and others refuted it \[[@CR15], [@CR60]\].

In recent years, many organisations and governments have promoted the use and access of ICT in the health sector \[[@CR24]\]. However, the findings of this study suggested that increasing ICT access and use in society has an adverse effect on health outcomes. A number of previous studies have discussed the potential health benefits of ICT development, access and use. Their findings concluded that ICT improves the diffusion and communication of health information and knowledge, therefore, facilitates increased health literacy and healthier lifestyles \[[@CR26], [@CR61]\]. As a result, increasing access and use of ICT can improve a country's public health outcomes \[[@CR62], [@CR63]\]. Conversely, Punamäki et al. \[[@CR31]\] found that ICT use is connected with lower perceived health status and yet others concluded that excessive usage of ICT can cause musculoskeletal problems \[[@CR64], [@CR65]\]. Other, empirical findings also associated ICT access and use with unhealthy lifestyles and obesity \[[@CR28]\], mental health disorders such as lack of sleep and depression \[[@CR58]\] and development of health compromising-behaviour in both men and women \[[@CR66]\]. Therefore, ICT use and access have both positive and negative impacts on health. On balance, the results of this study show that the negative impacts outweigh the benefits. Not surprisingly, improvement in ICT skills (measured by the literacy rate) has been found to promote better health.

The level of ethnic diversity in most OECD countries is relatively low and the study findings show that countries with lower levels of diversity enjoy higher health status. This conclusion is consistent with previous empirical works which showed that ethnic minorities are subject to unequal access to health care services, and poorer health status \[[@CR3], [@CR4], [@CR19], [@CR64], [@CR67], [@CR68]\]. Nevertheless, the detrimental consequences of higher ethnic diversity on health outcome can be compensated in part by devoting a larger share of publicly funded health care expenditure to minority groups.

A common point of argument among researchers and policymakers is the utility of public welfare expenditure, especially in the healthcare sector. In this context, the present study provides strong evidence that where health care is generously financed by government better health outcomes are evidenced. The study results support previous findings of \[[@CR5], [@CR69]--[@CR71]\] and contradicts Self and Grabowski \[[@CR72]\] who found no such relationship in developed countries.

Finally, growth in per capita GDP, advancement in technology and increasing levels of tertiary education are pivotal components of better health outcomes in OECD countries. On the contrary, obesity plays a major role in deteriorating health status. The percentage of the obese population (male) has increased by 31% (authors own calculation) over the period 2004 to 2015, which can be a considerable cause of concern for future public health policies in these countries.

Some policy implications can be derived from the findings of this study.

First, this study examined the relationship between the development of ICT access, use and skills with overall public health outcomes. This study supports the hypothesis that overuse of ICT (by the general population) can have a negative impact on public health. Governments should create public awareness regarding the ill-effects of ICT use on health and recommend necessary actions to avoid future health problems.

Second, policymakers should recognise that having a population of greater ethnic diversity makes it more challenging to promote better health outcomes for all. This may be part of the answer to the puzzle as to why some countries could achieve better health outcomes with lower per capita expenditure on health. Populous countries such as the USA with higher levels of ethnic diversity may find it difficult to achieve the health status of Japan (a country with lower ethnic diversity). Again, policymakers should realise that improving health care access for ethnic minorities with publicly funded health care can improve the overall health of a country.

Third, understanding the effectiveness of public health expenditure in producing a healthy nation is very important. The findings of the study may be an important indication for governments especially those with comparatively lower health status. Even with the burden of rising health expenditure, governments should continue to fund the health system adequately.

Last, the issue of rising obesity in OECD countries should be urgently addressed. Governments must actively promote healthy lifestyles and discourage or control social and individual practices that cause obesity.

Conclusions {#Sec12}
===========

A number of empirical studies have examined the factors affecting health outcome in OECD countries. These studies either used cross-section data or a singular variable to measure health outcome. This study constructed a health outcome index with indicators comprising both mortality and morbidity data for 30 OECD countries over the period of 2004 to 2015. The index represented overall health status in individual countries. Development of the index gave the opportunity to investigate the impact of income inequality, ethnic diversity and ICT development on health with a panel dataset. In addition, the estimated method controlled for heterogeneity and cross-section dependence which improved the reliability of the estimated outcomes.

The study results indicate that income inequality, ethnic diversity and developments in ICT access and use have a significant negative relationship with health outcome in the OECD countries. On the other hand, improvements in ICT skills level and a higher levels of publicly funded health expenditure show statistically significant and positive associations with better health outcomes. Rising national incomes, developments in technology and higher education levels also improve a country's health outcomes. On the contrary, rising obesity amongst the adult population has serious negative consequences for health.

There are some limitations of the study. The panel dataset is restricted to only 12 years due to unavailability of data for some of the variables such as perceived health status and income inequality. However, the time frame (2004--2015) is important because during this time ICT development accelerated rapidly (Table [2](#Tab2){ref-type="table"}) and ethnic diversity increased in the OECD countries. In addition, use of more morbidity indicators would have increased the credibility of the health outcome index. The panel data has shown significant heterogeneity and cross-section dependence. The results of this study, therefore, should be interpreted with caution.

Highlights {#Sec13}
----------

This study constructed a new health outcome index using both mortality and morbidity variables.Income inequality and higher level of ethnic diversity negatively impacts health outcomes in OECD countries.Improvement in ICT access and use negatively, and ICT skills positively effects health outcomes.OECD countries with larger portion of publicly funded health care show better health outcome compared to other countries.

Additional file
===============

 {#Sec14}

Additional file 1:The country list, variable definitions and additional statistical analysis \[[@CR73]--[@CR84]\]. (DOCX 126 kb)

**Electronic supplementary material**

The online version of this article (10.1186/s12992-018-0385-2) contains supplementary material, which is available to authorized users.

The paper is part of the first author's PhD studies at the University of Southern Queensland, Australia and funded through the USQ International Stipend Research Scholarship and USQ International Fees Research Scholarship. The authors would like to thank the editor and two anonymous reviewers for their constructive comments, which helped to improve the manuscript. Special thanks to Mr. Afshar Ali for his suggestions regarding constructing the index.

Availability of data and materials {#FPar1}
==================================

The sources of the data are:i.The OECD statistics. Published by Organization for Economic Cooperation and Development database 2016. Available at <https://data.oecd.org/>.ii.'World Development Indicators' database published by the 'World Bank'. The data are publicly available from the following link of the database,iii.[*http://databank.worldbank.org/data/reports.aspx?source=world-development-indicators*](http://databank.worldbank.org/data/reports.aspx?source=world-development-indicators)iv.Central Intelligence Agency (CIA). 2017. CIA The World Factbook. USA. <https://www.cia.gov/library/publications/resources/the-world-factbook/fields/2075.html>v.International Telecommunication Union (ITU): Measuring the Information Society Report 2017. ICT Facts and Figs. 2017. CH-1211, Geneva, Switzerland. <https://www.itu.int/en/ITU-D/Statistics/Documents/publications/misr2017/MISR2017_Volume1.pdf>.

In addition, the complete structured dataset used in this study are available from the corresponding author on reasonable request.

RHR designed the study, conducted statistical analysis and prepared the initial draft of the manuscript. KA assisted in designing the study, supervised data analysis and the writing of the manuscript. JG updated the research design and reviewed the article. All authors read and approved the final manuscript.

Ethics approval and consent to participate {#FPar2}
==========================================

Not applicable

Consent for publication {#FPar3}
=======================

Not applicable

Competing interests {#FPar4}
===================

The authors declare that they have no competing interest.

Publisher's Note {#FPar5}
================

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
