Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Estimating an individual-level deprivation index for HIV/HCV coinfected persons in Canada

  • Adam Palayew,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Epidemiology, Biostatistics & Occupational Health, McGill University, Montreal, Quebec, Canada

  • Alexandra M. Schmidt,

    Roles Conceptualization, Data curation, Investigation, Methodology, Software, Supervision, Writing – review & editing

    Affiliation Epidemiology, Biostatistics & Occupational Health, McGill University, Montreal, Quebec, Canada

  • Sahar Saeed,

    Roles Conceptualization, Data curation, Writing – review & editing

    Affiliation Epidemiology, Biostatistics & Occupational Health, McGill University, Montreal, Quebec, Canada

  • Curtis L. Cooper,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliations Ottawa Hospital Research Institute, Ottawa, Ontario, Canada, Canadian Institutes of Health Research, Canadian HIV Trials Network, Vancouver, British Columbia, Canada

  • Alexander Wong,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Medicine, University of Saskatchewan, Regina, Saskatchewan, Canada

  • Valérie Martel-Laferrière,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Microbiology and Infectious Diseases, Centre de recherche du Centre hospitalier de l’Université de Montréal, Montreal, Quebec, Canada

  • Sharon Walmsley,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliations Canadian Institutes of Health Research, Canadian HIV Trials Network, Vancouver, British Columbia, Canada, University Health Network, University of Toronto, Toronto, Canada

  • Joseph Cox,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliations Epidemiology, Biostatistics & Occupational Health, McGill University, Montreal, Quebec, Canada, Canadian Institutes of Health Research, Canadian HIV Trials Network, Vancouver, British Columbia, Canada, Division of Infectious Diseases and Chronic Viral Illness Service, Department of Medicine, McGill University Health Centre, Montreal, Quebec, Canada

  • Marina B. Klein ,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

    marina.klein@mcgill.ca

    Affiliations Epidemiology, Biostatistics & Occupational Health, McGill University, Montreal, Quebec, Canada, Canadian Institutes of Health Research, Canadian HIV Trials Network, Vancouver, British Columbia, Canada, Division of Infectious Diseases and Chronic Viral Illness Service, Department of Medicine, McGill University Health Centre, Montreal, Quebec, Canada

  • for the Canadian Coinfection Cohort Study Investigators

    Membership of the Canadian coinfection cohort study is provided in the Acknowledgments.

Abstract

Background

HIV-HCV coinfected individuals are often more deprived than the general population. However, deprivation is difficult to measure, often relying on aggregate data which does not capture individual heterogeneity. We developed an individual-level deprivation index for HIV-HCV co-infected persons that encapsulated social, material, and lifestyle factors.

Methods

We estimated an individual-level deprivation index with data from the Canadian Coinfection Cohort, a national prospective cohort study. We used a predetermined process to select 9 out of 19 dichotomous variables at baseline visit to include in the deprivation model: income >$1500/month; education >high school; employment; identifying as gay or bisexual; Indigenous status; injection drug use in last 6 months; injection drug use ever; past incarceration, and past psychiatric hospitalization. We fitted an item response theory model with: severity parameters (how likely an item was reported), discriminatory parameters, (how well a variable distinguished index levels), and an individual parameter (the index). We considered two models: a simple one with no provincial variation and a hierarchical model by province. The Widely Applicable Information Criterion (WAIC) was used to compare the fitted models. To showcase a potential utility of the proposed index, we evaluated with logistic regression the association of the index with non-attendance to a second clinic visit (as a proxy for disengagement) and using WAIC compared it to a model containing all the individual parameters that compose the index as covariates.

Results

We analyzed 1547 complete cases of 1842 enrolled participants. According to the WAIC the hierarchical model provided a better fit when compared to the model that does not consider the individual’s province. Values of the index were similarly distributed across the provinces. Overall, past incarceration, education, and unemployment had the highest discriminatory parameters. However, in each province different components of the index were associated with being deprived reflecting local epidemiology. For example, Saskatchewan had the highest severity parameter for Indigenous status while Quebec the lowest. For the secondary analysis, 457 (30%) failed to attend a second visit. A one-unit increase in the index was associated with 17% increased odds (95% credible interval, 2% to 34%) of not attending a second visit. The model with just the index performed better than the model with all the components as covariates in terms of WAIC.

