Abstract
There are many open questions concerning the development of calling, and longitudinal empirical evidence is limited. We know that a calling is associated with many beneficial outcomes, but we do not know how it changes through time and what predicts these changes. Previous studies have shown that calling is relatively stable at the sample level. We show that, at the individual level, calling shows huge variations through time. We identified nine developmental trajectories that are typical across facets of calling, and we found evidence that the development of a calling is fostered by the extent to which individuals have lived it out. We also observed that the more a calling has grown over a 2-year period, the more it is lived out during the third year. These results provide support for a developmental model of calling according to which having a calling and living it out reciprocally influence each other. The practical and theoretical implications of these results are discussed.
In the past, calling has been viewed as a stable construct that can be found or discovered (Bunderson & Thompson, 2009; Duffy & Sedlacek, 2007; Wrzesniewski, McCauley, Rozin, & Schwartz, 1997). Hence, researchers focused on the search for a calling (Duffy & Sedlacek, 2007) and investigated predictors and outcomes of having and living out a calling (Dobrow, 2013; Duffy, Allan, Autin, & Bott, 2013). One aspect for which there is almost no knowledge regards how a calling changes through time and how it develops at the individual level (Duffy, Douglass, Autin, & Allan, 2014). In this article, we investigate whether the results of previous longitudinal studies that investigated sample-average changes in calling are valid for single individuals, using a longitudinal cluster analysis method. We also explore predictors and outcomes of changes in calling, focusing on the relation between living out a calling and the latent growth curve of calling.
The Construct of Calling
Career calling is a multidimensional construct that has been defined in many ways (Dalla Rosa, Galliani & Vianello, 2014). In this article, we adopt a broad conceptualization of calling that includes the key dimensions of most previous conceptualizations (Bunderson & Thompson, 2009; Dik & Duffy, 2009; Dobrow & Tosti-Kharas, 2011; Hagmaier & Abele, 2012; Praskova, Creed, & Hood, 2015; Wrzesniewski et al., 1997). Indeed, it has been found that people who feel called to a career or domain are passionate about the calling domain (passion component: Bunderson & Thompson, 2009; Dobrow & Tosti-Kharas, 2011; Hagmaier & Abele, 2012; Wrzesniewski et al., 1997), have a purpose in life (purposefulness component: Dik & Duffy, 2009; Dobrow & Tosti-Kharas, 2011; Praskova et al., 2015), sacrifice other areas of their life in favor of their calling (sacrifice component: Bunderson & Thompson, 2009; Dobrow & Tosti-Kharas, 2011), feel that the calling’s domain is always on their mind (pervasiveness component: Dobrow & Tosti-Kharas, 2011; Hagmaier & Abele, 2012; Wrzesniewski et al., 1997), report that their calling comes from something greater than or beyond them (transcendent summons component: Dik & Duffy, 2009; Hagmaier & Abele, 2012; Hunter, Dik, & Banning, 2010; Zhang, Dik, Wei, & Zhang, 2015; Zhang, Herrmann, Hirschi, Wei, & Zhang, 2015), and that part of their calling is being useful to others (prosocial orientation component: Bunderson & Thompson, 2009; Dik & Duffy, 2009; Hagmaier & Abele, 2012; Praskova et al., 2015). In terms of outcomes, it has been repeatedly shown that having a calling fosters commitment and engagement with the calling’s domain and that individuals with a calling are more satisfied with their job and their life (see, e.g., Dalla Rosa, Galliani, & Vianello, 2017; Dobrow & Tosti-Kharas, 2011; Duffy et al., 2013; Hagmaier & Abele, 2012, for a review).
Given the importance that a calling has in life, it is normal to wonder how and why individuals develop a calling for a job, a life role, or a domain of activities. Drawing on the distinction between objective and subjective careers, Hall and Chandler (2005) theorized that the development of a calling is an ongoing, cyclical process of career self-exploration and discernment, in which self-confidence, goal setting, and efforts derived from the calling lead to both objective and subjective success that in turn reinforces the calling. Some of these hypotheses were later supported. Duffy, Douglass, Autin, and Allan (2014) found a reciprocal temporal relation between having a calling on the one hand and life meaning and personal growth initiative on the other hand. Duffy, Allan, Autin, and Douglass (2014) found a longitudinal reciprocal relation between living out a calling, career commitment, and finding meaning in one’s work. The distinction between the perception of having a calling and the actual experience of living out that calling in one’s job or life is critical if we assume that the development of calling is cyclical. “Perceiving a calling refers to the degree to which an individual believes she or he is called to a particular career. […] Living a calling refers to the degree to which an individual is currently engaging in activities or work that meet this calling” (Duffy & Autin, 2013, p. 220). The work as a calling theory (WCT; Duffy, Dik, Douglass, England, & Velez, 2018) recently formalized an assumption that has been alive for a long time in the literature: Having a calling motivates individuals to engage in activities related to the calling’s domain in order to live it out. Indeed, calling motivation is defined as the individual’s level of motivation to pursue her or his calling (Duffy, Autin, Allan, & Douglass, 2015). Hence, living out a calling may be seen as a form of subjective success that results from the interaction between individual efforts in pursuing one’s calling and favorable environmental opportunities.
Living out a calling may also be seen as a form of self-actualization. Indeed, the connection between living out a calling and self-actualization is very strong and dates back to the latest descriptions of Maslow’s theory of needs (1971, p. 40): Self-actualizing people are, without one single exception, involved in a cause outside their own skin, in something outside of themselves. They are devoted, working at something, something which is very precious to them—some calling or vocation in the old sense. They are working at something which fate has called them to somehow and which they work at and which they love, so that the work-joy dichotomy in them disappears […]. The devotion and dedication is so marked that one can fairly use the old world vocation, calling, or mission to describe their passionate, self-less, and profound feelings for their work. […] Such vocation-loving individuals tend to identify with their work and to make it into a defining characteristic of the self. […] The tasks are loved because they embody these (intrinsic) values. […] The self has enlarged to include aspects of the world and […] the distinction between self- and not-self (outside, other) has been transcended.
We argue that calling follows the same developmental process that Maslow theorized for self-actualization and self-transcendence: The more the individuals feel that they are living out a calling, the more they build on and understand their calling, which in turn strengthens their motivation and their efforts to live it out, in a cyclical process that ultimately leads to profoundly meaningful and self-transcendent activities. 1 This hypothesis is supported by the evidence that the relation between a calling and its positive outcomes is stronger when individuals are living out their calling: Living out a calling both mediates and moderates the impact of perceiving a calling’s meaning in work and life, work and life satisfaction, career commitment, and maturity (Dalla Rosa et al., 2017; Duffy et al., 2013; Duffy, Bott, Allan, Torrey, & Dik, 2012; Duffy & Dik, 2013). Comparing academics with an unanswered occupational calling to academics who were actually living out their calling and academics with no calling at all, Gazica and Spector (2015) found that the former group reported significantly poorer physical and psychological health. This result suggests that if an individual has a calling, living it out is critical for her or his well-being. Notably, Duffy and Autin (2013) suggested that perceiving a calling fosters living out a calling through the mediation of work volition: People with a calling are more able to live out their calling because of an increased sense of mastery and control in their careers. We expect that living out a calling will strengthen the perception of having a calling, which in turn will lead to increased efforts to live it out.
Evidence on the Development of a Calling
Empirical evidence on the development of a calling is limited to six longitudinal studies. Five of them followed participants for less than 2 years, and only two of them collected more than two waves of data. These studies focused on the temporal precedence between a calling and related variables, such as meaning in life, career planning, vocational identity, well-being, behavioral involvement, and social comfort in the calling’s domain. Hence, they did not specifically focus on intraindividual changes in calling. Rather, they studied sample-level means across time points. This is very common in the psychological literature, where most models are intended to describe processes that operate within each individual but are tested by looking at results averaged across many individuals. It has been shown that this approach may give a misleading picture of what is true for each one (Miller & Schwartz, 2018), especially with respect to the interpretation of on-average null effects. Previous longitudinal research observed that, on average, callings decrease very slightly through time. Yet Duffy, Manuel, Borges, and Bott (2011) observed a stronger decrease, and Duffy et al. (2014) observed a remarkable increase. Table 1 provides average annualized rates of change for all the longitudinal studies of calling currently published.
Average Standardized Rates of Change in Calling.
Note. When not reported in the original article, the annual rate of change (d) has been computed by taking the difference between mean levels of calling across observations and dividing it by the pooled standard deviation of observations across time points. This value has then been annualized by dividing the result by the time lag (in years) of each specific study. This effect size is a standardized mean difference, and can be interpreted using Cohen’s (1988) guidelines. Values smaller than .2 indicate small effects, values between .5 and .8 indicate moderate effects, and values greater than .8 indicate large effects. Participants’ age statistics were not provided in Duffy et al. (2011). We computed mean and SD of participants’ age using public data from the University of Cincinnati, College of Medicine (2010).
Although average rates of change in calling are small, variation across studies is not trivial and leads to a somewhat uninformative conclusion that, at the sample level, callings may decrease, increase, or even stay the same. Dobrow (2013) and Hirschi & Hermann (2013) suggested that this may happen for individuals as well, but empirical evidence on intraindividual changes in calling is very limited. Indeed, there is no study that provides a description of the individual trajectories of change or that investigates their shape at the individual level.
We do not know whether individual changes in calling across time are frequent or rare nor whether changes are typically small, moderate, or large. Current empirical evidence on changes in calling may be compatible with many different developmental patterns. Results of previous studies could have been generated by a wide range of individual growth curves. For instance, the results in Dobrow (2013) are compatible with both a scenario in which a large majority of people show a slight decline and a small minority show a steep decline, and one in which half of the participants decrease steeply and the other half increase somewhat less steeply. So far, only two studies have been published that collected at least three waves of data, which are necessary to model some simple shapes of developmental trajectories. One of them (Dobrow, 2013) did not investigate this question. The other one (Hirschi & Herrmann, 2013) found that, at the sample level, a nonlinear rate of change better fits the data. However, the authors did not investigate further the shapes of the individual trajectories. Specific knowledge on how and why calling changes would inform the theory on the dynamic nature of calling. For instance, a sample-level quadratic shape may suggest that calling is susceptible to change, but then elastically tends to return to the initial level. High intraindividual variability over time may suggest that calling is a state-like construct. On the other hand, high intraindividual stability would suggest that it is a dispositional or trait-like construct. Also, the lack of a general developmental trend across subjects would suggest that calling is not a dispositional construct but the result of continuous, idiosyncratic interactions with the environment. We decided to investigate the development of calling in a sample of college students. This choice has been informed by the fact that this phase of life is critical in the development of constructs that are similar to calling, such as vocational interests (Low, Yoon, Roberts, & Rounds, 2005).
Method
Participants
The data used in this article are part of a multiwave data collection that began in 2014 and is currently ongoing. A sample of 5,886 Italian college students of all grades was involved in the first wave, 1,700 students took part in the second wave of data collection, and 881 took part in the third wave of data collection, resulting in a total sample size of 6,370 students. The sample was mainly composed of women (63.8%, at T1). Participants’ mean age at T1 was 23.37 years (SD = 5.39, min age = 18 years, max age = 69 years). All participants were bachelor’s (68% of the total, min age = 18 years, max age = 64 years, M age = 22.55, SD = 5.38) or master’s (30%, min age = 21 years, max age = 69 years, M age = 25.12, SD = 4.90) students enrolled in 24 different study domains, 44% had work experience at T1. In this article, we analyze the responses of the total sample for analyses involving structural equation models (SEMs), estimating missing values by means of full information maximum likelihood algorithms. We adopted a listwise deletion procedure (Nunnally & Bernstein, 1994) for analyses involving longitudinal clustering procedures, resulting in the exclusion of participants who were totally absent from one wave of data collection. The subsample of students that participated in all three waves (N = 434) is composed of bachelor’s students (75%, min age = 18 years, max age = 57 years, M age = 21.12, SD = 3.98) and master’s students (25%, min age = 21 years, max age = 59 years, M age = 24.51, SD = 5.51). Sixty-nine percent were females, and 38% were employed.
A second article (Vianello, Dalla Rosa, Anselmi, & Galliani, 2018) has been published that reports the analyses that were run to validate the multidimensional measure of calling that we employ in the current article. Details on this work, including the open access article and the public data set, can be found here: https://osf.io/viwsu/. The validation study involved two samples: A sample of working adults (no overlap with the current article) and a sample of college students that partially overlaps with the data set used in the current study. All analyses reported in the validation study (exploratory factor analysis (EFA), confirmatory factor analysis (CFA) multigroup invariance, and external validity) were based on the responses collected at Time 2 (T2; N = 1,700), except for the longitudinal invariance test, for which data from all three waves were used (N = 434). The goals of these two manuscripts and the variables involved, with the exception of the multidimensional calling scale, are completely different. A third, unpublished, manuscript (Dalla Rosa, Vianello, & Anselmi, 2017) has been prepared for reporting the results of a longitudinal path analysis that identifies three predictors of the levels of having a calling 1 year apart. Details regarding this work, including the data set, the manuscript, and the analysis code, can be found here: https://osf.io/5jgrp/. The smaller data set used in the present article (N = 434) is the same that we used in Dalla Rosa et al. (2017), which reports that social support, engaged learning, and clarity of professional identity predict the average levels of calling 1 year later. The current article is completely different in scope, analyses, and variables involved. Anonymized data and analysis scripts for the current article are publicly available on the Open Science Framework (https://osf.io/he3zc/).
Procedure
Data were collected nonanonymously using an online survey. The second and third waves occurred, respectively, 12 and 24 months after the first data collection. The choice of the time lag was informed by previous research showing null or very small changes in calling or closely related constructs such as work passion (see, e.g., Tóth-Király, Bőthe, Jánvári, Rigó, & Orosz, 2019). Participants were invited to take part in the survey through an e-mail that briefly described the study and provided a link to the survey. The landing page of the web survey presented a more detailed description of the study and the predicted amount of time required to answer the questions. Then, participants were informed of their right to exit the study at any point and have their data deleted, and they were provided with a privacy statement and a contact e-mail. Finally, they were informed that going further in the survey would have been registered as a sign of their willingness to participate nonanonymously and to share their responses with the researchers.
Measures
Calling
There is currently no agreement on the core components of calling (Thompson & Bunderson, 2019). Hence, we employed the most comprehensive measure of calling that is currently available: The Unified Multidimensional Calling Scale (UMCS; Vianello et al., 2018). This 22-item Likert-type scale measures both classic and modern components of calling and draws heavily from previous work by Dik, Eldridge, Steger, and Duffy (2012) and by Dobrow and Tosti-Kharas (2011). The UMCS measures seven facets of calling: passion, purposefulness, sacrifice, pervasiveness, prosocial orientation, transcendent summons, and identity. The scale underwent an extensive validation study (Vianello et al., 2018), which provided evidence of satisfactory to excellent reliability and convergent, discriminant, and predictive validity.
The UMCS is the only multidimensional measure of calling for which strict measurement invariance over different observations has been demonstrated. Strict measurement invariance (equal factor loadings, equal indicator intercepts, and equal item uniquenesses; Meredith, 1993) is a key measurement property that should be tested before running any longitudinal studies, and a necessary requirement when scores on a scale are compared across repeated observations. Strict invariance across time is rarely achieved, and the interpretation of change in noninvariant measures can lead to heavily misleading conclusions because observed changes can be attributed to measurement artifacts rather than to true changes in the construct under investigation. The facets that are included in the UMCS have been extensively studied in the literature on calling and are included in the leading theoretical conceptualizations of the construct (Bunderson & Thompson, 2009; Dik & Duffy, 2009; Dik, Eldridge, Steger, & Duffy, 2012; Dobrow & Tosti-Kharas, 2011; Hagmaier & Abele, 2012; Praskova et al., 2015; Wrzesniewski et al., 1997). Adopting a conceptualization and measure of calling that employs all of them therefore maximizes the generalizability of results throughout the many different conceptualizations of calling that currently exist in the literature.
The passion component of the UMCS reflects the pleasure and satisfaction that individuals feel in carrying out activities related to the calling’s domain (example item: “I am passionate about what I am studying”). Purposefulness reflects how much the engagement in the calling’s domain provides the opportunity to find purpose and meaning in life (example item: “I see my academic and professional career as a path to purpose in life”). Sacrifice measures individuals’ willingness to sacrifice other areas of their lives to pursue their calling (example item: “I can give up many things to keep studying this subject”). Pervasiveness refers to a constant presence of the calling’s domain in individuals’ thoughts and consciousness, even when people are not directly involved in domain-related activities (example item: “My studies are always on my mind”). Prosocial orientation refers to other-oriented purposes that guide the involvement in the calling’s domain (example item: “The most important part of my future work is helping others to meet their needs”). Transcendent summons measures the presence of a spiritual push toward engaging in activities in the calling’s domain (example item: “I believe that I have been called to pursue my current line of study”). Finally, identity refers to the importance of the calling in defining who individuals are and the feeling of predetermination and inevitability in pursuing their calling (example item: “What I study is part of my destiny”). The identity scale of the UMCS was not used in this study because this scale was missing in the first wave of data collection. For all the analyses in this article, standardized factor scores are used unless otherwise noted. Having a calling was assessed at each time point, the reliabilities at T1 for the total sample ranged from .71 to .85. The estimation of factor scores was conducted in MPlus using full information maximum likelihood to handle partially missing data.
Living out a calling
The degree to which students felt they were currently living out their calling was assessed at T2 and T3 with the item “Are you living out your calling in the program you are enrolled in?” evaluated on a 4-point scale ranging from not at all to very much.
This single-item measure was validated on an independent sample (N = 205, M age = 46.8, SD = 11.18; 32.4% females), in which we measured living a calling with the scale developed by Duffy, Allan, and Bott (2012). Details on this study, including the codebook, the description of the study, and the public data set, can be found here: https://osf.io/ed9yh/. In line with previous studies (r xx = .94 in Duffy & Autin, 2013; r xx = .95 in Duffy et al., 2012 and in Duffy et al., 2013), the multi-item measure exhibited very high internal consistency (r xx = .95). The Spearman Brown prophecy formula (Nunnally & Bernstein, 1994, pp. 263–264) indicated that single-item measures extracted from this scale would have been satisfactorily reliable (r xx = .76). Hence, the item with the highest factor loading (.91) and the highest item-total correlation (r = .87) were selected for inclusion in this study.
We also checked whether in this study correlations with living a calling are attenuated by unreliability. Previous research showed that living a calling correlates with the presence of a calling (Duffy et al., 2013; Duffy et al., 2012; Duffy & Autin, 2013; r ranged between .32 and .49) and academic satisfaction (Duffy, Allan, & Bott, 2012; r = .29). In the current study, we found living a calling to relate in the expected directions and with similar or stronger intensity to presence of a calling (r t2 = .54; r t3 = .42) and academic satisfaction (r t2 = .43; r t3 = .44). Lastly, the single-item measure of living a calling that we employed shows a test–retest reliability of r xx = .54, which is in line with a previous result obtained using the multi-item measure (r xx = .53; Duffy, Allan, Autin, & Douglass, 2014).
Academic satisfaction
The degree to which participants were satisfied with their study program was measured at T2 and T3 with the question: “How much are you satisfied with your study program?,” evaluated on a 4-point scale ranging from not at all to very much. It has been shown that a global rating measure of satisfaction is likely to provide a more tailored picture of participants’ view (Nagy, 2002; Scarpello & Campbell, 1983). We observed a test–retest reliability of .58 after 1 year and correlations in the expected direction (Duffy, Allan, & Bott, 2012; Duffy & Autin, 2013; Duffy et al., 2012; Duffy et al., 2013; Nauta, 2007; Strahan & Credé, 2015) with living a calling (r t2 = .43; r t3 = .44), dropout intention (r t2 = −.19; r t3 = −.21), and grade point average (r t2 = .13; r t3 = .20).
Dropout intentions
Participants’ intention to drop out of college was evaluated at each time point with 1 item asking them to choose the option that best described their willingness to abandon their studies on a scale from 1 (I am going to finish my study program) to 4 (I am going to leave my study program). Correlations between scores on this single-item measure after 1 year (r t1−t2 = .31, r t2−t3 = .53, and r t1−t3 = .33) were in line with previous findings (Eicher, Staerklé, & Clémence, 2014; Perez, Cromley, & Kaplan, 2014). Test–retest reliability was r xx = .49 (change-free estimate: Heise, 1969, Equation 9). This measure correlates in the expected direction (Aitken, 1982; Strahan & Credé, 2015) with academic satisfaction (r t2 = .27; r t3 = .24) and grade point average (r t1 = .16; r t2 = .18; r t3 = .17).
Employment status and work experience
We asked participants to report if they had worked or were currently working at the time of the data collection. For those who had work experience, we asked to report how long they had worked or have been working. The response scale was from 1 to 7 (1 = some days; 2 = a few weeks; 3 = 1 month; 4 = a few months; 5 = 1 year; 6 = a few years; 7 = many years).
Age, gender, employment status, work experience, and mean grade were assessed at each time point. All single-item measures had to be employed to limit the number of items in the questionnaire and the time needed to answer. Too many questions would have had detrimental effects on both the partial nonresponse rate and the mortality of participants across the three waves of data collection (Drolet & Morrison, 2001).
Results
Attrition Analysis
We analyzed the nonrandom sampling effect of subject attrition following the strategies proposed by Goodman and Blum (1996). We analyzed whether attrition affected the probability of participating in all data collections based on measures of calling; work experience; dropout intention; mean grade; employment status, age, and gender at T1 and T2; and academic satisfaction and living out a calling at T2. Then, we analyzed whether reducing the sample to only those who participated in the three data collection points had an effect on means and variances.
Results suggest that the only nontrivial differences between those who remained and those who dropped out of the study regarded age (d = .28) and work experience (d = .26). Specifically, students who provided data at all three observations were younger (M age = 21.97) and had shorter work experience (M = 5.09) than students who dropped out earlier than T3 (M age = 23.49; mean work experience = 5.41 years). This probably occurred because students who graduated before T3 were less likely to read our invitation to the study.
Calling: State or Trait? Relative Versus Absolute Changes at the Sample Level
Our first aim was to investigate the sample-average changes in calling and compare them with previous longitudinal studies. In the current study, the standardized annual rate of change from T1 to T3 ranged from d = −.09 (transcendence) to d = .14 (purposefulness). Hence, at the sample level, our results confirm previous studies: Calling showed tiny differences across time. A different picture emerges if we analyze absolute changes, irrespective of their sign. In the aggregate, the mean absolute change from T1 to T2 and from T2 to T3 was Δ z 1 = .57 and Δ z 2 = .53, respectively. Opposite individual trajectories must compensate for each other, concealing the actual amount of intraindividual change. Disaggregated analyses showed that the most stable facet was passion, for which mean absolute changes of Δ z 1 = .49 and Δ z 2 = .46 were observed. The least stable facet was pervasiveness, for which we observed mean absolute changes of Δ z 1 = .60 and Δ z 2 = .63. In repeated-measures designs, these are moderate-to-large effect sizes. Very large changes (Δ z > .8) were present in 25% of participants.
Taken at face value, these results might seem incompatible with the notion that calling is a dispositional trait, which is relatively stable by definition. As a comparison parameter, a recent meta-analysis showed that, between ages 22 and 30, the Big Five traits change less than .28 SD (−.01 SD for Openness; Roberts, Walton, & Viechtbauer, 2006), which means less than .04 SD per year. Yet, clusters of individuals were identified (Morizot & Blanc, 2003) in which changes in personality were up to 2 SD between ages 17 and 30. Hence, in order to investigate whether individual trajectories follow similar patterns of change across participants, we conducted a longitudinal cluster analysis.
Changes in Calling at the Individual Level: Clusters of Change
To investigate how and how much calling changes at the intraindividual level, we preliminarily considered that, for each person and for each time lag, calling can increase, decrease, or stay the same. Hence, we identified nine theoretical trajectories, represented in the left panel of Figure 1. Then, we clustered individuals based on their trajectories of change over time using the R package “TSclust” (Montero & Vilar, 2014). This package implements a clustering procedure proposed by Chouakria and Nagabhushan (2007), which is based on the temporal correlation coefficient CORT (Chouakria, 2003). Being a correlation coefficient, the value of CORT ranges between −1 and 1. A value of 1 means that, over any time interval [ti , ti +1], the two trajectories increase or decrease simultaneously with the same growth rate. A value of −1 means that the two trajectories show a similar growth rate but opposite directions (one increases whereas the other decreases). Finally, a value of 0 means that there is no monotonicity between the two trajectories and that their growth rates are stochastically linearly independent. By taking into account the dependence between consecutive time points, CORT has been found to outperform classical correlation coefficients when comparing the behavior of two time series (Chouakria & Nagabhushan, 2007).

