Country-Level Factors Associated With Nurse Salaries: Empirical Evidence From Organisation for Economic Co-operation and Development Countries and Taiwan

Introduction

The current global shortage of health professionals is particularly severe with regard to nurses. According to World Health Organization (WHO, 2020) estimates, the global shortage of nurses was approximately 5.9 million in 2018 and was expected to be 5.7 million in 2030. A study of 23 member states of the Organisation for Economic Co-operation and Development (OECD) projected that these countries would have a shortage of 2.5 million nurses by 2030 (Scheffler & Arnold, 2019). Bridging this shortfall is both increasingly important because of global population growth and rapid population aging and critical to achieving health-related Sustainable Development Goals (WHO, 2020).

Nursing staff shortages are associated with adverse patient outcomes, with lower nursing staffing in hospitals contributing to higher hospital mortality (Aiken et al., 2018; Griffiths et al., 2019). A recent systematic review that evaluated 18 studies on the association between nurse staffing and missed care found that 14 of the 18 reported an association between lower nurse staffing and more reports of missed care (Griffiths et al., 2018). Furthermore, an umbrella review summarized the findings of 15 previous literature review studies and concluded there is strong evidence of an association between nurse staffing levels and length of stay, patient satisfaction, quality of nurse-delivered care, and likelihood of readmission (Blume et al., 2021). Meanwhile, a cross-country study found that higher nurse densities are associated with higher life expectancy (Amiri & Solankallio-Vahteri, 2019). Similarly, an analysis of data for 35 OECD countries reported that, in the long term, a 1% increase in nurse density is associated with 0.98%, 0.97%, and 0.96% reductions, respectively, in infant mortality, neonatal mortality, and perinatal mortality (Amiri et al., 2020). Finally, in addition to health outcomes for patients and populations, inadequate nurse staffing has been linked to negative outcomes for nurses. Lower rates of nursing staffing increase the risks of burnout, job dissatisfaction, and intention to leave (Dall'Ora et al., 2020; Shin et al., 2018).

Adequate salary significantly impacts nurses' motivation, performance, productivity, job satisfaction, and retention (McHugh & Ma, 2014; Muthmainnah & Hariyati, 2018). For health professionals, low pay reduces job satisfaction (Halcomb et al., 2018), and low rates of satisfaction negatively impact retention (Bimpong et al., 2020). In contrast, higher salaries and benefits for nurse practitioners have been associated with longer employment durations (Hagan & Curtis, 2018). However, raising salaries alone may have a limited effect on retention (Bimpong et al., 2020).

For nurses who are new to the profession, dissatisfaction with salary is the main factor cited in those considering a career change. Young nurses weigh their salary against the demands and responsibilities of their work (Flinkman et al., 2008). Attrition among young nurses reduces the number of nurses available to replace nurses who retire and those who leave the workforce before retirement age (Flinkman et al., 2008). In addition, the existence of salary compression (i.e., the lack of significant salary differentiation throughout one's career) in nursing can make the career unattractive to recruits and, among experienced nurses, reinforce dissatisfaction and intention to leave (Greipp, 2003).

International migration of the nursing workforce is another issue exacerbating the nursing shortage in certain regions, as nurses in lower-wage countries are enticed to seek nursing employment in higher-wage countries. Differences in nursing wage levels between countries are one of the main drivers of such migration (Awases et al., 2004; Dovlo, 2007). Nurses in developing countries leave their home countries to search for better professional development, salaries, working conditions, and living standards in developed countries that are experiencing nursing workforce shortages (Kline, 2003). The migration of health professionals particularly impacts developing countries with more fragile healthcare systems. For those who remain in their home countries, international emigration increases the workload and risk of overwork, making retention of healthcare professionals more challenging (Tomblin Murphy et al., 2016).

Over the past decade in Taiwan, it has been noted that the salary of nurses is low compared with those of physicians and other medical and nonmedical professionals (Fan & Chen, 2009). Moreover, amidst the 2019 coronavirus pandemic, high job stress and stagnant growth in nurses' wages were acknowledged by national legislators as essential factors contributing to nursing turnover (Huang, 2021). To determine whether and to what extent nurses in Taiwan are undervalued in terms of their salaries, cross-country comparisons may provide a useful reference. Few studies have adopted a cross-country approach to examine the determinants of nurse salary (NS) across different healthcare systems. Moreover, even fewer have examined this issue quantitatively.

This study was developed to examine country-level indicators associated with NS using OECD country data and to compare NS in Taiwan with those in OECD countries. An explanatory model for NS was developed based on OECD data and indicators commonly used in the literature to compare healthcare systems. The predictions of the model were then compared with actual NS rates in Taiwan.

Methods

A variable-oriented approach was used to identify health system indicators significantly associated with NS and to compare the NS level in Taiwan with those of counterpart OECD countries. The variable-oriented approach is mainly used to assess relationships between case (e.g., country) variables in large, “observational” samples with the goal usually being to specify general patterns that apply to the sample as a whole to make inferences and predictions (Ragin, 1999). The application of the variable-oriented approach to cross-country comparisons of health systems may be attributed to research in health economics (Cacace et al., 2013). This approach involves selecting data such as databases, comparators, and variables and using a quantitative method to identify links between variables. In the following subsections, we introduce the selection of data and variables and the statistical analysis methods used to build NS prediction models feasible for comparing NS levels between OECD countries and Taiwan.

Data

Data were extracted for this study from the OECD.Stat database (https://stats.oecd.org/) for the period of 2009–2018 for the following 28 OECD countries: Australia, Belgium, Canada, Czech Republic, Denmark, Estonia, Finland, France, Germany, Greece, Hungary, Iceland, Ireland, Israel, Italy, Japan, Korea, Mexico, Netherlands, New Zealand, Norway, Poland, Slovak Republic, Slovenia, Spain, Turkey, the United Kingdom, and the United States. Ten of the 38 OECD member countries were not included in this study because of missing data on crucial study variables. Most of the 28 countries selected were high-income countries and considered to be at an equivalent or higher level of development than Taiwan.

We collected indicators from the OECD.Stat database covering the themes of health, demography and population, national accounts, and social protection and well-being. These four themes contain most of the indicators related to the study topic. Although data collection methods and definitions used vary across these countries, the database provides high-quality data for comparing healthcare systems between countries (Reibling et al., 2019). The relevant statistics and indicators for Taiwan were obtained from the official websites of the Directorate General of Budget, Accounting and Statistics of the Executive Yuan, the Ministry of Health and Welfare, and the Ministry of Labor, Taiwan.

OECD.Stat was used as the primary data source to ensure consistency in terms of variable definitions, data collection methods, and data quality. There were missing values for some indicators in the OECD.Stat database. To obtain the data for these values, we consulted the Health at a Glance report, which is published every 2 years by the OECD (https://www.oecd.org/health/health-at-a-glance/). Additional missing data for the OECD countries in Europe were also obtained from Eurostat (https://ec.europa.eu/eurostat/web/health/overview) and from the official WHO website (https://www.who.int/data/gho) or official government statistics for non-European OECD countries. OECD.Stat and Eurostat cover a wide range of themes and indicators at the country level. We have compared them based on data of the common member countries and found consistency between them. The WHO database focuses on health-related themes. In this study, only a few missing values were replaced with the WHO data. After consulting the abovementioned sources, values for data still unavailable were imputed using the linear interpolation method and multivariate imputation using chained equations (Van Buuren, 2018).

The data in the abovementioned, publicly available databases were aggregated and cannot be used to identify individuals. Analyses based on these data do not involve human subject research.

Dependent Variable

The average annual NS for the OECD countries was treated as the dependent variable (DV) in this study. This salary was calculated as the average total annual income, including social insurance and income taxes. The calculation also included nonrecurring payments such as bonuses and nightshift and overtime compensation. The data included nurses who were public or private employees. However, in countries with predominantly publicly funded healthcare systems such as the United Kingdom, New Zealand, and Denmark, the data included nurses working in public healthcare institutions only. For most countries, the data only covered nurses working full-time in hospitals, although, for some countries and for certain years, the data also included nurses working in other settings and part-time nurses. Furthermore, for some countries, the data included only registered nurses, whereas others also included data on lower-level nurses. We did not further adjust the data for the above differences because of a lack of breakdown information. The average NS for each country was converted into U.S. dollars (USD) and adjusted for purchasing power parity (PPP).

The average NS level in Taiwan was calculated by consulting the Ministry of Labor's “Survey on Earnings by Occupation” (https://pswst.mol.gov.tw/psdn/). This survey is conducted each year in August. Thus, July was used as the reference period for each year's data. We extracted the data for nurses working in the healthcare system. Nurses' earnings comprised regular and nonregular earnings. We multiplied the monthly earnings by 12 to estimate annual salary and added 1.5 months of regular earnings to represent a year-end bonus. The annual NS level for Taiwan was also converted into USD and adjusted for PPP.

Independent Variables

To identify the indicators associated with NS, published research on comparative health systems and the health workforce was consulted. Salary level is generally determined by the supply of and demand for workers in the market, and the OECD proposed a general framework for assessing the main factors affecting the supply and demand of health workers (Ono et al., 2013). According to this framework, demand for health workers is derived from the demand for healthcare, affected by demographics, morbidity and associated health risks, the purchasing power for healthcare (gross domestic product [GDP] and health financing), and healthcare utilization, whereas supply of health workers is determined by the stock, inflow, and outflow of human resources for health (Ono et al., 2013), which are impacted by economic inequality (Squires et al., 2016). As the health system transforms inputs such as human resources, capital, and technology into outputs (i.e., population health; Frogner et al., 2015), the scope of the variable selection was extended in this study to encompass population health outcomes. Considering the relevance to NS based on the literature and theoretical grounds and the availability and completeness of data across countries, we selected 22 indicators from the OECD.Stat database and grouped them into seven dimensions, including demographics, socioeconomic status, health behaviors and risks, healthcare resources/input, health financing, healthcare utilization, population health outcomes, and economic inequality.

Demographics collected included the proportion of the population aged 65 years and over and the crude birth rate. Meanwhile, the socioeconomic context was represented by GDP per capita and the proportion of the population aged 15–64 years with a tertiary education. Health behaviors and risks comprised the prevalence of smoking, per capita alcohol consumption, and the prevalence of obesity.

Healthcare resources/inputs were categorized into health professional human resources, capital inputs, and the use of medical technology. Physician density and nurse density represented health professional human resources. Hospital bed density represented the capital inputs of a healthcare system. Computerized tomography scanners per million population represented the use of medical technology. Health financing was measured using health expenditure as a percentage of GDP and household out-of-pocket expenditures as a percentage of total health expenditures. The former represented the overall size of the healthcare market, whereas the latter represented the public/private-sector balance of the healthcare system.

This study took the number of doctor consultations per capita, the number of hospital discharges, and the average length of stay in hospital as indicators of health service utilization. The three indicators used to represent population health outcomes included the crude death rate, infant mortality rate, and life expectancy at birth. In addition, the two indicators used to represent healthcare quality were the 30-day in-hospital mortality rate for acute myocardial infarction (AMI) after hospital admission (indicating acute care quality in hospitals) and the hypertension hospital admission rate (representing the rate of avoidable hospital admissions and indicating primary care quality).

Economic inequality was indicated by the Gini coefficient and the gender wage gap. The Gini coefficient, ranging between 1 and 0, measures the degree to which incomes are concentrated or equally distributed across members of society, with values closer to 1 indicating incomes are more unequally distributed. The gender wage gap was defined as the difference between the median incomes of men and women relative to the median income of men.

The definitions and dimensions of the variables are summarized in Table 1. The means and standard deviations of these variables for the OECD countries and Taiwan are reported.

Table 1. - Description and Summary Statistics of Study Variables Across the Selected OECD Countries and Years, 2009–2018 Dimension and Variable Name Variable Description, Definition, and Relevance to Demand/Supply of Healthcare and Nurses OECD Countries
(N = 28; n = 280) Taiwan
(N = 1; n = 10) Mean SD Mean SD Dependent variable  NS Annual nurse salary, 1,000 USD, PPP 42.22 14.86 38.06 3.19 Demographic  AGE65 Proportion of the population aged 65+ years. An indicator of population aging indicating an increased demand for healthcare. 9.08 2.35 12.11 1.40  CRBIRTHR Crude birth rate, annual number of live births per 1,000 population. Infants and children have higher healthcare needs compared with adults. 11.72 3.18 8.53 0.75 Socioeconomic context  GDP GDP per capita, 1,000 USD, PPP. It indicates the purchasing power for services provided by health professionals. 38.57 11.76 43.78 5.05  HIEDU Proportion of the population aged 15–64 years with tertiary education. Higher education helps to produce health professionals. 33.97 10.47 40.60 3.32 Health behaviors and risks  SMOKER Proportion of population aged 15+ years who are daily smokers. A higher prevalence rate results in increased demand for healthcare. 19.26 5.02 16.93 2.32  ALCHCONS Liters of pure alcohol per capita of aged 15+ years. A greater consumption results in increased demand for healthcare. 8.74 2.57 3.19 0.47  OBESEP Self-reported, proportion of population aged 15+ years with a body mass index of 30 or over. A higher prevalence rate results in increased demand for healthcare. 16.68 5.58 15.59 1.22 Healthcare resources/input  NURSENO Nurse density, the number of practicing nurses providing services directly to patients, per 1,000 population. Increased nurse workforce supply may reduce nurses' salaries. 8.75 3.96 6.15 0.54  DOCTORNO Physician density, the number of practicing doctors providing services for individual patients, per 1,000 population. It possibly increases demand for nursing staff when there is a shortage of physicians. 3.21 0.86 1.82 0.12  HBED Hospital beds regularly maintained and staffed and immediately available for the care of admitted patients, per 1,000 population. It may increase the need for nursing staff according to medical regulations. 4.78 2.66 4.21 0.06  CTS Computed tomography scanners per million population. The use of medical technology represents the efficiency of healthcare. 24.56 19.90 15.43 1.15 Health financing  EXPHEAL Health expenditure as a percentage of GDP. The purchasing power for healthcare. 8.87 2.30 6.38 0.16  OOPP Household out-of-pocket expenditure as a percentage of health expenditure. Increased out-of-pocket may decrease the demand for healthcare. 18.97 7.62 34.28 0.76 Healthcare utilization  LOS Average length of stay, calculated by dividing the number of bed days by the number of discharges during the year. Intermediate output of nurses' labor inputs; it may reflect quality of care. 6.44 2.37 8.94 0.23  CONSULTD Number of physician consultations per capita. Intermediate output of nurses' labor inputs. 6.94 2.99 12.47 0.42  HDISCHAR Inpatient discharge, release of a patient who was formally admitted into a hospital for treatment and/or care and who stayed for a minimum of one night, per 1,000 people. Intermediate output of nurses' labor inputs. 151.67 41.40 134.43 3.88 Health outcomes  CRDTHR Crude death rate, annual number of deaths per 1,000 population. Final output of nurses' labor inputs. 8.61 2.05 6.83 0.41  INFMORTR Infant mortality, per 1,000 live births. Final output of nurses' labor inputs. 2.74 1.65 2.45 0.16  LE Life expectancy at birth. Final output of nurses' labor inputs. 80.17 2.49 79.80 0.57  AMI In-hospital mortality rates within 30 days after admission for acute myocardial infarction (AMI), per 100 patients. The quality indicator for acute care of hospitals. 7.42 a 4.55 8.00 c  ADMRHYPT Hypertension hospital admission, annual number of nonmaternal/nonneonatal hospital admissions with a principal diagnosis code of hypertension, per 100,000 population. The quality indicator for primary care. 83.19 b 84.50 9.40 d Economic inequality  GINI_PC The Gini coefficient. Economic inequalities can increase financial barriers to healthcare access and reduce demand for healthcare. 31.18 5.15 33.88 0.31  GWAGEGAP Gender wage gap, difference between the median income of men and women relative to the median income of men, 1,000 USD, PPP. Nurses are more likely to be female, and the lower gender wage gap means nurses are better paid. 13.76 7.34 15.88 1.23

Note. N = number of countries; n = number of observations; OECD = Organisation for Economic Co-operation and Development; USD = U.S. dollars; PPP = purchasing power parity; GDP = gross domestic product.

aBecause of missing data for some countries and periods, the number of countries is N = 27 and the number of observations is n = 243. bBecause of missing data for some countries and periods, the number of countries is N = 26 and the number of observations is n = 229. cData for year 2013 only. dData for year 2015 only.


Statistical Analysis

In this study, both the DV and independent variables (IVs) were repeatedly measured across multiple countries, and panel regression was used for the analysis. In this study, the random-effects model (REM) was used to consider both intercountry and intracountry (between and within) variations, whereas the fixed-effects model (FEM) was used to explore intracountry associations between the DV and IVs, controlling for fixed intercountry differences. In single-regression analysis, the REM and FEM were used to examine the relationship between individual IVs and the DV without adjusting for the other IVs. For the multivariate regression analysis, to avoid multicollinearity among the IVs and the overfitting problem, the backward elimination method was used to reduce the set of IVs. The Akaike information criterion (AIC) was used to measure the goodness of fit of the models with different sets of IVs (Yum, 2022), with lower AIC values indicating better goodness of fit. As the sample size was relatively small, a corrected AIC (AICc) was used in the analysis. The preceding analysis used complete data to select predictors related to nurses' remuneration. This analysis did not consider AMI in-hospital mortality or hypertension admission rates because of missing data for some countries or periods. However, using the incomplete data, we added these two variables to the models developed in the previous step and examined the relationships between the two quality-of-care indicators and NS level, respectively.

The 28 OECD countries included in this study consisted of two upper-middle-income countries (Mexico and Turkey) and 26 high-income countries, as defined using the economic classification method of the World Bank. We further conducted a sensitivity analysis by restricting the data to the 26 high-income countries to select variables and build multivariable regressions to compare the results with those of the 28 OECD countries.

Stata 14.0 (StataCorp LLC, College Station, TX, USA) was used to develop the REM and FEM. The AICc values for all models were generated using the “MuMIn” package (Version 1.46.0) in R (Version 4.1.2), available under the GNU General Public License V2.0 from the Comprehensive R Archive Network at https://cran.r-project.org/web/packages/MuMIn/index.html.

Ethical Considerations

The current study used publicly available databases only. The aggregated nature of these databases and the anonymization of all subjects mean that human subjects were not directly involved. Analysis of deidentified public data is not considered human subject research and thus does not require ethical approval by an institutional review board.

Results

Annual NS trends for the OECD countries and Taiwan for the period of 2009–2018 are shown in Figure 1. The annual NS level in the OECD countries rose from an average of 37,872 USD (PPP) to 47,886 USD (PPP) over this period, representing an increase of 26%. This level in Taiwan rose over the same period from an average of 33,030 USD (PPP) to 43,920 USD (PPP), representing an increase of 33%. By region, North America had the highest average NS level, followed by the Asia-Pacific region, with Europe reporting the lowest level. During this period, the average NS level in Taiwan was lower than the average for the Asia-Pacific region.

F1Figure 1.:

Annual Nurse Salary in 28 OECD Member Countries and Taiwan, 2009–2018Note. OECD = Organisation for Economic Co-operation and Development; USD, PPP = U.S. dollars (USD) corrected for purchasing power parity (PPP). The gray lines represent the annual salary of nurses for each of the 28 OECD member countries from 2009 to 2018. The filled black circles indicate the average nurse salary of the 28 OECD countries. The dotted lines marked by pink diamonds, green circles, and blue squares represent the average nurse salary of North American countries, Asia-Pacific countries, and European countries, respectively. The orange triangle indicates the salary of nurses in Taiwan.

The results of the simple regression model are shown in Table 2. In both the REM and FEM, NS was found to be positively associated with GDP per capita, percentage of high educational attainment, life expectancy at birth, nurse density, physician density, number of physician consultations, number of computerized tomography scanners per million population, and prevalence of obesity (p < .05) and to be negatively associated with crude birth rate, infant mortality, bed density, length of stay, household out-of-pocket expenditure as a percentage of health expenditure, prevalence of smoking, and gender wage gap (p < .05). Of the examined variables, only hospital discharges and the Gini coefficient were found to be not associated with NS in either model (p > .05).

Table 2. - Results of the Simple Regression Models on the Annual Salary of Nurses of OECD Countries Variable Random-Effects Model Fixed-Effects Model β t p β t p GDP 0.58 17.00 .002 0.57 16.17 <.001 AGE65 1.54 3.50 .036 0.87 1.78 .076 CRBIRTHR −2.43 −8.32 .007 −2.67 −8.77 <.001 HIEDU 1.04 15.73 .002 1.07 15.30 <.001 SMOKER −1.26 −11.29 .004 −1.25 −10.94 <.001 ALCHCONS −1.18 −2.11 .085 −1.62 −2.48 .014 OBESEP 1.98 12.47 .003 2.11 12.84 <.001 NURSENO 1.62 5.92 .014 1.51 4.89 <.001 DOCTORNO 9.59 9.84 .005 10.28 10.17 <.001 HBED −1.94 −3.60 .035 −2.21 −3.45 <.001 CTS 0.69 10.34 .005 0.79 10.87 <.001 EXPHEAL 1.65 3.41 .038 0.72 1.30 .196 OOPP −0.88 −5.03 .019 −0.88 −4.23 <.001 CONSULTD 1.15 3.11 .045 1.55 3.90 <.001 HDISCHAR −0.02 −0.73 .271 −0.01 −0.29 .769 LOS −1.30 −3.01 .048 −1.47 −3.18 .002 CRDTHR 2.04 2.70 .057 3.47 4.01 <.001 INFMORTR −5.93 −9.84 .005 −6.46 −10.16 <.001 LE 4.02 14.19 .003 4.09 13.72 <.001 GINI_PC −0.35 −1.48 .139 −0.43 −1.62 .106 GWAGEGAP −0.47 −3.33 .040 −0.53 −3.50 <.001

Note. Number of countries = 28, total number of observations = 280. β = coefficients of regressions; OECD = Organisation for Economic Co-operation and Development; GDP = gross domestic product per capita (expressed in 1,000 U.S. dollars, adjusted to the purchasing power parity, [PPP]); AGE65 = the proportion of the population aged 65 years and over (%); CRBIRTHR = crude birth rate (per 1,000 population); HIEDU = the proportion of the population aged 15–64 years who had tertiary education (%); SMOKER = the prevalence of smoking (%); ALCHCONS = alcohol consumption per capita (liters); OBESEP = the prevalence of obesity, self-reported (%); NURSENO = nurse density (per 1,000 population); DOCTORNO = physician density (per 1,000 population); HBED = hospital bed density (per 1,000 population); CTS = the number of computerized tomography scanners per million population; EXPHEAL = health expenditure as a percentage of GDP (%); OOPP = household out-of-pocket expenditure as a percentage of health expenditure (%); CONSULTD = the number of physician consultations per person; HDISCHAR = the number of hospital discharges (per 1,000 population); LOS = average length of stay in hospitals (days); CRDTHR = crude death rate (per 1,000 population); INFMORTR = infant mortality rate (per 1,000 live births); LE = life expectancy at birth (years); GINI_PC = Gini coefficient (%); GWAGEGAP = gender wage gap (expressed in 1,000 U.S. dollars, adjusted to PPP).

The results of the multivariable regression with backward elimination are shown in Table 3. In the REM, higher NS in the OECD countries was shown to be positively associated with GDP per capita (0.49, 95% confidence interval [CI] [0.41, 0.56]), proportion of population aged 65 years and over (2.72, 95% CI [2.17, 3.26]), crude birth rate (1.02, 95% CI [0.56, 1.49]), number of computerized tomography scanners per million population (0.26, 95% CI [0.17, 0.35]), alcohol consumption per person per year (0.94, 95% CI [0.26, 1.61]), and prevalence of obesity (0.64, 95% CI [0.40, 0.89]) and to be inversely associated with infant mortality (−3.13, 95% CI [−3.94, −2.32]), bed density (−0.99, 95% CI [−1.72, −0.25]), hospital discharges (−0.08, 95% CI [−0.11, −0.05]), household out-of-pocket expenditure as a percentage of health expenditure (−0.34, 95% CI [−0.56, −0.11]), and the Gini coefficient (−0.25, 95% CI [−0.50, −0.01]). Although nurse density was included in the REM to enhance the model's goodness of fit (as determined by the AICc), it was only weakly associated with NS (p = .054). For the FEM, the selected IVs' effects on NS were in the same direction and of a similar magnitude as those in the REM.

Table 3. - Results of the Multivariable Regression Models on the Annual Salary of Nurses of OECD Countries Using the Backward Elimination Method Variable Random-Effects Model

留言 (0)

沒有登入
gif