The Influence of Gender, Age, Matriline and Hierarchical Rank on Individual Social Position, Role and Interactional Patterns in Macaca sylvanus at ‘La Forêt des Singes’: A Multilevel Social Network Approach

A society is a complex system composed of individuals that can be characterized by their own attributes that influence their behaviors. In this study, a specific analytical protocol based on social network analysis was adopted to investigate the influence of four attributes (gender, age, matriline, and hierarchical rank) on affiliative (allogrooming) and agonistic networks in a non-human primate species, Macaca sylvanus, at the park La Forêt des Singes in France. The results show significant differences with respect to the position (i.e., centric, peripheral) and role (i.e., implication in the network cohesiveness) of an individual within a social network and hence interactional patterns. Females are more central, more active, and have a denser ego network in the affiliative social network tan males; thus, they contribute in a greater way to the cohesive structure of the network. High-ranking individuals are likely to receive fewer agonistic behaviors than low-ranking individuals, and high-ranking females receive more allogrooming. I also observe homophily for affiliative interactions regarding all attributes and homophily for agonistic interactions regarding gender and age. Revealing the positions, the roles, and the interactional behavioral patterns of individuals can help understand the mechanisms that shape the overall structure of a social network.


INTRODUCTION
Animal societies are complex systems in which individuals have non-random and complex interactions, and are likely to develop behavioral strategies (Dunbar, 1989). This leads to the formation of a multilayered and multi-behavioral structure. However, questions persist about the fundamental evolutionary process by which a society emerges, stabilizes, and adapts.
Age, or more precisely, ontogenesis (i.e., the development of an organism), influences the evolution and development of social relations and species-specific behaviors that are largely affected by interactional experiences with congeners (Harlow and Suomi, 1974;Hinde, 1974;Wilson, 1980;Shimada and Sueur, 2014). For example, hierarchical rank acquisition appears to be closely related to age (Borries et al., 1991) and the early experiences of juveniles (Mitchell et al., 1967;Olds et al., 1997). However, this influence differs according to species and gender (Sosa, 2015). Additionally, older individuals and females in particular are more likely to experience social exclusion (i.e., decrease in social interactions) (Hauser and Tyrrell, 1984).
One major kinship phenomenon among the animal kingdom is the matrilineal rank inheritance (MRI) (Kawamura, 1958) observed in macaques. It consists of the transmission of hierarchical rank from mother to daughter; the latter acquires the hierarchical rank directly below that of her mother. In addition, as according to the youngest ascendancy rule, young females outrank their older sisters (Thierry et al., 2004). The MRI process is made possible by nepotism, in that related females support each other during conflicts against non-kin females and help juvenile females outrank their older sisters (Cheney, 1977;Datta, 1983;Chapais and Gauthier, 1993). Furthermore, an adult female can outrank her mother when she is old and subsequently lacks kin support and has limited physical ability (Chapais and Berman, 2004). See Chapais and Berman (2004) and Hepper (2005) for an overview.
Social network analysis (SNA) is one approach used to analyze systems  as complex as animal societies. SNA was first applied in psychological studies and, for a few decades, in animal social research (see Prell, 2011 andBrent et al., 2011 for an overview of SNA epistemology). However, certain methodological precautions must be taken when using any of the various analytical techniques based on SNA (Wasserman and Faust, 1994;Krause et al., 2009;Brent et al., 2011). In this study, I describe an analytical protocol based on SNA tools that compensates for the intrinsic limitations of animal behavioral data (i.e., dependency of data) and allows the analysis of weighted networks (network with weighted links).
Several studies have used SNA tools to examine the position and role of group members in non-human primates and other animal species. Lusseau and Newman (2004) revealed that central individuals are key players in maintaining social cohesion and have greater knowledge of their environment. In some non-human primate societies, central individuals are highranking animals (Kanngiesser et al., 2011). Using an interspecific comparative approach, several studies have analyzed network metric variations and succeeded in linking them to variability in social structure and dominance style (Sade, 1972;Voelkl and Noe, 2008;Sueur et al., 2011). Previous studies have also found that individuals from the philopatric gender are more central within a network (Smuts, 1985;Matsuda et al., 2012). In this way, central individuals play an important role in group cohesion and their position depends on several individual characteristics. Thus, identifying these central individuals according to their attributes could allow us to better understand how a social structure is shaped.
SNA research also addresses the principles of homophily and heterophily that refer to preferential interactions between similar (homophily) or dissimilar (heterophily) individuals (Lazarsfeld and Merton, 1954). These phenomena have been observed in many animal species: cetaceans (Lusseau and Newman, 2004), fishes (Croft et al., 2005), marmots (Wey and Blumstein, 2010), and human (McPherson et al., 2001) and non-human primates (Silk, 2001;Cords, 2002;Carter et al., 2015). However, animal research has generally disclosed the existence of homophily for one behavior as related to a single attribute. In this study, I examine the existence and level of homophily as related to a variety of behaviors and attributes. Moreover, revealing such a phenomenon may help us understand how individuals build their networks depending on the attributes of other individuals.
Macaca (Macaca sp.) societies are characterized by their common social organization, but they are also known for their different social styles. Extensive research has shown that dominance hierarchies vary greatly in the macaque genus (i.e., dominance styles) (De Waal and Luttrell, 1989;Thierry et al., 2000;Sueur et al., 2011). Furthermore, the hierarchical structure of females in the Macaca taxon is a well-studied phenomenon that appears to be entirely dependent on the MRI (Thierry et al., 2000). In contrast, each Macaca species has stable multimale, multi-female, and multi-generational social groups in which females are philopatric and males migrate. These common characteristics allow the elucidation of the influence of individual attributes on the interactions between individuals and represents an excellent biological model for this study.
In this study, I use SNA tools to determine individual positions and interactional patterns according to four specific attributes (age, gender, matriline, and hierarchical rank) in affiliative and agonistic networks in M. sylvanus. Based on previous studies, several assumptions can be made in response to the following questions: (1) Who are the most central individuals? As in many cercopithecines, M. sylvanus females are the philopatric gender, which should increase their ability to form denser, stronger, and more perennial networks than males (Smuts, 1985). Thus, they can be expected to be the most active and central individuals in the affiliative network. Exploring such functions could reveal the significance of their role in facilitating group cohesion. Males are generally the more aggressive individuals (Gray, 1971), and therefore should be particularly active and central in the agonistic network.
(2) How age and gender influence the positions and roles of individuals? According to our extensive knowledge of the MRI process, we can expect to observe age-related behavioral variations in females that are highly correlated with their reproductive status (Borries et al., 1991;Chapais, 2004) and matriline. The social activity of young females would therefore be more intense (affiliatively and agonistically), with a decrease in activity during their latter ontogenesis, which in some cases may lead to social exclusion at an advanced age. In males, a minimum hierarchical level for older individuals may exist that enables them to maintain a certain ranking (Sosa, 2015). Such kinetics among males reduces their chances of experiencing social exclusion and thus they may face only a minor decrease in social activity, position, and role. (3) Do common interactional patterns exist among individuals according to their attributes? One sociological model predicts attractiveness to high-ranking females (Seyfarth, 1977). However, this model appears subject to variability. It is mainly observed in despotic societies, and attractiveness to low-ranking individuals has been reported in other species (Schino, 2001;Sueur and Petit, 2008). According to these findings, Seyfarth's model should not apply to M. sylvanus. I also expect to observe homophily related patterns such as allogrooming that target same-gender and same-age individuals and kin (Hirsch et al., 2012). As for agonistic behaviors, it is difficult to form any hypothesis, but heterophily can be suspected.
Responding to these questions allows to reveal how individual attributes and social structure (gender philopatry) produce behavioral divergences that lead to different positions and roles within the group (question 1), how these deviations evolve with the ontogenesis of an individual (question 2), and by which mechanisms individuals interact among themselves (question 3). This multilevel approach allows for a better understanding of how these different levels shape the overall structure of the society in M. sylvanus.

