Happy Cultures? A Multilevel Model of Well-Being with Individual and Contextual Human Values

Despite the abundant literature in Happiness Science, no paper to date has studied the joint effects of human values on subjective well-being at individual and contextual level. Using European Social Survey data for life satisfaction and Salomon H. Schwartz’s scale for human values with four and ten dimensions, this paper presents novel evidence on the direct effects of individual, regional, and national human values, utilizing two different ways of building cultural indicators of human values. We show that regional factors explain approximately 2% of the dispersion of individual life satisfaction, whereas national factors explain around 12%. The results on the effects of individual human values support Sortheix and Schwartz’s hypothesis, with a significant difference: Individual Conformity has a positive impact on well-being, not the negative sign Sortheix and Schwartz predict for Conservation values. We also find positive direct cultural effects for Benevolence and Conformity and negative effects for Tradition. Additionally, we propose a research agenda for human values and contextual effects on well-being studies.


Introduction
One contributor to differences among the social sciences is each discipline's approach to individuals in a social context. A number of subdisciplines that converge on the new science of well-being have appeared on the margins of social science disciplines. These subdisciplines include social, cultural, and political psychology; socioeconomics; and cultural sociology; among others. Despite expansion of this new science, major questions remain persistently elusive. What role do individual and cultural factors play in determining wellbeing? Although standard definitions of culture mention shared values, norms, and beliefs or expectations, the topic is challenging even if we reduce the focus to human values: Do specific cultures, or cultural attributes, promote individuals' well-being?
To limit the scope of this study, we start from Salomon H. Schwartz's theory of human values (Schwartz, 1992;Sagiv & Schwartz, 2022), which details ten basic human values. The proper empirical method for modelling the joint effects 1 of individual and cultural human values on individuals' well-being is to estimate multilevel models. Such models are not a panacea, but they are useful tools for analyzing the direct effects of contextual and individual attributes on human behavior or feelings (Duncan et al., 1998). Surprisingly, no previous paper has used multilevel models to evaluate the joint individual and cultural effects of Schwartz's human values on subjective well-being.
The goal of this paper is to study the direct effects of individual and cultural values on subjective life satisfaction, using Schwartz's scale of human values and multilevel models for European Social Survey (ESS) data. We study these effects using two levels of aggregation of human values and two levels of geographical aggregation in Europe, regions and countries. We also compare our results using two different ways of aggregating human values: the means for individual values in each geographical area and the percentage of people with high scores for each value in that territory.
We are trying to encourage debate on detection of cultural attributes that promote the subjective well-being of the average individual within a territory. As discussed below, investigating the positive effects of some cultural attributes on the average individual is compatible with acknowledging the existence of darker aspects of those cultural attributes. We build on Sortheix and Schwartz's (2017) hypotheses about the effects of individual human values on well-being. Although Morrison and Weckroth (2018) studied the effects of four Schwartz's basic human values on well-being at an ecological level, however, no previous research has analyzed individual-and-cultural-level effects jointly. Indeed, Veenhoven's (2015) survey on social conditions for human happiness ignores cultural human values. Therefore, this paper focusses more on detection of statistical cultural patterns than on modelling social well-being (Keyes, 1998), open societies (Krys et al., 2019), or modernization (Welzel & Inglehart, 2010). Our inquiry is related to Ye et al.'s (2015) analysis of the effects of Geert Hofstede's (1980) nine cultural dimensions, although their study uses ecological data. Evaluating the individual and aggregate effects of ten human values, we conclude that happy cultures are characterized by three of these values: high levels of Benevolence (help others and care for them) and Conformity (follow rules and behave properly), and low levels of Tradition (respect traditions and customs). We highlight these three cultural values to organize the paper's ideas.
Our main contributions are the following. Firstly, we estimate a joint multilevel model for Schwartz's ten basic human values rather than ten separate models for each human value (Sortheix & Lönnqvist, 2014;Sortheix & Schwartz, 2017). Secondly, we add two indicators of aggregate human values to this model to attempt to capture direct cultural effects. Thirdly, we compare cultural effects at regional and national level for the four higher-order value dimensions.
Our main results are as follows. Firstly, comparison of our regional and national effects shows that countries explain about 12% of variance in life satisfaction, whereas regional heterogeneity only explains an additional 2%. This result on the greater significance of the national level is not surprising, but its quantification is novel and has not been measured explicitly in previous multilevel regional papers on well-being (Pittau et al., 2010;Aslam & Corrado, 2012;Ballas & Tranmer, 2012;Neira et al., 2018).
Secondly, our findings on the effects of individual human values in individual wellbeing are consistent with Sortheix and Schwartz's (2017) hypotheses, with some significant differences. In the higher-order dimension Conservation, we find a negative effect of Tradition (as expected) but no significant effects of Security. We also find positive effects of Conformity, at both individual and national level. Thirdly, for first time in the literature, we show that individuals' well-being increases in cultures with high emphasis on Benevolence and Conformity and low emphasis on Tradition. The results on Conformity may be the most controversial, and they form part of our proposed research agenda.
The rest of the paper is organized as follows. Section 2 reviews the literature on human values, culture, and well-being, with emphasis on multilevel models. Section 3 describes the methodology and data. Section 4 discusses the results. Finally, Sect. 5 draws conclusions and proposes a research agenda. Fig. 1 The circular structure of 10 basic values, four higher order values and two underlying motivational sources. Source: Sortheix and Schwartz (2017) Sortheix and Schwartz (2017)