Conclusion

We estimated an individual-level deprivation index in the Canadian Coinfection cohort. The index identified deprivation profiles across different provinces. This index and the methodology used may be useful in studying health and treatment outcomes that are influenced by social disparities in co-infected Canadians. The methodological approach described can be used in other studies with similar characteristics.

Introduction

HIV/HCV co-infected individuals are often marginalized and of lower socio-economic status [1]. Routes of HIV and HCV transmission overlap and are exacerbated by up-stream social determinants of health such as injection drug use [2]. The concept of deprivation was defined by Townsend in 1987 as “a state of observable and demonstrable disadvantage relative to the local community or the wider society or nation to which that individual, family, or group belongs” [3]. However, deprivation is a latent construct and cannot be directly measured nor estimated using a single variable. Therefore, we must use several variables but these variables may be highly collinear. To measure deprivation, studies aim to isolate the common variation inherent to several variables and summarize it with numerical scores. The most common measure of deprivation in Canada is the Pampalon Index which was developed in Quebec [4]. The Pampalon index is a 6-item index based on an aggregate of a census area, described by the proportion of the population that are: employed, with high school education or more, separated or widowed or divorced, single parent families, live alone, and the average household income for a given census area [5]. The Pampalon index has been used to measure deprivation in other provinces in Canada, in HIV mono-infected, HCV mono-infected, and HIV-HCV co-infected individuals in British Columbia with administrative data [57].

While the Pampalon index is useful in describing the deprivation of a census area, it is limited in its ability to describe deprivation among individuals in observational studies. First, it is an area-based index. Area-based indices assume that deprivation is homogeneous across the entire geographic area of measurement, which can lead to a loss of individual-level information by grouping people together who have different individual needs [8, 9]. Second, the Pampalon index was primarily designed and validated to measure deprivation in the general population, which may not be applicable to more marginalized populations such as people living with HIV and HCV [10]. Variables that help provide good differentiation in the Pampalon index such as the proportion of people who are widowed, separated or divorced may be useful to distinguish between members of the general population, but may not differentiate well for more marginalized populations where marriage is less common. On the other hand, characteristics such as injection drug use and incarceration history may be more useful in distinguishing levels of deprivation within HIV-HCV co-infected populations.

The primary aim of this study was to develop and estimate an individual-level model-based deprivation index that considers social, material, and lifestyle variables for members of the Canadian co-infection cohort (CCC) using an item response theory (IRT) model. Secondly, we explored the differences in the index among four regions in Canada, and how the parameters of the index varied by the different regions. Additionally, to provide an example of how the proposed index could be used, we looked at an outcome potentially impacted by deprivation and compared the fitted model with only the index as a covariate against the model that had all the variables that make up the index included as covariates.

Methods

Study setting and participants

The CCC is a publicly funded open prospective cohort of approximately 2000 HIV/HCV co-infected individuals recruited from 18 centers across Canada. Since 2003, HIV+ adults with evidence of HCV infection (antibody positive) have been eligible to participate. Data on socio-demographic, behavioral, and clinical characteristics and laboratory data are collected bi-annually. Details on HIV and HCV treatments/responses, clinical diagnoses, and deaths are collected using standardized case report forms [11]. The data we used for this study spanned all first visits for the CCC from 2003 to the third quarter of 2018.

Development of the index

Variable pool development: We considered 19 variables potentially associated with material and social deprivation for inclusion in the model. For all analyses, all the variables were dichotomized. For our model, it was important that all the variables had a clear hierarchy, since the model assumes that all levels of every variable are monotonically increasing [12]. The variables considered for inclusion fell into two categories: 1) the 6 Pampalon variables and 2) social and behavioral variables associated with deprivation in the literature.

Pampalon variables.

