Hostname: page-component-586b7cd67f-t7czq Total loading time: 0 Render date: 2024-11-27T22:05:04.076Z Has data issue: false hasContentIssue false

Antiretroviral treatment attrition in Swaziland: a population-based study

Published online by Cambridge University Press:  22 August 2016

P. G. KHUMALO
Affiliation:
Antiretroviral Treatment Unit, Mbabane Government Hospital, Mbabane, Swaziland
Y. J. CHOU
Affiliation:
Institute of Public Health, National Yang-Ming University, Taipei, Taiwan
C. PU*
Affiliation:
Institute of Public Health, National Yang-Ming University, Taipei, Taiwan
*
*Author for correspondence: Dr C. Pu, 155 Linong St. Sec 2, Beitou, Taipei, Taiwan, ROC. (Email: [email protected])
Rights & Permissions [Opens in a new window]

Summary

Swaziland has the highest prevalence of human immunodeficiency virus (HIV) in the world. Attrition (loss to follow-up and mortality) in people living with HIV/AIDS (PLWHA) already on treatment is a major challenge, undermining achievements of the antiretroviral treatment (ART) programme in Swaziland. The contributing factors to attrition in the Swazi context are unclear. This study aims to (1) estimate attrition from the ART programme 12 months after ART initiation in Swaziland, and (2) determine the predictors of attrition in PLWHA treated with ART in Swaziland. A retrospective cohort study using national baseline data was conducted. A competing-risk Cox proportional hazard regression was used to determine the predictors of attrition. We estimated 10·3% (95% confidence interval 10·1–10·6) attrition in 16 423 participants that initiated ART in 2012. Attrition was significantly associated with sex, age, district, treatment supporter at initiation, co-infection of HIV and TB, functional status, WHO clinical stage, and ownership of facility. Our study can form a base of policies, plans, and service delivery strategies for preventing and controlling attrition in Swaziland.

Type
Original Papers
Copyright
Copyright © Cambridge University Press 2016 

INTRODUCTION

According to 2014 statistics [1], the national prevalence rate of human immunodeficiency virus (HIV) in Swaziland is the highest in the world (32% in adults aged 18–49 years), making HIV and acquired immune deficiency syndrome (AIDS) the most serious public health concern in the country. Prevalence has peaked at 54% for women aged 30–34 years and 47% for men aged 35–39 years [1]. Increased morbidity and mortality characterize the AIDS epidemic in Swaziland, and life expectancy has dropped from 56 years in 1995 to 49 years in 2013 [2]. The availability of antiretroviral treatment (ART) has transformed the previously fatal HIV infection into a manageable chronic illness [Reference Dlamini3], with over 80% of people living with HIV/AIDS (PLWHA) who are eligible for treatment actually benefiting from treatment in 2012 [Reference Mazibuko4]. To focus on countrywide HIV and AIDS service delivery, Swaziland established a national ART programme in 2004. At inception, one of the programme's goals was to achieve universal access to ART for eligible PLWHA. Strategies for achieving this goal have included making ART accessible free of charge, as well as expanding coverage through a three-phase process. In the first phase, ART service provision was extended from one hospital (Mbabane Government Hospital in Hhohho District) to five major hospitals distributed in all four districts of Swaziland. In the second phase, ART service provision was extended to five health centres and some private clinics. Finally, in the third phase, services were made available in private and public clinics [Reference van Schalkwyk5].

The number of facilities providing ART services increased from one in 2004 to 19 in 2006, reaching 110 sites by the end of 2011. Thirty-eight of these were fully capacitated roll-out sites (i.e. they met the minimum criteria for providing ART services), constituting 18 public or mission facilities and 20 private or non-governmental organization (NGO) facilities. The remaining 72 sites were outreach sites that mainly provided drug refills for ART patients. Twenty-nine of these outreach sites also initiated patients on ART, and were supported by the public or mission hospitals and health centres. The distribution of ART service sites in the four districts of Swaziland were 20 in Hhohho, 42 in Manzini, 25 in Lubombo, and 23 in Shiselweni [6].

