- Open Access
Measuring gender when you don’t have a gender measure: constructing a gender index using survey data
© The Author(s). 2016
Received: 7 October 2015
Accepted: 19 May 2016
Published: 28 May 2016
Disentangling the impacts of sex and gender in understanding male and female differences is increasingly recognised as an important aspect for advancing research and addressing knowledge gaps in the field of work-health. However, achieving this goal in secondary data analyses where direct measures of gender have not been collected is challenging. This study outlines the development of a gender index, focused on gender roles and institutionalised gender, using secondary survey data from the Canadian Labour Force survey. Using this index we then examined the distribution of gender index scores among men and women, and changes in gender roles among male and female labour force participants between 1997 and 2014.
We created our Labour Force Gender Index (LFGI) using information in four areas: responsibility for caring for children; occupation segregation; hours of work; and level of education. LFGI scores ranged from 0 to 10, with higher scores indicating more feminine gender roles. We examined correlations between each component in our measure and our total LFGI score. Using multivariable linear regression we examined change in LFGI score for male and female labour force participants between 1997 and 2014.
Although women had higher LFGI scores, indicating greater feminine gender roles, men and women were represented across the range of LFGI scores in both 1997 and 2014. Correlations indicated no redundancy between measures used to calculate LFGI scores. Between 1997 and 2014 LFGI scores increased marginally for men and decreased marginally for women. However, LFGI scores among women were still more than 1.5 points higher on average than for men in 2014.
We have described and applied a method to create a measure of gender roles using survey data, where no direct measure of gender (masculinity/femininity) was available. This measure showed good variation among both men and women, and was responsive to change over time. The article concludes by outlining an approach to use this measure to examine the relative contribution of gender and sex on differences in health status (or other outcomes) between men and women.
Better understanding and accounting for male and female differences has been gaining attention in many health-related research areas [1–3]. In the area of work and health, men and women differ in their work exposures and work-related health conditions. In addition, the relationships between work exposures and health outcomes may also differ for men and women. These male/female differences can be due to sex – referring to biological differences between men and women – or gender – referring to social differences between men and women .
It is increasingly recognised that both sex and gender matter in understanding the relationships between working conditions and health outcomes, and that research that fails to take sex and gender into account is limited in both quality and applicability [4, 5]. Stratifying analyses to examine the relationships between work and health separately for men and women has been proposed as one approach to better account for sex and gender . However, it is recognised that this approach does a better job of understanding “sex” differences than it does in understanding “gender” differences [7, 8]. Furthermore, sex and gender often interact, suggesting that differences between men and women might be due to a combination of both biological (sex) and social/cultural (gender) factors. To develop a better understanding of the relative contribution of each of these aspects requires measures of both sex and gender to be included in analyses [5, 8].
Measuring gender and sex can take different forms depending on the way data are being collected. When conducting primary data collection for quantitative studies researchers have the options to include measures of gender, such as the Bem-Sex-Role-Inventory (BSRI)  which asks participants to self-identify with personal traits, or the Masculine Gender Role Stress scale . However, are there options for researchers to measure gender if such scales are not present in existing data?
The concept of gender diagnosticity was first introduced by Lippa and Connelly  to estimate the probability of being male or female, based on some gender-related diagnostic indicator. In their original study, Lippa and Connelly used occupational preference ratings as a measure of gender-diagnosticity, finding that these preferences were distinct from responses to the Personal Attributes Questionnaire (PAQ)  and the BSRI . They also reported that occupational preference was more predictive of being a man or a woman than either the BRSI or PAQ indices . This suggests that a gender diagnostic approach offers an alternative method to measure social differences between men and women based on their roles and preferences, compared to indices such as the PAQ and BRSI that are based on gender stereotypes .
A decade later Lippa and colleagues used this same approach (occupational preferences) to examine the role of sex and gender on mortality . In this study they found masculinity (as measured by occupational preferences) was predictive of mortality among both men and women, resulting in the highest mortality rate being observed among the most masculine men, and the lowest mortality rate observed among the most feminine women . The objective in a gender-diagnostic approach is to identify indicators that best differentiate between different gender-based groups. As Lippa and Connelly noted in their original paper, multiple indicators would ideally be used to form a gender index, resulting in a more reliable scale .
Most recently a variation on this approach has been used to examine the impact of sex and gender on cardiovascular risk factors among individuals with premature acute coronary syndrome . In this study the gender index was comprised of information on whether the respondent was the primary earner in their household; their personal income; the number of hours and responsibility for housework; and level of stress at home – along with measures of masculinity and femininity from the BSRI . Similar to the previous gender diagnostic studies, this paper found that both sex and gender were important in predicting many cardiovascular risk factors, but that the gender score was generally more important than sex (male/female) in predicting risk in multivariable models .
The preceding studies relied on primary data collection where the concept of a gender index or gender-diagnostic approach was part of the study design. This is not always feasible in population-health or health services research studies that rely on existing surveys and administrative health records, despite the growing body of literature that indicated that gender differences matter to primary prevention and health care practices [15–17].
In this paper the aim is to develop a gender-index using the Canadian Labour Force Survey (LFS). The LFS was selected because it has a number of questions that are commonly available in other data sources, and a gender-index using this data may be readily applied (and modified or expanded upon) to other secondary data sources. We use data from the 1997 and 2014 Labour Force Surveys to accomplish three objectives: to develop a gender-index using existing population health survey data; to examine the distribution of our gender index across males and females (i.e. to ensure that it measured a separate concept to sex); and to examine if there have been changes in gender roles (as measured by the index) among male and female labour market participants between 1997 and 2014. We then discuss how this index (or a similarly constructed index) could be used in research that exploits secondary data to better understand the relative contribution of aspects of gender and sex in male/female differences in health outcomes.
This paper uses secondary data from Statistics Canada’s LFS. The LFS is a monthly survey carried out by Statistics Canada with the objective of providing information on trends in labour market participation and hours of work across major occupational and industrial sectors in Canada . The LFS surveys approximately 56,000 Canadian households per month. Households remain in the sample for six consecutive months, with one sixth of the sample rotated out, and replaced by a new group of households representing one sixth of the sample each month. The target population for the LFS is the civilian, non-institutionalised population 15 years of age and over residing in all of Canada’s provinces and territories. Persons living on Aboriginal reserves, full-time members of the Canadian Armed Forces, and the institutionalised population are excluded from coverage, as are households in extremely remote areas. Statistics Canada estimates these groups represent less than 2 % of the Canadian population aged 15 and over, and that the LFS is representative of its target population . For the purpose of this analysis, the Public Use files from the 1997 and 2014 Labour Force Surveys were used through Statistics Canada’s Data Liberation Initiative . The year 1997 was chosen as the start point for the analysis, as the questions asked in the LFS changed in this survey year. For each survey cycle, the analysis was restricted to respondents who were currently working for pay or profit in the past month, excluding unpaid family workers, regardless of the number of hours worked.
Labour Force Gender Index (LFGI)
Gender is a multidimensional construct that includes four dimensions: gender roles (behavioural norms applied to men and women); gender identity (how an individual sees themselves on the male/female continuum); gender relationships (how individuals are treated by others based on their ascribed gender); and institutionalized gender (how power and influence are distributed differently among men and women) . The LFGI constructed from the LFS focused primarily on the dimensions of gender roles and institutionalised gender among labour force participants. Given the data available, the LFGI was comprised of four main measures: responsibility for caring for children; occupational segregation; hours of work relative to partner/spouse; and education relative to partner/spouse. Differences in male and female participation rates in education in Canada and other developed countries have changed considerably since the early 1970’s,with women outnumbering men in university and post-secondary education completions [20, 21]. However, education was used in the construction of the LFGI as it is a measure of educational attainment relative to one’s partner/spouse, not a measure of absolute educational attainment. Each measure in the index is described in detail below.
Responsibility for caring for children
In each cycle of the LFS respondents are asked if they were away from work (either completely or partially) in the last week, and the reason for this absence, with one option being personal or family responsibilities. Respondents working less than 30 h per week are also asked the main reason they are not working more hours per week, with one option being caring for children and another being other personal or family responsibilities. Using responses to these questions, the following three category variable was created: 0 = no reduction in labour market participation due to personal or family responsibilities; 1 = part or full week absence due to personal or family responsibilities; 2 = working part-time due to personal or family responsibilities.
Self-reported occupation is coded into 47 major groups based on the National Occupational Classification system . For the LFGI responses to the 1997 LFS were used to classify each of these 47 occupations into one of four groups: 0 = occupations where less than 26 % of workers were women; 1 = occupations where 26 to 50 % of workers were women; 2 = occupations where 51 to 74 % of workers were women; and 3 = occupations where 75 % or more of workers were women. Occupations with the lowest participation of women are conceived as the most masculine occupations, while occupations with the highest participation of women are conceived as the most feminine occupations.
Hours of work relative to partner/spouse
Respondents are asked the usual number of hours they usually work each week. For respondents who are living with a spouse they are also asked the number of hours their spouse usually works per week. Using both these sources of hours worked each respondent was grouped into one of the following four categories: 0 = respondent working, but spouse not in the labour force; 1 = respondent working more hours than their spouse; 2 = respondent working the same number of hours as their spouse; and 3 = respondent working less hours than their spouse. If respondents did not have a spouse they were grouped with respondents working more hours than their spouse.
Education level relative to partner/spouse
Respondent’s and spouse’s highest level of education are reported in the following six categories: 0 to 8 years of education; some secondary education; graduated from high school; some post-secondary education; post-secondary certificate or diploma; and university degree. Using this information respondents were grouped into one of the following three categories: 0 = respondents with a higher level of education than their spouse; 1 = respondents with the same level of education as their spouse; 2 = respondents with a lower level of education than their spouse. Similar to work hours, respondents without a spouse were grouped with respondents with a higher level of education than their spouse.
To create the LFGI the values for the above four measures (caring for children, occupational segregation, hours or work and education level) were summed for each respondent providing a score ranging from 0 to 10, with higher scores indicating more traditionally feminine gender labour market roles of respondents and lower scores indicating more traditionally masculine gender labour market roles.
Correlations between the four measures of the LFGI were examined, and between each component and the final LFGI score. LFGI scores were then compared for men and women, and for the 1997 and 2014 LFS. Linear regression analyses then examined if the relationship between sex (male versus female) and LFGI scores changed between 1997 and 2014, after adjustment for differences in age, province, and month of survey participation between the 1997 and 2014 surveys. To examine if gender scores had changed for men and women between 1997 and 2014 a multiplicative interaction term between sex (male/female) and survey year was included in the model. The regression analysis was based on a 10 % random sample to avoid the possibility of a Type I error given the size of the LFS samples. All analyses were weighted to account for the initial probability of selection for each household, non-response and coverage errors, as specified by Statistics Canada . Analyses were conducted using SAS Version 9.3 .
Distribution of gender index components for Canadian men and women in 1997 and 2014
1997 LFS (N = 696,350)
2014 LFS (N = 729,132)
p-value for diff
p-value for diff
Responsibility for caring for children
No absence from work due to family or household responsibilities
Part or full-week absence due to family or household responsibilities
Works part-time due to family or household responsibilities
Occupation (based on 1997 LFS only)
Less than 26 % women
26 to 50 % women
51 % to 74 % women
75 % women
Hours of work
Respondent works spouse does not
Respondent works more than spouse/respondent does not have a spouse
Respondent works same amount as spouse
Respondent works less than spouse
Respondent higher level of education than spouse/respondent does not have a spouse
Respondent same education as spouse
Respondent lower education than spouse
Polychoric correlations between gender index and its components
1. Gender Index
2. Responsibility for caring for children
3. Occupation (based on 1997 LFS only)
4. Hours of work
Adjusted ordinary least squared (OLS) estimates for sex, survey year and their interaction in gender index score
Survey year/sex multiplicative interaction term
Disentangling the impacts of sex and gender in understanding male and female differences is increasingly recognised as an important aspect for advancing research and addressing knowledge gaps in the field of work-health . However, achieving this goal in secondary data analyses where direct measures of gender, such as the BRSI or PAQ, have not been collected is challenging. The objective of this paper was to demonstrate how a gender index – based primarily on gender roles – could be developed using routinely collected information from the Canadian Labour Force Survey. A second objective was to examine how gender scores were distributed among men and women (i.e. sex) and if there had been changes in gender roles among working Canadian men and women over the 17 year period between 1997 and 2014. Differences were observed between men and women in each component of the LFGI, with women generally having higher LFGI scores (indicating greater feminine gender roles) compared to men. While we found that women had higher LFGI scores in both 1997 and 2014 small increases in LFGI scores were observed for men between 1997 and 2014, and small decreases in LFGI scores for women over the same time period.
These study results should be interpreted taking the following strengths and limitations into account. The household-based sampling strategy employed by Statistics Canada in conducting the LFS resulted in a truly representative sample of the Canadian labour market, and findings can be generalised to Canadian labour market participants over the study time period. However, the large sample also increases the possibility of a Type I error and the inference of a meaningful difference that has no practical or meaningful importance. This may be the case for the observed differences in the LFGI score over time and the interpretation that men are taking on greater feminine gender roles while women are taking on greater masculine gender roles. In each of these cases the differences over time periods were less than 0.5 on an index that ranges from 0 to 10. To put this into context, if LFGI scores continue to increase among men and decrease among women at the same rate as observed over the 19-year study period (1997 to 2014), it will take until 2097 for men and women to have similar LFGI scores, indicating gender-equity in relation to labour market roles.
Three of the four measures that comprised the LFGI distinguished between social and occupational roles of men and women in the expected direction. However, a similar number of men and women had lower education than their partner/spouse in 1997, while men were more likely than women to have lower education than their partner/spouse in 2014. To some extent this result reflects the changing nature of characteristics previously thought of as masculine or feminine . For example, the BRSI has “ambitious” and “analytical” as masculine traits  while the PAQ has “likes math and science” and “intellectual” as masculine traits . Given changes in educational participation between men and women, along with the information presented in this paper, future work that creates indexes/measures that reflect gender roles and institutionalised gender may choose to exclude education (as both an absolute measure and in relation to the respondent’s partner/spouse) as a component of such measures.
Finally, the construction of the LFGI represents the sum of scores for the relative components. While this approach has the advantage of simplicity, making the approach easy to replicate, it does make assumptions about the relative contribution of each of the components of the index in relation to overall labour market gender roles, which may not be valid. Alternative approaches to constructing the index (e.g. factor analyses or cluster analyses) may be warranted, and researchers should weigh the advantages and disadvantages to each analytic approach if they choose to replicate the work in this paper.
How could the LFGI be used to better understand processes that create male/female differences in health?
Figure 2a presents a simple DAG where there is a difference in a health outcome for men compared to women (note this difference could be in either direction – i.e. more prevalent among men compared to women, or more prevalent among women compared to men). If a difference in the health outcome is present among men and women, path “A” in this DAG will be not equal to zero and will be statistically significant. For theoretical purposes, the relationship between male/female and the health outcome (path A) is assumed to be adjusted for all confounders, and that male/female and the health outcome (along with the additional variables included in Fig. 2b and c below) have been measured without error.
In order to understand why the risk of the health outcome is greater for men (or women), the DAG is extended to include intermediate or mediating variables, such as the LFGI, to understand the impact of “sex” and “gender” in male/female differences in the health outcome. This extended model is presented in Fig. 2b. Again, for theoretical purposes, we make no confounding for all paths and no measurement error assumptions for all variables. We now have a direct path (path A′) and an indirect path (paths B and C) linking male/female to the health outcome. The magnitude of the indirect path will be determined by the strength of the relationship between male/female at the LFGI (path B) and the relationship between the LFGI and the health outcome (path C). The estimate for path B is equivalent to the regression estimate for female (relative to male) presented in Table 3 previously. It is important to note that the estimate from Table 3 indicates that the LFGI score was strongly influenced, but not completely explained by whether the respondent was male or female. The estimate for the direct effect (path A′) can be interpreted as the difference in the health outcome for men compared to women that would remain if men and women had similar roles in relation to labour market status (i.e. if men and women had similar scores on the LFGI). The difference between path A and path A′ (which in an ordinary least-squared model will be equivalent for the product term of path B and C) , can be interpreted as the amount of the originally observed difference in the health outcome for men and women that can be explained by differences in labour market roles (as assessed by the LFGI) between men and women .
It is important to note that if path A′ is still associated with the health outcome then this indicates that male/female differences in the health outcome are not completely explained by differences in labour market roles only. The remaining differences between men and women (path A′) will likely be a combination of other biological (sex) differences between men and women that are relevant to the outcome, and other social (gender) differences between men and women that are both relevant to the outcome, and not captured in the LFGI. This has been explicitly described in Fig. 2c. If data was available to construct – either individually or as part of an index – all other sex and gender related factors that are relevant to the outcome of interest, then path A″ in Fig. 2c would approach zero, and one could examine the relative contribution of biological factors (paths F and G) and gender factors related to labour market roles (paths B and C) and non-labour market roles (paths D and E). The caveat for Fig. 2c is that each of the three pathways can be measured and estimated as distinct from each other. This highlights the need for the integration of “sex” and “gender” into the study design and data collection phase as part of a comprehensive research process , so that a more thorough examination of sex and gender into male/female differences in health status can be routinely undertaken using the approach outlined above.
In this paper we developed a measure of feminine and masculine gender roles, using self-reported survey data, where no direct measure of gender (masculinity/femininity) was available. This measure had face validity in terms of being related to, but distinct from sex (male/female), and was also responsive to change over time. Future research should examine the relative importance of including additional measures to an index such as the LFGI. For example, in the study by Pelletier and colleagues  primary earner status was the measure most strongly related to masculine BRSI scores, while the number of hours spent doing housework and responsibility for doing housework were the measures most strongly associated with high feminine BRSI scores. The LFGI in this paper included some indication of primary earner status and some indication of household responsibilities for respondents. However, a more detailed or direct measure of primary earner status and household responsibilities (in particular housework) may have allowed further refinement or distinction between masculine and feminine roles in the LFGI. In addition, it would be interesting to examine how an index with a reduced number of measures would perform in differentiating gender roles for men and women. As mentioned in the introduction to this paper, the early work by Lippa and colleagues focused only on differences in occupational preferences between men and women [11, 13]. Interestingly, occupational segregation, as well as hours worked, was the measure most strongly correlated with the LFGI in our study. We also suggest that research examining work-related health outcomes should (and in many cases can) integrate measures of sex and gender, using an approach similar to the one outlined in this paper.
BSRI, Bem-Sex-Role-Inventory; LFGI, Labour Force Gender Index; LFS, Canadian Labour Force Survey; PAQ, personal attributes questionnaire.
Peter Smith and Mieke Koehoorn are both supported by Research Chairs in Gender, Work & Health from the Canadian Institutes of Health Research.
PS and MK were both involved in the conceptual development of this paper. PS performed all data analyses and wrote the first draft of the manuscript. MK provided feedback on the manuscript. Both authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ristvedt SL. The evolution of gender. JAMA Psychiatry. 2014;71(1):13–4.View ArticlePubMedGoogle Scholar
- Phillips SP. Including gender in public health research. Public Health Rep. 2011;126(Supp 3):16–21.PubMedPubMed CentralGoogle Scholar
- Schiebinger L. Scientific research must take gender into account. Nature. 2014;507:9.View ArticlePubMedGoogle Scholar
- Johnson JL, Greaves L, Repta R. Better science with sex and gender: Facilitating the use of a sex and gender-based analysis in health research. Int J Equity Health. 2009;8:14.View ArticlePubMedPubMed CentralGoogle Scholar
- Doyal L. Sex and gender: the challenges for epidemiologists. Int J Health Serv. 2003;33(3):569–79.View ArticlePubMedGoogle Scholar
- Messing K, Punnett L, Bond M, Alexanderson K, Pyle J, Zahm S, et al. Be the fairest of them all: Challenges and recommendations for the treatment of gender in occupational health research. Am J Ind Med. 2003;43:648–29.View ArticleGoogle Scholar
- Nowatzki N, Grant KR. Sex is not enough: the need for gender-based analysis in health research. Health Care Women Int. 2011;32(4):263–77.View ArticlePubMedGoogle Scholar
- Springer KW, Mager Stellman J, Jordan-Young RM. Beyond a catalogue of differences: A theoretical frame and good practice guidelines for researching sex/gender in human health. Soc Sci Med. 2012;74(11):1817–24.View ArticlePubMedGoogle Scholar
- Bem SL. The measurement of psychological androgyny. J Consult Clin Psychol. 1974;42(2):155–62.View ArticlePubMedGoogle Scholar
- Eisler RM, Skidmore JR. Masculine gender role stress: Scale development and component factors in the appraisal of stressful situations. Behav Modif. 1987;11(2):123–36.View ArticlePubMedGoogle Scholar
- Lippa R, Connelly S. Gender Diagnosticity: a new Bayesian approach to gender-related individual differences. J Personal Soc Psychol. 1990;59(5):1051–65.View ArticleGoogle Scholar
- Spence JT, Helmreich R, Stapp J. The Personal Attributes Questionnaire: a measure of sex role stereotypes and musculinity-feminity. Cataloge Selected Doc Psychol. 1974;4(43):Ms No 617.Google Scholar
- Lippa RA, Martin LR, Friedman HS. Gender-related individual differences and mortality in the Terman Longitudinal Study: Is masculinity hazardous to your health? Pers Soc Psychol Bull. 2000;26(12):1560–70.View ArticleGoogle Scholar
- Pelletier R, Ditto B, Pilote L. A composite measure of gender and its association with risk factors in patients with premature acute coronary syndrome. Psychosom Med. 2015;77(5):517–26.View ArticlePubMedGoogle Scholar
- Borkhoff CM, Hawker GA, Kreder HJ, Glazier RH, Mahomed NN, Wright JG. The effect of patients’ sex on physicians’ recommendations for total knee arthroplasty. CMAJ. 2008;178(6):681–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Leung Yinko SS, Pelletier R, Behlouli H, Norris CM, Humphries KH, Pilote L, et al. Health-related quality of life in premature acute coronary syndrome: does patient sex or gender really matter? J Am Heart Assoc. 2014;3(4):e000901.View ArticlePubMedPubMed CentralGoogle Scholar
- Mosca L, Barrett-Conner E, Wegner NK. Sex/gender differences in cardiovascular disease prevention: what a difference a decade makes. Circulation. 2011;124(19):2145–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Statistics Canada. A guide to the Labour Force Survey. Ottawa: Statistics Canada; 2011. Report No.: 71-543-G.Google Scholar
- Statistics Canada. Labour Force Survey Public Use Microdata Files. 2013. Available through Statistics Canada’s Data Liberation Initiative.Google Scholar
- Frenette M, Zeman K. Why Are Most University Students Women? Evidence Based on Academic Performance, Study Habits and Parental Influences. Ottawa: Statistics Canada; 2007. Report No.: Catalogue no. 11F0019MIE, no.303.Google Scholar
- Statistics Canada. College graduates by program level, Classification of Instructional Programs, Primary Grouping (CIP_PG) and sex, annual (number). Postsecondary Student Information System (PSIS), CANSIM Table 477-0016. 2011.Google Scholar
- Human Resources & Skills Development Canada. National Occupational Classification Career Handbook. 2nd ed. Ottawa: Government of Canada; 2011.Google Scholar
- The SAS Institute. The SAS System for Windows, Release 9.3. Cary: SAS; 2011.Google Scholar
- Hoffman RM, Borders L. Twenty-five years after the Bem Sex-Role Inventory: a reassessment and new issues regarding classification variability. Meas Eval Couns Dev. 2001;34(1):39–55.Google Scholar
- Glymour MM. Using Causal Diagrams To Understand Common Problems in Social Epidemiology. In: Oakes JM, Kaufman JS, editors. Methods in Social Epidemiology. San Francisco: Jossey-Bass; 2006. p. 387-422.Google Scholar
- Schisterman EF, Cole SR, Plat RW. Overadjustment bias and unnecessary adjustment in epidemiolgical studies. Epidemiology. 2009;20(4):488–95.View ArticlePubMedPubMed CentralGoogle Scholar
- VanderWeele TJ, Robinson WR. On the causal interpretation of race in regressions adjusting for confounding and mediating variables. Epidemiology. 2014;25(4):474–84.Google Scholar
- MacKinnon DP, Fairchild AJ, Fritz MS. Mediation analysis. Annu Rev Psychol. 2007;58:593–614.View ArticlePubMedPubMed CentralGoogle Scholar