Evolution of family systems and resultant socio-economic structures

Families form the basis of society, and anthropologists have observed and characterised a wide range of family systems. This study developed a multi-level evolutionary model of pre-industrial agricultural societies to simulate the evolution of family systems and determine how each of them adapts to environmental conditions and forms a characteristic socio-economic structure. In the model, competing societies evolve, which themselves comprise multiple evolving families that grow through family labour. Each family has two strategy parameters: the time children leave the parental home and the distribution of inheritance among siblings. The evolution of these parameters demonstrates that four basic family systems emerge; families can become either nuclear or extended, and have either an equal or strongly biased inheritance distribution. Nuclear families in which children leave the parental home upon marriage emerge where land resources are sufficient, whereas extended families in which children staying at the parental home emerge where land resources are limited. Equal inheritance emerges where the amount of wealth required for a family to survive is large, whereas strongly biased inheritance emerges where the required wealth is small. Furthermore, the frequency of polygyny is low in the present model of agricultural societies, whereas it increases for the model of labour-extensive subsistence patterns other than agricultural societies. Analyses on the wealth distribution of families demonstrate a higher level of poverty among people in extended families, and that the accumulation of wealth is accelerated in families with strongly biased inheritance. By comparing wealth distributions in the model with historical data, family systems are associated with characteristic economic structures and then, modern social ideologies. Empirical data analyses using the cross-cultural ethnographic database verify the theoretical relationship between the environmental conditions, family systems, and socio-economic structures discussed in the model. The theoretical studies made possible by this simple constructive model, as presented here, will integrate the understandings of family systems in evolutionary anthropology, demography, and socioeconomic histories.


