Article Text

Original research
Contribution of economic and nutritional context to overweight/obesity dynamics in Indian women from 1998 to 2016: a multilevel analysis of national survey data
  1. Jeswin Baby1,2,
  2. Jithin Sam Varghese3,
  3. Shruthi Cyriac3,
  4. K M Venkat Narayan4,
  5. Anura V Kurpad5,
  6. Tinku Thomas6,
  7. Shivani Patel4
  1. 1Division of Epidemiology and Biostatistics, St John’s Research Institute, Bangalore, Karnataka, India
  2. 2Department of Statistical Sciences, Kannur University, Kannur, Kerala, India
  3. 3Nutrition and Health Sciences Doctoral Program, Laney School of Graduate Studies, Emory University, Atlanta, Georgia, USA
  4. 4Department of Global Health and Epidemiology, Rollins School of Public Health, Emory University, Atlanta, Georgia, USA
  5. 5Department of Physiology, St John's Medical College, Bangalore, Karnataka, India
  6. 6Department of Biostatistics, St John's Medical College, Bangalore, Karnataka, India
  1. Correspondence to Dr Tinku Thomas; tinku.sarah{at}sjri.res.in

Abstract

Background Overweight/obesity increased dramatically among Indian women since 2000. We evaluated the independent contributions of economic and nutrition context to the changing distribution of overweight/obesity among women from 1998 to 2016 across India.

Methods Individual-level data from 473 912 ever married Indian women aged 18–49 in the National Family Health Surveys (1998–1999, 2005–2006, 2015–2016) were merged with year-matched state-level economic and nutrition context indicators. Cross-classified generalised linear mixed models were estimated to quantify associations of contextual characteristics with overweight/obesity (body mass index ≥25 kg/m2) across survey rounds.

Results Between 1998 and 2016, age-standardised prevalence of overweight/obesity increased from 13.9% to 27.5% nationally at an annual growth rate of 0.8%. After accounting for a woman’s age, parity and social class, the adjusted OR (aOR) for overweight/obesity was 2.02 times higher for every unit of state log per capita gross domestic product (GDP) (95% credible interval (CrI) 2.00 to 2.03). Yet, the association of state GDP with overweight/obesity generally decreased over survey round. Women in states with higher per capita daily oil (aOR 1.02 per gram; 95% CrI 1.01 to 1.03) and sugar (aOR 1.05 per gram; 95% CrI 1.04 to 1.05) consumption were more likely to be overweight/obese, while women in states with higher cereal consumption were less likely to be overweight/obese (aOR 0.93 per 10 gram; 95% CrI 0.93 to 0.93).

Conclusions Indicators of state economic development and nutrition transition were independently associated with a woman’s likelihood of being overweight/obese. The impact of state wealth waned over survey round, suggesting that risks for overweight/obesity may be increasingly shaped by individual factors as economic development expands in India.

  • nutrition
  • statistics & research methods
  • epidemiology

Data availability statement

Data for NFHS surveys are available free on request at www.dhsprogram.com. Household-level data for NSSO surveys are freely downloadable from http://microdata.gov.in/nada43/index.php/dissemination. The code for the analysis is available at www.github.com/jeswinbaby/nfhs_overweight. Data is available at www.dhsprogram.com. The code for the analysis is available at www.github.com/jeswinbaby/nfhs_overweight.

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

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, appropriate credit is given, any changes made indicated, 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 explore the annual rate of change in age standardised prevalence between 1999 and 2015 across various sociodemographic categories.

  • We use a Bayesian cross-classified modelling framework to study the economic and nutritional contextual factors associated with overweight/obesity.

  • We consider contextual factors after accounting for individual-level and household-level characteristics from serial cross-sectional surveys using multilevel model.

  • Serial cross-sectional data precludes inference regarding development of overweight/obesity in women over their life course.

  • Lack of contextual data at district and subdistrict levels prohibit further granular analysis of impact of contextual factors.

Introduction

High body mass index (BMI) is the fastest growing contributor to death and disability in India.1 While India has lagged behind other nations in the nutrition transition from underweight to overweight/obesity, the country may experience the largest growth in overweight/obesity among the poor by 2040.2 3 Indian women in particular have been impacted by overweight/obesity. Overweight/obesity has nearly tripled since 1998 among rural women, and the prevalence of overweight/obesity is higher among women compared with men throughout adulthood in both urban and rural settings.4 5 The rise in overweight/obesity among women has also drawn concern because of its implications for potential metabolic risks to the children of overweight/obese pregnant women.6