In accordance with the Swaziland National ART and Care guidelines of 2014, people living with HIV/AIDS (PLWHA) with a cluster of differentiation 4 (CD4+) count of ⩽350 cells/mm3, or who exhibited stage III or IV symptoms according to WHO clinical staging definitions, were eligible for ART initiation. PLWHA with comorbid conditions such as HIV-related renal failure, hepatitis B virus co-infection, or active tuberculosis (TB) disease were eligible for ART regardless of their CD4+ cell count [1].

National efforts to expand the availability of ART have contributed to lower HIV-related morbidity and mortality. Nevertheless, attrition in PLWHA already on ART is a major challenge undermining ART programme achievements [7]. PLWHA who discontinue treatment are at a high risk of illness and death because of AIDS-related conditions [Reference Miller8].

By 2015, only two studies had addressed factors associated with retention in Swaziland, and neither has been published [Reference Mazibuko4, 9]. Both studies analysed data derived from smaller sample sizes. For example, in 2010, the Ministry of Health of Swaziland in conjunction with the WHO selected only four healthcare institutions (Mbabane Government Hospital, Hlatikulu Government Hospital, Raleigh Fitkin Memorial Hospital, and Good Shepherd Hospital) for their study; the cohort of 946 PLWHA yielded a loss to follow-up (LTFU) rate of 6·76% after a 12-month follow-up period [9]. However, the data on factors associated with LTFU were self-reported, which might have affected LTFU accuracy. According to Mazibuko's estimates [Reference Mazibuko4], attrition rates were 30·9% and 29·1% for doctor- and nurse-managed cohorts, respectively, for a sample size of 871 participants. Because the general population was not used in either study, their results are not generalizable.

Considering the variability of attrition predictors in sub-Saharan Africa and the shortage of attrition studies in Swaziland, we analysed a large population sample to estimate attrition for the ART programme 12 months after ART initiation and to determine the predictors of attrition in PLWHA treated with ART in Swaziland.

MATERIALS AND METHODS

Study design and sample

We conducted a retrospective cohort study using national routinely collected computerized baseline data sourced from the Institute of Health Management. The data were obtained from all medical institutes (both public and private) that provided ART services in Swaziland in 2012. Data were collected from each PLWHA at ART initiation and were aggregated into monthly and quarterly facility reports, which were electronically transmitted to the Institute of Health Management. A total of 16 423 persons (aged ⩾20 years) were newly initiated on ART treatment in 2012. Of these, 2383 were missing data in one or more variables and were therefore excluded from our regression analysis. Consequently, data on 14 040 PLWHA were included in our regression analysis.

Measurements

All participants were followed up until December 2013. Our primary outcome was attrition, which included two mutually exclusive groups: LTFU and mortality. LTFU was defined as a situation whereby a PLWHA has not attended a health facility within 3 months of their most recent scheduled appointment for medication collection, laboratory testing, and/or clinical review and has not been documented as having died. Mortality meant a participant was documented in the facility register as having died. The primary outcome for the survival analysis was retention, which was defined as participants who were alive and in active follow-up at the end of December 2013. Participants who were classified as LTFU during the study but eventually returned before the end of the study were counted as active because their return was a result of defaulter tracing interventions in the programme.

Information collected at baseline included: demographic and social factors (sex, age, district, the availability of a treatment supporter), clinical factors (TB status, functional status, WHO clinical stage), and facility factors (facility type, ownership). Facility types included primary healthcare facilities (clinics and public health units), secondary healthcare facilities (health centres) and tertiary healthcare facilities (hospitals, referral hospitals). Primary healthcare facilities are the first level of contact between individuals and families within the health system. Secondary healthcare facilities are the middle tier of health system, in which PLWHA from primary healthcare are referred to specialists. Tertiary facilities represents the highest level of the health system, in which specialized care is provided usually by referral from the lower levels.

Statistical analysis

PLWHA on ART were censored on 31 December 2013. We separated attrition into two mutually exclusive groups (LTFU and mortality) and calculated overall cumulative attrition, LTFU and mortality rates at 3, 6, 9 and 12 months. Rates were calculated by adding the number of participants who discontinued care at each interval and dividing by the total number of participants enrolled in ART in 2012. This was expressed as percentages with 95% confidence intervals (CIs).

