Article Text

Download PDFPDF

Adaptation and validation of the Distress Scale for Mexican patients with type 2 diabetes and hypertension: a cross-sectional survey
  1. Ingrid Patricia Martinez-Vega1,
  2. Svetlana V Doubova1,
  3. Rebeca Aguirre-Hernandez2,
  4. Claudia Infante-Castañeda3
  1. 1Epidemiology and Health Services Research Unit, CMN Siglo XXI, Mexican Institute of Social Security, Mexico City, Mexico
  2. 2Departamento de Farmacología de la Facultad de Medicina, Universidad Nacional Autónoma de México, México City, Mexico
  3. 3Instituto de Investigaciones Sociales, Universidad Nacional Autónoma de México, México City, Mexico
  1. Correspondence to Dr Svetlana V Doubova; svetlana.doubova{at}gmail.com

Abstract

Objectives The aim of this study was to adapt and validate the Distress Scale for Mexican patients with type 2 diabetes and hypertension (DSDH17M).

Setting Two family medicine clinics affiliated with the Mexican Institute of Social Security.

Participants 722 patients with type 2 diabetes and/or hypertension (235 patients with diabetes, 233 patients with hypertension and 254 patients with both diseases).

Design A cross-sectional survey.

Methods The validation procedures included: (1) content validity using a group of experts, (2) construct validity from exploratory factor analysis, (3) internal consistency using Cronbach's α, (4) convergent validity between DSDH17M and anxiety and depression using the Spearman correlation coefficient, (5) discriminative validity through the Wilcoxon rank-sum test and (6) test–retest reliability using intraclass correlation coefficient.

Results The DSDH17M has 17 items and three factors explaining 67% of the total variance. Cronbach α ranged from 0.83 to 0.91 among factors. The first factor of ‘Regime-related Distress and Emotional Burden’ moderately correlated with anxiety and depression scores. Discriminative validity revealed that patients with obesity, those with stressful events and those who did not adhere to pharmacological treatment had significantly higher distress scores in all DSDH17M domains. Test–retest intraclass correlation coefficient for DSDH17M ranged from 0.92 to 0.97 among factors.

Conclusions DSDH17M is a valid and reliable tool to identify distress of patients with type 2 diabetes and hypertension.

  • Instrumental study
  • Psychometric properties
  • Distress
  • Diabetes

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • We adapted and validated the Diabetes Distress Scale (DDS17) for patients with type 2 diabetes and hypertension, considering that both are frequent chronic diseases often occurring together that have common risk factors, severe long-term complications and require several common self-care strategies—circumstances that can cause distress.

  • The study provides the evidence that the Distress Scale for Mexican patients with type 2 diabetes and hypertension (DSDH17M) is valid and reliable and describes specific problems that contribute to distress patients in these populations.

  • The primary limitation of this study is that the DSDH17M validation was performed only in a group of patients treated at two family medicine clinics of Instituto Mexicano del Seguro Social (IMSS) in Mexico City, possibly limiting generalisability of the results.

Introduction

Distress is a non-psychiatric emotional burden in patients with chronic diseases.1 ,2 Emotional distress is non-specific to a chronic disease as observed in patients with diabetes, hypertension, cancer, multiple sclerosis and HIV, among others.2 Its prevalence among patients with chronic disease varies according to disease type, severity, treatment, availability of psychological support by healthcare providers, family and community.2–6 Emotional distress may occur in response to the disease manifestation, progression, complications, impact on self-concept and need of continuous non-pharmacological and pharmacological treatment, among others.2–5 Timely detection and management of emotional distress in patients with chronic disease is associated with better self-care, quality of life and health outcomes.7 Emotional distress needs to be differentiated from depression due to the differences in underlying treatment.