and on ESS documentation
Conservation: Values that emphasize order, self-restriction, and resistance to change Security-Safety, harmony, and stability of society, of relationships and of self v05-Important to live in secure and safe surroundings v14-Important that government is strong and ensures safety Conformity-Restraint of actions likely to upset others and violate social expectations or norms v07-Important to do what is told and follow rules v16-Important to behave properly Tradition-Respect, commitment and acceptance of the customs and ideas that traditional culture or religion provide v09-Important to be humble and modest, not draw attention v20-Important to follow traditions and customs Openness to Change: Values that emphasize independence of thought, action and feeling and readiness for change Self-Direction-Independent thought and action-choosing, creating, exploring v01-Important to think new ideas and being creative v11-Important to make own decisions and be free Stimulation-Excitement, novelty, and challenge in life v06-Important to try new and different things in life v15-Important to seek adventures and have an exciting life Hedonism-Pleasure and sensuous gratification for oneself v10-Important to have a good time v21-Important to seek fun and things that give pleasure Self-enhancement: Values that emphasize pursuit of one´s interests, relative success, and dominance Achievement-personal success through demonstrating competence according to social standards v04-Important to show abilities and be admired v13-Important to be successful and that people recognize achievements Power-Social status and prestige, control or dominance over people and resources v02-Important to be rich, have money and expensive things v17-Important to get respect from others Self-transcendence: Values that emphasize concern for the welfare and interests of others Benevolence-Preservation and enhancement of the welfare of people with whom one is in frequent personal contact v12-Important to help people and care for others well-being v18-Important to be loyal to friends and devote to people close Universalism-Understanding, and protection for the welfare of all and the environment v03-Important that people are treated equally and have equal opportunities v08-Important to understand different people v19-Important to care for nature and environment 2 Theoretical Framework Schwartz (1992) defines values as broad motivational goals that serve as guides to behavior and as criteria for judging people and events. Each value is defined by the goals towards which it is directed, that is, the motivation it expresses. Figure 1 presents the ten basic values and the goals that define them. Values form a circular structure, in which the more compatible any two values are, the closer they are on the circle, and the more they conflict, the farther apart. Table 1 below shows additional details on Schwartz's scale of human values and how they are measured in the ESS. The ten basic values may be grouped into four higher-order values that summarize the opposition between competing values: Conservation versus Openness to Change in the growth dimension and Self-Enhancement versus Self-transcendence in the self-protection dimension. The growth dimension involves relation to self-restriction and order versus novelty and independent thought. The selfprotection dimension captures motivation towards personal interests versus towards others. Moreover, Conservation and Self-transcendence are higher-order values with a social focus, whereas Openness to change and Self-enhancement have a personal focus.

Schwartz's Scale of Individual Human Values and Subjective Well-Being
Although the direction of causality is always debatable, 2 this study focuses on the effects of individual human values on subjective well-being, as surveyed by Schwartz and Sortheix (2018) and Sagiv and Schwartz (2022). Sortheix and Schwartz (2017) argue that growth-oriented values with a personal focus tend to contribute positively to well-being, whereas values oriented to self-protection that have a social focus tend to contribute negatively. Using the four higher-dimension values, therefore, Openness to change values (which combine growth orientation and personal focus) relate positively to subjective wellbeing, whereas Conservation values (which combine self-protection orientation and social focus) relate negatively. As to the effects of Self-enhancement and Self-transcendence on well-being, the opposing implications of growth-protection and social-personal orientation roughly balance each other. As to the four basic values underlying Self-enhancement and Self-transcendence, Sortheix and Schwartz (2017) expect either that one sign will dominate or that no association will occur.
Our main results for individual human values (presented below) are generally consistent with this theory, although we estimate the joint effects of all values, whereas Sortheix and Schwartz (2017) estimate a different model for each of the ten basic values.

