Research ArticleNETWORK SCIENCE

Cumulative effects of triadic closure and homophily in social networks

See allHide authors and affiliations

Science Advances  08 May 2020:
Vol. 6, no. 19, eaax7310
DOI: 10.1126/sciadv.aax7310
  • Fig. 1 Mechanisms of triadic closure and choice homophily.

    (A) A focal node selected uniformly at random (node at the center with green boundary) finds a candidate neighbor by either (B) selecting a node uniformly at random with probability 1 − c or by (D) closing a triangle with probability c. (C and E) For a focal node in group i, the candidate neighbor in group j is accepted with bias probability Sij (where i, j ∈ {a, b}). We parametrize Sij with tunable parameters sa, sb such that Saa = sa, Sbb = sb, Sab = 1 − sa, and Sba = 1 − sb. If the potential edge [dashed line in (B) and (D)] is accepted, an edge of the focal node (selected uniformly at random) is replaced by one between the focal and candidate neighbor nodes. Otherwise no edges are rewired. (F) Probability (T2)aa of choosing a candidate neighbor with triadic closure from the same group as the focal node (a) as a function of observed homophily oa = Taa for equally sized (na = nb) and equally connected groups (Taa = Tbb and Tab = Tba = 1 − Taa). If the network is not randomly mixed (Taa ≠ 1/2), the probability of triadic closure choosing two nodes of the same group is always larger than the same probability if the selection is done uniformly at random (1/2), implying that triadic closure amplifies existing observed homophily in the network (or suppresses heterophily). However, triadic closure without choice homophily is not enough to maintain the observed homophily in the network, which would make (T2)aa equal to observed homophily (dotted and solid lines cross). An exception is the case Taa = 1, where two completely separate components exist, and triadic closure cannot create edges between them.

  • Fig. 2 Interplay between triadic closure and choice homophily.

    (A) Schematics of the stationary states the network can converge to during the rewiring process (stable fixed points of the model dynamics based on a mean-field analysis; see the Supplementary Materials). Continuous (dashed) lines represent large (low) intra- and intergroup connectivity. The homophily amplification fixed point 0 has many edges within groups and a few between (high Taa and Tbb values). In the core-periphery fixed point + (−), the large (small) group forms the core and attracts edges from the small (large) periphery group: Tab and Tbb (Tba and Ta) are high, while Taa and Tba (Tbb and Tab) are low (see the Supplementary Materials for details on fixed point classification). (B) Observed homophily oa = Taa (fraction of neighbors in the same group of a focal node in group a) at the end of the dynamics as a function of choice homophily s = sa = sb. Group sizes are equal (na = nb), and triadic closure probability c is varied. Mean-field calculations (solid lines correspond to stable fixed points and dashed lines to unstable ones) and numerical simulations (crosses) agree very well (see the Supplementary Materials for systematic analysis). Inset: Case c = 0.9 where both stable/unstable points (continuous/dashed lines) exist. Since na = nb, fixed points +/− are equivalent, but for suitable choice homophily (0.56 < s < 0.70), one group becomes the core and the other periphery depending on the initial network and/or random chance.

  • Fig. 3 Phase diagrams of available fixed points.

    Areas of different shading correspond to parameter ranges in which the various fixed points exist (stable points 0, +, and − or the corresponding unstable points U0, U+, and U−) (see Fig. 2A). Note that unstable fixed points are saddle points with qualitatively similar relative connectivities than their stable variants. (Top) Choice homophily is the same for both groups (s = sa = sb), and fixed points are shown as a function of s and c. (Bottom) Choice homophily for the groups is varied, and fixed points are shown as a function of sa and sb. Here, the triadic closure probability is fixed to c = 0.95. (Left) Groups are of equal size (na = nb). (Right) Groups are of unequal size (na = 0.1). The dashed line corresponds to parameter values in Fig. 2B inset and the dot to Fig. 4A and 4B. Results in the main panel of Fig. 2B are in the white region with a single fixed point (0) (upper left panel).

  • Fig. 4 Model temporal evolution.

    Euclidean distance δ in (Taa, Tbb)-space between current state of the dynamics and fixed points (A) + and (B) U0 as a function of time t (in units of the average number of times an edge is selected for rewiring) for na = nb, c = 0.9, and s = 0.6 (dashed vertical line in Fig. 2B inset). Markers are simulation results, and lines are mean-field solutions. Colors correspond to three different initial network configurations. (B) The dynamics amplifies homophily first by quickly approaching the unstable point U0 and diverting away from it after ∼100 time steps. (A) The blue and orange curves slowly approach the core-periphery stable point +. The green curve approaches the other core-periphery stable point (−).

  • Fig. 5 Homophily amplification in real-world social networks.

    (Top) Estimated homophily amplification A (see Eq. 3) for empirical networks. Lines are mean-field fits as a function of triadic closure c. Points correspond to the estimated mean homophily amplification values using an ABC method, depicted at the estimated mean triadic closure 〈c〉. Edges, nodes, and attributes in each dataset are as follows: Facebook friendships for users by graduation year in the largest of the 100 networks (Facebook), links between political blogs by party (Polblogs), friendships between students by gender (School), Last.fm friendships for users by gender (Last.fm), and links between people in boards of directors in Norway in 2010 (Directors). The smaller group a is denoted in red. (Scatter plots) Estimated homophily amplification Aa and Ab, and mean choice homophily estimates saABC and sbABC (both from ABC method) in all 100 universities of the Facebook dataset, where the two groups i = a, b (na < nb) are graduation years. The largest network is marked with a red cross (see Materials and Methods and Table 1 for dataset details). (Bottom) Estimated distributions for choice homophily saABC and sbABC using the ABC method for the corresponding top panel dataset (colors denote group), where vertical lines depict mean estimates (saABC and sbABC).

  • Table 1 Properties of empirical datasets.

    Properties of empirical datasets.. List of real-world social networks used in this study, their main properties, and estimated model parameters. N is the number of nodes in the network, 〈k〉 is the average degree, na is the fraction of nodes in the smaller attribute group, HIa = (Taana)/(1 − na) and HIb = (Tbbnb)/(1 − nb) are the Coleman homophily indices (60) of the groups, sa(0) and sb(0) are the estimates of the bias parameters when c = 0 [equal to observed homophily, sa(0) = oa and sb(0) = ob], and sa(1) and sb(1) are the estimates when c = 1. saABC and sbABC are the mean ABC estimates for sa and sb.

    abNknaHIaHIbsa(0)/sb(0)sa(1)/sb(1)saABCsbABC
    Facebook2009200814,25142.4947.2%0.7680.7940.889/0.8930.666/0.6600.630.70
    PolblogsLiberalConservative1,22427.3148.0%0.8130.8100.909/0.9020.665/0.6640.8040.861
    SchoolFemaleMale2124.1249.1%0.3080.2130.656/0.6050.595/0.6020.5530.625
    Last.fmMaleFemale188,6725.4931.7%0.1680.1060.621/0.5400.588/0.5500.550.588
    DirectorsFemaleMale2,4126.1235.8%0.192−0.3820.658/0.3210.554/0.3720.5940.36

Supplementary Materials

  • Supplementary Materials

    Cumulative effects of triadic closure and homophily in social networks

    Aili Asikainen, Gerardo Iñiguez, Javier Ureña-Carrión, Kimmo Kaski, Mikko Kivelä

    Download Supplement

    This PDF file includes:

    • Sections S1 to S4
    • Figs. S1 to S10
    • Table S1
    • References

    Files in this Data Supplement:

Stay Connected to Science Advances

Navigate This Article