1. Introduction
This research proposes to explore and analyze the determinants of Romanian job satisfaction with respect to it’s seven development regions. We also intend to define these particular development regions by appealing to job satisfaction determinants, in order to better understand the introspective components that may or may not stimulate local organizational culture. Building job satisfaction may represent a desiderate for any organization seeking to be more efficient and effective, since it is clear that satisfied workers are more productive and more attached to their employer. Several scholars have pointed out that job satisfaction is an important determinant of institutional performance [
1]. Others emphasize the important role that job satisfaction has in any employment-oriented development policy, as it is a central component of individual well-being [
2]. Every employee has their own expectation, beliefs, values, and views, which makes more difficult the understanding of what makes peoples (dis)satisfied with their jobs. For instance, working condition, payments, or procedures implemented in these organizations can satisfy one employee in their work, but may dissatisfy another. Thus, job satisfaction is a complex phenomenon, and countless attempts have been made to define it [
3]. The definition of job satisfaction could be summarized, as it follows: “a positive (or negative) evaluative judgment one makes about one’s job or job situation” [
4]. As such, improvements in job satisfaction are beneficial for both workers and organizations. From the humanitarian perspective, job satisfaction is associated with the physical and psychological well-being of workers, being considered an evaluation instrument for employees’ good treatment within the organization. As such, job satisfaction can be seen as reliable construct which captures employees’ perceptions related to workplace conditions. According to the totalitarian perspective, dissatisfied workers’ behaviour can affect organizational functioning, productivity, and profitability [
5]. Better organizational performance can be achieved when staffed by highly satisfied employees [
6].
According to Judge [
7], the factors that influence job satisfaction are either environmental factors or personal characteristics/traits. Davies et al. [
8] and Bidewell et al. [
9] maintain that there is direct relation between life satisfaction and job satisfaction. Other studies emphasize that employees are motivated by unfulfilled needs, which include esteem (achievement and recognition) [
10], social needs (sense of belonging, appropriate support and protection) [
11], self-actualization (the possibility to develop new skills and reach their fullest potential) [
12], personality characteristics, and behaviors (Big Five personality traits, laziness, loneliness, reservation, person’s connectedness with self) [
13,
14], job/work characteristics (organizational commitment, work atmosphere, experienced stress, recognition of good performance, retiring framework conditions, pay, and work type) [
15,
16,
17]. Beside these factors, family well-being, education, number of previous jobs and residence were also found to influence attitudes about job satisfaction [
18,
19,
20].
Other researchers view job satisfaction as a bi-dimensional construct formed by “intrinsic” and “extrinsic” satisfaction dimensions [
21]. Intrinsic job satisfaction factors are centered on achievement, recognition, responsibility, advancement, growth, and work itself, while extrinsic factors are supervision, working conditions, co-workers, pay, policies and procedures, job security, status, and personal life [
22]. Several other scholars focused on the motivational effects of distributive justice, based on comparisons between the inputs and outcomes of oneself versus those of comparison others [
5]. Using equity theory to explain the concept of job satisfaction, Adams and Freeman [
23] emphasized the fact that an individual becomes satisfied when there is a balance between inputs and outputs, when he is compared with others who are doing similar work. Based on the theoretical and empirical findings from the literature, several job determinants were chosen as predictors in the analysis as they were related to job satisfaction. Also, we include those indicators for which we had data availability for each of the regions included in the study.
Our study investigates the main determinants of job satisfaction among aged citizens, who actively participate in the labor market. We chose this category of individuals in concordance with the National Strategy for Employment 2014–2020 main objectives and actions directed towards increasing labor market participation of the elderly. This category of citizens remains one of the most affected by recent economic crises and restructuring. The Romanian employment rate for older workers of 46.3 per cent in 2018 is still situated below the EU-28 average. According to ILO [
24], the increase in the employment of elder workers appears to be essential in diminishing the deficit of human and professional resources.
The novelties of this paper are two: first, most studies in the field analyze job satisfaction in the developed countries and are focused mostly on individual factors that improve job satisfaction. Our analysis explores job satisfaction in the context of a post-transition economy from the former communist bloc, a developing country which is struggling to attract and maintain its employees [
25], in the context of labor shortages. Moreover, the analysis encases a large variety of indicators (self-esteem, social needs, self-actualization, working characteristics, and personality traits among others), and focused on a specific age category which is currently lacking from most Romanian studies. Second, we apply logistic regressions with average marginal effects, using the latest SHARE-ERIC dataset (Wave 7) filtered for Romania (over 2000 records). This approach is supported using a well-known econometric model applicable to such a phenomenon, namely the binary logistic regression model, developed by Kwon and Remøy [
26], in the context of a developing and post-transition economy from the former communist bloc. Using a SHARE-ERIC survey dataset (2017) for a Romanian sample of people aged 50 and over, we explore and analyze the common and the specific peculiarities of job satisfaction (motivation) in the case of individuals whose reside in one of Romania’s seven development regions.
The paper is structured as follows: in the next section, we present the literature review on job satisfaction and, further, we formulate a series of hypotheses to be tested. In the following section, we detail the data and methods used in the research and then we present the main results. After these steps, we discuss them, and we present the most important conclusions.
3. Data and Methods
This research paper started from a set of questions formulated by SHARE-ERIC (Survey of Health, Ageing and Retirement in Europe—European Research Infrastructure Consortium). In 2017, SHARE-ERIC collected many observations used in this project, including more than 2000 of Romanian citizens aged 50 or more. From this data source, we were interested in eight categories, namely: AC—activities; CC—childhood circumstances; DN—demographics; GV_BIG5—personality traits; RA—retrospective accommodation; RE—retrospective employment; WQ—work quality, GV_ISCED—ISCED standards for classifying education [
64].
We observed several scientific principles that ensure robust research. First we use reliable data sources, which a methodology of transparent and supportive of replication. Secondly, named the triangulation principle, means we rely on many different but convergent approaches, techniques, and tools to properly investigate complex phenomena. Third, we employ the golden rule of cross-validation (using the cvlasso command in Stata) for identifying solid influences and avoid over-fitting (using the rlasso command in Stata using the rigorous LASSO or Least Absolute Shrinkage and Selection Operator procedure).
For merging multiple sources into a distinct data file for each of those eight aforementioned categories, finally we used 1:1 join statements in Stata 16. The next steps consisted in effective data processing, including renaming, numerical scale derivations, and missing values’ treatment, done mostly using functions and facilities available both in spreadsheet programs and in the Open Refine tool. All these steps preceded the use of automatic variable selection procedures, acting as a sort of data mining that was followed by statistical analysis. One of the most important derivations started from lists of declared residences for each respondent, by which we divided the data set by region, specifically seven large Romanian development regions, namely: C or the central region (counties of Mureș, Harghita, Covasna, Brașov, Sibiu, Alba, and Hunedoara), W or the western region (counties of Arad, Timișoara, Hunedoara, and Caraș-Severin), NW or the north-western region (counties of Bistrița-Năsăud, Cluj, Maramureș, Satu-Mare, Bihor, and Sălaj), SW or the south-western region (counties of Mehedinți, Dolj, Gorj, Vâlcea and Olt), S or the southern region (counties of Argeș, Dâmbovița, Prahova, Teleorman, Călărași, Ialomița, Ilfov, and the capital city, Bucharest), SE or the south-eastern region (counties of Vrancea, Galați, Brăila, Buzău, Tulcea, and Constanța), and NE or the north-eastern region (counties of Suceava, Botoșani, Neamț, Iași, Bacău, and Vaslui).
When processing the data for this region (Romania), we aimed for clear and trustful answers, and were also cognizant of traditional treatment procedures for missing values and their effect on classifier accuracy [
65]. Still, we did not assimilate missing values, responses of undecided, or unwillingness to answer to a given value of the original scale, but rather generated an extra grade (usually a middle one—
Table 1,
Table 2 and
Table 3, the value of 2 in the 0–4 derived scale); this came with the cost of artificially generated variance, but ensured a more balanced approach with more realistic values for the coefficient of determination (R squared), the accuracy of classification, and the ratios between the magnitudes of the most powerful resulting influences. The processed dataset served as input for further variable selection procedures and regression analysis in Stata 16.
To analyze the determinant factors that influence the probability of being fully satisfied with the job (compl_satisf_with_my_job was true (i.e., equal to one) when_satisf_with_my_job/wq727 was responded to at maximum Likert value (four) and otherwise false (i.e., equal to zero), as seen in
Table 1,
Table 2 and
Table 3) in our proposed models, we have started from a well-known econometric model (Equation (1)) applicable to such a phenomenon, namely the binary logistic regression model [
26]:
where:
p is the probability of being satisfied with the job;
k is the total number of independent variables, k = 2, …, m;
βk is the effect of a change in variable Xk on the probability of the analyzed state of the outcome (being satisfied with the job);
X
k is one explanatory variable (Equation (1)) from the array (∑) of the features in
Table 1; and
ε represents the error term.
Binary logistic regressions have been used to support robustness checks of the dual-core, and, to confirm particular regional influences, cases were filtered from among the remaining influences with respect to: lower p values corresponding to the size of errors when compared to that of coefficients; lower VIF (Variance Inflation Factor) values as proofs of a lack of collinearity—usually less than 10 [
66]; higher values resulting from goodness-of-fit (GOF) tests [
67] both for p values (to reject the null hypothesis) and chi square; higher values for the Area Under the Curve of Receiver Operating Characteristic, known as AUCROC, AUROC or shortly, ROC [
68] and indicating the accuracy of classification for a scenario/model; and larger R-square values which suggest better explanatory power for the resulting models.
The descriptive statistics, containing the list of variables selected for this study, are available as two subparts (
Table 2 and
Table 3), with four subsets each, which have been presented in descending order of the total number of observations for each. More details and explanations about these variables are available in
Table 1. All study sites (
Table 2 and
Table 3) reveal, from the very beginning, noticeable differences in terms of average intensity of the primary outcome and several possible predictors assumed to be most related with the phenomenon.
Table 2 and
Table 3 present the summary statistics for the entire dataset and the first three development regions and the other remaining four ones.
4. Results and Discussion
In terms of merging the original vertically partitioned data subsets, we finally chose a 1:1 merge statement in Stata 16. For cleaning the data and performing additional derivations, in most cases, the spreadsheets’ immediate visual feedback and insight, powerful built-in and user-defined functions, customizable filters for particular text patterns, and fast autofill and split-text-to-columns facilities, were deemed adequate for performing manual cleaning tasks.
We also filtered the resulting dataset, corresponding to Wave 7 (2017), by considering the original field W7-ac_country set with the value of “Romania”, for which we had 2144 unique records. Then, we conditioned the variable dn003_ (birth year) at less than or equal to 1967 to retreive only those responses of Romanian people aged 50 or more, of which there were 2056. Next, to identify the specific subsets corresponding to those seven development regions in terms of last declared residence, we started from those 2056 filtered records above and used a cascade of IFs to generate the last non-blank residence by considering the related fields (from ra025c_1 to ra025c_30). Thus, we obtained 2052 records corresponding to the observations mentioned in
Table 2,
Table 3,
Table 4 and
Table 5. The difference with the previous filtered amount of 2056 records consisted in only four observations with unspecified last residence, which were consequently dropped.
We performed binary logistic regressions for Romania overall and for its seven development regions and preserved only those influences satisfying selection rules depending on: significance (p), VIF, GOF, AUCROC, and R-sq. values. Next, we performed post-estimations and reported the average marginal effects (not raw coefficients) to ensure support for comparability when the magnitude was concerned for both intra- and inter-scenarios/models’ comparisons. These average marginal effects have been reported in two subparts (
Table 4 and
Table 5) in descending order of the total number of observations for each regional subset resulting from last residence and considering two scenarios for each: only the dual-core (a, c, e, g, i, k, m, o) and the dual-core plus particular influences (b, d, f, h, j, l, n, p).
Table 4 and
Table 5 present Romania overall and the first three regions and, also, the last four comparable regional models, respectively, in terms of average marginal effects on job satisfaction, with the specifications that: (1) the source was represented by its own calculations in Stata 16, and (2) standard errors were between round parentheses; (3) *, **, ***, **** indicate significance at 10%, 5%, 1%, and 1‰.
Assuming predictors’ potential high correlation and models’ overfitting as usual reasons for too-high R squared values, we performed additional post estimations, such as the maximum correlation coefficient in predictors’ matrices for each model in
Table 4 and
Table 5. All these values were well below the limit of 0.7, beyond which we usually observe a high correlation between predictors. Moreover, we assessed the computed VIFs in OLS (ordinary least squares regressions) against dynamic thresholds acting as maximum acceptable values, depending on models’ explanatory power (1/(1–R
2)). Our results indicated that all models in
Table 4 and
Table 5 met these threshold conditions. They reconfirmed the lack of multicollinearity.
Additionally, we applied the cvlasso and the rlasso commands for each already reported model. The latter, rlasso (or the rigorous LASSO variable selection procedure), is a well-known penalizing method to control overfitting and it removed none of the variables already reported in
Table 4 and
Table 5.
Moreover, we compared the values for the explanatory power (pseudo R^2), the Wald statistics (Wald chi^2(2)), and the maximum correlation coefficient between predictors for concurrent ordinal logit models when considering the reverse causality that often causes endogeneity problems. In these tests, we have first chose variables in the dual-core and a variable to analyze that corresponds with job satisfaction. Then we replaced the latter with each of the two (in their original form, meaning on a scale) by interchanging their roles. The results indicated better scores (largest Wald chi^2(2), pseudo R^2, and lowest correlation coefficient) when considering job satisfaction as the variable to analyze and the other two as predictors, rather than vice versa (e.g., work atmosphere/work gave recognition as outcome).
The results, after performing logistic regressions revealed interesting relationships. In the overall model (N = 2052,
Table 4), scenarios (a) and (b) underline the most powerful and significant dual-core of the models and rankings in terms of importance (both intrinsic and extrinsic factors; recognition and working condition; good atmosphere between colleagues). The model based only on this dual-core (
Table 4, scenario a) is powerful enough in terms of accuracy of classification (good to excellent for an AUCROC value of 0.8893), while the maximum VIF < 6 (way below 10) leaves enough room to further identify and consider other specific influences.
The value of Pseudo R square (0.3839) indicates good explanatory power for this simplified model. Therefore, these very promising results for the overall dataset encouraged us to explore further and test several hypotheses to identify particular patterns for each of the seven Romanian development regions. Hence, we additionally discovered the positive role exerted by existing opportunities to develop new skills, appropriate support in difficult situations at the workplace, protective measures from authorities in case of health hazards, how often the individual feels that life is full of opportunity, the manifestation of openness as a Big Five individual feature, family financial well-being, and the urbanity level of the respondent’s residence. The only negative influence discovered for this overall model (and most comprehensive set of specifications (
Table 4, scenario b)) is afferent to loneliness in childhood. The Pseudo R square (0.4266) increased significantly when adding these eight influences, as did the accuracy of classification (0.9059), which then becames excellent (AUCROC > 0.9), while collinearity remained acceptable (maximum VIF < 10).
First, in terms of regional models, when considering the north-west region of Romania, the dual-core is confirmed in the same order of its two components as in the overall model. In fact, this finding was valid for all regional models. Besides these basic findings, we emphasized other positive and negative influences on job satisfaction. Hence, the only negative findings corresponded to loneliness and an individual’s choice to perform paid work after retirement, suggesting that a low esteem behavior indeed negatively influences job satisfaction. Positive influences on job satisfaction were displayed by individuals who had the opportunity to develop new skills, who had retired from their first job, and those who considered themselves open persons. In this regional model, the explanatory power (pseudo-R square of 0.4528 and 0.5452), accuracy of classification (excellent for AUCROC of 0.9137 and 0.936), and accepted collinearity (maximum VIF of 6.26 and 9.82) indicated better values for both the simple and the most comprehensive sets of specifications (
Table 4, scenarios c and d) than the ones obtained for the overall model above.
Second, for the adults from SE, besides the accent put on the previous dual-core foundation, other positive influences for job satisfaction corresponded to just two other variables: interpersonal trust and opportunity. More precisely, based on these two additional findings, we can state that the higher the level of trust in other people, the greatest career satisfaction. In addition, the more often a person has the belief that life is full of opportunity (optimism and ambition as traits), the higher the level of career satisfaction. For this particular region, all three indicators (pseudo-R square, AUCROC, and maximum VIF) recorded good-to-excellent values (0.3793 and 0.4038; 0.8846 and 0.8952; 6.42 and 7.89), close to those of the overall model.
Third, in NE Romania, besides the already identified dual-core, which indicated the lowest magnitude for its first component, deserved recognition for, in comparison with the rest of regions. Some other predictors are underlined below. The only negative influence was laziness, suggesting that the job satisfaction becomes low in the case of a person who tends to be lazy and avoid high-pressure jobs. Accordingly, a rational individual who understands the negative consequences of idleness on the organizational status-quo and personal achievement tends to experience lower levels of career satisfaction. The additional positive determinants, except the dual-core, were: no previous jobs (number_of_jobs), family finance (family_financial_welfare), job protection and security (health_protect_from_authorities), and meaningfullness of life (how_often_feel_life_has_no_sense). Consequently, for this specific region, we can state that the more an individual has changed jobs, the more career satisfaction they find. The same positive relationship manifests when the well-being of one’s family is greater, when employers regularly take appropriate measures toward workplace safety, and the frequency with which one feels their life is meaningless. The influence of the last variable is quite strange, suggesting that meaning in life plays a negative role on career satisfaction, since much other research emphasized an opposite role of meaning in life in life satisfaction [
69]. We can understand this influence only in the following terms: sometimes, compromises in one’s personal life may lead to unexpected career success, but mostly because of a so-called work-life (im)balance. For this particular area, especially when considering the second, most comprehensive scenario (
Table 4, scenario h), all three indicators above (pseudo-R square, AUCROC, and maximum VIF) recorded better values than in the case of the overall model.
Fourth, SW Romania revealed other interesting influences. Here, we have not identified any negative predictors of job satisfaction (
Table 5, scenario j), only positive ones, for example, better familial economic circumstances impact career satisfaction. Also, personality traits, such as openness and being reserved, seem to be positive predictors. If a person holds the belief that life is full of opportunity more strongly, then their career satisfaction increases. Finally, feeling increasingly helpless over the circumstances of one’s life positively influence job satisfaction, revealing a particular attitude in such individuals: they put a greater emphasis on external loci of control when expressing career satisfaction, feeling that they are not agent in their. For this region (
Table 5, scenarios i and j), all three indicators above (pseudo-R square, AUCROC, and maximum VIF) recorded better values than in the case of the overall model.
Fifth, Romania’s southern region presents a series of elements that differ in comparison with other regions. For instance, feelings of diminishing agency in one’s life negatively influenced job satisfaction, suggesting a greater emphasis on an internal locus of control, a finding in contradiction with those in the SW region. Family welfare represented the factor of greatest magnitude, and was found to positively influence career satisfaction. Novel to the southern region model was the positive role played by urbanity, suggesting that the individuals who have their permanent residence in urban areas, here, are more satisfied with their career, we surmise due to better opportunities, infrastructure, and mobility. For this region and especially with consideration to the second, most comprehensive scenario (
Table 5, scenario l), all three indicators above (pseudo-R square, AUCROC, and maximum VIF) indicated better values than the same in the overall model.
Western Romania expressed different attitudes towards career satisfaction. Those who retired from the only job they had ever had were more satisfied with it (loyalty to the workplace, commitment, and orientation on a very long term). Moreover, those who had graduated from high school were more satisfied with their career, as were those who felt neglected more often, which is likely a function of balance between work and life. In addition, for this particular region, attitudes of laziness and loneliness betray low self-esteem and negatively affect job satisfaction. Though good-to-excellent in terms of accuracy (AUCROC = 0.8708,
Table 5, scenario n), with a decent explanatory power (pseudo-R square of 0.3415) and acceptable collinearity (maximum VIF of 6.01), this regional model had the lowest value of the second indicator (colleague_good_atmosphere).
Finally, central Romania emphasized other specific predictors and different additional influences. Our most intriguing finding, here, regards a negative relationship between life satisfaction and job satisfaction, a finding that contradicts much other research [
58], but which may be explained with reference to work-life balance [
70]. In addition, analysis of introspections about laziness in this population challenged established perspectives; those who regarded themselves as having become lazier also expressed more job satisfaction. As correlation, this was also found by Dalal [
71]. Several positive influences were related to the economic security of the family and to introversion (in our case, reserved personalities). For this region, especially when considering the second, most comprehensive scenario (
Table 5, scenario p), all three indicators (pseudo-R square, AUCROC, and maximum VIF) recorded better values than those of the overall model.
The results above (
Table 4 and
Table 5) clearly indicate that for respondents from the western region of Romani (W), receiving deserved recognition mattered most (it had almost the same magnitude of effect as workplace atmosphere) when compared with the other regions (confirmation of the second hypothesis, H2). We expected that, in light of our previous research, results from this region, which geographically included Transylvania under the former Habsburg occupation, showed patterned differences with respect to immigration and moral attitudes regarding generational differences and job satisfaction in those with considerable work experience [
60]. The additional results in
Table 5 (the most comprehensive scenarios of our regional models) clearly indicate further evidence for impact of financial well-being of the family (family_financial_welfare) on job satisfaction. We analyzed this further, together with the dual-core, by using a binary logistic model; family_financial_welfare passed all checks but not with the same significance as the dual-core (usually ** or ***), when considered together with it. Consequently, we consider familial financial well-being an interesting predictor of job satisfaction (and confirmation of the second part of the third hypothesis, H3); better familial economic security may positively influence career satisfaction due to increased psychological comfort, greater confidence in personal ability and aptitude, and less career uncertainty. Moreover, these results emphasized that life satisfaction is not always positively correlated with job satisfaction (a rejection of the first part of H3), in some cases acting as a negative predictor of job satisfaction (
Table 5, scenario p). In addition, the influence of laziness in W opposes C, which showed a positive influence, and NE and W, where it showed a negative influence on full job satisfaction; we take this to be a partial rejection of H4. The positive influence of interpersonal trust (although only in respondents from SE) partially confirmed H5. Further, degree of urbanity was positively correlated with job satisfaction (
Table 4—scenario b, the overall model;
Table 5—scenario l, region S). In addition, we must emphasize the positive influence associated with the respondent’s number of jobs held (
Table 4—scenario h, NE).
The results in
Table 6 and
Table 7 show the strongest influences, which persisted even when the regression model was changed (i.e., using ordinal logistic regression operating on the entire 0–4 scale of job satisfaction instead of the binary logistic ones).
Though not reported in this paper, additional regressions were performed for robustness checks of the two most powerful influences identified above (the dual-core), who’s strength we ascribe to having used holistic questions. These checks considered three criteria about the respondent (marital status, gender, and having graduated high school) and revealed a strong dual-core irrespective of how the overall dataset for Romania was divided (whether along marital status, sex, or education). All these findings have been obtained by assimilating missing, undecided, or unwilling responses to a added median grade (i.e., 2 in a final scale from 0 to 4). Still, when considering just the dual-core in the overall model, in each of the seven regional models, and in the splits (subsets) above meant for robustness proofs, the maximum VIF left room for exploring and adding new specific influences on this dual-core foundation. Therefore, we consider our first hypothesis (H1) fully confirmed.
In future research, we intend to explore the peculiarities of other European regions in terms of both job and life satisfaction, starting from the datasets provided by the SHARE-ERIC consortium.