At the individual level, aside from innate biological factors, excess weight is a result of energy imbalance from energy intake exceeding energy expenditure.7 While literature documents the higher prevalence of overweight/obesity among socioeconomically advantaged women, our understanding of whether and how individual-level risks are augmented by the environment in which women reside is limited.4 8–14 In particular, economic development and food environment are contextual forces that influence population-level energy balance. In low-to-middle income countries (LMICs), economic development is associated with decreased physical activity from mechanised labour and transportation, as well as higher purchasing power to acquire nutrient-dense foods.3 15 16 With sustained economic growth, families are able to afford diets that are more aspirational to conform to consumption behaviours of socioeconomically advantaged groups, introducing an assortment of more energy-dense foods to the food baskets, with high fat and sugar content.17 Understanding the role of the economic and nutrition context in the distribution of overweight/obesity among women may improve our ability to prevent overweight/obesity and address-related health burdens.

Drawing on nationally representative data collected at three time points between 1998 and 2016 in all states of India, we examined the contribution of multiple indicators of economic and nutrition context to overweight/obesity among women. We evaluated the hypothesis that economic and nutritional context influence overweight/obesity above and beyond individual socio-demographic factors, and that these contextual factors are drivers in changes in the distribution of overweight/obesity over time.

Methods

Data sources

We constructed a data set from multiple nationally representative data sources collected in or near 1998–1999, 2005–2006 and 2015–2016 with state and individual-level data to enable a multilevel analysis of population-level overweight/obesity among women of reproductive age across a 15-year period in India.

Individual-level data were from the National Family Health Surveys of India (NFHS) conducted in 1998–1999 (NFHS-2), 2005–2006 (NFHS-3) and 2015–2016 (NFHS-4), designed to be representative at the state and national levels.18–20 While the 2019–2020 NFHS survey has been fielded, data are currently not available. Hence this analysis was restricted to NFHS-2 to NFHS-4. These surveys employed a multi-stage stratified sampling scheme and were designed to be representative at the state and national levels. Strata were defined by urban–rural setting and the primary sampling units (PSU) were villages in the rural stratum and wards in the urban stratum. See online supplemental figure 1 and online supplemental note 1 for detailed description of sample exclusions.

Measures of economic and nutritional context were obtained from external data sources at the state level which allowed us to capture contextual measures that are not available in NFHS. The choice of state as the second level of analysis was pragmatically driven as district-level identifiers (for NFHS 3), and economic and nutrition indicators were not available at lower levels.

State-level contextual measures were derived from four data sources at time periods that were contemporaneous with NFHS rounds. State per-capita gross domestic product (GDP) data were obtained from the Ministry of Statistics and Programme Implementation for 1999–2000, 2004–2005 and 2015–2016. Literacy projections were obtained for 1998–1999, 2005–2006 and 2015–2016 from the Census of India.21 The National Survey Sample Office Employment-Unemployment surveys and Consumer Expenditure surveys provided information on occupation and food consumption in 1999–2000 (round 55), 2004–2005 (round 61) and 2011–2012 (round 68). Both surveys use multistage stratified sampling to select households and are representative at the state and national levels.

Overweight/obesity

Overweight/obesity was measured using directly assessed weight and height at the time of survey. BMI ≥25 kg/m2 was classified as overweight/obese based on WHO international classification.22 BMI<25 kg/m2 was treated as the reference category for all analysis. In analyses that define overweight/obesity using a BMI≥23 kg/m2, which has been documented as a threshold of excess diabetes risk in South Asians, we observed similar associations of contextual factors and overweight/obesity across NFHS rounds (data not shown). Therefore, we applied the internationally accepted definition of BMI ≥25 kg/m2 to define overweight/obesity in this study.

Individual sociodemographic characteristics