We also performed both univariable and multivariable competing-risk Cox proportional hazard regressions to determine cause-specific predictors of LTFU and mortality after 12 months. A competing-risk Cox regression accounts for the fact that LTFU and mortality are not independent [Reference Fine and Gray10]. To maximize the use of all available data in regression, we modelled missing data for each predictor as a separate category, i.e. ‘missing’. Only risk factors that were significantly associated with each outcome (LTFU and mortality) were included in the multivariable analysis. We reported hazard ratios with 95% CIs and corresponding P values. We used the variance inflation factor (VIF) to check for multicollinearity in predictors. No variable had a VIF value >10, with 1·28 being the mean VIF.

Ethical standards

All procedures performed in this study were in accordance with the ethical standards of the Swaziland Scientific and Ethics Committee, which approved the project.

RESULTS

Characteristics by attrition status

Table 1 shows that LTFU was higher in females (6·8%) compared to males (5·5%), and that mortality was higher in males (5·4% vs. 3·3% in females). Manzini District had the highest (9·2%) LTFU compared to other regions, while Shiselweni District had the highest mortality rate (5·7%). Individuals with advanced HIV/AIDS disease (at WHO stages III and IV) had greater percentages of both LTFU and mortality (7·6% and 7·7%, respectively) compared to those with mild HV/AIDS disease (WHO stages I and II). PLWHA initiating treatment in private facilities showed the least LTFU (2·2%), and those who initiated treatment at NGO facilities were least likely to die (2·5%). The level of missing baseline participant data in this study varied from 2·5% for TB status, 4·4% for functional status, 7·4% for sex, 9·6% for facility ownership, 14·4% for facility type to 14·5% for district.

Table 1. Patient and programme characteristics, stratified by primary outcome (retention and attrition)

LTFU, Lost to follow-up; NGO, non-governmental organization.

The overall attrition, LTFU, and mortality rates are presented in Table 2. After a 12-month follow-up period, overall attrition, LTFU, and mortality rates were 10·3% (95% CI 10·1–10·6), 6·2% (95% CI 5·9–6·5), and 4·1% (95% CI 3·8–4·4), respectively. LTFU was low (0·3%) during the first 3 months, increasing to 4% at 6 months. Mortality, on the other hand, was 2·6% at 3 months, steadily increasing to 4·1% at 12 months.

Table 2. Mortality and LTFU distribution by time on antiretroviral treatment (ART) (N = 16 423)

LTFU, Lost to follow-up.

Factors associated with attrition

In the multivariable competing-risk Cox hazards model (Table 3), age group, district, having a treatment supporter at baseline, TB history, WHO clinical stage, facility type, and facility ownership were independent predictors for LTFU. Sex, age group, TB history, functional status, WHO clinical stage, facility type, and facility ownership were independent predictors for mortality.

Table 3. Univariable and multivariable competing-risk Cox hazards of factors associated with programme attrition in patients who initiated antiretroviral treatmen in 2012 (N = 16 423)

LTFU, Lost to follow-up; HR, hazard ratio; CI, confidence interval; NGO, non-governmental organization.

* Statistically significant P value.

Compared to young adults, middle-aged and older PLWHA had a lower risk of LTFU, as shown by an adjusted hazard ratio (aHR) of 0·63 (95% CI 0·5–0·7) and 0·62 (95% CI 0·5–0·8), respectively. PLWHA who initiated treatment in facilities in Manzini and Shiselweni districts were at a higher risk of LTFU, with an aHR of 4·3 (95% CI 3·1–5·4) and 1·4 (95% CI 1·2–2·0), respectively, compared to those who initiated ART in Hhohho District. The risk of LTFU in PLWHA who had a treatment supporter on ART initiation was lower (aHR 0·5, 95% CI 0·4–0·7) than in those who did not have a treatment supporter. The hazard of LTFU was higher in PLWHA with a positive TB status (had pulmonary TB at initiation) (aHR 3, 95% CI 3·1–3·9) than in those with a negative TB status (did not have pulmonary TB at initiation). PLWHA with advanced HIV/AIDS had a significantly higher risk of LTFU (aHR 1·2, 95% CI 1·1–1·4) than PLWHA with a mild HIV/AIDS disease. Initiating ART in a secondary healthcare facility was associated with a higher risk of LTFU (aHR 4·1, 95% CI 2·9–5·6) than initiating ART in a primary healthcare facility. Moreover, initiating treatment in a tertiary healthcare facility was associated with a higher risk of LTFU (aHR 1·8, 95% CI 1·3–2·3) compared to initiating ART in a primary healthcare facility. PLWHA treated in NGO-owned facilities were 1·7-fold more likely to be classified as LTFU than PLWHA treated in government-owned facilities. Compared to government facility, private facility had a lower hazard of LTFU (aHR 0·4, 95% CI 0·3–0·6).