Study Site and Subjects
The current study was conducted over a period of 4 months (July to October 2011) in the park La Forêt des singes, in Rocamadour, France. The 141 M. sylvanus individuals in the park are divided into three groups and live in semi-free ranging conditions (Sugiyama, 2015) in a 20 hectare forest. They are fed in foraging areas twice per day and have water ad libitum. For more details on the management of the park, refer to de Turckheim and Merz (1984). The demographic data (gender, age, and matriline) were provided by the scientific director of the park, Ms. Ellen Merz. The study focused on one of the three groups. Four newborns were excluded from the observations (three males and one female), so that the number of individuals observed was N = 52. The group had a balanced gender ratio of 25 females and 27 males, with an age range between 1 and 25 years old. The individuals were previously identified during 1 month through their tattoos. The observations were conducted with the approval of the park management, an agreement that was subject to the specific condition that I would not directly contact nor handle individuals. As I performed simple observations without any type of intervention, I did not require authorization from the French National Advisory Ethics Committee.

Behavioral Observations
Observations were conducted by repeated focal samplings of 30 min per individual. Each individual was observed approximately 30 times (15 ± 2 h), for 786 observation hours of 52 individuals. Focal sampling time was determined after 2 months of pre-observation. To trade with bias of observation in time of day and feeding time, individuals were observed randomly from 8 am to 5 pm over the 4 months. During the observations, I registered allogrooming and agonistic behaviors (threatening face or growl, charge, avoidance, attack, chase, and aggressive slap, grab, or bite). A complete description of the ethogram of M. sylvanus can be found in Hesler and Fischer (2007). An iPad 1 tablet (Inc, 1976) computer and the WhatISee2.0 application (Inc, 2009) were used to register the individuals involved, and the direction, frequency, and duration of the behaviors. Directed and weighted agonistic and allogrooming matrices were built using the obtained behavioral frequencies (Figure 1). The overall observation yields a total of 5867 agonistic interactions and 1281 grooming interactions. The agonistic matrices allow us to calculate the hierarchical rank of each individual using David's Score (David, 1987) with R 3.0.1 (Ihaka and Gentleman, 1996) package steepness (de Vries et al., 2006).