Culture as Aggregation of Individual Human Values
Most of the literature cited in this paper comes from psychology and approaches wellbeing from the point of view of individuals' attributes. This approach is also standard in economics. As in other subfields related to socioeconomic development, however, part of the discussion of well-being has been oriented to distinguishing economic and cultural factors (Schyns, 1998). Political science points to institutional factors (Inglehart & Klingemann, 2000), whereas sociology emphasizes the role of social contexts and socialization (Hitlin & Piliavin, 2004;Veenhoven, 2015). One line of research on the value-well-being relationship focuses on the congruence between individuals' values and those of the environment (Sagiv & Schwartz, 2000;Morrison & Weckroth, 2018). This line connects with our cultural approach, although our study seeks to distinguish between the direct effects of individual and cultural human values.
Culture may be defined as the set of values, norms, and beliefs (or expectations) shared to some degree in a social group. Minkov and Hofstede (2014) view culture as values, beliefs, attitudes, and other patterns of thinking and feeling, which provide the software of the mind. Socio-spatial culture refers to the broader societal traits and relations that underlie the prevailing mind-sets and overall way of life in specific places (Huggins & Thompson, 2019). One preliminary question therefore concerns the interpretation of individual and shared values. A second involves how to operationalize this sharedness empirically.
Schwartz's theory of human values is especially suitable for studying the direct effects of individual and cultural values on well-being because it has adopted a cross-cultural perspective from the beginning (Schwartz, 1992;Goodwin et al., 2020;Sagiv & Schwartz, 2022). Debate exists, however, about possible differences in the meaning and structure of values at individual versus aggregate level (Fischer & Poortinga, 2012). A non-isomorphic relationship between the two levels arises if the individual-level value structure is primarily determined by individual-level variables and the country-level value structure primarily by factors that apply at country level. Hofstede states explicitly that his value scale apply only to countries. In a recent seminar, Schwartz declared this an unresolved question, 3 arguing that nationally representative samples enable us to measure what people are exposed to, that is, what we might call citizens' values rather than cultural values.
We are skeptical about these arguments. Generally speaking, well-being research is trying to estimate the linear effects of endogenous and correlated indicators of human values at individual and aggregate level, treating categorical variables as if they were continuous, ignoring multiple relevant interactions, ignoring many issues of culture (polarization, subcultures, etc.), on a topic with different approaches in different disciplines, under competing theories of human values, and using methodologies with different limitations. Theory-driven empirical inquiries are thus crucial. The search for statistical patterns using alternative methodologies helps, however, to define stylized facts useful to developing theories. Our paper follows the latter direction. For instance, we check the robustness of our results for culture comparing two different indicators of the sharedness of human values, as detailed in Sect. 3.2.
Moreover, Tov and Diener (2009) conclude both that there are pancultural experiences of subjective well-being that can be compared across cultures and that culture-specific patterns make cultures unique in their experience of well-being. This argument supports analysis of heterogeneity in multilevel models of well-being, as noted in the methodological section below. 4 Most studies of the effects of culture on well-being use ecological data, however (Suh et al., 1998;Ye et al., 2015;among others). Surprisingly, little multilevel research analyzes the individual and aggregate effects of human values on well-being, even though multilevel models are the proper tool for preventing or reducing the possibility of atomistic and ecological fallacies (Snijders & Bosker, 2012). Below, we summarize some studies involving multilevel research.