Age (18–25; 25–35; 35–49 years), parity (0, 1, 2 and 3 or higher) and highest educational qualification (no formal education; primary (up to fourth class); secondary (up to 10th class); or greater than secondary (11th class and higher)) were self-reported. Social caste was categorised based on constitutionally recognised groups: Scheduled Caste and Scheduled Tribe, Other Backward Castes and General.23–25 Relative socioeconomic position of each woman within her state and survey period was classified using quintiles of a household wealth index derived from 25 household assets common across the three NFHS rounds.4 Residence of each woman was classified as urban or rural by PSU.

State-level contextual variables

We considered three measures of socioeconomic environment (online supplemental note 2) that capture aggregate wealth, education and occupation—GDP per capita, literacy rate and employed adults engaged in sedentary occupation. As measures of food environment (online supplemental note 3), we computed the state-level median daily per capita intake (in grams) of cereals, sugars and oils from food expenditure data.

A composite economic and nutrition score was derived from state-level economic and nutrition contextual indicators, respectively. The composite scores were computed from round-level principal component analyses. The first principal component was used to study the association of collective economic and nutrition environments in overweight/obesity outcomes.

Statistical analysis

States were based on state boundaries defined in NFHS-2 (1998–1999). We computed the age-standardised prevalence of overweight/obesity across strata of sociodemographic characteristics at each survey round. Age-standardised prevalence was calculated as marginal mean estimated from logistic regression models of overweight/obesity with each subgroup as a covariate while adjusting for age. We calculate the annualised rate of change (as percentage points) at the state level as the slope of the linear regression of marginal mean on year. We used NFHS sample weights, scaled to represent the population of women 18–49 years for each survey, for prevalence estimation.26 Association between prevalence and state-level contextual indicators were examined using Pearson’s correlation coefficients.

Our primary outcome was individual-level overweight/obesity among women assessed in NFHS. Overweight/obesity was modelled as a function of individual-level sociodemographic characteristics and state-level contextual factors as higher-level variables. A four level (individual, PSU, state and survey year) multilevel model approach was adopted to deal with the hierarchical, clustered structure of the data, as well as to understand the individual and state level contextual factors on individual-level overweight/obesity status of women. A cross-classified Bayesian generalised linear mixed models (GLMM) with logistic link function were used to estimate the association of contextual variables with overweight/obesity with and without adjustment for contextual characteristics. We checked the assumptions for cross-classified nature of the data, namely homogeneity of variance and normality of residuals. The 4-level primary GLMM with individual covariates and random intercepts was specified as

Embedded Image

Embedded Image

where overweight/obesity status O for the ith individual in pth PSU cross-classified within sth state and rth NFHS survey year(online supplemental figure 2) was modelled as a function of M individual-level covariates (X).

We used a cross-classified nesting structure to appropriately reflect the hierarchical nature of the data, in which each woman resided in a PSU (rural village or urban ward) cross-classified under a given state and round (online supplemental note 4 and 5).27 We estimated the variance partition coefficient (VPC) across the different levels for the null model (no covariates) using a simulated Intra class correlation (ICC) (online supplemental note 6).28 We also estimated associations of contextual indicators with overweight/obesity separately by round.

All associations are reported as adjusted OR (aOR) and 95% credible intervals (CrI) conditional on random effects. These may be interpreted as the association of a given covariate for women in a typical village/ward nested within a typical state and round. Statistical analyses were conducted using R V.3.5.0, MLwiN V.3.03 using R2MLwiN V.0.8–V.5 and Stata V.16.1.

Results

The final analytical sample consisted of 473 912 ever-married women aged 18–49 years residing in 26 states in 1998–1999, 2005–2006 or 2015–2016. We noted sociodemographic changes among the sample over the study period (table 1). From 1998–1999 to 2015–2016, women with no formal education declined (51%–36%), women aged 35–49 years increased (46%–51%), women who were nulliparous increased (60%–70%), and women in constitutionally protected castes decreased (40%–25%). At the state level, inflation adjusted per-capita GDP doubled, and the literacy rate rose by over one-third (table 2). The population engaged in sedentary occupations was similar over the study period. There were no changes in median per capita cereal intake or sugar intake. However, median oil intake increased from 7.3 to 8.6 g. Despite large temporal variability in GDP, literacy, and oil consumption, the relative position of states remained consistent from 1998 to 2016 (rank correlation ≥0.9; Online supplemental table S1). Economic and nutrition indicators (except cereal intake) were positively correlated with overweight/obesity prevalence for all rounds (r=0.35–0.80; figure 1).

