It is well documented that the US food supply and the dietary intakes of the US population do not adhere to the federal government’s dietary recommendations for a healthy diet( Reference Wilson, Reedy and Krebs-Smith1, Reference Miller, Reedy and Kirkpatrick2). Compared with the Dietary Guidelines for Americans 2010, the average intakes of fruits, vegetables, whole grains, dairy and legumes are 20–80 % below the recommended levels while the amounts of sodium, added sugar, solid fats and refined grains in the diet are excessive(Reference Wilson, Reedy and Krebs-Smith1). Given that the majority of food dollars in the USA are spent on groceries (57·3 %), household food purchases may significantly influence diet quality(3).
The types of foods available in the home are strongly correlated with individual dietary intake and diet quality among adults and youth(Reference Appelhans, French and Tangney4–Reference Neumark-Sztainer, Wall and Perry12). Household grocery purchases may therefore serve as a proxy to study the home food environment and household members’ food consumption and diet quality. As food purchasing behaviours are modifiable by interventions(Reference Liberato, Bailie and Brimblecombe13, Reference Adam and Jensen14), grocery purchases may be an opportune target for public health efforts to enhance home food environments and improve the diet quality of all household members. In particular, some research groups have hypothesized that reducing the proportion of ultra-processed foods in the diet may be a simple and effective strategy to improve diet quality(Reference Ludwig15–Reference Monteiro, Moubarac and Cannon17).
Industrially processed foods and beverages constituted over 75 % of energy purchased by US households in 2012(Reference Poti, Mendez and Ng18). The role of industrially processed foods in a healthy diet is controversial(Reference Ludwig15–Reference Monteiro, Moubarac and Cannon17, Reference Weaver, Dwyer and Fulgoni19, Reference Dwyer, Fulgoni and Clemens20). Researchers have argued that so-called ‘ultra-processed foods’, defined as industrial formulations made with no or minimal whole foods, are inherently nutritionally unbalanced as these are generally high in total energy, sugar, sodium and/or saturated fat, while low in fibre, essential micronutrients and phytochemicals(Reference Ludwig15–Reference Monteiro, Moubarac and Cannon17). Conversely, others disagree that food processing level is useful predictor of nutritional quality and maintain that ultra-processed and processed foods are nutritionally important in the USA(Reference Weaver, Dwyer and Fulgoni19, Reference Dwyer, Fulgoni and Clemens20).
National and international studies support that diets higher in ultra-processed foods are of poorer nutritional quality, lower in fibre, protein and several essential micronutrients, and higher in saturated fat, carbohydrates and added sugars(Reference Martinez Steele, Popkin and Swinburn21–Reference Cediel, Reyes and da Costa Louzada27). Nevertheless, industrially processed foods provide a significant proportion of several essential micronutrients in the US diet(Reference Eicher-Miller, Fulgoni and Keast28). Foods of all processing levels have also been found to contribute to both ‘nutrients to increase’ and ‘food components to limit/reduce’ according to the Dietary Guidelines for Americans 2010(Reference Eicher-Miller, Fulgoni and Keast28).
To date, the contribution of ultra-processed foods to US household grocery purchases and its association with diet quality has not been empirically determined. It is unclear if interventions aimed at improving the healthfulness of food purchases should emphasize processing level or other metrics of diet quality, such as the Healthy Eating Index-2015 (HEI-2015). The current study therefore aims to: (i) describe the processing level defined by the NOVA framework and diet quality of US households’ grocery purchases for home consumption, as measured by NOVA and the HEI-2015, respectively; and (ii) examine the associations between food processing level and the diet quality of household grocery purchases.
Methods
Data source
The present study uses data from the US Department of Agriculture’s National Household Food Acquisition and Purchase Survey (FoodAPS) collected between April 2012 and January 2013. FoodAPS is a cross-sectional survey designed to collect comprehensive data regarding household food purchases and factors that may impact food purchase decisions from a nationally representative sample of US households(29). A multistage sample design was used to obtain the FoodAPS sample of households(29).
A total of 14 317 individuals from 4826 households participated in FoodAPS(29). The main food shopper or meal planner of each household was selected as the primary respondent for the survey and provided information regarding household demographic characteristics, income, employment status and food security status(29). All household members were asked to document detailed information regarding all food acquisitions over a 7 d period, regardless if foods were consumed or not(29). Each purchasing event was classified into two mutually exclusive groups, as: (i) grocery purchases, i.e. all foods acquired for at-home consumption; and (ii) away-from-home purchases, i.e. all foods acquired and consumed while outside the home(29). In cases where information on the quantities of foods acquired were missing, the Economic Research Service (ERS) imputed estimated quantities based on information about food items, the stores from which they were acquired and household characteristics(Reference Mancino, Todd and Scharadin30). The methodology for imputing missing quantities is described in detail elsewhere(Reference Mancino, Todd and Scharadin30).
Analytical sample
All households with data regarding grocery purchases were eligible to be included in the current study (n 4367). Households reporting grocery purchases <6 items (<10th percentile, n 375) or >150 items (>99·5th percentile, n 24) were excluded as purchases were considered unlikely to be representative of usual grocery purchases. Additionally, we excluded six households purchasing only unidentified items and one household identified as an outlier, as total energy purchased was >15 times greater compared with the household purchasing the second greatest amount of energy (11 367 204 v. 743 279 kJ (2 716 827 v. 177 648 kcal)). A total of 3961 households were included in the final analytical sample. A flowchart of the creation of the analytical sample is presented in Fig. 1.
Determination of level of processing of purchased items
All recorded food items were classified based on the extent and purpose of the industrial food processing applied, in accordance with the NOVA food classification(Reference Moubarac, Parra and Cannon31). NOVA distinguishes four mutually exclusive levels of processing: (i) ‘unprocessed or minimally processed foods’ including fresh, dry or frozen plant and animal foods; (ii) ‘processed culinary ingredients’ including table sugar, oils, fats, salt and other constituents used in kitchens to make culinary preparations; (iii) ‘processed foods’ including foods such as canned fish and vegetables, simple breads and artisanal cheeses which are manufactured by adding processed culinary ingredients to unprocessed or minimally processed foods; and (iv) ‘ultra-processed food and drink products’, which are industrial formulations made with no or minimal whole foods and produced with substances extracted from foods or synthesized in laboratories such as dyes, flavourings and preservatives, often using processing techniques with no domestic equivalent(Reference Moubarac, Parra and Cannon31). All items were further classified into thirty-four mutually exclusive subgroups. The online supplementary material, Supplemental Table S1, describes each NOVA processing level and lists the subgroups included in each level.
We classified foods by considering: (i) the US Department of Agriculture main food description and additional food description of food codes; and (ii) the ERS food group. The ERS food group classification captures level of processing and convenience of foods and provides information regarding the exact form in which each item was acquired (e.g. fresh, frozen, canned or further processed into a prepared dish)(32). ERS food group assignation is based on both the item description provided by the participant and the item information in the Information Resources, Inc. database(32). We used the ERS food groups to identify processing level of food items that could be bought in multiple forms and processed to various degrees. The ERS food groups were primarily used to distinguish between 100 % fruit and vegetable juices and juices with added sugar; to differentiate fresh/frozen from canned vegetables, fruits, legumes, meat and fish; and to discriminate between store-prepared ready-to-eat meals and industrially manufactured frozen, canned and packaged meals. In cases of conflicting information regarding processing level, we gave priority to the ERS food group over food code information, as food codes that accurately reflect all properties of the food item were not available for all reported items. Two authors independently reviewed the classification of each item. Utilizing data regarding all household grocery purchases, we calculated the relative contribution of each NOVA category to the total energy acquired by each household (percentage of energy; %E).
Healthy Eating Index-2015
Diet quality of grocery purchases was evaluated using the HEI-2015. The HEI was developed in 1995 by the US Department of Agriculture and is updated every five years to reflect current national dietary recommendations(Reference Guenther, Reedy and Krebs-Smith33, Reference Kennedy, Ohls, Carlson and Fleming34). The HEI-2015 is comprised of nine adequacy and four moderation component scores, the latter being reverse coded(35). The component scores are summed to determine the total HEI-2015 score, which has a maximum score of 100. Higher values indicate better diet quality.
Publicly available SAS code from the Division of Cancer Control and Population Sciences of the National Cancer Institute was used to calculate the HEI-2015 total score and its component scores for the grocery purchases of each household(36). Participating households were categorized into three groups of HEI-2015 total score: <40, 40–59 and ≥60. The cut-offs were established based on the sample distribution in order to ensure a sufficient number of households in each score group. The grading system for the HEI-2015 recommended in the literature categorizes scores from A (>90) to F (<60)(Reference Krebs-Smith, Pannucci and Subar37). As only 26·8 % of households in the analytical sample had a HEI-2015 total score >60, this grading system would not allow us to differentiate between varying diet quality in the majority of the sample. We therefore chose to use the data-driven approach outlined above.
Covariates
Potential covariates were selected based on the literature to be tested in univariate models. Covariates of interest included age of the primary respondent, race of primary respondent (categorized as White, Black, Hispanic and Other race), highest education level in household (categorized as less than high school, high-school degree or some college, bachelor’s degree and above), any smoker in the household (yes or no), participation in the Supplemental Nutrition Assistance Program (SNAP; yes or no), benefits from the Special Supplemental Nutrition Program for Woman, Infants, and Children (WIC; yes or no) and family income-to-poverty ratio (continuous)(Reference Hiza, Casavale and Guenther38, Reference Andreyeva, Tripp and Schwartz39). The number of people in the household (continuous), number of children under 18 years in the household (continuous), census region (categorized as Northeast, Midwest, South, West) and household location (categorized as rural or urban) were also considered potentially relevant covariates.
Statistical analysis
Descriptive statistics (means and frequencies) for household characteristics and food purchases were calculated for the full sample and according to HEI-2015 score group. Differences between household characteristics and food acquisitions according to HEI-2015 score group were determined by Pearson’s χ 2 test for categorical variables, and by unadjusted linear regression, treating HEI-2015 score group as an ordinal variable, for continuous variables.
Mean total and component HEI-2015 scores and associated standard errors were computed for tertiles of ultra-processed food purchases (%E) by unadjusted linear regression, treating tertile of ultra-processed food purchases as an ordinal variable. Unadjusted and multivariable linear regression analysis was used to determine the association between household purchases of ultra-processed foods and processed foods (%E; continuous), respectively, and HEI-2015 total score (continuous). Potential confounders were included separately in the unadjusted model and were retained in the final model if they changed the coefficients of the main exposures by >10 %, as previously published(Reference Makarem, Bandera and Lin40) or if theoretically important based on the previous literature. The final multivariable model was adjusted for race of the primary respondent, number of people in the household, family income-to-poverty ratio, highest education level in household, SNAP participation and any smoker in the household. The exposure variables (ultra-processed foods and processed foods) were included in the same model in order to determine the effect of each exposure independent of the other. Unadjusted linear regression was used to determine mean differences in the most commonly purchased foods (%E) within each NOVA processing level: (i) between households in the top tertile of ultra-processed food purchases (≥67·9 %E) with HEI-2015 total score <40 and ≥60, respectively; and (ii) between households in the bottom tertile of ultra-processed food purchases (<48·4 %E) with HEI-2015 total score <40 and ≥60, respectively.
All analyses were performed using the statistical software package Stata/SE version 14.0 and used the FoodAPS sample weights to account for oversampling of certain populations, non-response and population coverage. The Taylor series linearization variance approximation procedure was used to account for the complex sample design of FoodAPS in the variance estimation. Statistical significance was set to α < 0·05.
Results
Sample characteristics
Household characteristics of the analytical sample, overall and according to HEI-2015 total score groups, are presented in Table 1. Primary respondents were primarily non-Hispanic White (70·3 %) and had attended some college or obtained an associate’s degree (53·4 %). Mean HEI-2015 score was 54·7 (se 0·4). Households with HEI-2015 total score ≥60 were less likely to participate in SNAP (6·3 v. 21·0 %, P < 0·001) and WIC (19·4 v. 33·3 % of WIC-eligible households, P = 0·038), less likely to be located in a rural census tract (28·3 v. 40·8 %, P = 0·004), and the primary respondent was more likely to have a bachelor’s or higher degree (55·3 v. 23·0 %, P < 0·001), compared with households with HEI-2015 score <40.
HH, household; SNAP, Supplemental Nutrition Assistance Program; WIC, Special Supplemental Nutrition Program for Woman, Infants, and Children.
All values are means with their se unless otherwise noted. P values were estimated by unadjusted linear regression, treating HEI-2015 group as an ordinal variable, for continuous variables, and by Pearson’s χ 2 for categorical variables. Missing values: SNAP participation (n 1), anyone in HH receives benefits from WIC (n 3100), race of primary respondent (n 4).
* Of WIC-eligible households (n 861).
Household grocery purchases according to Healthy Eating Index-2015 score group
As shown in Table 2, ultra-processed foods provided 55·8 %E of grocery purchases, while minimally processed foods provided 28·4 %E. The energy share of ultra-processed foods was greater among households with a HEI-2015 total score <40 compared with ≥60 (67·3 v. 47·4 %E, P < 0·001).
%E, percentage of energy.
* Higher scores of the HEI-2015 total and component scores indicate higher diet quality.
Healthy Eating Index-2015 total and component scores according to ultra-processed food purchases
HEI-2015 total and component scores according household grocery purchases of ultra-processed foods (%E; tertiles) are presented in Table 3. On average, households purchasing the least ultra-processed foods had 10·7 points higher HEI-2015 total scores than those purchasing the most ultra-processed foods (HEI-2015 total score of 59·6 v. 48·9, P < 0·001). Households purchasing the highest proportion of ultra-processed foods were furthest from meeting the recommendations of the Dietary Guidelines for Americans 2015–2020 for all food groups and nutrients except whole grains, dairy, fatty acid ratio and saturated fats. HEI-2015 component scores for dairy, fatty acid ratio and saturated fat did not differ according to tertile of ultra-processed food purchases.
P values were calculated by unadjusted linear regression using tertile of ultra-processed food purchases as the independent variable (ordinal) and HEI-2015 score as the dependent variable.
Association between food processing level and Healthy Eating Index-2015 total score
A 10 % increase of purchased energy from ultra-processed foods was associated with a 1·8-point reduction in HEI-2015 total score, controlling for purchases of processed foods (%E), race of the primary respondent, number of people in the household, family income-to-poverty ratio, highest education level in household, SNAP participation and any smoker in the household (β = −1·8, 95 % CI −2·0, −1·5, data not shown). Purchases of processed foods were not associated with HEI-2015 total score in unadjusted (β = −0·3, 95 % CI −0·3, 1·0) and multivariable linear regression analyses (β = −0·1, 95 % CI −0·4, 0·7).
Commonly purchased foods within each NOVA processing level
Households with HEI-2015 total score ≥60 (n 1064) included households both in the lowest (<48·4 %E; n 515) and the highest tertile (≥67·9 %E; n 163) of ultra-processed food purchases, yet the average HEI-2015 total score was significantly higher among those in the lowest tertile of ultra-processed food purchases (71·3 v. 67·0, P < 0·001, data not shown). The energy contribution of the foods most commonly purchased by households in the top tertile of ultra-processed food purchases (≥67·9 %E) with HEI-2015 total score <40 and ≥60, respectively, and households in the bottom tertile of ultra-processed food purchases (<48·4 %E) with HEI-2015 total score <40 and ≥60, respectively, is presented in Fig. 2.
Households in the top tertile of ultra-processed food purchases with a HEI-2015 total score ≥60 purchased significantly more minimally processed fruits and vegetables, milk, plain yoghurt and ultra-processed breakfast cereals than households with equally high ultra-processed food purchases yet HEI-2015 total score <40 (n 385, P < 0·05). Households with high ultra-processed food purchases and a HEI-2015 total score ≥60 also acquired less animal fats; cheese; ultra-processed cakes, cookies and pies; prepared meals; as well as less ice cream (not shown) and carbonated soft drinks (not shown; P < 0·02). Conversely, households in the bottom tertile of ultra-processed food purchases with a HEI-2015 total score <40 (n 135) acquired significantly less minimally processed fruits and vegetables, ultra-processed breakfast cereals and salty snacks but more animal fats, sugar, ultra-processed desserts (not shown) and carbonated soft drinks (not shown) than households with equally low ultra-processed food purchases and HEI-2015 total score ≥60 (P < 0·05).
The relative energy contribution of all foods purchased by households in the top tertile of ultra-processed food purchases (≥67·9 %E) with HEI-2015 total score <40 and ≥60, respectively, and households in the bottom tertile of ultra-processed food purchases (<48·4 %E) with HEI-2015 total score <40 and ≥60, respectively, is presented in the online supplementary material, Supplemental Table S2.
Discussion
In 2012–2013, ultra-processed foods provided the majority of energy in the average US grocery cart (55·8 %E). Higher purchases of ultra-processed foods were associated with lower overall diet quality, defined by the HEI-2015. At similar levels of ultra-processed food purchases, households purchasing more minimally processed fruits, vegetables, roots and legumes had higher diet quality than households purchasing less of these foods.
To our knowledge, previous studies have not assessed the association between NOVA processing levels and diet quality of household purchases, defined by the HEI-2015. Poti et al. compared the content of saturated fat, sugar and sodium of consumer packaged goods across processing levels using data from the 2000–2012 Nielsen Homescan Panel(Reference Poti, Mendez and Ng18). It was observed that 94·7 % of ‘highly processed foods’ (defined as multi-ingredient formulated mixtures processed to the extent that they are no longer recognizable as their original plant or animal source) exceeded the recommendation for saturated fat (>10 %E) of the Dietary Guidelines for Americans 2010, while 94·5 % exceeded the recommendation for sugar (>15 %E) and 96·3 % exceeded the recommendation for sodium (>2400 mg/8368 kJ (2000 kcal)). A significantly higher proportion of highly processed foods (60 %) and ready-to-eat foods (27·1 %) exceeded the Dietary Guidelines for Americans 2010 limits for saturated fat, sugar and sodium, simultaneously, compared with less processed foods (5·6 %) and foods requiring cooking/preparation (4·9 %)(Reference Poti, Mendez and Ng18). These findings are in line with our results for the HEI-2015 added sugar and sodium component scores. It is noteworthy, however, that saturated fat component scores did not differ according to processing level of food purchases in the current study. This may be because we assessed the average saturated fat content of all foods purchased by households, i.e. the whole grocery basket, rather than individual products. In our sample, households in the top tertile of ultra-processed foods purchased the least animal fats, such as butter, lard and cream, which are high in saturated fats.
Our results are consistent with previous international and national studies reporting that diets higher in ultra-processed foods tend to be less nutritious(Reference Martinez Steele, Popkin and Swinburn21–Reference Cediel, Reyes and da Costa Louzada27). Among participants in the National Health and Nutrition Examination Survey (NHANES) 2009–2010, a linear inverse association was found between consumption of ultra-processed foods (%E) and a principal component analysis-derived nutrient-balanced diet pattern high in fibre, potassium, magnesium and vitamin C, and low in saturated fat and added sugars(Reference Martinez Steele, Popkin and Swinburn21). The average dietary content of protein, fibre, vitamins A, C, D and E, zinc, potassium, phosphorus, magnesium and calcium decreased significantly with greater intake of ultra-processed foods (%E), while carbohydrate, added sugar and saturated fat contents increased(Reference Martinez Steele, Popkin and Swinburn21). Studies assessing diet quality in relation to dietary intake of ultra-processed foods in Canada, Brazil, Chile and France report similar findings(Reference Moubarac, Batal and Louzada22–Reference Cediel, Reyes and da Costa Louzada27). Of note, all of the studies examining saturated fat found that intake of saturated fat increased with greater intake of ultra-processed foods(Reference Moubarac, Batal and Louzada22, Reference Adams and White23, Reference Julia, Martinez and Alles25, Reference Batal, Johnson-Down and Moubarac26). The discrepancy between these results and the findings of the current study may be because we considered grocery purchases only and did not account for ultra-processed foods consumed away from home that may be high in saturated fat, such as fast foods.
It is noteworthy that households purchasing the most ultra-processed foods achieved the highest component scores for whole grains, which is a source of dietary fibre. Previous studies have observed an inverse association between ultra-processed food consumption and fibre intake(Reference Martinez Steele, Popkin and Swinburn21–Reference Adams and White23, Reference Julia, Martinez and Alles25, Reference Batal, Johnson-Down and Moubarac26). The greater dietary content of whole grains with higher purchases of ultra-processed foods in our sample may reflect greater acquisition of ultra-processed wholegrain products such as breakfast cereals, granola bars and popcorn. These products are often also high in sugar and/or sodium(Reference Poti, Mendez and Ng18) for which the net effect on overall diet quality is uncertain. Conversely, the higher intake of fibre associated with diets lower in ultra-processed foods in previous studies may reflect greater intakes of other sources of dietary fibre, such as minimally processed legumes, fruits, vegetables, nuts and seeds.
Taken together, our findings, along with previous research using purchasing and consumption data, suggest that the energy share of minimally processed and ultra-processed foods in the diet may influence diet quality. Specifically, increasing the proportion of minimally processed plant foods and reducing the share of ultra-processed foods in the diet seems to be associated with better diet quality. Our results also highlight the importance of healthful food choices within each processing level. In particular, increasing minimally processed fruits, vegetables and legumes and reducing animal fats (processed culinary ingredients) as well as sugary ultra-processed foods, such as cakes, cookies, ice cream and carbonated soft drinks, in the diet may improve diet quality, regardless of the total energy contribution of ultra-processed foods.
Some limitations should be noted. FoodAPS does not have data regarding household consumption and waste of the purchased foods, and hence we cannot not make inferences about individual dietary intake. However, our estimates of food purchases according to processing level are similar to previously published estimates in NHANES 2009–2010 of the mean dietary share from each processing level(Reference Martinez Steele, Popkin and Swinburn21). For example, minimally processed foods provided 30·2 %E and ultra-processed foods provided 57·5 %E in NHANES 2009–2010 (v. 28·4 and 55·8 %E, respectively, in FoodAPS). Furthermore, a considerable proportion of purchased items lacked information on quantity, for which the ERS added imputed estimated quantities(Reference Mancino, Todd and Scharadin30). Nevertheless, the imputed values have been validated by the ERS(Reference Mancino, Todd and Scharadin30). Furthermore, we did not assess energy adequacy. The HEI-2015 uses a density approach and defines diet quality based on food groups and nutrients per 4184 kJ (1000 kcal). It is therefore unknown if participating households acquired an insufficient or excessive amount of energy in relation to their energy needs. Misclassification of foods may have occurred for items for which insufficient details were available to determine processing level with certainty. Nevertheless, the use of the ERS food groups enabled greater precision in the determination of processing level, minimizing the risk of misclassification. Finally, we only assessed household purchases for home consumption and do not know if households differed with regard to the amount and nutritional quality of foods consumed away from home.
Our study also has several strengths. We used comprehensive food purchasing data collected from a nationally representative sample of US households. As a result, the external validity of our findings is likely to be high. Processing level was determined by NOVA, a comprehensive, coherent and specific classification framework that has been recommended by the FAO(41). FoodAPS collects data regarding the exact form in which foods were purchased, thus reducing the risk of misclassification of processing level(32). Finally, the use of purchasing data constitutes another advantage of the current study, as it mitigates the reporting bias associated with assessment of food consumption. Yet, it is worth noting that our results are consistent with previously published NOVA contributions estimated using dietary intake data(Reference Martinez Steele, Popkin and Swinburn21).
In summary, our results indicate that processing level defined by the NOVA framework is useful to predict diet quality, defined according to the Dietary Guidelines for Americans 2015–2020. However, healthful food choices within each processing level are also of relevance for diet quality, e.g. choosing vegetable oils over animal fats and unsweetened over sugary breakfast cereals. Recommendations to consume more healthful unprocessed/minimally foods and less ultra-processed foods may enhance the current food-based Dietary Guidelines for Americans. Future studies should assess if dietary guidelines that combine information on processing level and type of foods facilitate healthier food choices and improve diet quality.
Acknowledgements
Financial support: This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors. B.S.S. was a research fellow of Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES; grant number 88881.131991/2016-01), as a visiting scholar at New York University School of Medicine, USA. CAPES had no role in the design, analysis or writing of this article. Conflict of interest: None. Authorship: F.J. and B.S.S. are joint first authors. N.P., F.J., B.S.S., M.V. and A.D. developed the overall research plan; F.J., E.M.-S. and J.L. provided essential materials; B.S.S. performed the statistical analysis; F.J. and B.S.S. wrote the paper; N.P., M.V., A.D., E.M.-S. and J.L. revised the manuscript for important intellectual content. All authors have read and approved the final manuscript and take full responsibility for the final content. Ethics of human subject participation: This study involved secondary data analysis and no personally identifiable information about individuals was available to the researchers. Ethical approval was therefore not required.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S1368980019001344