Happy cultures? Three Human Values as Cultural Attributes: Benevolence, Conformity, and Tradition
This study focuses on cultural direct effects that promote the well-being of individuals living in a specific social environment. We operationalize the concept of a happy cultural attribute as a positive, significant estimate of an aggregate indicator of human values in a multilevel model of individual well-being. This interpretation is compatible with the darker aspects of happy cultures, even for cultures that stress Benevolence. One positive side of cultural Benevolence may be related to enjoyment of social relationships, as the literature on bonding social capital argues (Neira et al., 2018;Bruna & Rungo, 2020). The details in Table 1 may, however, reveal complex effects. To give a dramatic example, Mafia members from family-oriented cultures will score high on item v18 but low on v12. 5 Moreover, a culture that values Conformity may encourage well-being because of clear and accepted collective rules. Huxley's 1932 Brave New World exemplifies one potential dark aspect of Conformity. 6 Further, Andreoni et al. (2021) argue that norms detrimental to social welfare sometimes persist when societies are locked into a conformity trap. Additionally, our models say nothing about how each culture treats minorities or people with deviant values.
We are not proposing a model of positive social health, such as Keyes' (1998) concept of social well-being. Our research question about happy cultures is atheoretical. We ask whether social aggregates of Schwartz's ten basic human values can affect individuals' well-being directly. Our motivation is similar to that of Ye et al. (2015), who studied the effects of nine dimensions of cultural characteristics on subjective well-being, although they based their research on Hofstede's (1980) cultural dimensions and did not use multilevel models. Ye et al. (2015) find that Power distance and Gender egalitarianism play the most important and stable role in determining well-being.
Another more theoretically oriented strand of literature tries to establish empirical links between a particular cultural trait and individual well-being. A first example of this approach is Krys et al.'s (2019) concept of open society, based on four values: trust, tolerance, civic engagement, and nonmaterialism. Using multilevel models, these authors conclude that open societies foster satisfaction and that the positive effect of Hofstede's individualistic cultural attribute disappears when controlling for the effects of open societies. A second example of this approach is Li et al.'s (2021) multilevel model, which reveals that cultural differences on the liberal-conservative continuum moderate the relation between gender inequality and well-being.
A final example of a theory-driven search for cultural effects is Ronald Inglehart's theory of social change. Welzel and Inglehart (2010) and Minkov et al. (2020) use multilevel modelling to analyze a specific causal channel in which modernization affects individualism and emancipative values. Inglehart and Klingemann (2000) show that well-being is 5 The possible conflict between v12 and v18 is analogous in Hofstede's human values scale, in which prioritizing family over work is considered a sign of individualism, increasing the indicator of individualism in family-oriented Latin countries, and reducing it in individual-oriented Anglo-Saxon countries. Morrison and Weckroth (2018) note that Benevolence may well strengthen the in-group, but at the expense of the outgroup. 6 Under an authoritarian or democratic regime, Conformity might make the average individual happier. Indeed, some degree of conformity may be necessary to preserve social equilibrium, since all forms of society have their social problems. This need is one reason why all socialization processes seek to promote adherence to social rules, which is related to the preservation of social structures of power. However, examination of this question is beyond the scope of this paper. related to economic development, Protestant societies, and democratic institutions. Though our approach is not based on modernization theory, the empirical part of this paper follows Sortheix and Lönnqvist's (2014) focus on the role of contextual socioeconomic development in discussing the robustness of our results when estimation is controlled by a general indicator of socioeconomic development.

Regional and National Heterogeneity and Multilevel Studies
Part of the debate about culture relates to the territory in which culture is evaluated. The interdependence of personality psychology and socio-spatial culture determines the outcomes of urban and regional development (Huggins & Thompson, 2019), which may affect subjective well-being. Weckroth and Kemppainen (2016), for example, use European regional data to show that Self-direction is a strong predictor of GDP. Morrison and Weckroth (2018) were the first to study the effects of Schwartz's basic human values on wellbeing at an ecological level. They study four human values and compare individuals from Finland's metropolitan region to those from the rest of the country. More generally, however, the degree of cultural heterogeneity across territorial levels is controversial.
Using the ESS and regional indicators of Hofstede's cultural dimensions, Kaasa et al. (2013) conclude that there are striking cultural differences between regions in most countries and that the deeper the subdivision, the larger the differences, which are masked when one indicator is used for the whole country. Conversely, Minkov and Hofstede (2014) also use ESS data to conclude that regions are acceptable but in no way imperative units of analysis. Nations can provide sufficient pictures of cross-cultural variation, at least in terms of values. This last conclusion supports Welzel and Inglehart's (2010) conclusion that cultural differences between nations are considerably greater than those within nations.
Focusing on multilevel studies of well-being, Schyns (2002) analyzes individual and national factors. A few papers have considered regional level: Pittau et al. (2010), Ballas and Tranmer (2012), Aslam and Corrado (2012), and Neira et al. (2018). This literature does not distinguish the proportion of well-being indicators explained by regional versus national factors. The first research question of the empirical part of our paper uses ESS data to detect the most significant level of geographical aggregation for studying cultural effects of well-being.