We adapted five out of six variables used to estimate the Pampalon index. The variables were: high school education or less, living alone (yes/no), relationship status (single/married), employed (yes/no), and income dichotomized into equal to or more than $1500 a month and less than $1500 (based on low income threshold for a single person in Canada in 2008, which was the last year an income poverty line was used in Canada) [13]. The proportion of single parent homes was the one variable from the Pampalon index that did not have corresponding individual level data collected by the CCC.

Variables selected from the literature.

These variables were selected after performing a literature search in Medline, EMBASE, and PsychInfo focusing on articles reporting on social/behavioral factors associated with HIV/HCV co-infection. Variables identified from the literature were then matched to individual level variables available in our cohort, the reference group for these studies were either HIV or HCV mono-infected [7, 1419]. We considered 14 dichotomous (yes/no) variables, which included: identifying as gay or bisexual, vulnerably housed (defined as being homeless or living in a shelter), snorted drugs in the last 6 months, injected drugs ever, injected drugs in last 6 months, ever incarcerated, sex client in the last 6 months, sex work in the last 6 months, sex client ever, sex work ever, sexually transmitted infection in the last 6 months, depression, hospitalized at a psychiatric institution, and diagnosed with schizophrenia.

Variables selection process.

Variables were selected for inclusion into the final model using a pre-defined variable selection process. For each group of variables, all possible combinations of two-way χ2 significance tests (set at an α of 0.05), and the grouping between variables in output from multiple correspondence analysis (MCA) were evaluated [20].

Proposed model

We fit an IRT model with all complete cases of the final selected variables [21]. We used an IRT model because it allowed us to estimate an individual-level deprivation index for every member of our cohort and evaluate the relative importance of every item in determining an individual’s score. We fit two models: one that ignored groupings and another that grouped observations hierarchically based on an individual’s province of recruitment. The purpose of allowing the items of the model to vary by province is that there are known differences in sociodemographic factors in participants residing from different provinces (such as the prevalence of Indigenous populations, income, drug use) that could influence an individual’s level of deprivation [22].

Our proposed model follows Quinn 2004 [23]. Let xijk be the response to the j-th item for person i from province k. The hierarchical IRT model assumes there is a latent continuous variable such that if and if xijk = 1 then , with , with ijk ~ N(0,1). The parameter αjk is the response parameter and the parameter βjk is the discrimination parameter. The parameter θik is the index associated with subject i, and is what we are considering as our individual-level deprivation index. Note that we allow αjk and βjk to vary with the province, a priori. This allows the data to drive the inference procedure and point out if there are differences at the provincial level in the way the variables influence the estimation of the index.

Item characteristic curves.

An item characteristic curve (ICC) is used to evaluate the probability of responding yes to a given item based on an individual’s latent score [12]. The ICC is defined for the two-parameter IRT model as . The response parameter (αjk) is the location of the inflection point of an ICC, and the discriminatory parameter (βjk) is the slope of an ICC around its point of inflection [24]. This allowed us to examine how the different items discriminate as a function of the subject parameter (θik), which is our index. We further examined how the ICCs varied in the different provinces.

Model selection.

We compared the model that considered provincial grouping (assuming αjk and βjk varied by province) to the one that did not (assuming αj and βj common across provinces). We calculated the widely-applicable information criteria (WAIC) to evaluate which model, between the fitted ones, best fits the data [25]. The model with the lowest WAIC was chosen as the index. A difference of 5 in WAIC was considered to be substantial [26].

Association between the index and not attending a second visit.

Secondary to estimating the index, we examined the association of the deprivation index at baseline and not attending a second cohort visit (scheduled every 6 months)—as a proxy for disengagement in care. This was to provide an example of a potential use of the estimated index, but not to accurately model the probability of attending a second cohort visit. We fitted a model that considered only our estimated index as a covariate, as it summarizes through a single number the characteristics of the individual, and another model that considered as covariates each of the variables that are considered in the estimation of the index. All models were fitted using a logistic regression. The WAIC was used to determine the model that best fit the data. We chose to use logistic regression because the time between the baseline and when the outcome is assessed is fixed for all individuals obviating the need for time to event analysis. This analysis was not to make predictions for future outcomes in the cohort.

