Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The role of gender in social network organization

  • Ioanna Psylla,

    Roles Conceptualization, Data curation, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliation DTU Compute, Technical University of Denmark, Kgs. Lygby, Denmark

  • Piotr Sapiezynski,

    Roles Conceptualization, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliations DTU Compute, Technical University of Denmark, Kgs. Lygby, Denmark, College of Computer and Information Science, Northeastern University, Boston, United States of America

  • Enys Mones,

    Roles Conceptualization, Visualization, Writing – original draft, Writing – review & editing

    Affiliation DTU Compute, Technical University of Denmark, Kgs. Lygby, Denmark

  • Sune Lehmann

    Roles Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing

    sljo@dtu.dk

    Affiliations DTU Compute, Technical University of Denmark, Kgs. Lygby, Denmark, The Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark

Abstract

The digital traces we leave behind when engaging with the modern world offer an interesting lens through which we study behavioral patterns as expression of gender. Although gender differentiation has been observed in a number of settings, the majority of studies focus on a single data stream in isolation. Here we use a dataset of high resolution data collected using mobile phones, as well as detailed questionnaires, to study gender differences in a large cohort. We consider mobility behavior and individual personality traits among a group of more than 800 university students. We also investigate interactions among them expressed via person-to-person contacts, interactions on online social networks, and telecommunication. Thus, we are able to study the differences between male and female behavior captured through a multitude of channels for a single cohort. We find that while the two genders are similar in a number of aspects, there are robust deviations that include multiple facets of social interactions, suggesting the existence of inherent behavioral differences. Finally, we quantify how aspects of an individual’s characteristics and social behavior reveals their gender by posing it as a classification problem. We ask: How well can we distinguish between male and female study participants based on behavior alone? Which behavioral features are most predictive?

Introduction

For many decades, gender differentiation has been studied as an interdisciplinary topic and within a variety of fields including psychology, social science, anthropology, history, and biology. Existing studies have explored the nature of the existing gender differences, their origin, and impact on individuals’ lives. How to interpret the observed deviation between women and men is subject to debate among scholars. It is, however, universally accepted that behavioral differences are rooted in the different biological roles, and are reinforced by a society’s values and cultural beliefs.

Previous research has shown that gender-specific inequalities might originate from biological predispositions (e.g. hormones [1], brain structure [2]), as well as the organization of the hunter-gatherer societies in which humans initially evolved [3]. This differentiation is subsequently aggravated by cultural/societal expectations [1, 4], which are likely to lead the two genders to develop and maintain their social ties in different ways [5]. The study of social networks is essential for understanding how gender role influences the nuances observed in the structure and evolution of these social interactions. Although it does not provide answers regarding the origins of the gender differences in social behavior, it can help identify and understand these discrepancies to a larger extent.

Below, we first explore individual-level characteristics, specifically the psychological traits and mobility behavior within the cohort, noting that both these aspects have been found to relate strongly to social behavior [6, 7]. Next, we focus on social traits. For the two gender-groups, we evaluate similarities and differences with respect to social network role. Our analysis of social networks is based on longitudinal data describing person-to-person interactions (physical proximity using Bluetooth scans), calls and text messages, and online friendships (based on Facebook communication activity). Finally, we use classification models to quantify the extent to which a person’s gender can be inferred from their observed characteristics and behavior.

Data

The basis of this paper is the Copenhagen Network Study (CNS), a study focusing on nearly 800 freshmen at the Technical University of Denmark [8] who volunteered to donate data via Nexus 4 smartphones. The bulk of data collection was behavioral data from from the smartphones, supplemented with data from online questionnaires and 3rd party APIs, such as the Facebook Graph API. The derived datasets include:

  • Friendship graph and interactions (comments on wall posts) from Facebook,
  • Person-to-person proximity events, measured using Bluetooth,
  • Telecommuncation data (call and text message logs; only metadata, no content),
  • Location records (based on GPS and WiFi),
  • Questionnaires (responses to personality questionnaires, described in detail below).

This work is based on data collected between September 2013 and May 2014. The number of active participants and the quality of their data varies over the duration of the observation. To eliminate the effect of missing data on statistics, we calculate all indicators and network properties on a weekly basis and average for each individual. Participants with three active weeks or less during the nine-month period are excluded from the analysis.

After this filtering, this dataset consists of 166 female and 601 male students. In order to avoid the difference in population sizes affecting the standard deviations, we apply subsampling over the male and female population separately and calculate the distribution over the mean of the random sub-samples.

Below, unless otherwise specified, we use the following strategy to compare the two (female/male) classes. We we draw 1000 random subsamples, each equal to the half of the original class size from each class. Then, we perform pairwise comparisons between subsamples. We test the null hypothesis that the means of the two sampling distributions are identical (a two-tailed test).

In order to compare results across domains (personality, mobility, social interactions, etc) we measure the differences in distributions between the two genders using effect size r, defined as the ratio between the means of each distribution x1, x2 and the pooled standard deviation σp: (1) where σp(x1, x2) is defined via (2) and σ2(x) is given by (3)

Personality