Diabetes and hypertension are two of the most common chronic diseases worldwide. Furthermore, the prevalence of hypertension in persons with diabetes varies between 50% and 90%.8 Up to 75% of patients with diabetes die from cardiovascular disease attributed to hypertension.9 Yet, most research on emotional distress focuses primarily on diabetes with little concern regarding to hypertension or their co-occurrence. Several Diabetes Distress Scales (DDS) have been validated, whereas according to our knowledge the General Health Questionnaire (12GHQ) is the only scale that has been used for measuring hypertension-related distress.5 Furthermore, 12GHQ and several specific distress scales (eg, Problem Areas in Diabetes (PAID) scale) lack some dimensions of the construct (eg, patients’ feelings about their healthcare provider) that are relevant for patients with chronic diseases.10 To overcome these problems, the DDS17 was developed with the aim of measuring patients distress levels in response to the emotional burden, care regimen, interpersonal factors and physician care.10 However, this scale focuses only on patients with diabetes.

Type 2 diabetes and hypertension are frequently coexistent chronic diseases with common socially determined lifestyle risk factors (eg, stress, sedentary lifestyle, unhealthy food). When poorly controlled, both diseases can cause severe long-term complications. Their control requires adherence to common self-care strategies (eg, diet, physical exercise and pharmacological treatment) and continuous monitoring (eg, visits to the healthcare professionals, laboratory tests, self-monitoring of blood pressure and blood sugar). All the factors mentioned above can provoke distress. The recognition of these common circumstances of chronic diseases has led healthcare researchers to the development of an Integrated Chronic Care Model that focus on integrative healthcare of chronic patients.11 Also, the magnitude of chronic diseases has prompted to build a measurement framework for the burden of treatment in patients with chronic conditions.12 All these efforts acknowledge that today patients have frequent multiple chronic comorbidity that requires general tasks of treatment and self-care.

In Mexico, the prevalence of type 2 diabetes and hypertension in adults is high. The prevalence of diabetes increased from 7% in 2006 to 9.2% in 2012, whereas the prevalence of hypertension has remained at 31.5%.13 Mexican studies on emotional disorders among diabetic and hypertensive patients focus primarily on depression and have reported a high prevalence (24–48%).14–16 Although it is recognised that emotional distress is frequent in patients with diabetes and hypertension, little is known about the magnitude of this problem in Mexico; hence, its study and clinical treatment are likely to be deficient. One of the reasons for the paucity of research on emotional distress in patients with diabetes and hypertension in Mexico is the lack of a rigorously validated instrument for assessment of this condition. To our knowledge, an attempt was made to validate the PAID scale.17 ,18 However, this Mexican adoption of PAID only has partial validation that includes construct validity through exploratory factor analysis, internal consistency with Cronbach’s α and convergent validity analysis and lacks the content validity; furthermore, it has the same deficiencies regarding the dimensions of the PAID construct as the original scale.17 ,18

To close these gaps in the literature, the objective of this study was to adapt and validate a Mexican version of the distress scale for patients with type 2 diabetes and hypertension.

Methods

In 2014, we conducted a cross-sectional survey in two family medicine clinics affiliated with the Mexican Institute of Social Security (IMSS). The IMSS extends coverage to 48.8% of the Mexican population, providing healthcare for employees in the formal labour sector and their families.19 Clinics were selected according to convenience from the list of IMSS Mexico City ambulatory clinics: one from the northern area and the other from the southern area of the city. The clinics had the same services and infrastructure.

Participants

The inclusion criteria considered patients 20 years of age or older, with type 2 diabetes, hypertension or both and with at least 3 months since their diagnosis. We focused on patients with type 2 diabetes, as at IMSS patients with type 1 diabetes usually attend follow-up visits with specialist at the secondary care hospital settings, while patients with type 2 diabetes (that accounts for more than 90% of all diagnosed diabetes cases at IMSS) receive care from family doctors at the family medicine clinics.

Two nurses received a 1-week training course to supervise the fieldwork including participant identification and recruitment and questionnaire application. Nurses interviewed all available patients with type 2 diabetes and/or hypertension in the morning and evening shifts. Patients fulfilled the inclusion criteria and agreed to participate in the study after signing informed consent. Nurses also verified the patients’ height, weight, diagnoses and treatment in the electronic health records.

Study variables

The primary study variable was emotional distress. We adopted the DDS1710 for patients with diabetes and hypertension and validated it in the sample of Mexican patients. The original DDS is a 17-item scale with four domains: distress associated with emotional burden, care regimen, interpersonal factors and physician care.10 Each item is rated on a 6-point scale ranging from ‘not a problem’ to ‘a very serious problem’. The mean score was obtained for each domain as the sum of the items in each domain divided by the number of items in the domain.