Statistical analyses.

All analyses, graphics, and proposed model fits were done following the Bayesian paradigm using the software R [27] and the Stan programming language [28]. Functions from the package RStan, provide samples from the posterior distribution of the parameters in the model via Hamiltonian Monte Carlo with a No U-Turn Sampler [29]. The 95% posterior credible intervals for each parameter is presented after 5000 iterations with a thinning interval of 1 and 2500 burn-in iterations for 4 chains. Convergence of the model was determined through the examination of trace plots and R-hat values. The WAIC for each model was calculated using the WAIC function from the loo Package [30].

Availability of data and materials

A fictitious sample generated based on the posterior summaries of the parameters is also available, but the original data is not publicly available as it requires approval from the CCC to be shared.

Ethics statement

The study has been approved by research ethics boards at each of the participating institutions as follows: community advisory committee of the CIHR-Canadian HIV Trials Network, the Biomedical B Research Ethics Board of the McGill University Health Centre (2006–1875 and BMB-06-006t), the UBC-Providence Health Care Research Ethics Board (H16-02853), the Institutional Review Board Services, Regina Qu’Appelle Health Region Research Ethics Board (REB-14-70), the Conjoint Health Research Ethics Board of the University of Calgary (REB15-0848 REN5), the Nova Scotia Health Research Ethics Board (1022541), the Windsor Regional Hospital Research Ethics Board (07-122-17), the Veritas Independent Review Board (7035–11:502-04-2020), the Hamilton Integrated Research Ethics Board (06–397), the Comité d’éthique de la recherche du CHUM (2003–1582, SL 03.008-BSP), the Comité d’éthique de la recherche du CHU de Québec-Université Laval (2012–876, C11-12-153), the Sunnybrook Health Sciences Centre Research Ethics Board (252–2008), the Research Ethics Board of Health Sciences North (605), the University Health Network Research Ethics Board (06-0629-BE), the Ottawa Health Science Network Research Ethics Board (2007229-01H) and the Biomedical Research Ethics Board (12–178). The study was conducted according to the Declaration of Helsinki. Patient records/information was anonymized and de-identified prior to analysis. All patients gave written informed consent before undergoing an initial evaluation and were followed at study visits every six months.

Results

Overall, 1842 individuals were enrolled in the CCC between its inception and the third quarter of 2018. At cohort entry, the median age was 54 years, 28% were female, 24% Indigenous and 82% reported ever using injection drugs. The median CD4 count was 410 cells/μL, 31 individuals had HIV RNA greater than or equal to 50 copies/mL, and 91% were receiving antiretroviral therapy; 46% had ever received HCV treatment.

Variable selection

We initially considered 19 variables (Table 1), 5 from the Pampalon group and 14 from the literature group. After applying the selection process, 9 variables, which are stared in Table 1, were selected and used to fit the final IRT model. The output for the χ2 tests and the MCA for both groups of variables are presented in S1 File along with the final list of candidate variables that were considered for inclusion into the index. We analyzed 1574 complete cases of a possible 1842 enrolled participants for the 9 variables included into the model. Participants from Nova Scotia (n = 13) and Alberta (n = 47) were excluded since there were not enough individuals to evaluate the group level parameters in the hierarchical model.

thumbnail
Table 1. Characteristics of the participants at baseline according to variables that were considered for inclusion into the model stratified by province.

https://doi.org/10.1371/journal.pone.0249836.t001

Model based index results

The simple model (which assumes αj and βj) yielded a WAIC value of 12413.7 (standard error, 116.0) compared to a WAIC of 11960.6 (standard error, 117.0) for the hierarchical model (which assumes αjk and βjk). Since the model that allowed parameters to differ across the provinces fit the data better, we only present the results from that model. Boxplots for the individual-level parameter had similar values (Fig 1) suggesting our index has a similar distribution across the four provinces.

thumbnail
Fig 1. Boxplot of index scores by province.

Posterior summaries of the individual scores for each participant presented as boxplots.