Building Matriline Categories
Kinship bonds among individuals were determined using two methods. First, data were provided by the park officials who, along with scientists, have been monitoring the population in the park. Second, matrilines were determined through genetic analyses of mitochondrial DNA using eight microsatellite markers. The collection and analysis of DNA samples were performed by the park authorities. The poor quality of DNA samples made some DNA results uncertain. For this reason, matriline groups were built only with individuals whose relatedness was confirmed based on direct observations and genetic analyses. To conserve only close kinship relationships, only the individuals with the same mother were considered related for each mitochondrial haplotype (Figure 2). Thus, individuals whose matrilines were uncertain did not belong to any matriline group (eight males and one female). In addition, matriline results must be carefully considered, as not all individuals were taken into account owing to a lack of information on their kinship bonds.

Data Consideration before Analysis
Collecting data from all members of the same social group led us to the construction of two social networks (agonistic and affiliative) through the existence of multiple interactions. The intrinsic nature of the collected data (interactions between samegroup individuals) underlies the non-independence of the data required by inferential statistical techniques (Wasserman and Faust, 1994;Krause et al., 2009;Brent et al., 2011). Several possibilities exist to deal with this fact. Link filtering is commonly FIGURE 1 | Social networks: (A) agonistic network, (B) allogrooming network. Yifan Hu layout. Gradient color of vertices (from yellow to red) represents individuals' age (from youngest to oldest). Shape of vertices represents individuals' gender (females: triangles; males: circles). Size of vertices represents individuals' degree. Size of edges represents the strength of interactions and the color is in accordance with the age of the individual that gives the behavior. used in animal SNA to delete interactions that can be attributed to random or "chance" events (Croft et al., 2008). However, at present, this filtering process has not been submitted to any formal methodology and has two major limitations: (1) the nonconsideration of weak ties (Granovetter, 1973) and thus the loss of important information (Croft et al., 2011); and (2) the sensitivity to data errors such as misidentification.
The approach adopted in this study is based on permutation tests and may help standardize the analysis of animal behavioral data obtained using SNA. In this study, I assumed an approach that allows analysis without the need to filter the links. To this end, the weight of the links must be taken into consideration, which can be done by using weighted social network metrics (Opsahl, 2009;Brent et al., 2011). Frequency-based data are less prone to sampling biases, yet by themselves, they do not solve the issue of data dependency. Therefore, I used weighted network metrics with Null Hypothesis Significant Tests (NHST) involving a permutation-based approach (Manly, 2006). This method generates a set of random values based on the real data set and creates the null hypothesis that the real structural measure X is not different from the random one. This hypothesis is accepted or rejected by comparing the observed value X to the random one. If the observed value is greater than the random one from 95%, then the null hypothesis is rejected. The use of permutation tests in the study of animal societies is discussed in details by Whitehead (2008) and Croft et al. (2011).
The following analyses were performed on both weighted allogrooming and agonistic matrices with 10000 permutations.