Figure 1

Prevalence of overweight/obesity by economic and nutrition indicators. X-axis represents state wise (A) log (per capita GDP in rupees) (B) literacy rate (in %) (C) employed adults engaged in a sedentary occupation (in %) (D) median per capita oil intake (in grams) (E) median per capita sugar intake (in grams) (F) median per capita cereal intake (in grams). each point represents a state value for different rounds (white squares—NFHS-2, black circles—NFHS-3 and plus sign—NFHS-4). Sedentary occupation, and nutritional intakes were computed using survey estimation procedures as per respective NSSO rounds. The values corresponding to rNFHS-2, rNFHS-3 and rNFHS-4 represent the Pearson’s correlation coefficient between the state economic and nutritional indicators and overweight/obesity prevalence by NFHS round. GDP, gross domestic product; NFHS, National Family Health Survey; NSSO, National Sample Survey Office.

Table 1

Sociodemographic characteristics of Indian women ages 18–49 years from 1998–1999 to 2015–2016

Table 2

Distribution of state-level contextual indicators

Between 1998–1999 and 2015–2016, age-standardised prevalence of overweight/obesity increased (13.9%–27.5%; table 3). The increase in prevalence was observed across all states with varying rates (online supplemental figure 3). The overall annualised increase in the prevalence of overweight/obesity was consistently higher among more socioeconomically advantaged groups (table 3).

Table 3

Age-standardised prevalence of overweight/obesity and annualised changes in Indian women from 1998–1999 to 2015–2016

In multilevel logistic regression analyses of the association between individual characteristics and overweight/obesity across the three time periods, we observed a strong positive association of individual-level education, wealth, and age with overweight/obesity status (online supplemental table S2). The aOR for overweight/obesity were 2.04 (95% CrI 1.99 to 2.09), 2.14 (95% CrI 2.09 to 2.21) and 3.95 (95% CrI 3.84 to 4.07) for women in the highest education, wealth and age groups, respectively, relative to reference levels of each domain. In contrast, there were inverse associations between higher parity and overweight/obesity as well as between belonging to a constitutionally protected caste and overweight/obesity. Similar associations were observed by urban and rural areas (Online supplemental table S2). The associations were of lower magnitude for rural areas. The association of individual-level covariates with overweight/obesity were not impacted by inclusion of state-level variables.

Across survey rounds, the variability (as measured by VPC) in overweight/obesity was largely at the individual level in both urban (85.8%) and rural (70.2%) areas (Online supplemental table S3; null model). Although geographical variability at the ward/village and state levels was comparable between urban and rural areas, there was more temporal heterogeneity in overweight/obesity between rounds in rural (19.5%) as compared with urban areas (5.1%). Similarly, state variability was higher in rural (6.1%) compared with urban areas (2.7%).

Independent of individual socio-demographic characteristics, composite state-level economic (aOR 1.15; 95% CrI 1.09 to 1.22) and nutrition (aOR 1.19; 95% CrI 1.16 to 1.22) scores were positively associated with prevalent overweight/obesity status in women (table 4). With respect to individual indicators, the adjusted relative odds of overweight/obesity were aOR=2.02 (95% CrI 2.00 to 2.03) for every unit of log state GDP per capita. Similarly, we observed positive associations between literacy rate (aOR 1.32; 95% CrI 1.31 to 1.34) and sedentary occupation (aOR 1.05; 95% CrI 1.02 to 1.06). Residing in a state with higher daily oil (aOR 1.02 per gram; 95% CrI 1.01 to 1.03) and sugar (aOR 1.05 per gram; 95% CrI 1.04 to 1.05) consumption was associated with greater odds of prevalent overweight/obesity. Women in states with higher cereal consumption were less likely to be overweight/obese (aOR 0.93 per 10 g; 95% CrI 0.93 to 0.93). All state-level economic and nutrition indicators were associated with overweight/obesity in both urban and rural areas, with stronger associations of log state GDP and literacy with overweight/obesity in rural areas.

Table 4

Association of state level characteristics with overweight/obesity in pooled data of Indian women in 1998–1999, 2005–2006 and 2015–2016