Methodology
Using the lme4 package for R (Bates et al., 2015), we estimate two-and three-level multilevel (mixed, hierarchical) models in which the one-level data are for individuals i. 7 Some of our models have a second level for regions ( j ) or countries ( k ), and the three-level models nest regions ( jk ) within countries ( k ). We estimate random intercept models, which capture between-context heterogeneity (Duncan et al., 1998;Neira et al., 2018) to account for aggregate differences in the dependent variable: the territory in which one is located tends to affect subjective well-being. We thus do not consider random slopes or cross-level interaction, which could model between-individual heterogeneity, in which similar types of people are behaving/feeling differently in different types of places.
Our models are linear. Despite the categorical dependent variable, prior literature summarized by Bruna and Rungo (2020) shows that ordered probit models have results that are qualitatively similar to linear models but harder to interpret. We thus maintain the linear assumption, standard procedure in this literature.
We estimate three types of equations. The empty or null model uses only a global intercept to explain the dependent variable, which is an estimate of the global mean of life satisfaction (see Table 3 below). Random effects enable estimation of region-and countrylevel deviations of life satisfaction from this global mean. 8 To represent the variance of the residuals at different levels, we use the following notation: 2 for individual level, 2 j for regional level, and 2 k for national level. 9 Estimating these variances enables us to calculate the proportion of variance of life satisfaction due to (measured or unmeasured) factors common to individuals in the same region or country. This proportion is called the intraclass correlation coefficient ( ICC ) and is defined as follows in our two-level models: for regional level; for national level: Our three-level models include random intercepts for regions nested in countries. We present two different ICC s. The first is the national coefficient ICC k , as noted above. The second is the joint ICC of regions nested in countries, as follows: 10 After the null models, we estimate three-level models including individual-level variables ( X ijk ). This method controls for compositional effects: regions or countries could have different propensities to life satisfaction merely because they have different types of individuals. After considering those individual-level variables, we add aggregate variables for regions ( X jk ) and countries ( X k ). Table 4 below shows the results of these models, considering the four higher-order human values from Schwartz's scale. The final models in the paper, presented in Table 5 below, are two-level models including national random effects and focusing on the ten basic human values in Schwartz's scale at individual and country level.
Estimated group effects are weighted averages that combine information from the group itself with information from the mean for all groups. Estimating with random effects is a conservative approach that gives less weight to less reliable data. The random effects are precision-weighted residuals called 'posterior residuals', 'empirical Bayes estimates', or 'shrunken residuals'. 9 The number of national units in the present study (23) should not bias the results (Elff et al., 2021). 10 Technically, the variance partition coefficients ( VPC ) are the proportion of response variance at each level of the model hierarchy. Therefore, although VPC k = ICC k , VPC j = 2 j ∕( 2 j + 2 k + 2 ) is not the same as the ICC jk defined in the text.
Although cross-sectional multilevel models cannot distinguish causality between characteristics of individuals and cultures, they provide a major advantage over traditional onelevel studies, which do not differentiate within-context from between-context effects. We therefore use multilevel models for general assessment of possible direct effects of individual and cultural human values on life satisfaction.
This study estimates unweighted multilevel models. "Few things are as confusing to applied researchers as the role of sample weights" (Angrist & Pischke, 2009, p. 91). The decision whether and how to weight depends on the research goal. Focusing on cultural effects, it is not useful to give more weight to more populated countries. Following the suggestion of Solon et al. (2015) for causal analysis, we compare some of our results to those derived from weighted estimation (Table 6). 11