Network Metrics
For each individual, I calculated the following weighted network metrics: indegree, outdegree, degree, eigenvector centrality, and clustering coefficient with Ucinet 6.375 (Borgatti et al., 2002). Briefly, the degree corresponds to the total number of individuals that directly interact with one given individual (Freeman, 1979). The weighted version takes into account the weight of the links. I also differentiated between indegree (incoming ties) and outdegree (outgoing ties). This metric is historically the first and conceptually the simplest centrality network metric, and in this case, can also be considered as the activity, or "involvement, " of an individual. The eigenvector centrality index is the sum of the connections to neighbors weighted by their degree. This index provides a metric that determines the individual centrality relative to the rest of the network and the "influence" of an individual on the network (i.e., connection to highdegree nodes) and thus, on the social structure. Additionally, it would appear to be a more pertinent centrality metric for nonhuman primate groups (Kasper and Voelkl, 2009). The weighted clustering coefficient gives weight to the neighborhood densities proportionate to their size and indicates the contribution of each individual in the connectivity and thus, in the cohesion of the network structure (Watts, 2003;Hanneman and Riddle, 2005). For an overview of the weighted network metrics and calculations, see: Wasserman and Faust (1994), Croft et al. (2008) and Whitehead (2008).

Individual Level
For the first analysis, I aimed to study gender, matriline, hierarchical rank, and age-related changes in each network metric. To this end, I used general linear mixed models (GLMM) in which weighted degrees, indegrees, outdegrees, eigenvectors, and clustering coefficients are the dependent variables in separate models. Exact ages, genders, hierarchical ranks, and matrilines are the independent variables.
To offset the non-independence of these data, I realized GLMM with permutation. The consequent biological null FIGURE 2 | Matrilines. This scheme represents individuals which were kept for kinship analysis. Gray cells are dead individuals. Numbers represent the different matrilines. Each subline represents an offspring of the corresponding mother.
hypothesis was that any individual could have any network metric value. Opting for this method has several advantages. First, it takes into account the non-independency of the data; second, it is a better option than multiple t-tests and ANOVA (which both need discrete variables and would increase the number of tests) with permutation or simple correlations; and finally, it facilitates the analysis of the interactions between factors.
Initially, I created two types of models: those with no interactions with the dependent variables and those with gender interacting with other individual attributes to examine whether age, matriline, and hierarchical rank dissimilarly influence individuals according to their gender. Only factors estimated higher than 0.009 were considered significant. This threshold is arbitrary and aims to consider only significant effects with sufficient weight. These analyses were performed using SPSS 17 (SPSS, 2008) GLM Procedure with Bootstrap option of p-value = 0.05.