One major difference between the risks of LTFU and mortality was that sex was a significant predictor of mortality; the risk of mortality for males was 1·4-fold higher than that of females. In addition, in contrast to the risk of LTFU, PLWHA enrolled in privately owned facilities had considerably elevated risks of mortality (aHR 5·2, 95% CI 3·4–8·2) compared to those enrolled in government-owned facilities.

DISCUSSION

We calculated the overall attrition rate of ART in Swaziland and separated attrition into LTFU and mortality. Based on our research, this is the first study to use national data incorporating all medical facilities in Swaziland. After a follow-up period of 12 months, we estimated an overall population attrition of 10·3%. This magnitude of attrition is lower than those recorded in other resource-limited treatment programmes, such as in Mozambique and the Democratic Republic of Congo [Reference Molfino11, Reference Ditekemena12]. It is also lower than those reported in the two studies conducted in Swaziland [Reference Mazibuko4, 9]. We attribute the latter difference to the higher accuracy of our data, because it covered the general population and was collected by health professionals for routine clinical management. Thus, our study is the first, most accurate, measure of attrition in Swaziland.

The difference between Swaziland's attrition and rates from other African settings could be caused by the decentralization of ART services that occurred prior to and during the study period. According to the Joint United Nations Programme on HIV/AIDS (UNAIDS), decentralization coupled with various forms of adherence support have proven effective in improving retention in care [7]. In Swaziland, the model of decentralization in HIV care involved point-of-care diagnostics, task shifting to nurse-led primary-care facilities, community support groups, and integration of services [Reference Turashvili13].

Consistent with our findings, several studies have shown that the risk of mortality is higher in males commencing ART, compared to their female counterparts [Reference Braitstein14, Reference Auld15]. A study of six African countries associated this higher risk to the late initiation of ART by male participants. It is hypothesized that men are more likely to delay access to healthcare because of stigma, male norms that discourage admitting ill-health, and employment responsibilities [16]. A higher background mortality in men in general, regardless of HIV status, could also explain the sex difference in mortality during ART follow-up [Reference Auld15].

Age was a major factor for both mortality and LTFU during the 12-month follow-up period. Similar to findings from previous studies [Reference Boulle17, Reference Nash18], a higher mortality hazard was observed in old-aged participants. One explanation for this association is that CD4 T-cell recovery after starting ART is generally less robust in older PLWHA than in younger participants [Reference Scholten19]. In addition, a greater risk of polypharmacy and its attendant adverse consequences (liver and kidney failure) in elderly participants increases their risk of death after initiating ART. Health promotion strategies, therefore, are vital for geriatric participants. For example, routine kidney and liver function tests are essential for this age group to ensure the early detection of drug toxicities and to facilitate early treatment.

Compared to young adults, both old and middle-aged participants had a lower risk of LTFU. This may be due to their relatively lower mobility [Reference Auld15]. In our analysis, participants who had a treatment supporter on ART initiation were protected from being classified as LTFU. Treatment supporters play a primary role in assisting with PLWHA adherence and retention. This may involve providing emotional support, accompanying the PLWHA for clinic appointments and drug collection, or providing assistance to ensure the participant takes drugs correctly [20]. Because not all PLWHA have relatives and friends available for supporting the continuum of care, an emphasis on community-based support systems can maximize programme retention.

Participants who joined the ART programme with advanced HIV/AIDS had higher risks of LTFU and mortality. Advanced HIV/AIDS is a known risk factor for attrition in PLWHA on ART and is typically a consequence of delaying the start of ART [Reference Molfino11, Reference Koole21]. As in our study, previous studies have found that PLWHA initiating care at PHCFs were more likely to be retained in ART care [Reference Mutasa-Apollo22Reference Fatti, Grimwood and Bock24]. In the context of Swaziland, this finding may be related to the PHCFs performing higher in-community sensitization through health education, the training of PLWHA, the monitoring of community-based activities, and psychological support compared to secondary and tertiary healthcare facilities. The value of sensitization and other community-based interventions in improving adherence has long been recognized [Reference Rueda25].