Data
We used the eighth round of the ESS to model life satisfaction, an indicator of well-being closely related to cognitive judgments about feeling fulfilled in life or living a good life. Our dependent variable is an indicator resulting from the survey item 'How satisfied are you with life as a whole?' The responses range on a scale from zero (extremely dissatisfied) to ten (extremely satisfied). 12 The regional classification follows the Eurostat Nomenclature of Territorial Units for Statistics (NUTS), which establishes four aggregation levels, from countries (NUTS 0) to the smallest harmonized territorial units (NUTS 3). 13 Israel has no regional division, so its data is defined at NUTS level 0). For the other countries, our analysis uses a mix of three NUTS aggregation levels (NUTS 1, 2, and 3), as do Aslam and Corrado (2012) and Neira et al. (2018).
Schwartz's "Computing Scores for the 10 Human values" 14 explains methodological details for computing human values scores in the ESS. Since this survey includes 21 items, 15 arithmetic means are used to transform the answers to those items into 10 indicators of basic human values, which can be aggregated to 4 higher-dimension human values. This method treats categories from Very much like me to Not like me at all as if they were values for continuous variables, with subjectively comparable significance across the 21 items of the 10 dimensions. Schwartz also recommends calculating each value relative to the sum of the individual's scores for all values. This procedure has pros and cons, 11 We used the ESS analysis weights (anweight), which correct for different selection probabilities within each country, as specified by sample design for nonresponse, noncoverage, and sampling error; and considering differences in population size across countries. The United Kingdom, Germany, France, and the Russian Federation constitute 21.9% of the observations in our sample and have a weight of 59.5% using anweight. 12 The dependent variable used by Sortheix and Lönnqvist (2014) and Sortheix and Schwartz (2017) is an average of indicators of life satisfaction and happiness. Neira et al. (2018) show some significant differences in the determinants of life satisfaction and happiness. 13 https:// ec. europa. eu/ euros tat/ web/ nuts/ backg round. 14 See also the following link: http:// essed unet. nsd. uib. no/ cms/ topics/ 1/4/ 4. html. 15 For instance, item v07 of the questionnaire on human values reads as follows: 'Now I will briefly describe some people. Please listen to each description and tell me how much each person is or is not like you. Use this card for your answer. She/he believes that people should do what they're told. She/he thinks people should follow rules at all times, even when no-one is watching'. The options are the following: Very much like me; Like me; Somewhat like me; A little like me; Not like me; Not like me at all; Refusal; Don't know; and No answer. however. The arguments against include the fact that we do not know why some individuals have extreme opinions and others have moderate ones. There is no reason to assume that personal priorities on any topic should be evaluated by centering values on the average score for all values of each individual. Moreover, multilevel models are usually estimated with data centered on the average of a group, which becomes confusing if the individual data are also centered on the average score of each individual for each value. We therefore use the original scores. Table 5 shows a robustness test when using data on relative values.
An additional contribution of this paper is its comparison of alternative methodologies to measure aggregate human values, in both the three-level models (Table 4) and the two-level models (Table 5) below. Aggregate human values are considered in two ways, as (regional or national) means and as shares of people (%) in the territory who select the options Very much like me or Like me in the questionnaire items about that human value. The indicators are globally centered and standardized by the standard deviation of the observations in each group, 273 regions or 23 countries (Fairbrother et al., 2019). To control for other determinants of individual well-being, we consider the following ten possible relevant attributes: Age, Gender, Marital status, Education, Subjective health, Domicile, Political orientations, Religiosity, Trust in others, and Income (see Sect. 4.3 below for details). 16 Following Sortheix and Lönnqvist (2014), we also control some of our models for a national-level variable, the UNDP's Human Development Index (HDI), as a measure of countries' socioeconomic development. Table 2 presents the descriptive statistics of the variables used in this paper prior to the transformations of variables explained in the notes on the subsequent tables. Table 3 shows the results for the null models, without explanatory variables. Checking the ICC j of the unconditional model in column (1) of Table 3, we see that 13% of the variance in life satisfaction is explained by regional factors. This amount is similar to the 12% in the ICC k of the model in column (3), which considers only country-level random intercepts. Column (4) presents the results for the three-level models. Comparison of the ICC of column (2) and the difference in ICC s for the equations in column (4) enables us to conclude that regional-level factors help to explain only 2.1% of the dispersion of life satisfaction when national-level communalities are controlled for, 17, 18 Europeans' well-being is given mainly by factors distributed at the individual level, but a significant 11.8% of the 16 Although some of these controls may be affected by underlying values, the results for the control variables when human values are excluded from the model (not shown) are very similar. 17 We use the expression fixed effects to indicate dummy variables capturing communalities, particularly geographical areas. Differences among countries might stem from the quality of democracy, freedom of the press, institutions, and/or other factors. See, for example, Veenhoven (2015). We do not use these effects in the following tables in the paper. Controlling contextual differences out through use of dummy variables removes contextual variance, losing important information (Bell and Jones 2015). 18 As mentioned in Sect. 3.2, Israel has no regional division, but the ICC s discussed in the text are roughly the same when Israel is excluded from the sample.  dispersion of life satisfaction is given by national factors and a marginal additional 2.1% by regional factors.