Group Level
The aim of the second analysis was to examine homophily and heterophily, for which I used NHST (Stephens et al., 2007) with permutation.
The principle of homophily, or the preferential interactions between same-attribute individuals, consequently determines whether the links within a same-attribute group have greater frequencies than the links between groups. Thus, to study homophily between genders, I used a simple t-test with permutation for comparing the mean of the links between and within the groups depending on gender.
To study homophily according to age (a continuous attribute), I used the Moran statistic which indexes the differences between the score of an actor and the mean, and then weights the cross products (Moran, 1950). Permutations are used to create a sampling distribution in which scores on the attribute are randomly assigned to actors. As for any permutation test, the real structural measure (the Moran statistic in this case) is compared to the random one (Hanneman and Riddle, 2005). The Moran "I" statistic of autocorrelation ranges from −1.0 (perfect negative correlation) through 0 (no correlation) to +1.0 (perfect positive correlation). These analyses were performed using Ucinet 6.375 (Borgatti et al., 2002).
The results of the agonistic social network model for gender interactions with other individual attributes are as follows (synthetized in Table    The allogrooming social network is primarily influenced by gender and age, with matriline having no significant effect. We also observe that the higher the hierarchical rank, the more an individual receives allogrooming (indegree: 1.366, p < 0.05). Furthermore, with age, individuals are less central (eigenvector: −0.014, p < 0.01), less active (degree: −5.625, p < 0.01), and give (outdegree: −3.082, p < 0.01) and receive (indegree: −2.064, p < 0.01) less allogrooming. Interestingly, whereas no significant difference was observed between males and females in the agonistic network, in the allogrooming network, females are more central (eigenvector: 0.231, p < 0.01) and more active (degree: 71.708, p < 0.01) than males. Allogrooming behaviors are mostly given (outdegree: 32.018, p < 0.01) and received by females (indegree: 39.691, p < 0.01). These results are synthetized in Table 3.
The results of the allogrooming social network model for gender interactions with other individual attributes are as follows (synthetized in

Group Level
With respect to agonistic behaviors, we obtain homophily for gender (difference in means: −1.651, p < 0.05) and for age (I = 0.273, p < 0.05). Testing genders separately, we obtain homophily by age for females (I = 0.336, p < 0.05) and for males (I = 0.211, p < 0.05). The results for matriline and individual hierarchical rank were non-significant. For allogrooming, homophily is observable for gender (difference in means: −1.942, p < 0.05) and for age (I = 0.318, p < 0.05). Testing genders separately to analyze if there are homophilic differences between genders according to age, we obtain homophily by age for females (I = 0.565, p < 0.05), but we do not obtain significant results for males according to age (I = 0.100, p = 0.106). Homophily is also observed by matriline (I = 0.321, p < 0.05) Testing genders separately, we FIGURE 3 | 3D histogram of indegree variation according to age and matriline in the agonistic network. In bold, significant attributes. observe significant homophily by matriline for females (0.458, p < 0.01), but not for males. Finally, we also observe homophily by hierarchical rank (I = 0.228, p < 0.05), yet testing genders separately, we do not observe significant results in either gender.

DISCUSSION
In this study, I established an analytical protocol that balances the inter-dependency of the data without filtering the links and that considers the weight of the links, and I analyzed the effects of several factors (gender, age, matriline, and hierarchical rank) at different levels of social organization in a non-human primate species, M. sylvanus. These findings reveal to what extent SNA facilitates the investigation of various aspects of animal societies by studying: (1) the position and influence of individuals according to their attributes; (2) the attributerelated network; and (3) the interactional dynamics reflected by homophily. In this way, I demonstrated that the sociogenesis process (rank acquisition) is intimately linked to ontogenesis (i.e., it is age-related), and differs between genders. Hence, individuals with common attributes have similar positions and roles in the group. I also stressed the existence of homophily in several behaviors, reflecting common individual behavioral patterns, including: (1) the acquisition of status within an age-related category, leading to intra-generational conflicts; (2) high-ranking individuals preferably groom similar-rank and opposite-gender individuals to secure better protection and support; and (3) the existence of homophily in grooming behaviors by gender, age, hierarchical rank, and matriline. The results suggest six main findings. First, we observe that variations in individual attributes have a greater impact on the position, role, and interactional patterns of females than on males. In most cases, these dissimilarities result from the social structure of females in M. sylvanus that is based on philopatry and MRI. Additionally, we observe significant disparities in activity and centrality between males and females.
Females are more central and active (for both received and given behaviors) in the allogrooming network. More specifically, they give and receive more allogrooming, mainly with individuals who have similar characteristics, namely females according to homophily results. These findings are in line with the literature that stresses that the philopatric gender plays a key role in affiliative behaviors (Aureli and de Waal, 2000;Silk, 2001). From a biological perspective, the philopatric gender has more time to develop a denser, stronger, and perennial network than the non-philopatric gender. In addition, female matriline homophily results emphasize the relevance of kinship bonds among females in affiliative behaviors. Individuals with high centrality and activity thus preferentially contribute to the establishment of the global network structure (Lusseau and Newman, 2004;Sosa, 2014) and cohesion of the group. In M. sylvanus, these key individuals are unquestionably the females.
Second, we observe that for female M. sylvanus, network metrics decrease with age in the agonistic and allogrooming networks. In the agonistic network, older females are less active (degree, indegree, and outdegree) and less involved in the cohesion of the network (clustering coefficient). In the allogrooming network, older females are less active (degree, indegree, and outdegree) and less central (eigenvector). During the early years, high centrality and activity for allogrooming behaviors is likely related to a long period of mother-infant and kin-related preferential interactions (which is supported by the results of matriline homophily for allogrooming behaviors) that generates kin recognition and later, kinbiased affiliative interactions (Pereira and Fairbanks, 1993). Furthermore, juveniles learn how to interact by relating to their close relatives. The observed decrease of allogrooming centrality and activity with age is likely related to the progressive stabilization of the affiliative networks of females. In addition, the decline with progressing age in the agonistic network (as related to activity and the role in the cohesion of the network) is probably a result of the stabilization of the hierarchical ranks of females when sexual maturity is attained (Chapais, 2004). These results show that the sociogenesis process, or rank acquisition, is intimately linked to ontogenesis (i.e., it is agerelated), with the latter being closely related to the reproductive status of females (menarche and postmenopause) (Borries et al., 1991). This ontogenetic process can be characterized into three stages. The first occurs before sexual maturity, when the female has numerous social interactions in order to establish her position within the group. Second, once the female is mature, she has fewer social interactions, which indicates a period of stabilization of her position. The final stage corresponds to the postmenopausal period, which can lead to even fewer social interactions resulting from exclusion (Borries et al., 1991;Sosa, 2015). Unlike females, males are not subject to the phenomenon of declining social interactions, as none of their network metrics significantly decreases with age.
Third, the frequency of given agonistic behaviors increases with the hierarchical rank of a female. More specifically, the higher her hierarchical position, the greater number of submissive individuals with whom to ensure her rank a female has (Tokuda and Jensen, 1969) and the more she intervenes in conflicts to provide support (De Waal and Roosmalen, 1979;Seyfarth and Cheney, 1984;De Waal, 1997). This does not appear to be a response to received agonistic behaviors as the indegree would also increase, which is not the case. Instead, the agonistic indegree, together with the eigenvector, decreases with the hierarchical rank of a female. This decline in agonistic indegree is also observed in males, stressing that high-ranking individuals, regardless of gender, receive fewer agonistic behaviors than lowranking ones. This reveals the benefits of dominant positions, with the reduction of associated risks (Gartlan, 1968;Bernstein, 1976;Chapais, 1991). Additionally, a significant relationship between the frequency of received allogrooming behaviors and hierarchical rank is observed in females. The absence of such phenomenon among males can be attributed to the fact that attractiveness to high-ranking individuals in allogrooming for males seems species-specific (Watts, 2000). Nonetheless, this phenomenon is observed among females according to GLMM results, which is in accordance with the theory advanced by Seyfarth (1977) in which dominant females should be preferred allogrooming partners as they can provide better protection (Watts, 2000;Cheney and Seyfarth, 2008) and support (Cheney and Seyfarth, 1990).
The fourth finding is that, similar to the hierarchy results, females within the same matriline have similar centralities and activities, and the higher the matriline, the more central and active the female is in the agonistic network. This trend can be attributed to the fact that in the genus Macaca, the hierarchical ranks of females are intimately linked to their matrilines. Furthermore, a closer examination of the previously discussed high agonistic activity of immature females shows that the indegree is more intense for high-born ones (Figure 3). Before sexual maturity, females must settle their dominance relations with lower matriline-ranking females through agonistic interactions (Chapais (2004). Thus, in their early years, females compete to establish their hierarchical rank on multiple fronts: (1) within their own matriline (as supported by the MRI phenomenon); (2) within their age category (in accordance with age homophily results for agonistic behaviors); and (3) toward older lower-ranking females.
Fifth point is that homophily results for the agonistic network show that agonistic interactions are mainly directed within the same age and gender. These findings, combined with GLMM results, yield interesting biological interpretations. Higher activity and connections in the agonistic ego network among young individuals (GLMM results) can be interpreted as a phenomenon of hierarchical rank acquisition. Homophily results highlight that this rank acquisition occurs mainly between sameage and same-gender individuals (Chapais, 1988;Holekamp and Smale, 1991), stressing the existence of a particular phenomenon that we could call the intra-generational conflict. Furthermore, the fact that affiliative behaviors are also primarily directed toward same-age individuals (age homophily results for allogrooming behaviors) underlines the trend of individuals to build their affiliative network within their age category.
Sixth, we observe homophily in allogrooming behaviors according to hierarchical rank, but only when both genders are taken into account. In other words, opposite-gender individuals with similar hierarchical ranks have preferential affiliative interactions. Allogrooming behaviors facilitate the creation of affiliative bonds and potential support in future conflicts (De Waal and Roosmalen, 1979;Seyfarth and Cheney, 1984;De Waal, 1997). Subsequently, this behavioral pattern can explain the decrease of agonistic behaviors received by high-ranking individuals (observed in GLMM results) owing to the highranking support of a third party.
Homophily has previously been reported in many species (McPherson et al., 2001;Lusseau and Newman, 2004;Massen and Koski, 2013). Extensive research on human homophily stressed that it is a major mechanism in stranger cooperation (Haun and Over, 2015), social learning (Buttelmann et al., 2013), and cultural and norms transmission (Chudek and Henrich, 2011). Recent studies argue that homophilic preferences may explain the gap between animals and humans regarding these abilities (Haun and Over, 2015). Revealing homophily in several behaviors and as it is influenced by different attributes highlights the importance of these mechanisms in a non-human primate species. However, many methods exist to evaluate the presence or absence of homophily (E-I index, ERGM, assortativity, Moran I statistic), each one of them with inherent pros and cons that would need to be evaluated before determining which of these approaches is more relevant for studying animal societies.
This analytical protocol can be used to study other animal societies and might enable interspecific comparisons. I believe that the important findings of this study might help understand the global patterning of a non-human primate society, and likely other animal societies, from an evolutionary perspective.

AUTHOR CONTRIBUTIONS
The author confirms being the sole contributor of this work and approved it for publication.

ACKNOWLEDGMENTS
First of all, I would like to formulate special thanks to Marine Huser for her appreciated help during this study. I am grateful to Francesc Salvador Beltran and Vicenç Quera for guiding and advising me throughout this investigation. Many thanks to Ellen Merz who allowed me to lead the research field work at 'La Forêt des Singes' and provided me the demographic data; I sincerely thank her for her advice. Salvador Herrando, David Leiva Ureña and Roland Hilgartner provided helpful advice at different stages of the analytical protocol. I would also like to thank Frontiers in Psychology reviewers for their precious suggestions, which largely participated to the improvement of this work. Finally, thanks to all the anonymous who have encouraged me from near or far and allowed the existence of this work.