In this section, we investigate how gender differences are expressed through personality metrics. Data from responses to personality questionnaires show that although there are considerable variations within a gender, differences between males and females exist in a number of traits and at every age. As part of the CNS study, we consider the following dimensions of personality, which are listed below along with the central gender-related results pertaining to that measure.

  1. Big Five. The Big Five Inventory (BFI) is a widely used method for assessing human personality using five broad factors: openness, extraversion, neuroticism, agreeableness, and conscientiousness [9]. To measure big five, we use the questionnaire developed in Ref. [9]. Previous work has consistently found women to be more neurotic and agreeable than men [3, 1014]. There is less of a consensus with respect to gender differences in the remaining BFI attributes. For instance, some studies report higher conscientiousness and openness among women, while others find men as more conscientious [11, 12]. Detailed description of each personality trait and reference to additional literature are provided in Table 2.
  2. Self-esteem. We use the definition that self-esteem is a feeling of self-worth [15] and use Rosenberg’s 10-item instrument to measure it [15]. Feingold [14] found that males have slightly higher self-esteem than females, and Kling et al. [16] showed that this effect increases considerably in late adolescence. However, other studies exist that show no significant difference between males and females with respect to self-esteem [4].
  3. Narcissism. Narcissism has been previously found to be positively correlated with self-esteem [17, 18]. Here we assess Narcissism using the Narcissistic Admiration and Rivalry Questionnaire (NAR-Q), which integrates two distinct cognitive and behavioral aspects of narcissism: the tendency to approach social admiration through self-enhancement, and the tendency for an antagonistic self-defense (rivalry) [19]. The literature is consistent here: men tend to be more narcissistic than women, regardless the age and income [2023].
  4. Stress. Several studies have been conducted to measure stress levels among students in higher education, reporting that female students tend to have more stress (and more stressors) than male students, regardless of the instrument used for measurement [2427]. In this study, we measure stress with the widely used Perceived Stress Scale (PSS) [28].
  5. Locus of control. Locus of control reflects the extent to which a person perceives a reward or reinforcement as contingent on his own behavior (internal locus) or as dependent on chance or environmental control (external locus) [29]. We measure locus of control using a simplified, 13 item scale proposed by Goolkasian (see: http://www.psych.uncc.edu/pagoolka/LC.html). A lower score corresponds to internal locus, whereas a higher score indicates external locus. In general, the two genders have not been found to differ with respect to this psychological trait [14, 2931]. Lefcourt [32] argued that those who are classified as having an internal locus of control not only perceive but also desire more personal control than individuals with an external locus and found a that females desired greater internal control than males. However, women have been found to favor external control in items related to academic achievements [33].
  6. Satisfaction with life. Satisfaction with life constitutes a judgment of one’s life in which the criteria for judgment are up to the person [34]. We use the satisfaction with Life Scale (SWLS) instrument [35], which has been widely used to assess subjective well-being within various groups of population. The SWLS includes five generic statements, in which a subject must respond with a 1-7 scale, indicating the degree of agreement or disagreement. Results regarding gender have been shown to be highly dependent on age. Specifically, adolescent and elderly males have higher life satisfaction than females, while no observable difference is found among young adults [36, 37].
  7. Loneliness. We measure loneliness using the UCLA Loneliness Scale, a 20-item scale, in which a subject must indicate how often they feel an item characterizes them [38]. Male college students have been found to be more lonely than female students [38, 39]. It has also been shown that men are less willing to acknowledge feelings of loneliness, due to their more pronounced negative consequences of admitting to this feeling [40, 41].

In summary, women tend to score higher with respect to negative emotionality (such as neuroticism and stress) than men, but it has been argued that this may be due to females more readily admitting to or perceiving such intense feelings. Individualism also plays an important role in personality differences between the two genders. In the present study, we analyze the gender effect on the aforementioned personality traits in an environment where females are the numerical minority, and within a highly specific group of individuals (students at a technical university). The diverse dataset, however, allows us to combine the results from the questionnaires with the participants’ behavior in a natural setting.

Results

We test the null hypothesis that the two samples have equal means. We start with the Big Five Inventory, and measure effect sizes between the two genders. Fig 1 shows the normalized difference (i.e., the effect size) observed between males and females with respect to neuroticism, conscientiousness, agreeableness, extraversion, and openness. Each histogram represents the distribution of the difference in means normalized by the pooled standard deviation, and the mean in a subsample of females is subtracted from the mean in a subsample of males (for details on the effect size, see Eq (1)). The horizontal bars denote 5 and 95% percentiles. Neuroticism exhibits the largest deviation, positioned far to the left from the zero baseline (with a mean of dneu = −0.635), indicating that women score significantly higher in this personality characteristic than men. We also find significant, albeit less pronounced, differences with respect to conscientiousness (dcon = −0.436) and agreeableness (dagr = −0.259). Finally, we do not find statistical significance in the average values of extraversion (dext = −0.118) and openness (dope = 0.143). Fig 2 depicts the results describing the remaining personality measures. Stress is significantly higher among women (dstr = −0.451), while it is clear that men score higher in self-esteem (dse = 0.423). Overall, narcissism is higher among male students (dnar = 0.349), but mainly due to rivalry (driv = 0.334), which is its antagonistic aspect and less because of admiration (dadm = 0.241), which constitutes the assertive aspect of narcissism. We find that women score higher in I-E Rotter Scale (dloc = −0.157), indicating a greater average sense of external locus of control. Women also score higher with respect to satisfaction with life (dsat = −0.149). Finally, we find no statistically significant difference for loneliness (dlon = 0.095).

thumbnail
Fig 1. Signed effect sizes of the BFI measured between men and women.

We find female participants to be more neurotic and agreeable, in line with previous research [3, 1014]. Women in our study tend to be also more conscientious, and we identify no significant differences in scores for extraversion and openness. Negative values indicate that women achieve higher scores. Histograms show the distribution of effect sizes defined by Eq (1), horizontal bars denote 5 and 95% percentiles.

https://doi.org/10.1371/journal.pone.0189873.g001

thumbnail
Fig 2. Signed effect sizes of personality traits measured between men and women.

As indicated in previous studies on college population, women tend to feel more stress than men [2427] and have a more external locus of control in items related to academic achievement [33]. Men prove to be more narcissistic [2023] and to have a higher feeling of self-worth [14, 16]. Despite these results, and in contrast with previous literature, women in our study tend to report a higher satisfaction with life. Negative values indicate that women achieve higher scores. Histograms show the distribution of effect sizes defined by Eq (1), horizontal bars denote 5 and 95% percentiles.

https://doi.org/10.1371/journal.pone.0189873.g002

Mobility

In this section we verify whether there are observable differences in mobility traces between the participants of the two genders. We begin by providing a brief overview on the literature discussing differences between male/female mobility patterns, then discuss findings from our cohort.

There is a general consensus that mobility patterns are not gender neutral and womens’ mobility through urban space is distinguishable from mens’. Differences between men and women in their mobility have been ascribed to various components of the gender role, such as gender-related tasks, distinct family roles, and labor market position [42]. Men and women are assumed to perform a similar number of trips, but with distance traveled and the mode of transportation differing between them. Specifically, surveys conducted in Western countries in the ’90s have demonstrated that women travel fewer kilometers than men and make more trips as pedestrians and using public transportation [43, 44]. Moreover, the purpose of travel tends to differ, with women traveling most frequently for household errands and men making a majority of trips to work. Other studies explain the shorter commuting distances of women as a result of their weaker position in the labor market [45]. Interestingly, females have been observed to travel longer distances and explore larger areas in foraging tribes, the reason for this difference is argued to originate from the fact that women are expected to return home more frequently while gathering than men are while hunting [46].

Recent studies based on mobile phone records, however, have not observed substantial differences in travel distances [47], regularity, and predictability of movements [48] between male and female commuters. However, a study using travel diary data collected in Portland reports higher levels of activity among part-time employed women than those of part-time employed men throughout the day [49].

In conclusion, despite of recent advances in studying mobility behavior in detail based on high resolution observational data, gender-based differences are rarely observed.

Results

We follow the same procedure as in personality-related measurements: we apply subsampling to obtain equal sample sizes, calculate the effect size, and test the null hypothesis that the means of the two distributions are equal.

A common quantitative description of mobility behavior is given by the distribution of unique locations visited by an individual over some time period, e.g., using Pu(l), which is the relative frequency of visiting location l by individual u [48, 50]. Relative frequency is given by the relative time the individual spent at some location on a weekly basis. We analyze location data obtained by periodically collecting the position estimate from the location sensor of the students’ phones. The list of unique locations that characterize an individual is extracted as a list of clusters of location measurements a DBSCAN-based algorithm developed in Ref. [51] and validated in Ref. [52].

To further quantify individual mobility patterns, we measure the heterogeneity of the locations visited over time using entropy. Entropy is a measure of uncertainty or predictability of a distribution. Here we use entropy to capture the heterogeneity of an individual’s time spent across unique locations. Using Pu(l), it is defined as (4) where L(u) denotes the set of locations for user u. Individuals distributing their time more evenly are characterized by higher entropy. The effect sizes measured in the location related metrics are shown in Fig 3. We find that women both visit more unique locations over time, and they have more homogeneous time distribution over their visited locations than men, indicating that time commitment of females is more widely spread across places.

thumbnail
Fig 3. Effect sizes of the mobility indicators.

We show that women are characterized by higher uncertainty in their travel patterns and visit more unique places. This deviation has not been observed in recent literature. Negative values indicate that women score higher than men. Histograms show the distribution of effect sizes defined by Eq (1), horizontal bars denote 5 and 95% percentiles.

https://doi.org/10.1371/journal.pone.0189873.g003

Networks and interactions

Now we turn our attention towards social interactions among the students. We begin by providing a brief overview on the literature discussing differences between male/female network structures. We then discuss findings based on our cohort.

Previous work suggests that the sizes of real-world ego networks of the two genders are drawn from similar distributions [5355]. In contrast, women tend to have more friends online, as seen in multiplayer games [56] or social networking services [5759]. A study based on Facebook data describing around 1800 U.S. college students found that females show higher social activity and have greater betweenness centrality in their Facebook network compared to males [60].

Social networks display high gender homophily, both offline [61, 62] and online [59, 63]. The extent of preferring same gender friends varies with age, with e.g. girls forming smaller, more homogeneous groups than boys at young ages [64]. As soon as adolescents begin forming romantic ties during puberty, women start to invest more heavily in opposite-sex relationships; but they shift preferences to younger women (presumably daughters) as they age [65]. Men, on the other hand, are shown to increase their female contacts as they get older and particularly at the end of their life cycle [66]. Interestingly, heterophily between genders is prominent among the strongest ties. For instance, calls and text conversations are both more frequent and longer among mixed-sex pairs of individuals [66, 67].

Homophily has been studied as a function of transitivity (a measure of the probability of two individuals being connected provided they are both connected to the same alter) [68]. In this case, structural factors, such as network proximity, have been found to have a stronger effect on triadic closure compared to homophily [69]: a high number of shared contacts is a better indicator of triadic closure than sharing an attribute. A study based on data from several U.S. elementary schools reports that females form more triads than males and that dyads consisting of females are more likely to be in triangles [70]. Kovanen et al. [71] studied temporal gender homophily in 3-motifs using a large dataset of mobile phone records. They find that female-only motifs are over-represented compared to a reference model, whereas male-only motifs are under-represented. Contradicting the aforementioned findings, however, a study based on data from the Spanish social networking site Tuenti, found high levels of homophily in females’ dyads but a higher tendency of male users to form same-sex triangles [63].

Women have not only been found to be more actively engaged in online interactions, but also to spend more time engaged in phone conversation [72]. In a review paper, Smoreda and Licoppe [73] report that women tend to disclose more information to correspondents (especially about intimate topics) and are more expressive than males, which results in longer conversations, whereas men communicate mainly for instrumental purposes. In addition, other studies have shown that calls to a woman are longer than calls to a man regardless the gender of the caller [66, 73]. Circadian rhythms in call patterns have revealed further differences between men and women, with women making longer phone calls in the evenings and during the night, and mainly when the recipient is a man (which indicates an emphasis to romantic relationships) [74]. Likewise, young women have been reported to send a greater number of text messages, especially if the receiver belongs to the opposite gender [66].

In summary, previous studies found clear differences in the way men and women engage with their social networks. However, most of the studies focus on a single channel of interactions (e.g. online communications, behavior in an organization), failing to capture a potential persistence or deviation of the characteristics across different settings. Here we use the CNS data to compare communication across a number of different channels.

Results

We consider three types of communication: physical proximity (i.e., person-to-person) interactions, Facebook activity, and mobile phone communication (calls and text messages). Previous studies have shown that each channels may describe different aspects of social ties and potentially corresponding to different levels of connection intensity [7579]. To illustrate these differences, in Fig 4 we report the fraction of active links in each communication channel over time. Note that the vertical scale is logarithmic, indicating an increased presence of proximity links (purple pentagons), a moderate level of active Facebook connections (red circles) and a comparably low level of active telecommunication links (bottom curves).

thumbnail
Fig 4. Fraction of active dyads in the different networks.

Curves show the ratio between existing and potential links between participants in each network. All students attend classes at the same campus and eat at the same cafeterias, so their proximity network is very dense (with 40% of dyads active). Only about 2-3% of them actually connect as friends on Facebook, and less 1% communicate using calls, text messages, or Facebook interactions. Each network is unweighted and aggregated over a month—a link exists between two nodes if they had any interactions in a given channel during that month.

https://doi.org/10.1371/journal.pone.0189873.g004

Fig 5 visualizes a single-day snapshot of the three networks in this study. The community structure revealed in the person-to-person network is a fingerprint of the classes the students attend together. Less structure is observed in the Facebook feed network, which shows a single large component; the call network is the most fragmented of all of the networks. The high level of homophily in the proximity network is evident from the frequency of orange and blue lines, which represent the female-female and male-male connections respectively, as well as cliques that contain nodes of the same color. The call network shows the highest fraction of mixed gender connections, a possible indicator of couples.

thumbnail
Fig 5. Snapshot of the three interaction networks: (a) person-to-person aggregated over one weekday, (b) Facebook interactions (aggregated over a month) (c) phone call networks of the students (aggregated over a month).

Blue and orange nodes correspond to male and female students respectively, link color denotes whether the connection is between the same genders; orange for female-female, gray for mixed, and blue for male-male connections. The size of each node is proportional to their degree and width of the links represents the frequency of interactions. The person-to-person network shows clear separation into study groups, while this structure is no longer visible in communication networks.

https://doi.org/10.1371/journal.pone.0189873.g005

When investigating the networks between the study participants, we apply a different approach to accounting for the imbalance male/female subjects than in case of personality and mobility. Here, subsampling would alter the network structure and, thus, render e.g. measurements of homophily and other network metrics meaningless. Instead, we use the following reference model: we randomly permute genders between participants with uniform probability and then perform the calculations. Overall, we produce 2E network realizations, where E is the number of edges in the network.

To approve or reject the null hypothesis that the network is independent of gender homophily, we calculate the z-score, given by: (5) Here, x is an indicator, and are the mean and standard deviation of the indicator in the reference model. The z-score is expected to be zero if the null hypothesis is true.

To test the null hypothesis of no difference between the two gender groups, we draw the permutation distribution of the differences between the two genders and measure where this distribution falls relative to the mean difference of the empirical data. The p-values then are calculated by dividing the number of permuted mean differences that are larger/smaller than the one observed in the empirical data, by the number of items (2E) in the permutation distribution.

We explore the influence of gender homophily on formation of friendships in the various networks among the participants. To do this, we first identify the fraction of same-gender friends out of all friends an individual has. Fig 6 shows the z-scores of various network connections obtained by comparison with the permutation model (see also Methods). Women have remarkably more same gender friends than the ones measured with the reference model in online interactions and person-to-person interactions (z-score is 13.10 and 12.24 respectively). On the other hand, men also show a preference for forming homophilous ties through mobile communications, though to a less extent. To study whether men and women tend to form closed triangles with same-gender alters, we count the various motifs in each network. Results are shown in Fig 6 (color bars): male only (blue), female only (orange) and mixed (brown). Furthermore, we compare the results with the respective distribution of the expected motifs found in the reference model for the Facebook network (Fig 7). We find that male-only triads are insignificantly underrepresented compared to the reference model, and that there are more female-only motifs than what we would observe by chance (z = 13.101, p < .0001). Whereas a similar pattern is observed in person-to-person interactions, same-gender motifs are overrepresented for both genders in mobile communications. We conclude that women prefer other women for both their dyadic and triadic relationships in every form of interaction, while homophily is noticeable among males only in their trusted interactions through the phone.

thumbnail
Fig 6. Homophilous connection patterns in the three interaction networks.

Bars represent the z-scores of connection motifs (containing two or three individuals) obtained from the comparison with the permutation model. Grey color corresponds to dyadic connections, color bars denote triangles: male only (blue), female only (orange) and mixed (brown). Women in the study are more likely to form same-gender triangles in all types of interaction networks.

https://doi.org/10.1371/journal.pone.0189873.g006

thumbnail
Fig 7. Distributions of same-gender triangles in the Facebook friendship network.

Empirical observations (black dashed lines) are compared to the distribution of triangles measured in the ensemble of null models. Values on the x-axis correspond to the number of same-gender triangles in the network. It is noticeable that females form triangles with other females to a higher extent compared to the reference model, whereas males show the reverse pattern.

https://doi.org/10.1371/journal.pone.0189873.g007

We find that women tend to have a significantly higher number of contacts than men in both online and mobile networks, whereas the size of the person-to-person networks are similar. Fig 8 shows how degree varies over time in terms of mobile communication (calls); females have more contacts during nearly the entire period of interest. We measure betweenness centrality (see Methods for the definition) of each individual to investigate whether one gender tends more prominently positioned in a network than the other. We find that women consistently show higher betweenness indices, regardless of the mode of interaction.

thumbnail
Fig 8. Degree of males and females per week in calls.

Curves show the mean number of call contacts for males (blue) and females (orange). Women tend to have more call contact on average than men throughout the duration of the observation.

https://doi.org/10.1371/journal.pone.0189873.g008

Next we study the entropy of interactions. Similarly what we did for mobility distribution in Eq (4), we calculate the entropy of the distribution of interactions over the contacts: (6) where Pu(i) is the probability that user u interacts with their i-th contact in his ego-network N(u). The value of Su is estimated by the corresponding number of interactions relative to all interactions performed by user u. Individuals who interact equal amounts with many friends will have high entropy (and therefore can be characterized by lower predictability [80], whereas those who limit the vast majority of their interactions to a small set of others are expected to have low entropy (more predictable). In Fig 9, we plot the distribution of entropy effect sizes measured between males and females for the three interaction networks. We observe a significant difference for Facebook and calls, with women displaying higher entropy than men, indicating that females distribute their interactions with friends considerably more homogeneously. In addition, females exchange remarkably more text messages than males (p < .001). With respect to time spent on social interactions, we find that in our study, women are described by significantly longer conversation times during phone calls than men, regardless of the initiator of the call (p < .0001). The longest calls (on average) are observed on ties where a male initiates contact to a female (with an average duration of 117 seconds), with second longest average call-durations observed between females (an average duration of 114.56 seconds). The shortest average duration (71.52 seconds) is measured between pairs of males.

thumbnail
Fig 9. Effect size distribution of network metrics.

Women tend to exchange more text messages than men and distribute their communication more evenly across their call and Facebook contacts. We do not identify a significant difference in proximity entropy. Person-to-person proximity is, to a large extent, driven by class schedule and participants have less choice regarding who they interact with. Negative values indicate that women achieve higher scores. Histograms show the effect sizes of network related indicators (number of text messages and entropy) between males and females. Horizontal bars denote the interval of 5 and 95% percentiles.

https://doi.org/10.1371/journal.pone.0189873.g009

Gender prediction

Based on the findings presented above, we consider the classification problem of predicting gender based individual and social characteristics. In the literature, there have been attempts to predict gender based on Call Detail Record (CDR) data using semi-unsupervised techniques and deep-learning algorithms [81, 82]. De Montjoye et al. [83] found that gender is a strong predictor of neuroticism, a trait that is seen in the literature to be consistently higher among women. In this study, we combine the questionnaire data, mobility patterns, as well as social interaction habits of each participant, to build a dataset that offers adequate complexity for achieving a good performance in the gender-inference problem. Additionally, the machine learning process also provides insight into the question: What are the most predictive behavioral indicators of gender.

Results

We use the behavioral measures calculated above as features to train four different models: logistic regressor, AdaBoost, support vector classifier (SVC), random forest, and gradient boosting classifier implemented in the scikit-learn Python package [84]. Each models is evaluated using 10-fold stratified cross-validation. Each of the models underwent a hyper-parameter fine tuning procedure described in detail in the Methods section.

Men constitute 78% of the study participants. This poses a significant imbalance in the data and therefore, we measure performance using the area under receiver operating characteristic curve (ROC-AUC) which is robust against imbalance, as well as using the F1 score that is sensitive to the imbalance. The value of ROC-AUC can be interpreted as the probability that the classifier is able to identify the female in of a male/female pair. The F1 score is the harmonic mean of precision (what fraction of people identified as women are actually women) and recall (what fraction of women are identified as women) at a selected threshold.

Results are summarized in Fig 10 for each classifier along with the corresponding values of the random classifier based on the imbalance present in the data. All classifiers after the hyper-parameter fine tuning procedure perform similarly well, with ROC-AUC values of 0.86 and F1 scores of 0.5 and higher (compared to a random classifier with ROC-AUC of 0.5 and F1 of 0.22).

thumbnail
Fig 10. Classifier performance in gender prediction.

Reported are the ROC-AUC (black) and F1 (gray) scores for the five classifiers considered in the gender prediction. Dashed lines mark the baseline obtained from the random classifier that takes the imbalance into account. High values of area under receiver operating characteristic (ROC-AUC) indicate a strong separation between men and women in the feature space: when presented with a random male and a random female, the classifier will identify the female correctly 87% of times.

https://doi.org/10.1371/journal.pone.0189873.g010

Next, we investigate the question of which behavioral features are most informative regarding gender. To do this, we use the feature importance obtained by fitting a random forest model to the data (Fig 11). We find that that a tendency towards gender homophily in the social networks is the most important behavioral feature; this is true for all three types of interactions that we consider. Some aspects of personality are also important. Within the big five traits, neuroticism and conscientiousness are most predictive, while narcissism and self-esteem are most powerful among the remaining personality tests. High on the list, we also find various communication related network characteristics. With respect to feature types, network indicators are the most important ones, occupying five of the top six indicators.

thumbnail
Fig 11. Feature importance obtained by random forest.

Color shows the corresponding feature types: BFI (red), personality traits (green), mobility (yellow) and network (purple). The most important features are those who help separate men from women most: women, on average, reveal higher degree of homophily, neuroticism, and conscientiousness, while men tend to score higher in tests on narcissism, and self-esteem.

https://doi.org/10.1371/journal.pone.0189873.g011

Discussion

Conclusion

In this work we have studied gender differences within a population of freshmen. We have been able to identify a number of gender differences in personality traits (measured via questionnaires), mobility patterns, as well as social network behavior based on person-to-person, telecommunication, and online social networks.

Personality. Starting from gender differences with respect to personality, our findings are in accordance with observations in the psychology literature on gender differences. Discrepancies (or differences that are not significant) correspond to personality traits that, according to previous research, display ambiguous behavior over genders.

Mobility. With respect to mobility behavior, our results are not consistent with findings in the literature. Previous work has found a restricted travel space for women [42], but we find that women travel more than men on average.

Networks. Humans use multiple channels when we communicate: real-world conversations unfold when we meet person-to-person, we call each other and send text messages, we engage with Facebook posts and write comments, send email, and use other instant messaging platforms. Based on the communication channels we have access to here, we find that each of these channels plays a slightly different role with respect to gender similarities and differences. First, within all networks, we observe differences with respect to gender homophily, specifically an over-representation of female-only dyadic and triadic connections. This over-representation is mostly emphasized in weak links, which is consistent with the literature [56, 59, 71, 85]. Males show a lower level of homophily. For stronger social ties (that is, contacts require more effort to maintain), both genders show similar level of homophily. Second, we find that women tend to simply communicate more. On average, females maintain more contacts than males in the population, they exchange significantly more text messages, and talk longer on the phone. This is expressed both via a larger number of contacts as well as a higher entropy of neighbors. Furthermore, in agreement with a previous study on Facebook [60], we find that, in general, women have more central positions in the network, as expressed through a higher average betweenness centrality.

Predicting gender. Finally, we use the features discussed above to predict gender based on personality and behavioral patters. The prediction task is based on combined personality, mobility and network features for each individual study participant, allowing us reveal the relative importance of each feature in predicting gender. We find that personal characteristics and social behavior can be used to identify the gender of an individual with high performance (AUC = 0.87). We find that network features are highly revealing, followed by personality test scores.

Limitations

Our results point out significant differences in various aspects of social behavior between males and females, based on a population of nearly 800 freshmen at a large European university. However, in order to have a clear understanding of the results, it is important to note the limitations of the dataset as well as the methods applied, which we address in the following paragraphs.

Population sizes.

In our dataset, male population is around four times larger than the female population. This skewed female/male ratio presents may present certain biases simply because, in addition to gender differences, women may behave like a minority in some cases (e.g. with respect to homophily in the social networks). The female/male ratio also presents methodological challenges. We describe our approach to mitigate these issues, both at the individual level as well as for the network analysis in the Methods section.

Demographics.

The cohort of the CNS experiment consists of Danish and international freshmen at a technical university. Population impose strong constraints on demographics, with respect to age and social embeddedness. Furthermore, we do not have detailed demographic information regarding the contacts the participants made outside the experiment. Although demographic information is necessary to understand the results, and individuals located in Denmark indeed display different personal behavior (for example, a low overall level of neuroticism), our results regarding the comparison of males and females are in agreement with existing literature on personality traits.

Non-binary gender identification.

In this work, as in Copenhagen Networks Study in general, the participants reported their gender through a questionnaire that only offered two options: female and male. This limiting distinction might have contributed to additional noise in the measurement of differences as well as to lowering the performance of the models in the gender-inference task.

Materials and methods

Ethics statement

The Danish Data Protection Agency has approved the overall project structure (data collection, anonymization, and storage), as well as the content of the current study, cf Journal: 2012-41-066. The project complies with both local and EU regulations. All participants in the study have provided written informed consent.

The data obtained from Facebook was collected in accordance with the Terms of Service.

Network metrics

We construct three types of networks representing the various interactions among the participants: physical proximity, Facebook, and call networks. We then aggregate them over time windows of one week. Only consenting participants of the CNS are included in the networks, since we do not possess complete information (e.g, gender and social activity) about the other contacts of the students. However, extending the ego networks of the students with individuals outside the experiment (for instance, in the call networks), we can extract additional descriptors about the participating students, such as their total number of contacts or distribution of mobile phone conversation times.

In the present study we show statistics over different network metrics that are related to the local structure of the graphs and the position and role of participants in the global network. In the following, we provide a detailed description of the applied network metrics.

  1. Degree. The number of contacts ki an individual has in their respective social network. We calculate the degree in two different settings. First, by considering all contacts of a student, we can infer the total degree (without referring to the gender of the contacts). Second, by limiting the interactions to the participants, we calculate the degree describing same-gender contacts.
  2. Betweenness centrality. This measure quantifies the importance of an individual with respect to information flow on the network, when the shortest paths are taken into account. It is defined as: (7) Here, denotes the number of shortest paths between individuals j and k among which number of paths go through individual i. Therefore, betweenness effectively measures the fraction of shortest paths that pass an individual, which is a precursor of their relevance in case of any flow on the network (rumor propagation, spread of an infectious disease, etc).

Imputation of missing values

Due to the method of data collection some fraction of students has missing data in various channels. Overall, 21.5% of the participants exhibit missing data in at least one channel. To address the problem, we first remove participants with missing features in more than two of the five feature categories (personality, location, call, Facebook, and person-to-person interactions). We then apply a KNN based imputation to the remaining data [86], described as follows. For each user we find their k-nearest neighbors (with k = 7) by calculating the average difference of non-missing features with other users. We only use features that are present in the potential neighbor’s feature set, that is, if Luv = FuFv denotes the set of overlapping features of users u and v, the distance between the users is given as: (8) where and are the values of the i-th feature for users u and v respectively, and |Luv| is the size of the overlap set. Once the k-nearest neighbors of all users are determined, for each student we impute their missing values by the average of the corresponding non-missing feature values of their neighbors. If there is a single neighbor, their value is assigned.

Fine tuning of the machine learning models

We fine-tuned each of the models used in the gender prediction task. Through a grid search with cross validation we found the set of hyper parameters for which each model achieved the highest harmonic mean between F1 score and ROC-AUC on previously unseen data. Table 1 lists the parameter values in the grid. Optimal values are bold.

thumbnail
Table 1. Each model used in the gender prediction task was fine-tuned.

Optimal values for each parameter are marked in bold.

https://doi.org/10.1371/journal.pone.0189873.t001

Personality traits

We consider eleven personality traits in the main paper, with five traits forming the Big Five Inventory. Table 2 lists all the personality traits along with their definition and references to the corresponding literature.

thumbnail
Table 2. Definition of the personality traits considered in this paper for the male and female freshmen students.

Traits marked with an asterisk (*) are part of the Big Five Inventory. References provide further reading on the various personality traits.

https://doi.org/10.1371/journal.pone.0189873.t002

References

  1. 1. Wood J. Gendered lives. Nelson Education; 2012.
  2. 2. Reisberg D. Cognition: Exploring the science of the mind. WW Norton & Co; 1997.
  3. 3. Schmitt DP, Realo A, Voracek M, Allik J. Why can’t a man be more like a woman? Sex differences in Big Five personality traits across 55 cultures. Journal of personality and social psychology. 2008;94(1):168. pmid:18179326
  4. 4. Maccoby EE, Jacklin CN. The psychology of sex differences. vol. 1. Stanford University Press; 1974.
  5. 5. Deaux K, Major B. Putting gender into context: An interactive model of gender-related behavior. Psychological review. 1987;94(3):369.
  6. 6. Kalish Y, Robins G. Psychological predispositions and network structure: The relationship between individual predispositions, structural holes and network closure. Social Networks. 2006;28(1):56–84.
  7. 7. Cho E, Myers SA, Leskovec J. Friendship and mobility: user movement in location-based social networks. In: Proceedings of the 17th ACM SIGKDD international conference on Knowledge discovery and data mining. ACM; 2011. p. 1082–1090.
  8. 8. Stopczynski A, Sekara V, Sapiezynski P, Cuttone A, Madsen MM, Larsen JE, et al. Measuring large-scale social networks with high resolution. PloS one. 2014;9(4):e95978. pmid:24770359
  9. 9. John OP, Srivastava S. The Big Five trait taxonomy: History, measurement, and theoretical perspectives. Handbook of personality: Theory and research. 1999;2(1999):102–138.
  10. 10. Costa P Jr, Terracciano A, McCrae RR. Gender differences in personality traits across cultures: robust and surprising findings. Journal of personality and social psychology. 2001;81(2):322. pmid:11519935
  11. 11. Williams JE, Satterwhite RC, Best DL. Pancultural gender stereotypes revisited: The five factor model. Sex roles. 1999;40(7–8):513–525.
  12. 12. Rubinstein G. The big five among male and female students of different faculties. Personality and Individual Differences. 2005;38(7):1495–1503.
  13. 13. Schmitt DP, Allik J, McCrae RR, Benet-Martínez V. The geographic distribution of Big Five personality traits patterns and profiles of human self-description across 56 nations. Journal of cross-cultural psychology. 2007;38(2):173–212.
  14. 14. Feingold A. Gender differences in personality: a meta-analysis. Psychological bulletin. 1994;116(3):429. pmid:7809307
  15. 15. Rosenberg M. Society and the adolescent self-image. Princeton university press Princeton, NJ; 1965.
  16. 16. Kling KC, Hyde JS, Showers CJ, Buswell BN. Gender differences in self-esteem: a meta-analysis. Psychological bulletin. 1999;125(4):470. pmid:10414226
  17. 17. Twenge JM, Campbell WK. “Isn’t it fun to get the respect that we’re going to deserve?” Narcissism, social rejection, and aggression. Personality and Social Psychology Bulletin. 2003;29(2):261–272. pmid:15272953
  18. 18. Emmons RA. Factor analysis and construct validity of the narcissistic personality inventory. Journal of personality assessment. 1984;48(3):291–300. pmid:16367528
  19. 19. Back MD, Küfner AC, Dufner M, Gerlach TM, Rauthmann JF, Denissen JJ. Narcissistic admiration and rivalry: Disentangling the bright and dark sides of narcissism. Journal of Personality and Social Psychology. 2013;105(6):1013. pmid:24128186
  20. 20. Grijalva E, Newman DA, Tay L, Donnellan MB, Harms PD, Robins RW, et al. Gender differences in narcissism: A meta-analytic review. Psychological bulletin. 2015;141(2):261. pmid:25546498
  21. 21. Foster JD, Campbell WK, Twenge JM. Individual differences in narcissism: Inflated self-views across the lifespan and around the world. Journal of Research in Personality. 2003;37(6):469–486.
  22. 22. Bushman BJ, Baumeister RF. Threatened egotism, narcissism, self-esteem, and direct and displaced aggression: does self-love or self-hate lead to violence? Journal of personality and social psychology. 1998;75(1):219. pmid:9686460
  23. 23. Farwell L, Wohlwend-Lloyd R. Narcissistic Processes: Optimistic Expectations, Favorable Self-Evaluations, and Self-Enhancing Attributions. Journal of personality. 1998;66(1):65–83. pmid:9457770
  24. 24. Clark EJ, Rieker PP. Gender differences in relationships and stress of medical and law students. Academic Medicine. 1986;61(1):32–40.
  25. 25. Dahlin M, Joneborg N, Runeson B. Stress and depression among medical students: A cross-sectional study. Medical education. 2005;39(6):594–604. pmid:15910436
  26. 26. Misra R, McKean M, West S, Russo T. Academic stress of college students: Comparison of student and faculty perceptions. College Student Journal. 2000;34(2):236–246.
  27. 27. Misra R, McKean M. College students’ academic stress and its relation to their anxiety, time management, and leisure satisfaction. American Journal of Health Studies. 2000;16(1):41.
  28. 28. Cohen S, Kamarck T, Mermelstein R. A global measure of perceived stress. Journal of health and social behavior. 1983; p. 385–396. pmid:6668417
  29. 29. Rotter JB. Generalized expectancies for internal versus external control of reinforcement. Psychological monographs: General and applied. 1966;80(1):1.
  30. 30. Schultz D, Schultz S. Theories of personality. Cengage Learning; 2016.
  31. 31. De Brabander B, Boone C. Sex differences in perceived locus of control. The Journal of social psychology. 1990;130(2):271–272. pmid:2342352
  32. 32. Lefcourt HM. Research with the locus of control construct: extensions and limitations. Elsevier; 2013.
  33. 33. Strickland BR, Haley WE. Sex differences on the Rotter I–E scale. Journal of Personality and Social Psychology. 1980;39(5):930.
  34. 34. Pavot W, Diener E. Review of the satisfaction with life scale. Psychological Assessment. 1993;5(2):164.
  35. 35. Diener E, Emmons RA, Larsen RJ, Griffin S. The satisfaction with life scale. Journal of Personality Assessment. 1985;49(1):71–75. pmid:16367493
  36. 36. Neto F. The satisfaction with life scale: Psychometrics properties in an adolescent sample. Journal of youth and adolescence. 1993;22(2):125–134.
  37. 37. Pavot W, Diener E, Colvin CR, Sandvik E. Further validation of the satisfaction with life scale: Evidence for the cross-method convergence of well-being measures. Journal of personality assessment. 1991;57(1):149–161. pmid:1920028
  38. 38. Russell DW. UCLA Loneliness Scale (Version 3): Reliability, validity, and factor structure. Journal of Personality Assessment. 1996;66(1):20–40. pmid:8576833
  39. 39. Wiseman H, Guttfreund DG, Lurie I. Gender differences in loneliness and depression of university students seeking counselling. British Journal of Guidance and Counselling. 1995;23(2):231–243.
  40. 40. Borys S, Perlman D. Gender differences in loneliness. Personality and Social Psychology Bulletin. 1985;11(1):63–74.
  41. 41. Cramer KM, Neyedley KA. Sex differences in loneliness: The role of masculinity and femininity. Sex Roles. 1998;38(7–8):645–653.
  42. 42. Law R. Beyond ‘women and transport’: towards new geographies of gender and daily mobility. Progress in Human Geography. 1999;23(4):567–588.
  43. 43. Polk M. Swedish men and women’s mobility patterns: Issues of social equity and ecological sustainability. In: Proceedings from the Second national conference on women’s travel issues, Baltimore, October; 1996. p. 23–26.
  44. 44. Rosenbloom S. Trends in women’s travel patterns. In: Women’s Travel Issues Second National Conference; 2000.
  45. 45. Camstra R. Commuting and gender in a lifestyle perspective. Urban studies. 1996;33(2):283–300.
  46. 46. Hanson S, Hanson P. Gender and urban activity patterns in Uppsala, Sweden. Geographical Review. 1980; p. 291–299.
  47. 47. Kang C, Gao S, Lin X, Xiao Y, Yuan Y, Liu Y, et al. Analyzing and geo-visualizing individual human mobility patterns using mobile call records. In: Geoinformatics, 2010 18th International Conference on. IEEE; 2010. p. 1–7.
  48. 48. Song C, Qu Z, Blumm N, Barabási AL. Limits of Predictability in Human Mobility. Science. 2010;327(5968):1018–1021. pmid:20167789
  49. 49. Kwan MP. Interactive geovisualization of activity-travel patterns using three-dimensional geographical information systems: a methodological exploration with a large data set. Transportation Research Part C: Emerging Technologies. 2000;8(1):185–203.
  50. 50. Sekara V, Stopczynski A, Lehmann S. The fundamental structure of dynamic social networks. arxiv:150604704. 2015;.
  51. 51. Cuttone A, Lehmann S, Larsen JE. Inferring Human Mobility from Sparse Low Accuracy Mobile Sensing Data. In: Proceedings of the 2014 ACM International Joint Conference on Pervasive and Ubiquitous Computing: Adjunct Publication. UbiComp’14 Adjunct. New York, NY, USA: ACM; 2014. p. 995–1004.
  52. 52. Wind DK, Sapiezynski P, Furman MA, Lehmann S. Inferring stop-locations from wifi. PloS one. 2016;11(2):e0149105. pmid:26901663
  53. 53. Caldwell MA, Peplau LA. Sex differences in same-sex friendship. Sex Roles. 1982;8(7):721–732.
  54. 54. Dunbar RI, Spoors M. Social networks, support cliques, and kinship. Human Nature. 1995;6(3):273–290. pmid:24203093
  55. 55. Miritello G, Lara R, Cebrian M, Moro E. Limited communication capacity unveils strategies for human interaction. Scientific reports. 2013;3. pmid:23739519
  56. 56. Szell M, Thurner S. How women organize social networks different from men. Scientific reports. 2013;3. pmid:23393616
  57. 57. Acar A. Antecedents and consequences of online social networking behavior: The case of Facebook. Journal of website promotion. 2008;3(1–2):62–83.
  58. 58. Sheldon P. The relationship between unwillingness-to-communicate and students’ Facebook use. Journal of Media Psychology. 2008;20(2):67–75.
  59. 59. Thelwall M. Social networks, gender, and friending: An analysis of MySpace member profiles. Journal of the American Society for Information Science and Technology. 2008;59(8):1321–1330.
  60. 60. Lewis K, Kaufman J, Gonzalez M, Wimmer A, Christakis N. Tastes, ties, and time: A new social network dataset using Facebook. com. Social networks. 2008;30(4):330–342.
  61. 61. Brass DJ. Men’s and women’s networks: A study of interaction patterns and influence in an organization. Academy of Management journal. 1985;28(2):327–343.
  62. 62. Onnela JP, Waber BN, Pentland A, Schnorf S, Lazer D. Using sociometers to quantify social interaction patterns. Scientific Reports. 2014;4.
  63. 63. Laniado D, Volkovich Y, Kappler K, Kaltenbrunner A. Gender homophily in online dyadic and triadic relationships. EPJ Data Science. 2016;5(1):1–23.
  64. 64. McPherson M, Smith-Lovin L, Cook JM. Birds of a feather: Homophily in social networks. Annual review of sociology. 2001; p. 415–444.
  65. 65. Palchykov V, Kaski K, Kertész J, Barabási AL, Dunbar RI. Sex differences in intimate relationships. Scientific reports. 2012;2. pmid:22518274
  66. 66. Stoicaa A, Smoredaa Z, Prieurb C, Guillaumec JL. Age, gender and communication networks. 2010;.
  67. 67. Leskovec J, Horvitz E. Planetary-scale views on a large instant-messaging network. In: Proceedings of the 17th international conference on World Wide Web. ACM; 2008. p. 915–924.
  68. 68. Wasserman S, Faust K. Social network analysis: Methods and applications. vol. 8. Cambridge university press; 1994.
  69. 69. Kossinets G, Watts DJ. Origins of homophily in an evolving social network1. American journal of sociology. 2009;115(2):405–450.
  70. 70. Goodreau SM, Kitts JA, Morris M. Birds of a feather, or friend of a friend? using exponential random graph models to investigate adolescent social networks*. Demography. 2009;46(1):103–125. pmid:19348111
  71. 71. Kovanen L, Kaski K, Kertész J, Saramäki J. Temporal motifs reveal homophily, gender-specific patterns, and group talk in call sequences. Proceedings of the National Academy of Sciences. 2013;110(45):18070–18075.
  72. 72. Rakow L. Gender on the line: Women, the telephone, and community life. University of Illinois Press; 1992.
  73. 73. Smoreda Z, Licoppe C. Gender-specific use of the domestic telephone. Social psychology quarterly. 2000; p. 238–252.
  74. 74. Aledavood T, López E, Roberts SG, Reed-Tsochas F, Moro E, Dunbar RI, et al. Daily rhythms in mobile telephone communication. PloS one. 2015;10(9):e0138098. pmid:26390215
  75. 75. Sekara V, Lehmann S. The Strength of Friendship Ties in Proximity Sensor Data. PLoS ONE. 2014;9(7):1–8.
  76. 76. Onnela JP, Saramäki J, Hyvönen J, Szabó G, Lazer D, Kaski K, et al. Structure and tie strengths in mobile communication networks. Proceedings of the National Academy of Sciences. 2007;104(18):7332–7336.
  77. 77. Eagle N, Pentland AS, Lazer D. Inferring friendship network structure by using mobile phone data. Proceedings of the National Academy of Sciences. 2009;106(36):15274–15278.
  78. 78. Wilson RE, Gosling SD, Graham LT. A Review of Facebook Research in the Social Sciences. Perspectives on Psychological Science. 2012;7(3):203–220. pmid:26168459
  79. 79. Wang SS, Stefanone MA. Showing Off? Human Mobility and the Interplay of Traits, Self-Dosclosure, and Facebook Check-Ins. Social Science Computer Review. 2013;31(4):437–457.
  80. 80. Fano RM, Hawkins D. Transmission of information: A statistical theory of communications. American Journal of Physics. 1961;29(11):793–794.
  81. 81. Frias-Martinez V, Frias-Martinez E, Oliver N. A Gender-Centric Analysis of Calling Behavior in a Developing Economy Using Call Detail Records. In: AAAI spring symposium: artificial intelligence for development; 2010.
  82. 82. Felbo B, Sundsøy P, Pentland A, Lehmann S, de Montjoye YA. Using Deep Learning to Predict Demographics from Mobile Phone Metadata. arXiv preprint arXiv:151106660. 2015;.
  83. 83. de Montjoye YA, Quoidbach J, Robic F, Pentland AS. Predicting personality using novel mobile phone-based metrics. In: International Conference on Social Computing, Behavioral-Cultural Modeling, and Prediction. Springer; 2013. p. 48–55.
  84. 84. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: Machine Learning in Python. Journal of Machine Learning Research. 2011;12:2825–2830.
  85. 85. Reeder HM. The effect of gender role orientation on same-and cross-sex friendship formation. Sex Roles. 2003;49(3–4):143–152.
  86. 86. Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, et al. Missing value estimation methods for DNA microarrays. Bioinformatics. 2001;17(6):520–525. pmid:11395428
  87. 87. Goldberg LR. Language and individual differences: the search for universals in personality lexicons. Review of Personality and Social Psychology. 1981;2:141–165.
  88. 88. Digman JM. Personality structure: emergence of the five-factor model. Annu Rev Psychol. 1990;41:417–440.
  89. 89. Russel D, Peplau LA, Ferguson ML. Developing a measure of loneliness. Journal of Personality Assessment. 1978;42(3):290–294.