We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
Population versus hospital controls for case-control studies on cancers in Chinese hospitals
BMC Medical Research Methodologyvolume 11, Article number: 167 (2011)
Correct control selection is crucial to the internal validity of case-control studies. Little information exists on differences between population and hospital controls in case-control studies on cancers in Chinese hospital setting.
We conducted three parallel case-control studies on leukemia, breast and colorectal cancers in China between 2009 and 2010, using population and hospital controls to separately match 540 incident cases by age, gender and residency at a 1:1 ratio. Demographic and lifestyle factors were measured using a validated questionnaire in face-to-face interview. Odds ratios (ORs) and 95% confidence intervals (CIs) were obtained using conditional logistic regression analyses.
The two control groups had closely similar exposure distributions of 15 out of 16 factors, with the only exception being that hospital controls were less likely to have a BMI ≥ 25 (OR = 0.71, 95% CI: 0.54, 0.93). For exposure of green tea drinking, the adjusted ORs (95% CIs) comparing green tealeaves intake ≥ 1000 grams annually with non-drinkers were 0.51 (0.31, 0.83) and 0.21 (0.27, 0.74) for three cancers combined, 0.06 (0.01, 0.61) and 0.07 (0.01, 0.47) for breast cancer, 0.52 (0.29, 0.94) and 0.45 (0.25, 0.82) for colorectal cancer, 0.65 (0.08, 5.63) and 0.57 (0.07, 4.79) for leukemia using hospital and population controls respectively.
The study found that hospital controls were comparable with population controls for most demographic characteristics and lifestyle factors measured, but there was a slight difference between the two control groups. Hospital outpatients provide a satisfactory control group in hospital-based case-control study in the Chinese hospital setting.
Correct control selection is crucial to the internal validity of case-control studies [1, 2]. The function of controls is to provide valid information on the distribution of exposure within the population at risk of becoming a case . When cases are ascertained from hospitals serving a defined geographic area, a probability sample of unaffected individuals from the population of that geographic region can be used to enhance the likelihood that cases and controls come from the same source population. A concern with this approach is that population controls sampled from a convenient population register may not be representative of the true population at risk of being a member of a hospital case series; population controls may also lack comparability in recruitment fraction and recall of information [1–5]. But there are also concerns that hospital controls, especially those with other diseases, may fail to provide an unbiased sample of the population at risk with respect to exposure status [1–5].
When we applied successfully to Australia's National Health and Medical Research Council to perform case-control studies in China on the effects of green tea on the incidence rates of colorectal cancer, breast cancer and adult leukemia, it became a condition of funding that we recruit both hospital and population controls for around one fifth of the case series to determine if there was any difference of practical importance.
A handful of studies have compared different control groups in western countries where non-emergency access to hospital care generally depends on referral from medical practitioners in non-institutional settings [6–14]. When researchers have enrolled both hospital and population controls, differences in effect estimates using the two groups have varied from inconsequential [6, 7] to problematic [8–10], prompting discussions about characteristics of the exposure, selection protocols and clinical dynamics. Some studies have specifically reported differences between hospital and population controls in key exposure distributions [11–14], but elsewhere the two control groups were found to have very similar exposure distributions .
It was important to understand whether exposures measured in the study were in agreement when ascertained from population and hospital control groups, especially because our previous case-control studies of green tea and cancer in China had used hospital controls. This study aimed, therefore, to compare key exposure distributions between hospital controls and alternative population controls to determine if there was any difference of practical importance, and to examine if inferences drawn from green tea effect estimates for different cancer case groups using hospital controls would be different from those using population controls in the Chinese hospital setting.
Study design and participants
The study, as a validation component of three large case-control studies of malignancies, was conducted in Shenyang, the capital city of Liaoning Province, Northeast China. In August 2009 to July 2010, 540 incident cases with a primary diagnosis of leukemia, or breast cancer, or colorectal cancer were identified from histopathology and haematology records at the First Hospital of China Medical University, a public teaching hospital with 2,249 beds, around 32,000 inpatients annually and 3,000 outpatients daily. The eligible cases were permanent residents of urban Shenyang aged 18 to 85 years. During the study period, 540 population controls and 540 hospital controls were selected from Shenyang residents and from outpatients at the same hospital to separately match cases in a 1:1 ratio. Most of controls were interviewed within three months after cases were interviewed.
The methods of recruiting population controls were similar to those used in case-control studies in Shanghai, China [16, 17]. Population household registries, which kept records of all permanent residents in urban Shenyang, were used to select controls from the five metropolitan districts in Shenyang, namely Heping, Shenhe, Dadong, Huanggu, and Tiexi. With the assistance of the local community councils, residents who lived at their registered address during the study period were randomly selected from household registry rolls. Residents were eligible as population controls if they matched with individual cases of an updated list by gender and year-of-birth quinquennium on a given selection day.
Hospital controls were drawn from the population of patients in the hospital . A systematic selection process used in our previous studies was adopted for hospital control recruitment . They were selected from outpatients who attended the Medical Examination Centre at the First Hospital of China Medical University and were permanent residents of urban Shenyang. The eligible hospital controls were those without any malignancy after they had consulted their doctors. Each hospital control was selected as the first attendee on a given selection day to match the next case on a daily updated list of cases by sex and 5-year age group. Hospital outpatients were excluded as a control if they were not matched their corresponding cases by gender, year-of-birth quinquennium and living areas, and if they had a diagnosis of any malignancy before or after recruitment. The project protocol had received ethics approval from both the Human Research Ethics Committee of The University of Western Australia and the First Hospital of China Medical University authority.
Questionnaire and interview
Subjects were briefed regarding confidentiality and anonymity issues and the general aims of the study to investigate lifestyle factors. An appointment for an interview was made after obtaining the respondent's consent via an initial contact. A face-to-face interview was then conducted by the first author, using a structured questionnaire and usually took 30-40 minutes. The validated questionnaire, available from the authors upon request, was used to collect the information on: (i) demographic and lifestyle characteristics, e.g., area of residence, education, smoking, alcohol and tea consumption and physical activity; (ii) dietary intake assessed by a food frequency questionnaire (FFQ); and (iii) factors relevant to hormonal status and family history of cancer. The questionnaire was adapted from that used in our previous studies on cancers . This instrument was originally modified from one used for studying cancers in Shanghai in order to ensure cultural relevance . The questionnaire was translated into Chinese and checked using back-translation by professional Chinese translators. The internal consistency and reliability of the questionnaire was assessed in a preliminary study and then evaluated by a test-retest. The intraclass correlation coefficients for mean daily intake of tea and alcohol were 0.83 and 0.88 . Thus high coefficients for test-retest reliability suggested that the questionnaire may be relied upon in assessing selected demographic characteristics and lifestyle factors .
All data were checked at the end of each interview for completeness and were coded and analysed using SPSS version 18.0. Participants' self-reported current height in meters and weight in kilograms were used to calculate body mass index (BMI) (weight/height2). Daily energy intake and alcohol consumption were assessed using the FFQ. The frequencies of 100 food items, including beer, wine, and liquor intake, were assigned into nine categories: never or hardly ever; once a month; 2-3 times a month; once a week; 2-3 times a week; 4-6 times a week; once a day; twice a day; and ≥ 3 times a day. Food and alcohol consumption was based on habitual diet and a 'reference' recall period was set as one year prior to diagnosis for cases or interview for controls. If there was any recent change in habits, only information on the habits before the change was used in data analysis. The frequency and quantity variables derived from the FFQ were converted into daily food consumption, adjusted for the edible portions of foods, cooking methods, seasonal factors, and market availability . Total energy intake was estimated using Chinese Food Composition Tables . The frequency and quantity variables for beer, wine, and liquor were converted into daily intake in ml. Amounts of ethanol ingested were calculated by assuming 10 g of ethanol per 285 ml of beer, per 100 ml of wine, and per 30 ml of liquor based on a method used in a previous study . Physical activity was expressed in terms of weekly metabolic equivalent task hours (MET hrs/week) . MET scores of 6, 4.5, and 2.5 were assigned respectively for vigorous, moderate, and walking activity based on a compendium of physical activities . The quantitative variables for physical activity (MET) were divided into quartiles based on the corresponding empirical distribution of population controls, with the lowest quartile being the reference category.
Selected demographic characteristics and lifestyle factors in both control groups were first compared using a t-test for continuous variables and chi-square test for categorical variables. Associations of hospital/population control status with exposure variable were then assessed using adjusted odds ratios (ORs), 95% confidence intervals (95% CI) and p-values for trend estimated from conditional logistic regression adjusted for education, income, household size, BMI, smoking, alcohol, tea consumption, energy intake (kcal), and physical activity. The ORs were for the odds of being a hospital rather than population control. Each quantitative or ordinal quantitative measure was subjected to a linear trend test.
All cancer case groups combined, breast cancer, colorectal cancer, and leukemia only subgroups were analysed separately using the two types of matched control groups with regard to green tea intake using conditional logistic regression, adjusted for education, BMI 5 years ago, smoking, passive smoking, alcohol consumption, energy intake (kcal), physical activity, and cancer in first degree relative. For breast cancer only subgroup, the regressions analyses also adjusted for menopausal status, oral contraceptive use, and number of children breastfed.
Data were collected on a total of 1080 controls. Interviews were completed for 540 (89.7%) of the 602 eligible population controls who were approached to participate. The only reason for non-participation was refusal to provide a blood sample, which was an essential part of the study. Of 563 eligible hospital controls who were approached to participate, interviews were completed also for 540 (95.9%). For the 540 recruited cancer cases, the response fraction was 98.5%.
Comparisons of distributions of demographic characteristics and lifestyle factors between hospital and population controls are shown in Table 1. The two control groups were remarkably similar in their exposure distributions for age (forced by matching), marital status, household size, education, income, body mass index, smoking, alcohol intake, physical activity, tea drinking and family history of malignancy.
Table 2 presents the adjusted ORs for demographic characteristics and lifestyle factors, representing the tendency for each factor to associate independently with status as a hospital rather than population control. The 29 points estimates of non-baseline ORs ranged from 0.46 to 1.55, whilst 19 were within the range of 0.80-1.25. Two of the 29 confidence intervals excluded the null value, being associated with low ORs for a household size of 3 members compared with 1-2 members (0.46) and a BMI ≥ 25 compared with < 25 (0.71). However, there was no consistent trend in the association of household size with hospital/population control status, the OR for a household size of ≥ 4 compared with 1-2 members being 0.83 (p-value for trend = 0.21). For the key important study exposure of green tea drinking frequency, the ORs revealed no evidence of association with control type, being 0.94 at ≤ 6 times a week and 1.08 at ≥ 1 time per day compared with never or seldom drinking green tea (p-value for trend = 0.58).
Table 3 presents ORs according to green tea intake for the cancers combined and for each cancer case subgroups, using the two types of matched control groups. All ORs for the frequency of green tea drinking and amount of dried green tea leaves consumed (g/year) were on the same side of the null value and quite close in all the comparisons estimated using the 2 types of control groups. In the final models, compared with nondrinkers vs. dried green tealeaves intake ≥ 1000 grams annually, the adjusted ORs (95% CI) were 0.51 (0.31, 0.83) and 0.21 (0.27, 0.74) for all cancer case groups combined; 0.06 (0.01, 0.61) and 0.07 (0.01, 0.47) for breast cancer; 0.52 (0.29, 0.94) and 0.45 (0.25, 0.82) for colorectal cancer; 0.65 (0.08, 5.63) and 0.57 (0.07, 4.79) for leukemia using hospital outpatient and population controls respectively in each instance.
This validation study, as a component of three large case-control studies of malignancies, used separate hospital outpatient and population control groups to assess differences in the distributions of demographic characteristics and lifestyle factors between the two control groups. Furthermore, the study evaluated if different conclusions were reached when comparing cases with matched hospital controls or population controls in the effect of green tea intake. The study found that the distributions of study exposures in hospital outpatient controls were similar to those in population controls. In addition, we obtained similar results in assessing the effect of green tea drinking on three cancers together and separated, using hospital outpatient controls and population controls respectively. Therefore, our results suggest that in the context of a case-control study conducted at a Chinese hospital using urban cases, regardless of whether controls were selected from hospital outpatient attendees without malignancy or drawn from population household registries covering populations in the catchment area of the participating hospital, the two control groups appeared to have similar exposure distributions. Our results identified one exception to this general conclusion. It is possible that population controls may tend to report higher levels of being overweight than hospital outpatient controls.
The findings of this investigation were consistent with a smaller pilot study conducted in 1999-2000 by our research team , and with other studies which suggested that population and hospital controls were generally similar in their studies [8–10], as well as with a research by other investigators, who compared non-cancer outpatients with population controls in a case-control study in Japan .
A key exposure in our work was green tea drinking measured in three parallel case-control studies. For three cancers together, breast cancer, colorectal cancer, and leukemia, the inversed associations were similar in their direction and size regardless of which control series was employed to represent the underlying exposure to higher green tea consumption. The results, using either type of control series, were consistent with our previous finding that increasing the quantity of green tea consumption reduces the risk of breast cancer in Chinese women , as well as the results from four studies pooled in a meta-analysis, which reported a reduced risk of breast cancer for highest versus non/lowest intake of green tea . The association between colorectal cancer and green tea consumption assessed in this study were consistent with another meta-analysis combining results from eight studies for high green tea intake . Although we recruited less leukemia cases, an inverse association was also observed, although it fell short of statistical significance, which was similar to our previous study . The effects of green tea consumption estimated in this study suggest that hospital outpatient controls performed little different from population controls in green tea exposures.
Whilst hospital outpatient and population controls selected in this study may function in a similar manner for a hospital case-control design in China, this does not entirely resolve the underlying theoretical debate as to which of the two types of control group is the more defensible. For BMI, for example, where our study did find evidence of some practical and statistical differences, one is left with the question of which control exposure distribution is more representative of that in the population at risk of becoming a case. The result may have simply been due to chance, although the possibility that population controls were more likely to be overweight cannot be ignored. On the other hand, the lower response fraction of population controls (89.7% vs 95.9% for hospital controls) may have provided the circumstances where normal weight individuals in the general population were less inclined to participate in a research study on health and lifestyle.
In China, public hospitals provide treatment services for acute, severe, and critical outpatients and inpatients, and they fulfill education, research, prevention, and rehabilitation objectives by providing special medical services. All public hospitals have preventive health care branches, which are responsible for reporting infectious diseases, health checks, health counseling, community prevention services, health education, disease screening, family planning and birthing guidance, and care of staff . Patients living in cities readily visit hospitals as non-referred outpatients for check-ups. Cultural factors and insurance arrangements lead patients to maintain a strong relationship with one particular hospital, where they receive a complete range of health care [32, 33]. Survey data in 2008 showed low levels (14%) of community health service utilization, suggesting that community health services are not yet the first point of contact with the health system in China . Therefore, hospital outpatients in China are somewhat similar to ambulatory patients visiting GP clinics in Western countries.
As there is almost never one ideal control group, some researchers believe that population registers provide the most valid way of sampling controls in a hospital case-control study [5, 11–14], because the main theoretical strength of population controls is the potential to provide information of exposure that is unaltered by associations with illness . Others suggest that generally, hospital controls should be preferred in a hospital case-control study regarding practicability, cost and travel time for face to face interviews [9, 10, 15]. There may also be differences in the capacity to recall and report exposures between randomly selected population members and those who are actively engaged in the health system . In addition, only hospital controls have shown some evidence that in the event of developing the cancer, they would be likely to attend the hospital and become a case in the study [1–3].
With respect to the demographic characteristics and lifestyle factors included in the present study, our results suggest that hospital outpatient controls performed little different from population controls for most exposures. For the key exposure, the effect estimates for different cancer case groups using hospital outpatient controls were similar to those using population controls. Therefore, even though some points of concern exist, hospital outpatients provide a satisfactory control group in hospital-based case-control study in the Chinese hospital setting.
metabolic equivalent task hour
Body mass index
- 95% CI:
95 percent confidence interval.
Wacholder S, Silverman DT, McLaughlin JK, Mandel JS: Selection of controls in case-control studies. I principles; II types of controls; III design options. Am J Epidemiol. 1992, 135: 1019-1050.
Rothman KJ, Greenland S, Lash TL: Case-control studies. Modern epidemiology. Edited by: Rothman KJ, Greenland S, Lash TL. 2008, Philadelphia: Wolters Kluwer Health/Lippincott Williams & Wilkins, 111-127. 3
Miettinen OS: Theoretical Epidemiology: Principles of Occurrence Research in Medicine. 1985, New York: John Wiley & Sons, 46-83.
Miettinen OS: The concept of secondary base. J Clin Epidemiol. 1990, 43: 1017-1020. 10.1016/0895-4356(90)90092-4.
Grimes DA, Schulz KF: Compared to what? Finding controls for case-control studies. Lancet. 2005, 365: 1429-1433. 10.1016/S0140-6736(05)66379-9.
Olson SH, Kelsey JL, Pearson TA, Levin B: Characteristics of a hypothetical group of hospital controls for a case-control study. Am J Epidemiol. 1994, 139: 302-311.
Moritz DJ, Kelsey JL, Grisso JA: Hospital controls versus community controls: differences in inferences regarding risk factors for hip fracture. Am J Epidemiol. 1997, 145: 653-660.
Cogswell ME, Bitsko RH, Anderka M, Caton AR, Feldkamp ML, Sherlock SMH, Meyer RE, Ramadhani T, Robbins JM, Shaw GM: Control selection and participation in an ongoing, population-based, case-control study of birth defects. Am J Epidemiol. 2009, 170: 975-985. 10.1093/aje/kwp226.
Infante-Rivard C: Hospital or population controls for case-control studies of severe childhood disease?. Am J Epidemiol. 2003, 157: 176-182. 10.1093/aje/kwf174.
Laerlev L, Lynge E, Sabroe S, Olsen J: Colon cancer controls versus population controls in case-controls studies of occupational risk factors. BMC Cancer. 2004, 4: 15-20. 10.1186/1471-2407-4-15.
Malagoli C, Vinceti M, Pellacani G, Sieri S, Krogh V, Seidenari S, Bergomi M: Diet and melanoma risk: effects of choice of hospital versus population controls. Tumori. 2008, 94: 669-673.
Ruano-Ravina PA, Pérez-Ríos M, Barros-Dios JM: Population-based versus hospital-based controls: are they comparable?. Gac Sanit. 2008, 22: 609-613. 10.1016/S0213-9111(08)75363-9.
Neupane B, Walter SD, Krueger P, Loeb M: Community controls were preferred to hospital controls in a case-control study where the cases are derived from the hospital. J Clin Epidemiol. 2010, 63: 926-931. 10.1016/j.jclinepi.2009.11.006.
Dantas O, Ximenes RAA, de Albuquerque Mfpm, Montarroyos UR, de Souza WV, Varejão P, Rodrigues LC: Selection bias: neighbourhood controls and controls selected from those presenting to a Health Unit in a case control study of efficacy of BCG revaccination. BMC Med Res Methodol. 2007, 7: 11-17. 10.1186/1471-2288-7-11.
Inoue M, Tajima K, Hirose K, Hamajima N, Takezaki T, Kuroishi T, Tominaga S: Epidemiological features of first-visit outpatients in Japan: comparison with general population and variation by sex, age, and season. J Clin Epidemiol. 1997, 50: 69-77. 10.1016/S0895-4356(96)00297-1.
Chang SC, Rashid A, Gao YT, Andreotti G, Shen MC, Wang BS, Han TQ, Zhang BH, Sakoda LC, Leitzmann MF: Polymorphism of genes related to insulin sensitivity and the risk of biliary tract cancer and biliary stone: a population-based case-control study in Shanghai, China. Carcinogenesis. 2008, 29: 944-948. 10.1093/carcin/bgn025.
Hsing AW, Sakoda LC, Rashid A, Andreotti G, Chen J, Wang BS, Shen MC, Chen BE, Rosenberg PS, Zhang M: Variants in inflammation genes and the risk of biliary tract cancers and stones: a population-based study in China. Cancer Res. 2008, 68: 6442-6452. 10.1158/0008-5472.CAN-08-0444.
Porta M, Greenland S, Last JM: A dictionary of epidemiology. 2008, New York: Oxford University Press, 53-5
Zhang M, Holman CDJ, Huang JP, Xie X: Green tea and the prevention of breast cancer: a case-control study in Southeast China. Carcinogenesis. 2007, 28: 1074-1078.
Zhang M, Yang ZY, Binns CW, Lee AH: Diet and ovarian cancer risk: a case-control study in China. Brit J Cancer. 2002, 86: 712-717. 10.1038/sj.bjc.6600085.
Ji BT, Chow WH, Yang G, Mclaughlin JK, Zheng W, Shu XO, Jin F, Gao RN, Gao YT, Fraumeni JF: Dietary Habits and stomach cancer in Shanghai, China. Int J Cancer. 1998, 76: 659-664. 10.1002/(SICI)1097-0215(19980529)76:5<659::AID-IJC8>3.0.CO;2-P.
Zhang M, Binns CW, Lee AH: A quantitative food frequency questionnaire for women in southeast China: Development and reproducibility. Asia Pac J Public Health. 2005, 17: 29-35. 10.1177/101053950501700108.
Whitemore AS, Wu-Willians AH, Lee M, Zheng S, Gallagher RP, Jiao DA, Zhou L, Wang XH, Chen K, Jung D: Diet, physical activity and colorectal cancer among Chinese in North American and China. J Natl Cancer Inst. 1990, 82: 915-926. 10.1093/jnci/82.11.915.
Institute of Nutrition and Food Hygiene, Chinese Academic of Preventive Medicine: Food Composition Table (National Representative Values). 1999, Beijing: People's Health Press, 1
Kropp S, Becher H, Nieters A, Chang-Claude J: Low-to-moderate alcohol consumption and breast cancer risk by age 50 years among women in Germany. Am J Epidemiol. 2001, 154: 624-634. 10.1093/aje/154.7.624.
Zhang M, Lee AH, Binns CW: Physical activity and epithelial ovarian cancer risk: a case-control study in China. Int J Cancer. 2003, 105: 838-843. 10.1002/ijc.11165.
Ainsworth BE, Haskell WL, Whitt MC, Irwin ML, Swartz AM, Strath SJ, O'Brien WL, Bassett DR, Schmitz KH, Emplaincourt PO, Jacobs DR, Leon AS: Compendium of physical activities: an update of activity codes and MET intensities. Med Sci Sports Exerc. 2000, 32: S498-S504. 10.1097/00005768-200009001-00009.
Sun CL, Yuan JM, Koh WP, Yu MC: Green tea, black tea and breast cancer risk: a meta-analysis of epidemiological studies. Carcinogenesis. 2006, 27: 1310-1315. 10.1093/carcin/bgi276.
Sun CL, Yuan JM, Koh WP, Yu MC: Green tea, black tea and colorectal cancer risk: a meta-analysis of epidemiologic studies. Carcinogenesis. 2006, 27: 1301-1309. 10.1093/carcin/bgl024.
Zhang M, Zhao X, Zhang X, Holman CDJ: Possible protective effect of green tea intake on risk of adult leukaemia. Brit J Cancer. 2007, 98: 168-170.
Tao D, Hawkins L, Wang H, Langenbrunner J, Zhang S, Dredge R: Fixing the public hospital system in China. 2010, Washington: World Bank, 2: 8-20. Main report
Wang J, Kushner K, Frey JJ, Du XP, Qian N: Primary care reform in the Peoples' Republic of China: implications for training family physicians for the world's largest country. Fam Med. 2007, 39: 639-643.
Hu S, Tang S, Liu Y, Zhao Y, Escobar ML, Ferranti D: Reform of how health care is paid for in China: challenges and opportunities. Lancet. 2008, 372: 1846-1853. 10.1016/S0140-6736(08)61368-9.
Bhattacharyya O, Delu Y, Wong ST, Bowen C: Evolution of primary care in China 1997-2009. Health policy. 2011, 100: 174-180. 10.1016/j.healthpol.2010.11.005.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2288/11/167/prepub
Acknowledgements and funding
The authors acknowledge with gratitude the participation of patients, outpatients and community residents in Shenyang. We are grateful for the collaboration received from the participating hospital and the staff, and from community councils. In particular, we would like to thank Professor Liu Yun-peng and Dr Shi Jing of the First Hospital of China Medical University, for their kind assistance.
This work was supported by the National Health and Medical Research Council of Australia project grant (APP ID 572542). L Li was supported by the Australian Postgraduate Award and the University of Western Australia Establishment Award.
The authors declare that they have no competing interests.
LL collected the data, conducted data analysis, interpretation, and drafted the manuscript. MZ designed the study, directed all aspects of its implementation, including quality assurance and control, assisted data analysis, and reviewed drafts of the article. DH designed the study, supervised all aspects of its implementation, and reviewed drafts of the article. All authors have read and approved the final manuscript.