Clusters of individuals with similar trajectories of change for each facet of calling. Gray lines represent observed individual trajectories. The thick black line represents the mean trajectory across subjects. Theoretical trajectories are provided in the first column.
The clustering procedure allows the researcher to modulate the contribution of the proximity of trajectories with respect to the behavior (as expressed by coefficient CORT) and their proximity with respect to values (e.g., the distance between scores at any time point). We maximized the contribution of proximity with respect to behavior (i.e., the shape of trajectories). The clustering procedure might have difficulty in producing a cluster of individuals with a stable trajectory (a stability cluster). The coefficient CORT is undefined (CORT = 0/0) if at least one of the two trajectories is invariant over time (i.e., the same measure is observed in all time points). An undesired result occurs when both the trajectories are invariant over time, that is, the value of CORT is undefined instead of being 1. In the present work, the clustering procedure divided trajectories with the same monotonicity (i.e., that increased or decreased simultaneously) but different growth rates into different clusters before producing a stability cluster. Therefore, such a stability cluster was not obtained before a large number of clusters was extracted. Given the difficulties of the clustering procedure in producing the stability cluster, we decided to manually assign to this cluster the individuals who showed differences from T1 to T2 and from T2 to T3 that were smaller in size than .1 standard deviations, which are usually denoted as trivial in size (Cohen, 1988). The clustering procedure was implemented to distinguish the remaining individuals. Solutions with nine clusters were estimated and accepted unless similar trajectories were divided into different clusters. In this case, we estimated solutions with one fewer cluster until the more parsimonious solution best approximated the theoretical clusters. This clustering procedure was applied on each facet of calling separately.
The right panels of Figure 1 represent the individual trajectories of change for each cluster and for each facet. The vertical axis of each panel shows the standardized scores of the calling facets under investigation. The black line represents the average trajectory across subjects, and the gray lines represent the trajectory of each individual in that cluster. The number of subjects in each cluster goes from 23 to 102 (mean size = 49, SD = 18). The median CORT coefficient between trajectories in the same cluster goes from .91 to .99, indicating that the solution categorized individuals with very similar trajectories. All of the theoretical clusters have been observed except for the “increase-stable” cluster in the pervasiveness facet. As can be seen in Figure 1, individual trajectories can be extremely variable. The same subject could go three standard deviations down from T1 to T2, and then three standard deviations up from T2 to T3 (see, e.g., the “decrease-increase” cluster for transcendence).
Are Individual Changes in Calling Consistent Across Components?
Consistently with our multidimensional conceptualization and operationalization of calling, we expect to observe similar developmental trajectories in different components of calling within the same individuals. Hence, we analyzed the degree of consistency between individual changes across facets of calling, computing for each individual the CORT coefficients between the trajectories of the facets of calling. Results show that individuals who show a trajectory in a facet tended to show a similar trajectory in other facets. Yet, the degree of consistency depends on facets. Table 2 shows the median CORT coefficients across individuals, disaggregated by facets of calling. We observed a moderate-to-strong correlation between the trajectories of different facets. Overall, prosocial orientation and transcendent summons changed more independently than the other facets. Passion, sacrifice, and pervasiveness tended to change very consistently with one another.
Median CORT Coefficients Between Trajectories of the Same Individuals by Facets of Calling.
Note. N = 434.
Predictors and Outcomes of Change in Calling
So far, we have observed that calling is extremely variable and that developmental trajectories can take any possible shape. In this section, we investigate whether age, gender, mean grade, work experience, and living a calling have an impact on individual change trajectories. To this end, we specified and estimated an unconditional latent growth curve model (LGCM; Bollen & Curran, 2006; Little, 2013). Being a multilevel SEM, this model estimates the trajectories of change in terms of slopes and intercepts. In LGCM, latent variables for both the initial level (intercept) and the rate of change across time (slope) are specified. Hence, change trajectories were modeled as a distribution across subjects whose mean and standard deviation are estimated. As per standard practice in latent growth curve modeling, the loadings of the intercept factor were all constrained to 1. This set of constraints gives meaning to the intercept factor, which can be interpreted as the distribution of the initial level of individual trajectories. The shape of the trajectories was defined by constraining at 0 and 1 the loadings from the slope factor to the observed measure of calling at T1 and T2, respectively. The unconstrained estimation of the loading at T3 gives the trajectory a shape and informs the amount of change that happens from T2 to T3 in terms of the proportion of change that occurred from T1 to T2. For instance, an estimated loading of 2 at T3 would indicate a positive linear growth, whereas an estimate of 4 would indicate a positive quadratic growth. Hence, the slope factor can be interpreted as a distribution of the shapes of the growth curves across individuals.
An LGCM was estimated for each facet of calling. The differences in parameter estimates across models were trivial. Hence, for conciseness and the purposes of presentation, we computed a composite score of calling, averaging the scores of the six facets, and used this variable in the model depicted in Figure 2. 2 In this model, the unstandardized T3 loading was 1.15, meaning that the change that occurred at T3 was 85% less than the amount of change that we would have expected, assuming linear growth. In other words, 78% of the total change occurred before T2. Although fluctuations are in general very large, they decreased during the second year. Hence, calling tends to stability, and this result is compatible with the hypothesis according to which we observe calling in a phase of life in which it is highly malleable. This developmental pattern has been previously observed for vocational interests (Low et al., 2005), so we reasoned that it would have been worth testing whether changes in calling are predicted by the age of respondents or by their work experience. It may be that older students and students in employment change less than younger, full-time students. Hence, we added participants’ age and work experience to the model. We also included gender and mean grade as predictors, and mean grade, dropout intentions, and academic satisfaction as outcomes. For all outcomes and all predictors, we specified relations with both the intercept and the slope of growth trajectories. The final model has five predictors (age, gender, work experience, mean grade at T1, and living out a calling at T2) and four outcomes of change (living out a calling, mean grade, academic satisfaction, and dropout intentions at T3). Results are reported in Figure 2, which only displays significant paths. The left side of Figure 2 displays the relations between the predictors and the growth trajectories (intercept and slope). We can see that all predictors influence the initial level of having a calling (intercept factor). Yet the size of the relations, which is provided by the estimate of the regression path between predictors and the latent factors, is trivial for age, gender, and mean grade. Female, experienced, and older students tend to exhibit a slightly higher initial level of calling. The relations between age, gender, mean grade, and work experience are not significant for the slope factor; hence, the rate of change in calling is not influenced by these variables. The strongest predictor of the intercept and the slope of the growth trajectory is living out a calling. The more the students felt that they were living out their calling, the higher they started and the more their calling increased across 2 years.