https://doi.org/10.1371/journal.pone.0249836.g001

However, we found differences in how the response and discriminatory parameters varied by province (ICC curves; Fig 2). The ICC curves show the probability of responding positively to the question (y-axis) based on the value of the deprivation index (x-axis) and the different colored lines show the variation in response by province (Fig 2).

thumbnail
Fig 2. Variation of IRT by variable for the four provinces.

Variation of the item characteristic curves for each province by variable. The black line represents the mean value for all the indices, the blue line for Quebec, red line for Ontario, yellow line for British Columbia, and the green line for Saskatchewan. The x−axis shows the value of the estimated deprivation index and the y−axis the probability of responding positively to the question.

https://doi.org/10.1371/journal.pone.0249836.g002

For example, we found that history of incarceration and whether someone identified as gay or bisexual were more important in distinguishing between individuals in the province of Ontario than in the other three provinces. However, there were also some variables that showed little regional variation in the discriminatory parameter like unemployment, which had similar ICCs across the four provinces (Fig 2). The province that differed most from the others was Saskatchewan which had the highest response parameters for people who were Indigenous, do not identify as gay or bisexual, and for education high school or less due to a higher prevalence of those subgroups relative to the other provinces. This reflects that the population of coinfected participants from Saskatchewan differs demographically from other provinces, having more Indigenous people, fewer gay and bisexual males, and a lower proportion of individuals with post-secondary education compared to the rest of the provinces.

Association of the index and non-attendance at a second visit

Of 1537 eligible participants, 457 (30%) failed to attend a second visit. In the univariable logistic regression, a one-unit increase in the index score was associated with a 17% odds (95% posterior CrI: 2%, 34%), of not attending a second visit. Many of the covariates that make up the index (e.g. non-gay bisexual men who have sex with men (non-GBMSM), and unemployed) were also associated with not attending a visit in univariable analyses, but these estimates were not very precise (Table 2). The multiple regression analysis demonstrates that there is potential for individual components contained in the index to be confounded by, or correlated with, each other if they are included as independent covariates to adjust for deprivation. As expected, including all the covariates in a single multiple regression model reduced precision of individual estimates as the variables are all highly correlated to one another.

thumbnail
Table 2. Univariable and multivariable regressions.

The point estimates are the mean of the posterior summary and the 95% credible interval of the exponentiated odds ratios.

https://doi.org/10.1371/journal.pone.0249836.t002

Discussion

We developed an individual-level index incorporating social, material, and lifestyle components for Canadians who are coinfected with HIV/HCV. We examined how the characteristics of individuals differed for different values of the estimated index. Individual components of the model varied by province—reflecting local HIV-HCV epidemiology–while values of the index were distributed similarly across the country. High values of this index identified participants with more vulnerable profiles. We illustrated how the index could be used to evaluate an outcome that should be affected by deprivation–higher values of the index were associated with lack of attendance at a second cohort visit.

The index has several advantages over the most commonly used deprivation index, the Pampalon index [6]. Unlike the Pampalon index, an area based index, that values of our index represents an individual-level measure of deprivation compared with assigning everyone in the same census area the same value of deprivation [31]. This allows us to assign an individual a deprivation value that is reflective of their personal characteristics [32]. Additionally, our index included variables which are particularly relevant for deprivation in HIV-HCV coinfected individuals and have not been used in any previous deprivation indices to our knowledge such as incarceration history [3, 4, 10, 3338].

As a single summary measure of deprivation, the index offers value in that it captures the effects of multiple variables that contribute to deprivation simultaneously, preserving degrees of freedom and thus offering improved precision. As a tool for adjusting for deprivation in future studies, our index summarizes the shared variation among the variables eliminating the need for including multiple covariates to control for deprivation and therefore is more parsimonious, which is beneficial. We illustrate this when evaluating a potentially relevant outcome, disengagement from care. The model with the index alone provided a much lower WAIC compared to the multivariable model with all variables contained in the index assuming a difference in WAIC of greater than 5 between two models is considered as a substantial difference [26]. Additionally, the individual estimates from the model that contained all the covariates led to wider ranges of the 95% posterior credible intervals, indicating that the posterior variance of the coefficients was greater when compared to the coefficients from the univariate model counterparts. This is probably because these covariates are correlated (Table 2). Our index naturally accounts for this correlation among these covariates, summarizing the characteristics of an individual into a single number. This allows for us to be able to compare individuals on a continuous range of numbers based off a mathematical transformation of 9 dichotomous variables.