The highest risk observed in our study was that of dying after initiating treatment at a privately owned facility. Some studies have associated private facilities with suboptimal adherence, lower viral suppression rates, treatment failure, drug resistance, and high mortality [Reference Sulzbach, Wang and Ohanlon26, Reference Ekstrand27]. According to Sulzbach et al. [Reference Sulzbach, Wang and Ohanlon26], a lack of quality HIV services caused by limited training on HIV issues in private sector providers can explain the higher risk of mortality observed in this study.

The following limitations should be considered in interpreting our findings. The retrospective design of the study makes it vulnerable to missing information. Missing data is most likely related to incomplete documentation in facility reports. The problem of missing values underscores the emphasis on attentiveness and complete data collection during treatment initiation in the ART programme in Swaziland. The method we used to deal with missing data may render our results susceptible to bias.

The second limitation is that the true mortality rate may be slightly higher than that reported in this study, because we relied on health professionals to record deaths rather than linking our data to the death registry; the accuracy of our mortality estimates depends on the quality of records provided by these personnel. However, we are confident that estimations from the data collected by health personnel are correct; because the patient tracing programme and patient home visits would have revealed deaths that might not have been otherwise reported by relatives when they occurred.

Our results demonstrate current attrition rates and identify some key factors associated with the risk of attrition for the overall population. Having studied all 110 facilities providing ART services in Swaziland, we can confidently extrapolate our findings to the entire national ART programme. In addition, our findings are likely to yield estimations that are more accurate than previous studies in Swaziland because we analysed secondary data, and the data are not subjected to recall bias. With secondary data, the data collection process is often informed by expertise and professionalism that may not be available to research projects of smaller scale. Our study can form a base of policies, plans, and service delivery strategies for preventing and controlling attrition in Swaziland.

ACKNOWLEDGEMENTS

This study was supported by the Ministry of Science and Technology, Taiwan (grant no. 102-2314-B-010-015-MY3).

DECLARATION OF INTEREST

None.

References

REFERENCES