Predictors and outcomes of changes in calling. Only significant paths involving the intercept and slope of the growth curve are reported. All estimates are standardized except for the loading of Calling (T3) to the slope factor, which is more interpretable if not standardized. All possible relationships among predictors, outcomes, intercept, and slope, including paths, errors, residuals, and correlations have been estimated but are omitted for clarity from the figure. Italicized numbers are fixed values that define the meaning of the constructs, identify the model, and set the scale. χ2(71) = 61.52, p < .001, CFI = .99, RMSEA = .019, 95% CI [.014, .024].
The right side of Figure 2 displays outcomes of changes in the perception of having a calling. We can see that both the initial level (intercept) and the rate of change (slope) positively predicted the extent to which students reported that they were living out a calling at T3. Finally, we observed that both the intercept and the slope of the growth trajectory significantly predict participants’ dropout intentions, academic satisfaction, and mean grade. The higher participants’ calling was at T1 (intercept), and the faster their calling grew over the years (slope), the more participants were satisfied with their program, received higher grades, were less willing to leave the program, and lived out their calling at T3. The strongest effects were observed for academic satisfaction and living a calling. The intercept and slope factors account for 38% of the variance in academic satisfaction and living a calling at T3. The variance accounted for by the intercept and slope factors is much smaller for dropout intentions (R 2 = .13) and mean grade (R 2 = .08).
These results clearly suggest that the development of a calling is fostered by the extent to which individuals live it out. Individuals who are living out a calling are more likely to have higher and steeper growth trajectories. In parallel, the more a calling grows, the more individuals feel that they are living it out. Having and living a calling are reciprocally related. In addition, we observed that both the intercept and the slope of the change trajectories strongly and positively influenced important academic outcomes such as dropout intentions, mean grade, and academic satisfaction at T3.
Discussion
Results of previous studies on the development of calling are inconclusive. Different authors observed different temporal changes: At the sample level, calling increased (Duffy, Douglass, et al., 2014), decreased (Duffy, Manuel, Borges, & Bott, 2011), and remained substantially stable (Bott & Duffy, 2015; Dobrow, 2013; Hirschi & Herrmann, 2013; Praskova, Hood, & Creed, 2014). Our findings account for these discrepancies: Interindividual differences in change trajectories can explain the inconsistency of results obtained at the sample level of analysis. We observed that intraindividual changes in calling can be huge and that individual trajectories are very different from one another. These individual variations would not be surprising if we assumed that a calling changes as a function of the interaction between an individual and environmental conditions, events, and situations. Pending new evidence that contradicts or limits the generalizability of these results, we suggest that the calling construct at least moderately depends on the interaction with the environment rather than on stable individual factors such as genes, values, or personality traits. The notion that a calling is largely influenced by state-like components is also supported by the moderate intensity of autocorrelations (r = .5 to .7) that have been reported by Bott and Duffy (2015), Hirschi and Herrmann (2013), Dobrow and Heller (2015), Dobrow (2013), Duffy, Manuel, Borges, and Bott (2011), and Duffy, Douglass, et al. (2014). However, we should only tentatively accept the conclusion that calling is a state (as opposed to a trait) concept. Indeed, we did not include in our design a comparison parameter such as extroversion or conscientiousness against which the developmental patterns of calling could have been more meaningfully interpreted.
The implications of our results for the conceptualization of calling depend on the generalizability of the observations. If these results are validated and generalized to the whole life span, then we will be able to conclude that calling is not a stable disposition, but a state that heavily depends on interactions with the environment, such as experiences or feedback in the calling’s domain. However, we might have observed calling during a phase of life in which it is more flexible, perhaps because participants were in the process of building it. Indeed, results of the LGCM showed that the growth curve of calling tends to stability, as changes between T1 and T2 are more intense than those from T2 to T3. Although the effects of age or employment status on calling changes are tiny and practically trivial, a pattern of stability is more likely to be present in older students. The mean age of subjects in the stability cluster (M age = 22.96, SD = 5.47, N = 147) was slightly higher than that of the other clusters (M age = 21.62, SD = 4.46, N = 287). These results are compatible with the notion that we observed calling in a phase of life in which it is highly malleable. Calling could behave like vocational interests: It may change, even dramatically, during college and then stabilize in adulthood or during the first experiences of work socialization (Low et al., 2005). In light of these considerations, we can hypothesize that calling is a dispositional trait whose development occurs during the first experiences in the calling domain, such as college or any other vocational development program. We look forward to longitudinal research on working adults for whom stability of calling is expected. In the absence of such research, any conclusion regarding the stability of calling, or its state or trait nature, would be premature.
In this study, we used a model of calling (Vianello et al., 2018) that merges together many different facets (UMCS). We observed that these facets tend to covary, that is, share the same trajectory of change. This result can be interpreted as further evidence that the six dimensions are facets of the same construct. Yet the large range of correlations between .36 and .89 suggests that the dimensions of a calling can develop independently from one another and that some facets are more independent than others. Specifically, we observed that prosocial orientation and transcendent summons tend to change more independently than the other components of calling. This result may raise questions concerning the factor structure of the UMCS because it could be interpreted as supporting the notion that prosocial orientation and transcendent summons are not part of the same construct. This interpretation is debatable. On the one hand, we think that it is normal for the components of a multidimensional construct to develop independently from one another because they are likely to be influenced by different experiences. For instance, personality traits develop differently throughout life (Srivastava, John, Gosling, & Potter, 2003), yet they were consistently shown to form two higher order factors (Digman, 1997, Schermer & Goffin, 2018). On the other hand, Vianello et al. (2018) reported that prosocial orientation and transcendent summons showed small second order loadings in college students and adequate loadings in a sample of workers. Altogether, these results may support the hypothesis that prosocial orientation and transcendent summons become part of a calling during a later stage of its development.
One more important result of this study concerns the role of academic performance. We found that mean grade is an outcome of calling rather than an information used by individuals to understand whether they have a calling. The tiny or null relation that we observed between mean grade at T1 and the rate of change in calling suggests that specific skills in the calling’s domain are not a necessary condition for the development of a calling. This finding is partially consistent with Dobrow (2013), who found that an individual’s ability in the calling’s domain was neither a predictor of initial levels of the calling nor of temporal changes in the calling; this adds to the literature the demonstration that ability is a consequence, rather than a predictor, of levels and changes in calling. In our model, mean grade at T3 is influenced by both the intercept and the slope of the growth curve: Higher levels and increases in calling improve academic performance. We might speculate that this relation is mediated by commitment and dedication to the calling’s domain. We hope that future research will test this hypothesis and, more generally, investigate the process by which calling fosters success. We also observed that both the intercept and the slope of the growth curve predict dropout intentions and academic satisfaction. Not only do students with higher levels of calling show higher academic satisfaction and lower dropout intentions at T3, but also so do those that, regardless of their initial level, increased their calling. This result is of great practical importance since it represents the first empirical evidence that the process of developing a calling in itself is beneficial to the individual and predicts academic achievement.
Finally, and perhaps most importantly, we showed that living out a calling is the favored path toward the growth of a calling. We found support for Hall and Chandler’s (2005) cyclical model, and we identified a virtuous cycle between having and living out a calling, such that they reciprocally foster each other. These results support some of the predictions made by the WCT (Duffy et al., 2018), according to which a reciprocal relation exists between perceiving and living a calling. We found that higher levels of a calling and changes in a calling are both outcomes and predictors of higher levels of living out a calling. On the other hand, higher levels of living out a calling are both predictors and outcomes of levels of and changes in calling. We suggest that environmental opportunities to live out a calling are key, not only for experiencing meaningfulness at work, career commitment, and the general well-being associated with living out a calling but also for the development of a calling itself. In agreement with the WCT (Duffy et al., 2018), we believe that a calling develops through experiences in the calling’s domain, and the perception of living out a calling through these experiences is used by individuals to strengthen their beliefs of actually having a calling. On the other hand, the lack of opportunities to live out a calling will inevitably weaken individuals’ beliefs in actually having a calling. This result is also in agreement with the role attributed to the person–environment fit by Hagmaier and Abele (2012), and it is consistent with Praskova, Creed, and Hood’s (2015) definition of calling as “a developmental and dynamic construct” (p. 102).
Limitations of the Present Study
Our results clearly represent a starting point rather than an end point: We accounted for 10% of all changes in calling and look forward to research that can help us understand which other factors account for the remaining 90%. Second, our sample was composed of college students, so we could not test our hypothesis regarding the stabilization of calling during later adulthood. We hope future research will investigate how and how much calling changes in working adults and that it will thoroughly examine the virtuous cycle between having a calling and the feeling of living it out during professional activities. Furthermore, the generalizability of our results is still to be proven. Although it is likely that replications of our study on college students would observe the same pattern of results, we urge readers not to generalize these results to workers or ages other than young adulthood. Third, our measure of living out a calling was composed of only 1 item and was specifically referred to students’ academic experience. This is likely to exclude students who live a calling outside of their studies. As a consequence, we might have underestimated the number of students living out a calling. Also, living a calling during college might be a complex experience. Students may be aware that they are actually learning something that they feel called to, but they might lack the engagement in activities that help them to live out their calling. Hence, our measure of living a calling might not differentiate well how much participants were living a calling from their perception of having one, thus inflating the relation between them. 3 Although we believe that attending classes, reading course materials, discussing topics, work placements, and internships are activities that inform students about how much they are living their calling toward a domain, we cannot rule out the possibility that our estimates were bigger than they would have been if participants had been working in, rather than learning, the calling domain. Future studies might investigate this point and consider the development of a scale for measuring the construct of living a calling in college students. Finally, a limit of this research is imposed by the time lag that we employed. Our choice was informed by the implicit assumption that calling is a dispositional trait and by previous studies showing basically no short-term changes in calling. In light of the results discussed in this article, we would now choose a shorter time lag. Although we suggest interpreting the huge intraindividual variability as a sign of its development, it is necessary to exclude that changes in calling are exceedingly rapid and frequent. If high intraindividual variability is observed daily or weekly, then conceiving calling as a state-like construct would probably be wiser.
Conclusion and Future Prospects
This study answers some questions and opens many others. In general, our results emphasize how little we know about the development of a calling. The approach that previous research followed in conducting the analysis focused on sample-level averages. We wonder whether this approach might prevent us from analyzing the dynamic nature of calling. Hopefully, future research will consider adopting a person-specific approach to the study of calling development. Indeed, we observed that, across individuals, calling develops in very different ways. Also, we observed that these change trajectories are at least partially influenced by the extent to which people can live out their calling. Finally, we showed that positive outcomes are not only related to high levels of calling but also to changes. When controlling for the initial level of the growth trajectory, an increase in calling determines a longitudinal increase in important academic outcomes such as dropout, mean grade, and satisfaction with the program.
We still can’t provide a definitive answer on the state-or-trait nature of a calling nor on its origin, but we paved the way to future studies by showing that this answer may be found by investigating person–environment interactions.
Footnotes
Authors’ Note
Michelangelo Vianello coordinated the project, designed the study, developed materials, collected and analyzed the data, and wrote the report. Elisa Maria Galliani developed materials, analyzed the data, and wrote the report. Anna Dalla Rosa developed materials, collected and analyzed the data, and wrote the report. Pasquale Anselmi analyzed the data and wrote the report. Elisa Maria Galliani is now affiliated with Department of Political and Juridical Sciences and International Studies.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by an MIUR grant [grant number PRIN2012 LATR9N] and by a Departmental grant [SID16_01] awarded to the first author. The funding source had no involvement in study design, in the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the article for publication.