Limitations

The main limitation of our study was that the items selected for inclusion in the index were partly subjective. Experts may differ on what variables they believe should be included in the model. In our model, we included the 9 variables using a variable selection process, which considered the peer-reviewed literature, statistical associations amongst the candidate variables and a priori expert opinion. Another limitation is that the index is based exclusively on the data collected in our cohort, therefore, our score may not generalize well to coinfected individuals in other settings. However, the CCC is a very diverse cohort and is representative of a wide range of individuals across Canada in multiple provinces [11]. Additionally, this index cannot be used directly in another cohort. The process would need to be replicated by fitting the IRT model to the other cohort’s data. However, our index identified distinct characteristics that differed amongst the four provinces reflecting local epidemiology, despite their diverse socio-demographics, it is reasonable to assume that if we were to refit this model in another high-income country similar to Canada, like Australia, that it will be able to be flexible in accounting for group differences that may exist, provided that the same variables are available and data are robust. We did not impute data for missing variables as missing data were rare. We acknowledge that using complete cases could potentially introduce bias if data were not missing at random. From a Bayesian point of view, a model for the missing mechanism could be jointly fitted with the IRT model and the missing values could have been estimated in situations with higher rates of missing data.

Furthermore, we dichotomized the variables that were included in the model and this could lead to a loss of information. Although there are more levels in many of our variables we decided to dichotomize them due to concerns about sample size in small sub categories. Moreover, a strict monotonically increasing order is required for all variables in the model. When adding more categories, we need to make sure that there is a clear hierarchy among them, which is not possible for some variables such as Indigenous status, or non-GBMSM, which would need to have more than 2 categories. Finally, the outcome that we used was non-attendance of a second cohort visit, and acknowledge that this may not accurately capture disengagement from care as individuals could attend a subsequently scheduled appointment. Gaps in attendance in care is common in more marginalized populations who have difficulty accessing the healthcare system for a variety of reasons and may only attend erratically.

Conclusions

We developed and estimated a novel individual-level index that measures deprivation in a national representative cohort of HIV/HCV coinfected Canadians. Due to the hierarchical nature of the model, variables were able to contribute differently depending on the province each individual was from and reflected the local epidemiology. The proposed index is a valuable tool to understand the profile of any individual that is part of the CCC study. In particular, our proposed approach is able to provide an estimate for the index of new participants, potentially allowing health workers to better understand the profile of a new individual with respect to others in the study. If validated in other settings, it may be a valuable tool for studying health outcomes in HIV-HCV coinfection in observational studies.

Supporting information

S1 File. Selection criteria output.

This contains supplemental output for the selection criteria that were used to select variables into the final model.

https://doi.org/10.1371/journal.pone.0249836.s001

(PDF)

Acknowledgments

We would like to acknowledge the participants of the Canadian Co-Infection Cohort (CTN222), as well as the cohort co-investigators, Lisa Barrett (Dalhousie Univeristy), Jeff Cohen (Windsor Regional Hospital Metropolitan Campus), Brian Conway (Vancouver Infectious Disease Research and Care Centre), Curtis Cooper (University of Ottawa), Pierre Côté (Clinique Médiciale du Quartier Latin), Joseph Cox (Montreal University Health Centre), M. John Gill (Southern Alberta HIV Clinic), Shariq Haider (McMaster University), Mark Hull (BC Centre for Excellence in HIV/AIDS), Valérie Martel-Laferrière (Université de Montréal), Julio Montaner (BC Centre for Excellence in HIV/AIDS), Erica E. M. Moodie (McGill University), Neora Pick (Oak Tree Clinic), Danielle Rouleau (Université de Montréal), Aida Sadr(BC Centre for Excellence in HIV/AIDS), Steve Sanche (University of Saskatchewan), Roger Sandre (Hôpital Régional de Sudbury Regional Hospital), Mark Tyndall (University of British Columbia), Marie-Louise Vachon (Centre Hospitalier Universitaire de Québec), Sharon Walmsley (University Health Network), Alexander Wong (University of Saskatchewan), and the study coordinators/nurses.