Introduction
F amilies are basic components of society. Studies of family and kinship lay at the core of anthropology (Fox, 1983;Harrell, 1997;Lévi-Strauss, 1958;Shenk and Mattison, 2011). Families comprise people who are connected by three basic relationships, i.e. husband-wife, parent-child, and inter-sibling relationships (White, 1963), and a variety of rules (or patterns) are observed concerning them (Harrell, 1997;Laslett, 1988;Todd, 1999). Cultural traits pertaining to family relationships are slow to change, because they tend to be inherited vertically from parent to child, and are regulated by social norms (Cavalli-Sforza and Feldman, 1981). Indeed, empirical studies have revealed such slow changes based on population history (Guglielmino et al., 1995;Minocher et al., 2019;Mulder et al., 2001). Family traits have also attracted significant attention as relatively stable basic factors of social characteristics in the study of history (Braudel, 1992a, b, c) and historical demography (Macfarlane, 2002), among others.
The appearance of these traits was previously explained by cultural transmission or adaptation to social and ecological conditions (Goldschmidt and Kunkel, 1971;Laslett, 2015;Todd, 2011). Social studies unveiled a correlation between the period of exposure to the Western Church and the emergence of nuclear families (Schulz et al., 2019). In evolutionary anthropology, relationships between family members are explained based on parental investment theory and intra-family competition for reproductive resources (Ji et al., 2014;Trivers, 1974). Subsistence patterns and other socio-ecological conditions that affect family traits have recently been revealed quantitatively (Colleran, 2014;Gibson and Gurmu, 2011;Macfarlan et al., 2019;Ross et al., 2018). Phylogenetic comparative analyses have been performed to infer the origin and historical change of family traits by controlling the statistical non-independence due to the shared ancestry (Fortunato and Jordan, 2010;Fortunato et al., 2006;Minocher et al., 2019;Mulder et al., 2001). In particular, they revealed that the presence of heritable resources that are typically observed in agricultural/ pastoralist society leads to sibling competition over inheritance (Gibson and Gurmu, 2011) and to patriliny . Differentiation of very rich elite from the majority in agricultural society leads to a lower frequency of polygyny (Ross et al., 2018).
Indeed, correlations have been observed between family traits and several socio-ecological conditions. However, it is unclear whether social factors determine family traits or vice versa (Mace and Jordan, 2011). The reverse effect from family traits to social conditions is also reported in historical demography (Macfarlane, 2002).
As two characteristics in the family system, we consider parent-child and inter-sibling relationships. Although a variety of characteristics can represent these relationships, we focus on residence and inheritance patterns. The residence pattern can refer to nuclear families that comprise a pair of parents and their unmarried children, or to extended families that can involve parents and their married children. The inheritance pattern can refer to either the equal or strongly biased distribution of inheritance among siblings. On this basis, we can suppose four ideal types. (1) Absolute nuclear families, which are nuclear families with unequal inheritance. (2) Egalitarian nuclear families, which are nuclear families with equal inheritance. (3) Stem families, which are extended families with unequal inheritance. (4) Community families, which are extended families with equal inheritance. Indeed, other characteristics of family systems, such as marriage patterns, would be necessary to classify family systems comprehensively, which needs future works. Nevertheless, there is notable variation in the residence and inheritance patterns in pre-industrial agricultural societies in Europe, Northern Africa, and Asia, among others (Harrell, 1997). The influence of agricultural societies on political systems in the modernising era has been described (Rösener, 1993;Wallerstein, 2011). The link between the above four family systems and modern social ideologies has been discussed (Laslett, 1988); Liberalism, liberal egalitarianism, social democracy, and communism are dominant in regions with absolute nuclear, egalitarian nuclear, stem, and community families, respectively (Todd, 1990(Todd, , 1999. However, the discussion regarding the relationship between family systems and ideologies remains largely psychoanalytical. Hence, theoretical studies to unveil the conditions of the evolution of each family system and to connect family systems with socio-economic structures need to be conducted. To consider the interaction of ecological conditions, family systems, and social structures, a constructive approach to reveal the relationships between them is required. In addition, given that family traits are inherited from parents with slight changes over generations, it is appropriate to model their long-term evolution through the accumulation of small variations, as represented by mutations. By modelling the evolution of family systems to adapt to socio-ecological conditions, which in turn form the societylevel economic structure, we aim to integrate the understandings in evolutionary anthropology, demography, and socioeconomic histories. To model the economic consequence of family behaviour, we focus on agricultural society, where family systems determine residence and inheritance patterns in land usage (Todd, 1990(Todd, , 2011. Children may cultivate lands of their own or work together on their parents' land. One heir may inherit the land and most of the property exclusively, or the land and property may be divided equally among family members (Berkner, 1976;Kaser, 2002). Characteristics of the pre-industrial agricultural society include the importance of human labour, land and property in production (Colleran, 2014), a positive correlation of wealth and the number of offspring (Gibson and Gurmu, 2011), and the diminishing returns to labour input (Evenson and Mwabu, 2001;Ricardo, 1891). Hereby, we build the minimal model that is appropriate as long as these conditions are satisfied.
Notably, families constitute society, whereas society, as well as ecology, provides the environment for families. Hence, our model adopts a framework involving multi-level evolution for a hierarchical system. The multi-level evolution was originally introduced to explain the evolution of cooperative behaviour among eusocial insects by examining the conflict between the fitness of an individual and that of a group (Wilson, 1997;Wilson and Wilson, 2007). This framework is generally applied to the evolution of group-level structure in hierarchical systems (Spencer and Redmond, 2001;Takeuchi et al., 2017;Traulsen and Nowak, 2006;Turchin and Gavrilets, 2009). In the previous study, the framework was applied to construct a mathematical model for the evolution of kinship structures in clan societies, which revealed the environmental dependencies of diverse kinship structures Kaneko, 2020, 2021a). The variety of family traits regarding cousin marriage preferences and clan exogamy, as well as descent systems, were investigated in detail therein. Here, in contrast, we mainly focus on parent-child residence patterns and inter-sibling inheritance patterns, and briefly mention the conditions for polygyny.
In this study, we investigate the evolution of family systems and social structures by introducing an agent-based multi-level evolutionary model of pre-industrial agricultural societies. Competition is considered at two levels: that of family, which is an individual agent of the model, and society, which is a group of families. Families produce wealth through family labour and reproduce their population. They possess two strategy parameters concerning the time children leave their parents' home and the distribution of inheritance among them. These parameters are transmitted with slight mutations in each generation. Evolutionary simulations show that four family systems emerge depending on environmental parameters that characterise the land scarcity and perturbations that damage society. Then, the model is extended by adding the marriage process. We show that this extension affects the above result only minimally, whereas it facilitates the discussion of conditions for son-biased investment and polygyny. We then describe the characteristics of social structure in terms of the distribution of wealth in society and relate them to family systems.
Finally, the theoretical results are verified through a data analysis using the standard cross-cultural sample (SCCS), a global ethnographic database of premodern societies (Kirby et al., 2016;Murdock and White, 1969). SCCS contains data from 186 societies, which are thought to be culturally and linguistically independent of each other. In the discussion section, we show the relationships between family systems, socioeconomic structures, and the development of political ideology in the modernising era, by referring to socioeconomic histories.

Model
Overview of the model. The model is described below in general terms in this subsection (see the following subsection for the mathematical formulation). A schematic of the model is shown in Fig. 1a, b. Society consists of families. In each family, individual members live and work together. Children build families of their own by inheriting their parents' wealth. Wealth w is accumulated through production, which in turn increases the level of production and the population, as empirically reported (Gibson and Gurmu, 2011;Macfarlane, 2002). Each society splits into two when the number of families therein doubles its initial value N f , and each family is randomly assigned to one of the two daughter societies. At this time, another society is removed at random so that the number of societies in the entire system remains fixed to N s . This process can be interpreted as invasion, imitation, or the coarse-grained description of a growing system. Therefore, societies that grow faster replace others, resulting in society-level evolution. This multi-level evolution of families and societies follows the hierarchical Moran process (Itao and Kaneko, 2020;Takeuchi et al., 2017;Traulsen and Nowak, 2006).
In each simulation, the environmental parameters are given. Of particular importance are the following two: the capacity c × N f , which represents the amount of available land resources in society (hereafter called land capacity), and ϵ, which is wealth required for a family to survive the perturbations that damage society (hereafter called wealth required for survival). Agricultural production depends on the amount of available land resources, whereas insufficient capacity limits land resources per family. In this model, when the number of families in a society exceeds the capacity, the land resources and the production rate for each family decrease inversely with the number of families at that time. As for the wealth required for survival, ϵ must be paid by a family at the moment of its building. If wealth w is less than ϵ, its members die without reproducing.
Families have population and wealth, as well as the two strategy parameters, i.e. λ, which represents the inequality in the inheritance of wealth among siblings, and s, which represents the probability of children staying at their parents' home to produce together. The mth child's share of the inheritance is proportional to expðÀλmÞ. Therefore, λ = 0 represents an equal division of inheritance, whereas a larger λ represents the eldest child inheriting more. In some societies, the youngest child inherits the most instead of the eldest (Todd, 2011). If necessary, the order of children could be arranged in reverse to include such a case. Each family forms an extended family with probability s or otherwise forms a nuclear family (explained in detail below). These parameters may be determined by intra-family competition among parents and siblings (Ji et al., 2014;Trivers, 1974). Here, we do not model the competition explicitly but do so implicitly by tracing the evolution of λ and s.
The life cycle of the families depends on the strategy s. With the probability of 1−s, siblings are separated before agricultural production, i.e. they inherit some property determined by λ, lose ϵ of wealth, and build families of their own to produce independently. In contrast, with a probability of s, siblings remain in their parents' family to produce together, after which they are separated. According to the law of diminishing returns, productivity increases sub-linearly with labour force input (Bacci, 2017;Malthus, 1798;Ricardo, 1891). Following a study on pre-industrial farming (Evenson and Mwabu, 2001), we assume that production increases in proportion to the logarithm of labour input, and is perturbed by Gaussian noise with a mean of 0 and a variance of σ 2 , resulting from internal and environmental fluctuations.
Then, as long as the law of diminishing returns is satisfied, the total output of siblings is always higher when each sibling produces independently to form a nuclear family than when they concentrate their labour in an extended family. In this model, the productivities of N members are Nlog 2 and log ðN þ 1Þ for nuclear and extended families, respectively. However, under a limited capacity of available land resources, the total output of the society consisting of nuclear families will be lower than that of extended families, as a result of inefficient land usage (the productivities are 1 N Nlog 2 and log ðN þ 1Þ for nuclear and extended families, respectively). In other words, there is a conflict between family-and society-level preferences for nuclear versus extended families under the conditions of Societies (green) comprise families (blue) in which family members (black) live together. The grey frame represents a single generation. When the population of a society reaches twice the initial population, a society splits and another society is removed from the system at random to keep the number of societies fixed-i.e. we adopt the hierarchical Moran process.
b The life cycle of families. The black squares, circles, and triangles represent members of different generations. The size of the blue circle reflects the amount of wealth that each family possesses. The dashed arrows show the separation of siblings, and the arrows show the production of wealth through family labour. Each family has a strategy parameter s.
With the probability of 1−s, siblings are separated before production, and then they produce independently and reproduce the next generation, as in a nuclear family. On the other hand, with the probability of s, siblings produce together, after which they build their own home and reproduce the next generation, as in an extended family. When we extend the model to consider marriage, its process is added just before production.
limited capacity. As long as one considers pre-industrial farming, it is expected that family members work together on their land, and the law of diminishing returns is satisfied. However, this formulation will be inappropriate for modern farming or other subsistence patterns, which is beyond the scope of this model. The results for different formulations of labour-extensive subsistence patterns are briefly discussed. After sibling separation and the production of wealth, each family reproduces. The number of children in families is positively correlated with their wealth in pre-industrial societies (Gibson and Gurmu, 2011). Here, we assume that it follows the Poisson distribution with a mean of b + fw, where b and f represent the minimal birth rate and the increment of birth rate by wealth w, respectively. Children culturally inherit s and λ from their parents, with a slight variation through 'mutation' at the rate of μ, according to previous studies on cultural evolution (Cavalli-Sforza and Feldman, 1981;Creanza et al., 2017). At the time of altering generations, families lose dw of wealth, where d represents the decay rate of wealth due to ageing equipment, disaster, or taxation, for example. Additionally, each society splits if the number of families reaches 2N f at that time. The parameters are summarised in Table 1. The above is a minimal model to discuss the diversification of family systems concerning parent-child and inter-sibling relationships. To consider husband-wife relationships, we extend the model by assigning each family male and female populations, and the strategy for inheritance distribution for sons and daughters. In the extended model, one can have multiple spouses by paying sufficient bridewealth. Polygyny increases both production and reproduction. This model is explained in the Supplementary Text in detail.
Algorithm of the model. In this subsection, we show the mathematical formulation of our model. We adopted the following algorithm for changes in the wealth and population of families. For the parent family i and its jth child's family i, j, the population N and the amount of wealth w at time t are expressed as follows: With probability 1 À s t i , otherwise, Then, where In each simulation step (generation), families lose a proportion d of their wealth (Eq. (1)). With probability 1−s, siblings leave their parents' home (Eq. (2)), distribute the inheritance with a wealth loss of ϵ (Eq. (3)), and produce independently (Eq. (4)). Otherwise (with probability s), siblings produce together at their parents' home (Eq. (6)) and build their own families after production (Eq. (5)) instead of Eqs.
(2)-(4). Here, the production of wealth is inversely proportional to the number of families in society if the capacity is exceeded (Eq. (9)), is proportional to the logarithm of labour and increases with linear feedback from wealth, and is perturbed by noise resulting from internal and environmental fluctuations following a normal distribution (Eq. (10)). Finally, families produce offspring (Eq. (7)) and strategy parameters are transmitted with slight mutation (Eq. (8)). The birth rate increases linearly with wealth. Here, families reproduce without considering marriage explicitly. The extended model with the marriage process is explained in Supplementary Text.
In the simulation, the initial values of strategies are s = 0.5 and λ ¼ log 2 for all families. Hence, at the initial state, families can form nuclear or extended families with equal probability, and the inheritance is moderately biased so that the share of the inheritance received by subsequent children is half of that received by the preceding children. In other words, the families are not differentiated as nuclear or extended, or as equal or strongly biased inheritance providers. However, no qualitative changes are observed under other initial conditions. The source code has been made publicly available in the Dataverse repository (Itao and Kaneko, 2021b), https://doi.org/10.7910/DVN/3ZGCQI.

Results
Evolution of family systems. The simulations are performed for 2000 time steps. The time series of the evolution of family strategy parameters are shown in Fig. 2. Family strategies do not diverge within each society, but are concentrated around a specific value adapted to each given environmental condition. The probability to form extended families s (plotted in red) and the inequality in inheritance λ (plotted in blue) evolve depending on the values of land capacity c and wealth required for survival ϵ, respectively. The evolution of the strategy converges within~1000 steps in every parameter region.
We conducted multi-level evolutionary simulation 100 times for each condition and averaged the strategy parameters of families in the last 1000 steps. Figure 3a, b show the dependence of s and λ on c and ϵ. Increasing land capacity c makes nuclear family strategies more preferable at both the family-and society-levels. Then, the probability of forming an extended family s decreases, implying the evolution of nuclear families. Increasing ϵ, the wealth required for survival, increases the demand for inheritance by younger siblings. This results in a smaller inequality in inheritance λ, indicating the evolution of the equal inheritance. Figure 3c shows the phase diagram of family systems. Here, we classify the family systems as extended if s ≥ 0.5, and as nuclear if s < 0.5. Similarly, we classify them as unequal inheritance sharers if λ ≥ log 2 and as equal if λ<log 2. This criterion is determined according to whether the share of the inheritance awarded to the subsequent children is smaller (or larger) than half of that awarded to the preceding children. Therefore, the four basic family systems evolve depending on the two environmental parameters c and ϵ. A stem family (plotted in green) evolves if c and ϵ are small. A community family (plotted in orange) evolves if c is small and ϵ is large. An absolute nuclear family (plotted in blue) evolves if c is large and ϵ is small, and an egalitarian nuclear family (plotted in yellow) evolves if both c and ϵ are large.
We show the dependence of the phase diagrams of family systems upon other parameters in Figs. 4 and S1. The phase diagrams plotted against c and ϵ are qualitatively robust and independent of the other parameter values. However, quantitative trends exist. Generally, N s and N f determine the strength of selection pressures at the family-and society-levels (Takeuchi et al., 2017;Traulsen and Nowak, 2006). For a large N f or small N s (ultimately, if N s = 1), family-level competition becomes dominant rather than society-level competition, which leads to the evolution of selfish behaviour. Conversely, when society-level competition is dominant because of small N f or large N s , cooperative behaviour evolves. Figure 4a, b show that, if N f is large or N s is small, nuclear families evolve even when c is small. Recall that the total production of siblings increases if they work independently, but that of society decreases because of the inefficiency in land usage if the capacity is limited. Therefore, choosing a nuclear family under small c is a selfish strategy that evolves for small N s and large N f . As wealth w accumulates, ϵ would become relatively small for the wealth, and accumulation would be accelerated. However, when the minimal birth rate b is high, the number of offspring increases, and each share decreases. Therefore, less wealth is accumulated, and equal distribution evolves in larger parameter regions because of the relatively large ϵ, as shown in Fig. 4c. The dependence of the phase diagram on the parameters mutation rate μ, decay of wealth d, and increment of birth rate by wealth f are shown in Fig. S1.
Evolution of husband-wife relationships in the extended model. The simulation results are shown in Fig. S2 for the extended model considering the marriage process. Parental investment is biased for sons almost four times as much as daughters for most of the parameter regions. Investment for sons is advantageous because wealthy men can have many wives and increase their fitness. Evolutionary anthropologists have reported that daughter-biased investment evolves under paternity uncertainty , which is beyond the scope of our model. Hence, it is reasonable that only son-biased investment evolves in our model.
Results also show that the frequency of polygyny is <20%. The bias of parental investment and the frequency of polygyny are almost independent of land capacity c or wealth required for survival ϵ. Furthermore, consideration of the marriage process only minimally affects the results for parent-child and intersibling relationships. Hence, we will analyse the economic structures of evolved societies by using the previous minimal model in the following section.
Additionally, we studied a model in which the diminishing returns of family labour in production are relaxed, to consider labour-extensive subsistence patterns other than agriculture. Figure S3 shows that if production increases linearly to the number of family labourers, the frequency of polygyny increases to more than 20% almost independently of c and ϵ, even though the parental investment bias is almost the same as the above model. In this model, the increase in polygyny results from a larger fraction of relatively wealthy people. This scenario is consistent with empirical reports for foraging, horticultural and agropastoral societies (Ross et al., 2018). Extended families are dominant even when land capacity c is sufficient because they are no less preferable than nuclear families even at the family-level in this model. However, both nuclear and extended families are observed in most of the subsistence patterns (Murdock and White, 1969). It suggests that nuclear families can evolve owing to some reasons not covered by our model. To discuss the variation of family systems depending on subsistence patterns, it will be necessary to consider the difference in productivity and lifestyle.
Wealth distribution and evolution of social structure. Subsequently, we investigated the wealth distribution of families for each society after evolution. Note that data from the wealth distribution in modern society suggest an exponential-type tail for the rich side (Chakrabarti et al., 2013; Tao et al., 2019) (say a log-normal (Gibrat, 1931) or gamma distribution (Chakraborti and Patriarca, 2008)), and a power distribution for the poor side (Reed, 2003) (say a gamma distribution (Chakraborti and Patriarca, 2008)). The gamma distribution is obtained by assuming the wealth growth with positive feedback and nonlinear saturation, as well as a multiplicative stochastic process, which are included in our model (see Supplementary Text). Figure 5a shows the frequency distribution of the wealth of families within each society. In every parameter region, the distribution approximately follows a power-law on the poor side and has an exponential tail on the rich side, which is consistent with the above data. Because inherited wealth depends on birth order, we also plotted the distributions of wealth by distinguishing the eldest siblings from the others in Fig. 5b. With decreasing wealth required for survival ϵ, and consequently increasing inheritance inequality, the distributions of the wealth of siblings separate further. As a result, the accumulation of wealth by heirs is accelerated. Decreasing land capacity c and the evolution of extended families result in poorer younger siblings, whereas greater land capacity c and the evolution of nuclear families give rise to wealthier younger siblings.
Although the wealth distribution follows the power-law w α on the poor side and the exponential distribution expðÀβwÞ on the rich side universally, the heaviness of tails depends on the environmental parameters. We fitted values for the lightness of the poor tail α and those for the rich tail β and averaged them over 100 trials for each environmental parameter c and ϵ in Fig. 6a, b, respectively. Smaller c results in smaller α, i.e. the heavier tail is on the poor side, while smaller ϵ results in smaller β, i.e. the heavier tail is on the rich side. These results suggest the characteristics of the wealth distribution in the four corresponding family systems. However, it remains unclear whether they result from environmental conditions or family systems. To confirm the relevance of family systems, we computed the dependence of wealth distributions on family systems by sampling each family system using fixed environmental parameters near the boundary of the four phases of family systems, where the values of s and λ are distributed to cover all four family systems.
The average values of α and β for each family system, which were sampled for the fixed environmental parameters (plotted in black), are shown in Fig. 6c, d. They demonstrate the trend that the poor tail is heavier for extended families with poor younger siblings, whereas the rich tail is heavier for unequal inheritance with rich heirs. By comparing these results with the results  Table 1. ARTICLE HUMANITIES AND SOCIAL SCIENCES COMMUNICATIONS | https://doi.org/10.1057/s41599-021-00919-2 averaged over several environmental parameters around the phase boundary (plotted in red), it is shown that the above trend depends on each family system, and is further intensified by environmental parameter values.
Empirical data analyses. Next, we verify our results on the relationship between environmental conditions, family systems, and economic structures. Using the global ethnographic database of 186 premodern societies, called SCCS (Kirby et al., 2016;Murdock and White, 1969), empirical data analyses were conducted.
First, we classified family systems of pre-industrial agricultural societies. We then identified pre-industrial agricultural societies by using the Subsistence economy: dominant activity variables (5-7 correspond to agriculture). Then, we identified the family systems by using Domestic organisation (1-5 correspond to nuclear families and 6-8 correspond to extended families) and Inheritance distribution for movable property (1 corresponds to equal and 2-4 correspond to strongly biased inheritance) (see Supplementary Tables for a detailed explanation of these variables). Out of 186 societies in SCCS, 91 societies conducted agriculture and inheritance of movable properties. Among them, 14 societies were  classified as stem families, 30 as community families, 17 as absolute nuclear families, and 30 as egalitarian nuclear families. Figure S4 shows their geographic distribution. Here, we used the data on the inheritance of movable properties to identify intersibling relationships. However, similar trends on the following variables were achieved, even when we used those pertaining to the inheritance of real property, as shown in Table S2. Next, we conducted Spearman's rank correlation analyses and calculated the correlation between SCCS variables and parent-child (nuclear or extended) or inter-sibling (strongly biased or equal) relationships. The database contains various variables of socio-ecological factors. By calculating the correlation for each variable and listing the variables in descending order in the absolute value of the correlation, we found those variables related to the parameters in our model in the top of the list. The variables that are highly correlated with parent-child and intersibling relationships are listed in Table S1 and Table S2, respectively. Among them, we show the variables that can be related to the model parameters in Table 2. Table 2 shows the dependence of family systems on environmental conditions. Communality of land and Land Shortage, suggesting larger and smaller land capacity c, respectively, are correlated with extended families (Corr. −0.31 (P = 0.03) and Corr. 0.26 (P = 0.09), respectively). This is consistent with the theoretical results showing the evolution of extended families for smaller c. On the other hand, Frequency of internal warfare and Acceptability of violence within society suggest that violent conflict is more frequently observed in societies with equal inheritance (Corr. 0.36 (P = 0.08) and Corr. 0.27 (P = 0.13), respectively). Such violence will damage goods and require families to have more wealth to survive; as a result, the wealth required for survival ϵ increases in our model. Accordingly, equal inheritance is more frequent for larger ϵ, as is consistent with our results.
Furthermore, the data suggest the correlation between family systems and economic structures. Number of poor implying smaller α is positively correlated with extended families (Corr. 0.29 (P = 0.06)), whereas Number of rich people implying smaller β is negatively correlated with equal inheritance (Corr. −0.37 (P = 0.01)). Thus, the empirical data are consistent with our simulation results, concerning the relationship between environmental conditions, family systems, and society-level economic structures. See Supplementary Tables for the explanation on values of SCCS variables.

Discussion
By simulating the multi-level evolution model of family systems, we demonstrated the evolution of family systems depending on the environmental parameters for the capacity of available land resources c and the amount of wealth required for a family to survive ϵ. As for parent-child relationships, if there is sufficient land capacity, nuclear families evolve, whereas extended families evolve under a land shortage. As for inter-sibling relationships, if the wealth required for survival is large, equal inheritance evolves, whereas strongly biased inheritance evolves when that is small. Therefore, the four basic family systems characterised by both relationships above are represented as 'phases' depending on c and ϵ. By considering marriage, we then confirmed son-biased investment and infrequent polygyny. Additionally, we clarified the characteristics of wealth distribution determined by the dominant family systems within societies. The tail of the poor side is heavier (that is, many poor people) for extended families, and that of the rich side is heavier (that is, many rich people) for families with unequal inheritance. Empirical data analyses of premodern societies in SCCS supported our results. Now, we refer to demographics and socioeconomic histories in the premodern and modernising era, especially in Western Europe and East Asia. The land capacity c in our model can be measured approximately by the period since the onset of agriculture. In the areas where agriculture started early, population growth resulted in the exhaustion of available land, and labourintensive farming developed, as observed in China (Pomeranz, 2000;Wallerstein, 2011), Russia (Hizen, 1994), and Japan (Hayami, 2015). In Western Europe, especially Holland, the Paris Basin, Southern England, and Central Spain, the capacity was large until industrialisation, and labour-saving farming was developed (Pomeranz, 2000) as a result of the following reasons: agricultural progress in medieval times enabled virgin land cultivation by gathering the children not inheriting the parental lands (Bacci, 2017;Cameron et al., 1993;Grigg, 1980;Pirenne, 1956); the population stagnated in premodern times because of religious wars and plagues (Bacci, 2017); and colonies were established early on (Wallerstein, 2011). Accordingly, the model result concerning c implies that extended families evolve in the areas where agriculture started early. Table 2 also suggests that the exhaustion of land leads to the evolution of extended families.
As for wealth required for survival ϵ, demographics report that the frequency of violent conflict decreased in the following order in Eurasia: the centre of the continent, peripheral, and island regions (Khazanov and Wink, 2012;Macfarlane, 2002;Umesao, 2003). The regions close to the pole of civilisation and/or those frequently attacked by foreign people would have a large ϵ. Hence, the model result concerning ϵ implies that equal inheritance is dominant in the centre of the Eurasia continent and other regions that are vulnerable to warfare (see Fig. S4). The results of the empirical data analysis in Table 2 also support the correlation between such violent conflict with the evolution of equal inheritance.
From geohistorical reports discussed above, the family systems in each region can be explained: absolute nuclear families (nuclear family, unequal inheritance) in England and the Netherlands, where available land resources were sufficient and wealth required for survival was small; egalitarian nuclear families (nuclear family, equal inheritance) in France, Spain, and Italy, where both land capacity and necessity of wealth were large; stem families (extended family, unequal inheritance) in Japan, Germany and many parts of rural Western Europe, where both of them were small; and community families (extended family, equal inheritance) in China, Russia, and Northern India, where land capacity was small and necessity of wealth was large (Berkner, 1972;Todd, 2011).
Apart from these environmental conditions, the number of families within a society N f , the number of competing societies N s , and birth rate b are also relevant parameters for determining the family system. N f is large for large-scale land management as seen in England, the Netherlands, France, and Spain, whereas N f is small and N s is large in family farm management as observed in China, Russia, Japan, and Germany (Bacci, 2017;Cameron et al., 1993;Hizen, 1994;Pomeranz, 2000;Wallerstein, 2011). The trends in Fig. 4a, b are consistent with the observation of nuclear families in the former regions and extended families in the latter. The birth rates were low in Japan and Western Europe, especially in England (Macfarlane, 2002), and higher in Russia (Hizen, 1994). The observation of unequal inheritance in the former regions and equal inheritance in the latter demonstrates a similar tendency to Fig. 4c. Next, we examine the validity of our results regarding the relationships between family systems, socio-economic structures, and modern social ideologies. Figure 6 suggests that, in England and the Netherlands involving absolute nuclear families, the tail of wealth distribution is heavy on the rich side and light on the poor side. Indeed, wealthy farmers prospered and employed a majority as labourers who had better living standards than those in poorer regions (Laslett, 2015;Macfarlane, 2002;Shaw-Taylor, 2012;Tawney, 1912;Todd, 1990;Wallerstein, 2011). The accumulation of capital and independent labour forces explains the development of individual liberty and capitalism in England (Braudel, 1992a, b, c;Todd, 1990;Wallerstein, 2011). Wealth distribution in France and Spain, involving egalitarian nuclear families, was suggested to have light tails on both the rich and poor sides. That is, agricultural societies were less differentiated and weakly stratified (Dupeux, 1972;Rösener, 1993;Wallerstein, 2011), which forms the basis of the values of freedom and equality. Our results suggest that wealth distribution in Germany, Sweden, and Japan, involving stem families, had heavy tails on both rich and poor sides. Wealthy farmers prospered by exploiting others and the stratification of society advanced in accordance with the order and class distinctions (Hayami, 2015;Hayami and Kurosu, 2001;Kastner, 1978;Mager, 1981;Rösener, 1993;Todd, 1990). Wealth distribution in Russia and China, involving community families, was suggested to have a light rich tail and a heavy poor tail. Indeed, the middle class was significantly sparse, and people were uniformly poor (Rösener, 1993), which led to the adoption of communism (Thaxton, 1997;Weber, 1995). In this manner, the wealth distribution obtained in our model connects family systems with society-level characteristics observed in socio-economic history. A study of political ideology showed that people supported authoritarianism in the presence of many individuals being exposed to threats, and egalitarianism in the absence of strong inequality or power imbalance (Claessens et al., 2020). Our results are consistent with this, because the heavier poor and rich tails imply the presence of vulnerable and privileged people, respectively.
Note that our results regarding the family systems and the socio-economic structures are expected to be rather general. The conclusion here is independent of the details of the present model, as long as the production increases sub-linearly with labour input, and multi-level selection of families and societies is considered.
One can also discuss long-term changes in land capacity c and wealth required for survival ϵ, and their influence on family systems. At times of cultivation, a nuclear family evolves because of sufficient capacity. As the population increases and capacity becomes limited, an extended family would replace it. Additionally, because of the dense population, the risks of invasion from surrounding areas and conflict within societies increase, and accordingly, owing to the loss of wealth by violence, ϵ would increase gradually. This scenario explains the historical change of family systems from a nuclear family to a stem family, and then to a community family (Todd, 2011).
Environmental factors change gradually owing to the interaction between society and the environment. However, the change in environmental factors, in turn, will alter family systems and social structures. Such historical dynamics have been discussed as the interaction of factors on different time scales (Braudel, 1992a, b, c). To discuss such interaction of factors at different levels, the present constructive model will give a basic explanation.
The present model has some limitations. First, the differentiation of people between the elite and the majority was not considered. As society becomes stratified, people start to rent land from the elite. This results in the divergence of environmental factors and family systems between them (Todd, 2011). A model is needed for handling social stratification and the interaction of classes to discuss broader issues. Second, we did not model intrafamily competition for resources explicitly and the relationships between family systems and such competition remain unsolved. Intra-family competition has attracted attention in evolutionary anthropology (Ji et al., 2014;Trivers, 1974). In fact, the previous studies have revealed the high status of the elderly in extended families (Lee and Kezis, 1979), and sibling competition over the land inheritance (Gibson and Gurmu, 2011). The model, then, needs to have three levels, i.e. individual, family, and society. Finally, the current model focuses only on pre-industrial agricultural societies. Models focusing on other subsistence patterns are needed to discuss the diversity of family systems widely. For example, Fig. S3 shows that if the diminishing returns of family labour are relaxed, the frequency of polygyny increases. Furthermore, in the modern world, agricultural societies should no longer be regarded as isolated systems, but constitute components of a world-system (Wallerstein, 2011). A new model needs to be developed to consider the interaction between towns and agricultural societies, as well as international, political, and commercial networks.
Furthermore, the present empirical data analyses have some limitations. First, we could only analyse the correlation between cultural variables and family systems. Although it is desirable to conduct better analyses (such as classification learning) to reveal features relevant to family systems, it was infeasible due to data insufficiency. In this study, we used SCCS in which the data for a variety of socio-ecological factors are available. SCCS enabled us to test the correlation between model parameters and family systems empirically, although the sample size was limited. Future works should make use of other databases that include more societies (but fewer variables per society), such as the Ethnographic Atlas (Murdock, 1967). Second, we could not analyse causal relationships owing to a lack of chronological data. Our correlation analyses are insufficient to examine whether history progresses as our model predicted.
Here, the collaboration of field studies, historical analyses, and theoretical modelling is necessary to further elucidate the historical dynamics of societies. Field studies describe the individualor family-level behaviour and society-level structures synchronically. Historical or phylogenetic analyses unveil the diachronic change of such factors. In addition, Leach has emphasised the importance of generalising ethnographic findings by using mathematical formulation to unveil universal structural patterns that may appear in any type of society (Leach, 1961). The constructive model, as presented here, provides a simplified mathematical expression of family behaviours, gives a general framework that allows comparison of various societies, explains the universal patterns between family-and society-level factors, and reveals their historical evolution.
To summarise, we presented a multi-level evolution model to account for the emergence of the four basic family systems and the resultant socio-economic structures depending on environmental conditions, as is consistent with family-level anthropological studies and society-level economic histories. Here, the microscopic characteristics of families determine the macroscopic economic structures, which forms the basis for the development of societies. This study allows an explanation of the universal evolutionary constraint that human societies satisfy.