To describe the characteristics of the study population and validate the distress scale for patients with diabetes and hypertension (DSDH17M), we collected information regarding patient's age; sex; education and marital, employment and nutritional status. Nutritional status was classified into four groups according to body mass index (BMI): low weight (BMI<18.5 kg/m2), normal weight (BMI 18.5–24.9 kg/m2), overweight (BMI 25.0–29.9 kg/m2) and obese (BMI≥30.0 kg/m2). A stressful event was considered if the response was affirmative to the question of whether he/she had at least one negative stressful event during the previous month. To facilitate the answer to this question we provided a list of 28 stressful events selected from the scale of Holmes and Rahe.20 Also, at the end of this list, the patient was asked if he/she had any other negative stressful event during the previous month.

We included information on the clinical history such as duration of diabetes and/or hypertension after diagnosis and presence of any diabetes and/or hypertension-related complications. The Hospital Anxiety and Depression Scale (HADS) assessed depression (seven items) and anxiety (seven items)21 using a 4-point Likert response format with each subscale ranging from 0 to 21. Scores ≥11 identify the possible presence of anxiety or depression. Adherence to pharmacological treatment was ascertained when a patient reported complying with the indicated doses and schedules for all prescribed medications during the prior month. Blood pressure control was defined when, during the medical visit, blood pressure was <130/80 mm Hg for patients with comorbidity of diabetes and hypertension and <140/90 mm Hg for those who only had hypertension.22 Glucose control was defined when the patient had glycated haemoglobin (HbA1C) <7 or in the absence of this test when their fasting plasma glucose was <130 mg/dL in the last measurement.23

Procedures

The sample size was based on the practice of ensuring a person-to-item ratio of 10:1.24 It also considered 20% of losses during the survey (<80% of the information obtained). Therefore, minimum sample size calculated was 221 patients for each diagnosis (diabetes, hypertension or both). Patients from two clinics were included to obtain the necessary sample size and greater generalisability of the study results.

Validation of the DSDH17M was performed in two stages as follows:

Stage 1: DDS17 was translated into (Mexican) Spanish by a bilingual translator and then translated back to English (back translation) by another bilingual translator. Adoption of the scale for both populations (patients with diabetes and/or hypertension) was achieved by including the phrase ‘diabetes and/or hypertension’ in each item. Then, a group of experts evaluated the content validity of the scale for both diseases. This group included two health systems researchers specialised in family medicine, one psychologist, one medical sociologist and one nurse. All members of the group had experience in the care of patients with diabetes and hypertension and validation of health assessment tools. All members were invited to assess the appropriateness of the language and whether each item of the adapted DSDH17M was representative of the concept that the instrument purports to measure. Finally, experts rated each item as 1=’not relevant’, 2=’useful but not relevant’ or 3=’relevant’. Expert responses were grouped and those who rated the items as ‘relevant’ were counted. Content validity index (CVI) for each item was calculated as (ne–N/2)/(N/2), where ‘ne’ was the number of experts who indicated the item as ‘relevant’ and ‘N’ was the total number of experts.25 Content analysis showed that all DDS17 items can be applied for both diseases since all the items are formulated in general terms and avoid treatment and self-care particularities. There were no indications that items needed to be replaced or eliminated as all items obtained a CVI >0.7. Also, a pilot pretest of the scale in 25 patients was performed to assure clarity of items from a linguistic and cultural view. Online supplementary appendix 1 presents the adopted DSDH17M questionnaire.

Supplemental material

Stage 2: Descriptive analysis served to depict participant's characteristics, calculate mean, SD, skewness, kurtosis and frequency of distribution of each item of the adopted DSDH17M questionnaire. Corrected item-total correlations were computed to determine if the behaviour of each item was consistent with the other items in the same factor of the original scale.26 Exploratory factor analysis was carried out based on the polychoric correlation matrix of the whole sample. Three and four factors were extracted using generalised least squares. In all cases, the solution was analysed without being rotated and after varimax, and oblimin rotations were applied.27 Factor loadings >0.4 were considered significant.28 Percentage of variance explained by each factor and communality of each item were computed.