References

  1. 1. Klein MB, Rollet KC, Saeed S, Cox J, Potter M, Cohen J, et al. HIV and hepatitis C virus coinfection in Canada: challenges and opportunities for reducing preventable morbidity and mortality. HIV Med. 2013;14(1):10–20. pmid:22639840
  2. 2. Meissner EG. Update in HIV-hepatitis C virus coinfection in the direct acting antiviral era. Curr Opin Gastroenterol. 2017;33(3):120–7. pmid:28234770
  3. 3. Townsend P. Deprivation. Journal of Social Policy. 1987;16(2):125–46.
  4. 4. Pampalon R, Raymond G. A Deprivation Index for Health and Welfare Planning in Quebec. Chronic Diseases in Canada. 2000;21(3):10. pmid:11082346
  5. 5. Pampalon R, Hamel D, Gamache P, Philibert MD, Raymond G, Simpson A. An Area-based Material and Social Deprivation Index for Public Health in Québec and Canada. Can J Public Health. 2012;103(2). pmid:23618066
  6. 6. Pampalon R, Hamel D, Gamache P. Health inequalities in urban and rural Canada: comparing inequalities in survival according to an individual and area-based deprivation index. Health Place. 2010;16(2):416–20. pmid:20022551
  7. 7. Butt ZA, Shrestha N, Wong S, Kuo M, Gesink D, Gilbert M, et al. A syndemic approach to assess the effect of substance use and social disparities on the evolution of HIV/HCV infections in British Columbia. PLoS One. 2017;12(8):e0183609. pmid:28829824
  8. 8. Pampalon R, Hamel D, Gamache P. A comparison of individual and area-based socio-economic data for monitoring social inequalities in health. In: Canada S, editor. Health Reports: Statistics Canada; 2009.
  9. 9. Piantadosi S, Byar DP, Green SB. The ecological fallacy. American journal of epidemiology. 1988;127(5):893–904. pmid:3282433
  10. 10. Labbe E, Blanquet M, Gerbaud L, Poirier G, Sass C, Vendittelli F, et al. A new reliable index to measure individual deprivation: the EPICES score. Eur J Public Health. 2015;25(4):604–9. pmid:25624273
  11. 11. Klein MB, Saeed S, Yang H, Cohen J, Conway B, Cooper C, et al. Cohort profile: the Canadian HIV-hepatitis C co-infection cohort study. Int J Epidemiol. 2010;39(5):1162–9. pmid:19786463
  12. 12. Baker FB, Kim S. The Basics Of Item Response Theory: Springer; 2017.
  13. 13. Income Statistics Division. Low income cut-offs for 2008 and low income measures for 2007. In: Division IS, editor.: Statistics Canada; 2008.
  14. 14. Abadie R, Welch-Lazoritz M, Bilal K, Dombrowski K. Social Determinants of HIV/HCV Co-Infection: A case Study from People Who Inject Drugs in Rural Puerto Rico. Addict Behav Rep. 2017;5:29–32. pmid:28983502
  15. 15. Kakchapati S, Maharjan M, Rawal BB, Dixit SM. Social determinants and risk behaviors associated with prevalent Hepatitis C and HIV/HCV co-infection among male injection drug users in Nepal. Arch Public Health. 2017;75:39. pmid:28878895
  16. 16. March JC, Oveido-Joekes E, Romero M. Factors associated with reported hepatitis C and HIV among injecting drug users in ten European cities. Enferm Infecc Microbiol Clin. 2007;25(2):7. pmid:17288906
  17. 17. Pantalone DW, Hessler DM, Bankoff SM, Shah BJ. Psychosocial correlates of HIV-monoinfection and HIV/HCV-coinfection among men who have sex with men. J Behav Med. 2012;35(5):520–8. pmid:21947975
  18. 18. Pereira M, Fialho R, Canavarro MC. Prevalence and correlates of emotional distress in HIV/HCV coinfection. AIDS Care. 2014;26 Suppl 1:S56–64. pmid:24720733
  19. 19. Rourke SB, Sobota M, Tucker R, Bekele T, Gibson K, Greene S, et al. Social determinants of health associated with hepatitis C co-infection among people living with HIV: results from the Positive Spces, Health Places study. Open Medicine. 2011;5(3):12.
  20. 20. Greenacre M. Correspondence analysis in practice. 3 ed. Boca Raton Florida: CRC Press; 2017.
  21. 21. Johnson VE, Albert JH. Ordinal data modeling: Springer Science & Business Media; 2006.
  22. 22. Martel L. Recent changes in demographic trends in Canada: Statistics Canada Ottawa; 2015.
  23. 23. Quinn KM. Bayesian Factor Analysis for Mixed Ordinal and Continuous Responses. Political Analysis. 2004;12(4):338–53.
  24. 24. Reise SP. Item Response Theory. The Encyclopedia of Clinical Psychology 2014. p. 1–10.
  25. 25. Watanabe S. Asymptotic Equivalence of Bayes Cross Validation and Widely Applicable Information Criterion in Singular Learning Theory. Journal of Machine Learning Research. 2010;11:24.
  26. 26. Reich B, Ghosh S. Bayesian statistical methods: CRC Press, Taylor & Francis Group; 2019.
  27. 27. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2018.
  28. 28. Gelman A, Lee D, Guo J. Stan: A Probabilistic Programming Language for Bayesian Inference and Optimization. Journal of Educational and Behavioral Statistics. 2015;40(5):530–43.
  29. 29. Stan Development Team. RStan: the R interface to Stan. 2018.
  30. 30. Vehtari A, Gelman A, Gabry J. Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Statistics and Computing. 2016;27(5):1413–32.
  31. 31. Braveman PA, Cubbin C, Egerter S, Chideya S, Marchi KS, Metzler M, et al. Socioeconomic Status in Health Research One Size Does Not Fit All. Jama. 2005;294(22):10.
  32. 32. Rehkopf DH, Haughton LT, Chen JT, Waterman PD, Subramanian SV, Krieger N. Monitoring Socioeconomic Disparities in Death: Comparing Individual-Level Education and Area-Based Socioeconomic Measures. American Journal of Public Health. 2006;96(12):4. pmid:16809582
  33. 33. Ribeiro AI, Mayer A, Miranda A, Pina MF. The Portuguese Version of the European Deprivation Index: An Instrument to Study Health Inequalities. Acta Med Port. 2017;30(1):17–25. pmid:28501033
  34. 34. Guillaume E, Pornet C, Dejardin O, Launay L, Lillini R, Vercelli M, et al. Development of a cross-cultural deprivation index in five European countries. J Epidemiol Community Health. 2016;70(5):493–9. pmid:26659762
  35. 35. Krieger N, Waterman PD, Spasojevic J, Li W, Maduro G, Van Wye G. Public Health Monitoring of Privilege and Deprivation With the Index of Concentration at the Extremes. Am J Public Health. 2016;106(2):256–63. pmid:26691119
  36. 36. Bell N, Hayes MV. The Vancouver Area Neighbourhood Deprivation Index (VANDIX): A Census-based Tool for Assessing Small-area Variations in Health Status. Canadian Journal of Public Health. 2012;103(2):6. pmid:23618068
  37. 37. Havard S, Deguen S, Bodin J, Louis K, Laurent O, Bard D. A small-area index of socioeconomic deprivation to capture health inequalities in France. Soc Sci Med. 2008;67(12):2007–16. pmid:18950926
  38. 38. Ebert U, Moyes P. An axiomatic characterization of Yitzhaki’s index of individual deprivation. Economic Letters. 1999;68(1):8.