Regional and National Factors Explaining Subjective Well-Being
To explore this result further, Table 4 shows the estimation of three-level models after controlling for compositional effects that might explain different types of people in different regions or countries (see Sect. 4.3). Column (1) shows that country-level factors explain 8.3% of the unexplained part (residuals) of life satisfaction after controlling for these individual characteristics. Regional random effects add 1.6-1.7% additional explanation in all columns of Table 4.
Equations in columns (2) to (6) of Table 4 introduce standardized variables of human values, whose estimates are comparable. Columns (3) and (4) include our two indicators of aggregate human values defined at regional level. The two regional indicators of Self-transcendence seem initially to be statistically significant. However, in columns (5) and (6), The table shows descriptive statistics of the variables before possible transformations done for the equations in later tables. Further details about the individual control variables are provided in Table 7. The six categories of each of the 21 ESS items detailed in Table 1 have been inverted, so higher figures in the resulting variables indicate more importance given to that topic. Despite being categorical, individual human values are considered as continuous variables, as it is common in this literature. For the estimations in this paper, they are group mean centered (by the regional or national averages, depending on the estimated model) and standardized. Aggregated human values are (regional or national) means, or shares of people in that territory selecting high identification. They have been globally centered and standardized, the same than the national Human Development Index this cultural trait becomes significant at national but at not regional level, when nationallevel indicators are included in the model. Country-level between-group differences thus seem to be much more significant than regional differences in explaining cultural influences of human values on life satisfaction. For simplicity, Table 5 shows an analysis of the effects of the ten basic human values in a two-level model including only national random intercepts. Table 4 shows that two dimensions of the four-dimensional scale of individual human values have positive effects on life satisfaction, Self-transcendence and Openness to change. Self-enhancement, in contrast, seems to have negative effects. Conservation is not statistically significant. Our results are consistent with Sortheix and Schwartz's (2017) theory except in this last variable, as explained in Sect. 2.1. Table 5 (columns 1-3) presents more detailed results for the ten-dimensional scale of individual human values. 19 In Table 5, individual Hedonism is the basic human value with the largest (positive) estimate, explaining why Openness to change shows the highest estimate in Table 4. As predicted by Sortheix and Schwartz (2017), Openness to change has positive significant effects in Table 4, although Stimulation is not significant in Table 5. In Table 4, individual Self-transcendence has a positive sign due the positive role of Benevolence in Table 5, since Universalism is not significant in Table 5. This finding is consistent with Sortheix and Schwartz's (2017) hypothesis. Although the second-largest estimate in Table 5 is for Benevolence, the combined negative significant effect of Achievement and Power makes Table 3 Two-and-three-level multilevel null models for life satisfaction of 31,386 individuals ( i ) in 273 regions ( j ) of 23 countries ( k) p values, symbolized through stars, are calculated using the Satterthwaite approximation to the degrees of freedom. Equations in columns (1) and (2) include regional random intercepts. Equation in column (3) includes national random intercepts. Equation in column (4)  Self-enhancement the second-largest higher-order estimate in Table 4, even though Sortheix and Schwartz's (2017) model yields an ambiguous expectation for Achievement. We will discuss the results of individual Conservation values below.

Basic Human Values at Individual and National Level
Columns (2) and (3) of Table 5 show the results when our globally standardized country-level indicators of human values are added to the equation in column (1): national means of the individual human values and national percentages of citizens selecting high Table 4 Three-level multilevel models of life satisfaction for four dimensions of human values *p < 0.05 **p < 0.01 ***p < 0.001 Tables 2 and 7 for the individuals' control variables, omitted here. Individual human values are centered on the regional means and standardized. All the models include random intercepts by regions nested in countries (1) (2) (3) (4) (5) (6)
As mentioned in Sect. 3.2, Table 5 also includes a robustness analysis of the measurement of human values. The estimated model in columns (4) and (5) is the same, but column (4) shows the results when human values are centered on the mean of all human values for each individual before subsequent transformation of variables for our multilevel models. National Tradition becomes non-significant in column (4), but its Satterthwaite's p-value is low (0.058) (not shown in the table). The main difference thus involves national Conformity, not statistically significant in column (4).
Additionally, the last two columns of Table 5 check the robustness of our results on inclusion of the HDI. National Tradition is no longer statistically significant at 5% level, though its p-value is still relatively low, 0.09 in column (7) and 0.11 in column (8) (not shown in Table 5). As discussed in Sect. 2.3, this result is consistent with modernization theory and does not contradict the conclusion that cultures with high presence of traditional values constrain subjective well-being. Moreover, the standardized estimates of national Benevolence are about 40% lower when the HDI is included in the model. This result is also consistent with modernization theory and confirms an independent role of the cultural channel in explaining individual well-being (Inglehart & Baker, 2000). Our preferred models are thus those in columns (5) and (6), which stress the direct link between culture and individual subjective well-being.
Following the discussion at the end of Sect. 3.1, Table 6 shows the results for the models in columns (5) and (6) of Table 5 under weighted estimation, focusing on our cultural indicators. Examination of the t-values shows that weighted estimation produces similar results and reinforces the positive role of aggregate Conformity and Tradition. Table 7 shows the results of our control variables only, for the equation in column (6) of Table 5, as these results are very similar in all estimations. They are consistent with previous literature cited in this paper. Age, for example, displays the U-shaped relationship found Table 6 Robustness analysis of Eqs. (5) and (6) in Table 5 for the three national variables using weighted estimation See the notes in previous tables. Results in Table 5 are estimated using lme4 package (Bates et al., 2015), which does not calculate p values. The p values shown in this table or symbolized with stars in Table 5 have been calculated using other R packages under Satterthwaite's approximation to the degrees of freedom. The columns of weighted estimation have been estimated using WeMix package (Bailey et al., 2021), which produces clustered standard errors as in Stata mixed function, but currently does not produce p values. Level-one weights are the ESS analytical weights, reweighted to sum to the sample size after deletion of missing data. Level-two weights are 1, given that the included countries are not the result of a stratified process of selection of countries. For the weighed estimation, each standardized variable was recalculated using weighted means and standard deviations.