Figure 2 shows the differing associations of contextual indicators with overweight/obesity by round. The association of state GDP with overweight/obesity was positive across all rounds but declining over time. Nationally, the aOR of log per capita GDP for overweight/obesity was attenuated from 2.08 (95% CrI 1.99 to 2.18) in 1998–1999 to 1.20 (95% CrI 1.18 to 1.21) in 2015–2016. We found a similar pattern in urban and rural models. Similarly, the association of literacy, sedentary occupation, oils and sugar tended to reduce over time.

Figure 2

Associations of state-level characteristics with overweight/obesity in Indian women in 1998–1999, 2005–20006, 2015–2016. Separate models were fit for each state-level variable. All values are aOR (95% CrI) estimated from logistic cross-classified random intercept model adjusted for woman’s age category, educational attainment, wealth index, parity and caste. aOR, adjusted OR, CrI, credible interval; GDP, gross domestic product.

Comparing the differential impact of contextual measures, the composite state-lvel nutrition score yielded the best model fit based on deviance information criterion (DIC) in urban and rural areas separately (table 4). State-level literacy rate and per capita sugar yielded the best model fit with respect to individual contextual indicators.

Discussion

Between 1999 and 2015, the age-standardised prevalence of overweight/obesity more than doubled in Indian women, increasing by 0.81 percentage points each year. Independent of individual and household socio-demographic characteristics, state-level economic and nutritional contexts was associated with being overweight/obese across the 20-year period of analysis. Across the study period, we observed robust positive associations of state-level log per GDP capita, literacy rate, sedentary occupation, oil consumption and sugar consumption with individual-level overweight/obesity nationally and within urban and rural settings. Residing in states with higher cereal consumption was associated with lower odds of overweight/obesity in both urban and rural women.

There is consensus that cross-national, sub-national, and temporal variations in overweight and obesity are driven by differences in economic development and nutritional environment, as described in the literature around the nutrition and obesity transition.16 29 30 Within LMIC settings, the risk of overweight and obesity are generally higher among wealthier individuals.16 29 30 Across countries, obesity is more prevalent in high income countries (HIC), and over time, the difference in obesity prevalence between HICs and LMIC has diminished. Yet there is relatively limited empirical research quantifying the direct impact of specific contextual characteristics on prevalent overweight/obesity in LMICs at the national or individual level.31 Our descriptive analyses begin to close this gap and lend support to the socioecological paradigm underpinning the rise of overweight and obesity in India.

Stratifying by urbanicity and survey round revealed important findings. With respect to urbanicity, we found that the impact of state wealth (measured as log per capita GDP) on overweight/obesity was larger in rural compared with urban areas. While the relative odds of overweight/obesity were 1.39 times for every unit increase in urban areas, the relative odds of overweight/obesity were 2.16 times in rural areas across the study period. If urban is to be compared with HIC and rural with low income countries, the above finding is comparable to the higher prevalence in wealthy observed in the low income countries and among the disadvantaged populations in the higher-income countries.31 With respect to survey round, our analysis indicates that the impact of contextual characteristics has generally waned over time. Nationally, the association between state-level log per capita GDP and overweight/obesity waned over time. Taken together, the larger associations of GDP and overweight/obesity in rural areas and the declining role of GDP over time may indicate that the role of contextual factors in shaping obesity risk attenuates as the prevalence of overweight/obesity increases, household wealth and purchasing power rises, and environments become more homogeneous with economic development. This hypothesis must be further evaluated using well-characterised longitudinal data with measures of environment and weight status.

The positive though modest association of oil and sugar and the negative association of cereals with overweight/obesity ring an alarm. The country should be cautious of rise in consumption of energy-dense foods which are associated with overweight/obesity that are replacing cereals.32–35 Cereal consumption may have decreased over time (table 2) to maintain balance with reduced energy requirement of the population.36 The working group on addressing consumption of high fat, salt and sugar and promotion of healthy snacks in schools of India (constituted by the Government of India) recognised the rising intake of processed foods in both rural and urban areas. The ease of availability, taste, low cost, peer pressure and aggressive marketing and advertisements were identified as reasons for this rise among children and adolescents.37