We also performed separate factor analysis for patients with diabetes and hypertension and for those patients with both diseases, and carried out a multiple group confirmatory factor analysis29 to test whether the three groups of patients differed with regard to the underlying constructs. When successive invariance models (configural, weak, strong and strict) were compared, the incremental comparative fit index was <0.01, that is, the measurement model is invariant across the studied populations. Therefore, we present only the analyses for the whole sample of patients.

Cronbach's a and the average interitem correlation were used to measure the internal consistency of items associated with the same factor; a coefficients around 0.80 were considered acceptable.26 A score for each factor was calculated by summing all item responses within each factor and dividing by the total number of factor items. Convergent validity between the score factors and other measures of psychosocial well-being—anxiety and depression—were assessed by calculating Spearman's rank correlation coefficient; r<0.35 was considered to represent weak correlations, r ranging from 0.36 to 0.67 represent moderate correlations and from 0.68 to 1.0 high correlations.30 Wilcoxon rank-sum test served to evaluate the ability of the score factors to discriminate between subpopulations defined by different clinical and personal characteristics of the participants. We hypothesised that females with less education, those who reported stressful events and emotional problems, had obesity and complications of diabetes and/or hypertension, disease duration <3 years, did not adhere to treatment and did not have blood pressure and glucose control may present higher distress in all domains.4–6 ,31

Test–retest reliability was assessed by applying the DSDH17M twice with a lag time of 5 days to a subgroup of 35 participants. For each participant, two global mean scores and two mean scores for each factor were calculated and compared using the intraclass correlation coefficient (ICC); ICC >0.75 was considered excellent, 0.40–0.75 fair to good and <0.40 poor.32 All analyses were performed with R Core Team 2015 V.3.1.3 and STATA V.12.

Results

Of 856 eligible patients, 722 (84.3%) agreed to participate. Reasons for not accepting to participate were the lack of time (61.9%), no interest in responding to the survey (27.6%) and fatigue or feeling poorly (10.5%).

Table 1 describes the general characteristics, nutritional status, medical history, adherence to pharmacological treatment and disease control of study participants. The study included 722 patients, the majority being female (73%). Mean age was 58.2±10.4 years. Almost half (47%) completed only elementary school or less and were not employed. Most participants lived with a life partner. A total of 0.3% had a low weight, 16.0% had a normal weight, 36.1 were overweight and 47.6% had obesity; 32.5% had diabetes, 32.3% had hypertension and 35.2% had both diseases. Median time elapsed since diagnosis was 8 years; 61.5% had diabetes and/or hypertension complications, 12.3% reported some stressful event in the last month, 22% had anxiety and 9% had depression according to the HAD scale, whereas 54.8% attained good treatment adherence, 54.4% had blood pressure control and 46% of patients with diabetes had glucose control. Furthermore, online supplementary appendix 2 describes patients’ characteristics according to the disease group (diabetes, hypertension and both diseases).

Table 1

General characteristics, nutritional status, medical history, adherence to treatment and disease control (n=722)

Table 2 presents summary statistics of the DSDH17M. ‘Not a problem’ was the most frequent response for all items of the questionnaire. The three items in which most participants responded not having a problem were feeling that ‘the physician does not have sufficient knowledge about diabetes and/or hypertension’, ‘friends or family do not provide emotional support’ and ‘friends or family do not appreciate how difficult living with diabetes and/or hypertension can be’. The three items with the highest percentages of responses of having somewhat serious, serious or very serious problems were feeling that ‘he/she will end up with serious long-term complications regardless of what they do’ (28.7%), ‘diabetes and/or hypertension are consuming too much of their mental and physical energy’ (28.6%) and ‘he/she is not closely adhering to a good meal plan’ (28.3%).

Table 2

Descriptive statistics and item-total correlations of DSDH17M (n=722)

Corrected item-total correlations were >0.47, therefore the response to any given item was consistent with those of other items in the same dimension as in the original DDS17. SDs were >1, thus an adequate variability among responses was observed. Nevertheless, means ranged from 1.5 to 3 so the distribution of the answers had a positive skew for every item. Most responses were in the lower end of the Likert response format; specifically, more than 70% of the patients chose one of the following responses: 1=not a problem, 2=a slight problem or 3=a moderate problem.