Individual-Level Control Variables
Detail of results in Table 5 Weighted in other studies, meaning that the young and the old tend to report higher well-being. 21 Our selection of indicators and parsimonious way of defining them also make minor contributions to well-being science.

Conclusions
This is the first paper to analyze the joint effects of human values on subjective well-being at individual and contextual levels. We use multilevel modelling to study the direct effects of cultural human values on Europeans' life satisfaction. We study these effects using two levels of geographical aggregation, two methods to aggregate the data from individuals to territories, and two scales of human values. One of our main results is that geographical scale matters when analyzing subjective well-being in European countries. Using the European Social Survey, we conclude that around 12% of the dispersion of life satisfaction is due to national factors, whereas only 2% is due to regional factors. Indeed, regional direct cultural effects are non-significant if national effects are considered in our models.
At individual level, we find that well-being is affected by the following human values, ordered by size of estimate: Hedonism (+), Benevolence (+), Achievement and Power (−), Tradition (−), and Conformity and Self-direction (+). This result is generally consistent with Sortheix and Schwartz's (2017) hypotheses on the effects of human values on wellbeing. The main difference is that we find a positive effect of individual Conformity, even though Conformity is a Conservation value, whose sign should be negative in Sortheix and Schwartz's (2017) theory.
Our main research question addresses the possible existence of happy cultures, or cultural attributes that promote the average individual's subjective well-being. We conclude that happy European cultures are those of countries with high levels of Benevolence and Conformity, and low levels of Tradition. When the Human Development Index is added to the model, Tradition ceases to be statistically significant, and the estimates of Benevolence and Conformity are lower. This finding is consistent with modernization theory. We are not studying cultural change, however, but the direct effects of cultural values on individuals' well-being, which are confirmed by our results.
Our research has several limitations. Firstly, our conclusions about the relative importance of regional and national factors are data dependent. As we analyzed regional effects using a mix of NUTS levels in a European sample, our results must be compared to those of other studies using different data. Secondly, despite the advantages of multilevel modelling for studying the direct effects of individual and aggregate human values, our results say nothing about causality among individual variables or among individual and aggregate variables. We indicate relevant statistical patterns to be considered in future research. Neira et al. (2018) propose a research agenda for multilevel modelling as possible cuttingedge research in the science of well-being. Some of their proposals were to include contextual economic and cultural variables, and to perform different levels of geographical analysis. Our work makes contributions in these directions and can be continued in several ways. We have studied cultural human values, but these values only form part of other possible cultural attributes, such as beliefs or attitudes. Theoretical contributions from economics, sociology, and political science may suggest other contextual variables for consideration. Multilevel models with random slopes and cross-level interactions (mediations) are useful for researching why the same types of people may feel different in different cultures. The cause-effect relationship in happiness equations is not entirely clear (Ballas & Tranmer, 2012) and should be further studied with longitudinal data. Further research is needed on the specific cultural dimension of Conformity: What distinguishing mechanisms are associated with Conformity and Tradition that affect subjective well-being in opposing directions? Can our results be confirmed using other data? Should Schwartz's approach to Conservation be reformulated?
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.