1. Ministry of Health – Swaziland. Swaziland HIV Incidence Measurement Survey (SHIMS) First Findings Report. 2014 (https://www.k4health.org/sites/default/files/SHIMS_Report.pdf). Accessed 10 February 2015.Google Scholar
2. United Nations Development Programme (UNDP). Sustaining human progress: reducing vulnerabilities and building resilience. HDR. 2014 (http://www.unaids.org/sites/default/files/en/dataanalysis/knowyourresponse/countryprogressreports/2012countries/ce_SZ_Narrative_Report[1].pdf). Accessed 4 January 2015.Google Scholar
3. Dlamini, PS, et al. HIV stigma and missed medications in HIV-positive people in five African countries. AIDS patient care and STDs 2009; 23: 377387.Google Scholar
4. Mazibuko, S. Decentralization of antiretroviral treatment in Swaziland: outcome of nurse initiated versus doctor initiated treatment: University of Pretoria, 2006.Google Scholar
5. van Schalkwyk, C, et al. Outcomes and impact of HIV prevention, ART and TB programs in Swaziland–early evidence from public health triangulation. PLoS ONE 2013; 8: e69437.Google Scholar
6. Ministry of Health – Swaziland. Round of the HIV Sentinel Serosurveillance Survey. 2011 (http://s3.amazonaws.com/zanran_storage/www.fhi.org/ContentPages/2304502.pdf). Accessed September 2014.Google Scholar
7. UNAIDS. Swaziland country report on monitoring the political declaration on HIV and AIDS. 2012.Google Scholar
8. Miller, CM, et al. Why are antiretroviral treatment patients lost to follow-up? A qualitative study from South Africa. Tropical Medicine & International Health 2010; 15: 4854.Google Scholar
9. Ministry of Health-Swaziland, World Health Organization. Assessment of loss-to-follow-up and associated factors among ART clients in Swaziland. 2010 (http://www.stopaidsnow.org/sites/stopaidsnow.org/files/filemanager/TasP_M-Health_MaxART-LTFU_Report_2010_Dec.pdf). Accessed August 2014.Google Scholar
10. Fine, JP, Gray, RJ. A proportional hazards model for the subdistribution of a competing risk. Journal of the American Statistical Association 1999; 94: 496509.CrossRefGoogle Scholar
11. Molfino, L, et al. High attrition among HIV-infected patients with advanced disease treated in an intermediary referral center in Maputo, Mozambique. Global Health Action 2013; 7: 23758.Google Scholar
12. Ditekemena, J, et al. Antiretroviral treatment program retention among HIV-infected children in the Democratic Republic of Congo. PLoS ONE 2014; 9: e113877.Google Scholar
13. Turashvili, M, et al. Decentralisation of HIV\TB care in Shiselweni region of Swaziland: making a difference. 2014. Vienna Evaluation Unit (http://fieldresearch.msf.org/msf/bitstream/10144/325599/1/49_Swaziland%20HIV_TB_decentralization_combined%20report_July%202014.pdf). Accesed October 2014.Google Scholar
14. Braitstein, P, et al. Mortality of HIV-1-infected patients in the first year of antiretroviral therapy: comparison between low-income and high-income countries. Lancet 2006; 367: 817824.Google Scholar
15. Auld, AF, et al. Temporal trends in mortality and loss to follow-up among children enrolled in Côte d'Ivoire's national antiretroviral therapy Programme. Pediatric Infectious Disease Journal 2014; 33: 11341140.Google Scholar
16. Centers for Disease Control and Prevention. Differences between HIV-infected men and women in antiretroviral therapy outcomes-six African countries, 2004–2012. Morbidity and Mortality Weekly Report 2013; 62: 945.Google Scholar
17. Boulle, A, et al. Seven-year experience of a primary care antiretroviral treatment programme in Khayelitsha, South Africa. AIDS 2010; 24: 563–72.Google Scholar
18. Nash, D, et al. Long-term immunologic response to antiretroviral therapy in low-income countries: collaborative analysis of prospective studies: The Antiretroviral Therapy in Lower Income Countries (ART-LINC) Collaboration of the International epidemiological Databases to Evaluate AIDS. AIDS (London, England) 2008; 22: 2291.CrossRefGoogle Scholar
19. Scholten, F, et al. Health and functional status among older people with HIV/AIDS in Uganda. BMC Public Health 2011; 11: 1.Google Scholar
20. World Health Organization. Guidelines on co-trimoxazile prophylaxis for HIV-Related infections among children, adolescents and adults. 2006. (http://www.who.int/hiv/pub/guidelines/ctxguidelines.pdf). Accessed January 2015.Google Scholar
21. Koole, O, et al. Retention in a NGO supported antiretroviral program in the Democratic Republic of Congo. PLoS ONE 2012; 7: e40971.CrossRefGoogle Scholar
22. Mutasa-Apollo, T, et al. Patient retention, clinical outcomes and attrition-associated factors of HIV-infected patients enrolled in Zimbabwe's National Antiretroviral Therapy Programme, 2007–2010. PLoS ONE 2014; 9: e86305.Google Scholar
23. Decroo, T, et al. Lessons learned during down referral of antiretroviral treatment in Tete, Mozambique. Journal of the International AIDS Society 2009; 12: 6.Google Scholar
24. Fatti, G, Grimwood, A, Bock, P. Better antiretroviral therapy outcomes at primary healthcare facilities: an evaluation of three tiers of ART services in four South African provinces. PLoS ONE 2010; 5: e12888.CrossRefGoogle ScholarPubMed
25. Rueda, S, et al. Patient support and education for promoting adherence to highly active antiretroviral therapy for HIV/AIDS. Cochrane Database of Systematic Reviews 2006; 3 (3).Google Scholar
26. Sulzbach, S, Wang, W, Ohanlon, B. Assessing the role of the private health sector in HIV/AIDS service delivery in Ethiopia, 2009.Google Scholar
27. Ekstrand, ML, et al. Suboptimal adherence associated with virological failure and resistance mutations to first-line highly active antiretroviral therapy (HAART) in Bangalore, India. International Health 2011; 3: 2734.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Patient and programme characteristics, stratified by primary outcome (retention and attrition)

Figure 1

Table 2. Mortality and LTFU distribution by time on antiretroviral treatment (ART) (N = 16 423)

Figure 2

Table 3. Univariable and multivariable competing-risk Cox hazards of factors associated with programme attrition in patients who initiated antiretroviral treatmen in 2012 (N = 16 423)