Online supplementary appendix 3 provides information on descriptive statistics and item-total correlations of DSDH17M according to the disease group (diabetes, hypertension and both diseases).

Table 3 presents the three-factor model with a varimax rotation. This model was similar to the three-factor model with an oblimin rotation. The four-factor models (with a varimax rotation and with an oblimin rotation) did not produce a loading matrix with a simple structure and a Heywood case was reported in one case. The model in table 3 explained 67% of the total variance. Ten items were associated with factor 1 labelled as ‘regimen-related distress and emotional burden’. This factor explained 33% of variance. Items associated with factor 1 had communality values between 0.37 and 0.78. Three items defined factor 2 referred to as ‘interpersonal distress’. This factor explained 16% of the total variance and had communalities between 0.80 and 0.93. Finally, four items were related to factor 3 called ‘physician-related distress’. Factor 3 explained 18% of the total variance and communalities ranged from 0.51 to 0.88.

Table 3

Factor analysis of DSDH17M (n=722)

Table 4 presents the results of reliability and convergent validity analysis. The average interitem correlation was >0.5 for the three factors. Internal consistency as measured by Cronbach a was 0.91 for factors 1 and 2 with 95% CIs (95% CI 0.89 to 0.93) and (95% CI 0.85 to 0.96), respectively. Cronbach a for factor 3 was 0.83 (95% CI 0.78 to 0.88). For the whole set of items, Cronbach a was 0.91 (95% CI 0.89 to 0.92) and the average interitem correlation was 0.37.

Table 4

Reliability and convergent validity of DHDS17M (n=722)

The scores of the three domains of the DSDH17M had a positive and significant correlation with the anxiety and depression HADS subscales. Factor 1 of the DSDH17M was moderately correlated with anxiety (r=0.533, p<0.0001) and depression subscales (r=0.525, p<0.0001) of HADS, whereas the other two factors had a weak correlation with anxiety and depression subscales.

Discriminative validity analysis is summarised in table 5. Wilcoxon rank-sum test indicated that women, patients with obesity and those with stressful events and non-adherence to pharmacological treatment had significantly higher distress scores in all DSDH17M domains.

Table 5

Discriminative validity of DSDH17M through analysis of Wilcoxon rank-sum test (n=722)

The test–retest ICCs were 0.92 for physician-related distress, 0.96 for regimen-related distress and emotional burden, 0.97 for interpersonal distress, and for the average of all items.

Discussion

This study has two important contributions to the analysis of patients with chronic distress. It adapts DDS17 for patients with diabetes and/or hypertension and provides the evidence that the Distress Scale for Mexican patients with type 2 diabetes and hypertension (DSDH17M) is valid and reliable.

Adaptation of psychological instruments is a complex process requiring evidence of the semantic equivalence of the items, cultural fit of the instrument and adequate psychometric properties.33 The original DDS17 is a well-known scale for measuring levels of diabetes distress and areas of concern that cause the distress.10 This scale was developed in English10 and later translated, adapted and validated in several cultures, including Danish,34 Norwegian,35 Chinese36 and Chilean.37 The adaptation of a DDS17 for diabetes and hypertension is justified by the frequent coexistence of both diseases, by the presence of general tasks of treatment and self-care that these diseases share and that can cause distress. Also, content validity of the DDS17 scale was achieved by the confirmation of the relevance of all the items of the DDS17 for both diseases by the group of experts.

During the content analysis of the DSDH17M, a concern was raised about the difficulties in interpreting the results when both diseases co-occur. The conclusion was that in this case, the scale allows evaluating the level of distress and the area(s) of greater distress for both conditions and can be interpreted as the impact/burden of comorbidity on the patient. Because one cannot be sure that some dimensions of the distress caused by one chronic disease do not overlap with dimensions of the distress that other chronic disease causes in the same patient, if measured with two scales that measure very similar domains of what causes the diseases distress. In fact, one could expect to find a joint effect that not even the patient can separate.