While the annual rate of change in overweight/obesity was faster among advantaged groups, our results reiterate the finding that that overweight/obesity increased substantially across all education, wealth and caste groups.4 Socioeconomically vulnerable groups may require new vigilance for emergence of obesity-related chronic disease from a public health perspective. At 65%, the lifetime risk of diabetes in women in Indian cities is among the highest in the world, and the rise of overweight/obesity in rural women may usher in similar risks.38

Our study has numerous methodological strengths. Using a Bayesian GLMM has numerous advantages from providing stable estimates incorporating the hierarchy of space and time, better convergence to point estimates and modelling a cross-classified framework. Partitioning of variance using a simulated ICC allows for fewer approximations of nesting. Additionally, our study used the best available data sources which were nationally representative avoiding the pitfall of correlated errors when computing contextual effects from the same dataset.

Our study has several limitations. The use of serial cross-sectional data does not allow us to make inferences regarding the development of incident overweight/obesity in women over their lifetimes. Furthermore, we lacked individual-level data on modifiable risk factors for obesity such as diet and physical activity, sleep, which may have changed over time. Many of these health behaviours, however, are considered mediators that link context to weight status. BMI as measure of overweight and obesity, and not waist circumference or body fat percentage may be an inadequate measure of adiposity and related health risks in South Asians.3 11 Data on state-to-state migration of women were not available, and thefore, the analysis assumed that the current state of residence appropriately measured state-level contextual exposures for each woman. While contextual sources were matched closest to date of survey, lack of such data at lower geographical units (district or village) prohibit more granular analysis of the impact of contextual factors. We did not adjust for survey weights in the multilevel models due to software limitations. We analysed data on women 18–49 years and restricted our analysis to those who were non-pregnant/non-lactating, thus limiting generalisability.

Our analysis of overweight/obesity dynamics in India was constrained by several data limitations which motivate comment on several gaps and opportunities in the data infrastructure. For example, district-level identifiers across all survey around would have enabled a more granular evaluation of contextual determinants of excess weight, high quality administrative data on contextual characteristics would enable rigorous evaluation of specific features of context, and longitudinal data tracking individual weight status over time would provide a basis for understanding the risk (rather than prevalence) of overweight/obesity over the past two decades. This research infrastructure gap presents policymakers, private foundations, and academic institutions with an opportunity to establish—and make public—data sources that allow investigators to prospectively study the temporal association of regional and local environment on incidence of cardiovascular disease and diabetes for decentralised decision making. For example, placing integrated data reserves such as those at TATA-NIN (https://www.nintata.res.in) in the public domain would substantially enhance the ability of scientists to study the complex interplays of various factors over time.

Conclusion

This study identified a significant rise in overweight/obesity over the past two decades among Indian women under the age of 50, an age group that makes up over 75% of all women in India.21 Based on lessons from other settings, the rise in overweight/obesity in woman signals a likely future rise in overweight/obesity among children and men.16 While residing in a state with higher per capita GDP, oil and cereal consumption was associated with higher odds of overweight/obesity, the impact of these contextual risk factors may be diminishing over time. Further investigation into what may explain the apparent changing associations between contextual characteristics and overweight/obesity may provide insight on environmental characteristics that can be leveraged to halt the rise of excess weight in India. Future studies may also consider mechanisms linking environment with weight outcomes, in particular on the relative importance of food availability and social norms around diet are driving these results.

Data availability statement

Data for NFHS surveys are available free on request at www.dhsprogram.com. Household-level data for NSSO surveys are freely downloadable from http://microdata.gov.in/nada43/index.php/dissemination. The code for the analysis is available at www.github.com/jeswinbaby/nfhs_overweight. Data is available at www.dhsprogram.com. The code for the analysis is available at www.github.com/jeswinbaby/nfhs_overweight.

Ethics statements

Patient consent for publication

Ethics approval

The analysis was reviewed by the St. John’s Medical College Institutional Review Board; it was considered exempt from full review because it is based on an anonymous public use data set with no identifiable information on survey participants.

References

Supplementary materials

Footnotes

  • TT and SP contributed equally.

  • Contributors SP and TT conceptualised the manuscript and statistical analysis plan. JB led the statistical analysis and JSV contributed to the analyses. JB, JSV, SP and TT interpreted the results and drafted the manuscript with inputs from SC, AVK and KMN. All authors reviewed the manuscript. TT is the guarantor.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

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

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.