The DSDH17M adapted for Mexican patients with diabetes and/or hypertension has acceptable psychometric characteristics. Consistent with the original DDS17, the DSDH17M preserved 17 items. This is different from the results of the DDS validation in China36 and Chile37 where researchers eliminated several items (two and four, respectively). The factorial structure of DSDH17M has three factors. The first factor combined distress related to the treatment regimen and emotional burden, whereas the other two factors (interpersonal and physician-related distress) remained separate. This factor structure is different from the original scale that has four factors but is similar to the DDS validated in China (CDDS).36 The DSDH17M explained 67% of the total variance that was greater, comparing with the Chinese version (62.4%).

Internal consistency of the DSDH17M was acceptable with Cronbach's α values ranging from 0.83 to 0.91, similar to the original, Danish, Norwegian and Chinese scales and different from the results of the DDS validation in Chile,37 where Cronbach's a ranged from 0.62 to 0.75. DSDH17M also has test–retest reliability higher than the Chinese version (0.97 vs 0.74).36 Regarding the convergent validity, the factor ‘regimen-related distress and emotional burden’ of DSDH17M was moderately correlated with anxiety and depression found during validation of Chinese and Norwegian versions of the DDS.

Discriminative validity analyses of the DSDH17M revealed higher distress for females, patients with obesity, those with stressful events and non-adherence to pharmacological treatment. Chinese researchers also found a correlation of DDS with obesity and quality of life. Contrary to results from other studies, DSDH17M score did not differ significantly according to the complications of diabetes35 or hypertension, glucose1 ,34–37 or blood pressure control.

Multiple group confirmatory factor analysis indicated that the measurement model is invariant across studied populations. This finding is congruent with previous research identifying the presence of distress in patients with diabetes or hypertension. It also highlights the common problems related to distress in these patients. Congruent with the findings in Danish patients with type 1 diabetes, we found that for Mexican patients with diabetes 2 and/or hypertension the least frequent concerns were those related to physician-related distress. A high percentage of patients in both cultures reported that there was no problem regarding ‘the physician does not have sufficient information about diabetes and/or hypertension’. The items in which most participants responded having serious problems were from the domain of regimen-related distress and emotional burdens: ‘feeling that they will end up with serious long-term complications regardless of what they do’ and that ‘diabetes and/or hypertension are consuming too much of their mental and physical energy’ were major concerns followed by ‘feeling that they are not closely adhering to a good meal plan’. These problems were similar to those reported by the Danish study. These findings highlight the need for comprehensive team-based care as proposed by the Chronic Care Model where the physician, nurse, psychologist, nutritionist and other health providers work collaboratively to meet the needs of patients with chronic illnesses and build their practical skills and confidence in self-care and improve health outcomes.11 ,38

The primary limitation of this study is that the DSDH17M validation was performed only in a group of patients treated at two family medicine clinics of IMSS in Mexico City, possibly limiting generalisability of the results. Therefore, we recognise the need to validate DSDH17M in other settings in Mexico to corroborate the factor structure of this scale. Furthermore, the DSDH17M was validated only in a sample of patients with hypertension and type 2 diabetes; however, the original DDS17 was validated for patients with type 2 and type 1 diabetes. Thus, further validation of this scale in patients with type 1 diabetes in the Mexican context is advisable. Additionally, it would be worth to assess whether the instrument has sensitivity for detecting changes in distress over time.

In conclusion, DSDH17M is a valid and reliable tool to identify distress in patients with type 2 diabetes and/or hypertension and describe specific problems that contribute to patients with distress in these populations.

Acknowledgments

The authors would like to thank medical directors of the IMSS Family Medicine Clinics Number 20 and 15. They are also grateful to nurse Leticia Rodriguez and Ana Jesica Maya García, their participation was invaluable to the successful completion of their field work.

References

Footnotes

  • Contributors IPM-V and SVD conceptualised the study and conducted the literature review. IPM-V, SVD and CI-C participated in the process of content validity and in the field work. RA-H, IPM-V and SVD performed statistical analysis. IPM-V and SVD wrote the manuscript. RA-H and CI-C critically reviewed the manuscript for significant intellectual content. All authors approved the final manuscript.

  • Funding This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Patient consent Obtained.

  • Ethics approval IMSS National Research and Ethics Committee approved the project. Authorisation number R-2013-785-034.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement This article has the supplementary files (see online supplementary appendixes 1–3). The dataset is available on request (svetlana.doubova@gmail.com).