Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Stochastic block model

From Wikipedia, the free encyclopedia
Concept in network science
Part ofa series on
Network science
Network types
Graphs
Features
Types
Models
Topology
Dynamics
  • Lists
  • Categories

Thestochasticblock model is agenerative model for randomgraphs. This model tends to produce graphs containingcommunities, subsets of nodes characterized by being connected with one another with particular edge densities. For example, edges may be more common within communities than between communities. Its mathematical formulation was first introduced in 1983 in the field of social network analysis byPaul W. Holland et al.[1] The stochastic block model is important instatistics,machine learning, andnetwork science, where it serves as a useful benchmark for the task of recoveringcommunity structure in graph data.

Definition

[edit]

The stochastic block model takes the following parameters:

The edge set is then sampled at random as follows: any two verticesuCi{\displaystyle u\in C_{i}} andvCj{\displaystyle v\in C_{j}} are connected by an edge with probabilityPij{\displaystyle P_{ij}}. An example problem is: given a graph withn{\displaystyle n} vertices, where the edges are sampled as described, recover the groupsC1,,Cr{\displaystyle C_{1},\ldots ,C_{r}}.

Special cases

[edit]
An example of the assortative case for the stochastic block model.

If the probability matrix is a constant, in the sense thatPij=p{\displaystyle P_{ij}=p} for alli,j{\displaystyle i,j}, then the result is theErdős–Rényi modelG(n,p){\displaystyle G(n,p)}. This case is degenerate—the partition into communities becomes irrelevant—but it illustrates a close relationship to the Erdős–Rényi model.

Theplanted partition model is the special case that the values of the probability matrixP{\displaystyle P} are a constantp{\displaystyle p} on the diagonal and another constantq{\displaystyle q} off the diagonal. Thus two vertices within the same community share an edge with probabilityp{\displaystyle p}, while two vertices in different communities share an edge with probabilityq{\displaystyle q}. Sometimes it is this restricted model that is called the stochastic block model. The case wherep>q{\displaystyle p>q} is called anassortative model, while the casep<q{\displaystyle p<q} is calleddisassortative.

Returning to the general stochastic block model, a model is calledstrongly assortative ifPii>Pjk{\displaystyle P_{ii}>P_{jk}} wheneverjk{\displaystyle j\neq k}: all diagonal entries dominate all off-diagonal entries. A model is calledweakly assortative ifPii>Pij{\displaystyle P_{ii}>P_{ij}} wheneverij{\displaystyle i\neq j}: each diagonal entry is only required to dominate the rest of its own row and column.[2]Disassortative forms of this terminology exist, by reversing all inequalities. For some algorithms, recovery might be easier for block models with assortative or disassortative conditions of this form.[2]

Typical statistical tasks

[edit]

Much of the literature on algorithmic community detection addresses three statistical tasks: detection, partial recovery, and exact recovery.

Detection

[edit]

The goal of detection algorithms is simply to determine, given a sampled graph, whether the graph has latent community structure. More precisely, a graph might be generated, with some known prior probability, from a known stochastic block model, and otherwise from a similarErdos-Renyi model. The algorithmic task is to correctly identify which of these two underlying models generated the graph.[3]

Partial recovery

[edit]

In partial recovery, the goal is to approximately determine the latent partition into communities, in the sense of finding a partition that is correlated with the true partition significantly better than a random guess.[4]

Exact recovery

[edit]

In exact recovery, the goal is to recover the latent partition into communities exactly. The community sizes and probability matrix may be known[5] or unknown.[6]

Statistical lower bounds and threshold behavior

[edit]

Stochastic block models exhibit a sharp threshold effect reminiscent ofpercolation thresholds.[7][3][8] Suppose that we allow the sizen{\displaystyle n} of the graph to grow, keeping the community sizes in fixed proportions. If the probability matrix remains fixed, tasks such as partial and exact recovery become feasible for all non-degenerate parameter settings. However, if we scale down the probability matrix at a suitable rate asn{\displaystyle n} increases, we observe a sharp phase transition: for certain settings of the parameters, it will become possible to achieve recovery with probability tending to 1, whereas on the opposite side of the parameter threshold, the probability of recovery tends to 0 no matter what algorithm is used.

For partial recovery, the appropriate scaling is to takePij=P~ij/n{\displaystyle P_{ij}={\tilde {P}}_{ij}/n} for fixedP~{\displaystyle {\tilde {P}}}, resulting in graphs of constant average degree. In the case of two equal-sized communities, in the assortative planted partition model with probability matrixP=(p~/nq~/nq~/np~/n),{\displaystyle P=\left({\begin{array}{cc}{\tilde {p}}/n&{\tilde {q}}/n\\{\tilde {q}}/n&{\tilde {p}}/n\end{array}}\right),}partial recovery is feasible[4] with probability1o(1){\displaystyle 1-o(1)} whenever(p~q~)2>2(p~+q~){\displaystyle ({\tilde {p}}-{\tilde {q}})^{2}>2({\tilde {p}}+{\tilde {q}})}, whereas anyestimator fails[3] partial recovery with probability1o(1){\displaystyle 1-o(1)} whenever(p~q~)2<2(p~+q~){\displaystyle ({\tilde {p}}-{\tilde {q}})^{2}<2({\tilde {p}}+{\tilde {q}})}.

For exact recovery, the appropriate scaling is to takePij=P~ijlogn/n{\displaystyle P_{ij}={\tilde {P}}_{ij}\log n/n}, resulting in graphs of logarithmic average degree. Here a similar threshold exists: for the assortative planted partition model withr{\displaystyle r} equal-sized communities, the threshold lies atp~q~=r{\displaystyle {\sqrt {\tilde {p}}}-{\sqrt {\tilde {q}}}={\sqrt {r}}}. In fact, the exact recovery threshold is known for the fully general stochastic block model.[5]

Algorithms

[edit]

In principle, exact recovery can be solved in its feasible range usingmaximum likelihood, but this amounts to solving a constrained orregularized cut problem such as minimum bisection that is typicallyNP-complete. Hence, no known efficient algorithms will correctly compute the maximum-likelihood estimate in the worst case.

However, a wide variety of algorithms perform well in the average case, and many high-probability performance guarantees have been proven for algorithms in both the partial and exact recovery settings. Successful algorithms includespectral clustering of the vertices,[9][4][5][10]semidefinite programming,[2][8] forms ofbelief propagation,[7][11] and community detection[12] among others.

Variants

[edit]

Several variants of the model exist. One minor tweak allocates vertices to communities randomly, according to acategorical distribution, rather than in a fixed partition.[5] More significant variants include the degree-corrected stochastic block model,[13] the hierarchical stochastic block model,[14] the geometric block model,[15] censored block model and the mixed-membership block model.[16]

Topic models

[edit]

Stochastic block model have been recognised to be atopic model on bipartite networks.[17] In a network of documents and words, Stochastic block model can identify topics: group of words with a similar meaning.

Extensions to signed graphs

[edit]

Signed graphs allow for both favorable and adverse relationships and serve as a common model choice for various data analysis applications, e.g., correlation clustering. The stochastic block model can be trivially extended to signed graphs by assigning both positive and negative edge weights or equivalently using a difference of adjacency matrices of two stochastic block models.[18]

DARPA/MIT/AWS Graph Challenge: streaming stochastic block partition

[edit]

GraphChallenge[19] encourages community approaches to developing new solutions for analyzing graphs and sparse data derived from social media, sensor feeds, and scientific data to enable relationships between events to be discovered as they unfold in the field. Streaming stochastic block partition is one of the challenges since 2017.[20]Spectral clustering has demonstrated outstanding performance compared to the original and even improved[21]base algorithm, matching its quality of clusters while being multiple orders of magnitude faster.[22][23]

See also

[edit]

References

[edit]
  1. ^Holland, Paul W; Laskey, Kathryn Blackmond; Leinhardt, Samuel (1983). "Stochastic blockmodels: First steps".Social Networks.5 (2):109–137.doi:10.1016/0378-8733(83)90021-7.ISSN 0378-8733.S2CID 34098453.
  2. ^abcAmini, Arash A.;Levina, Elizaveta (June 2014). "On semidefinite relaxations for the block model".arXiv:1406.5647 [cs.LG].
  3. ^abcMossel, Elchanan; Neeman, Joe; Sly, Allan (February 2012). "Stochastic Block Models and Reconstruction".arXiv:1202.1499 [math.PR].
  4. ^abcMassoulie, Laurent (November 2013). "Community detection thresholds and the weak Ramanujan property".arXiv:1311.3085 [cs.SI].
  5. ^abcdAbbe, Emmanuel; Sandon, Colin (March 2015). "Community detection in general stochastic block models: fundamental limits and efficient recovery algorithms".arXiv:1503.00609 [math.PR].
  6. ^Abbe, Emmanuel; Sandon, Colin (June 2015). "Recovering communities in the general stochastic block model without knowing the parameters".arXiv:1506.03729 [math.PR].
  7. ^abDecelle, Aurelien; Krzakala, Florent; Moore, Cristopher;Zdeborová, Lenka (September 2011). "Asymptotic analysis of the stochastic block model for modular networks and its algorithmic applications".Physical Review E.84 (6) 066106.arXiv:1109.3041.Bibcode:2011PhRvE..84f6106D.doi:10.1103/PhysRevE.84.066106.PMID 22304154.S2CID 15788070.
  8. ^abAbbe, Emmanuel; Bandeira, Afonso S.; Hall, Georgina (May 2014). "Exact Recovery in the Stochastic Block Model".arXiv:1405.3267 [cs.SI].
  9. ^Krzakala, Florent; Moore, Cristopher; Mossel, Elchanan; Neeman, Joe; Sly, Allan; Lenka, Lenka; Zhang, Pan (October 2013)."Spectral redemption in clustering sparse networks".Proceedings of the National Academy of Sciences.110 (52):20935–20940.arXiv:1306.5550.Bibcode:2013PNAS..11020935K.doi:10.1073/pnas.1312486110.PMC 3876200.PMID 24277835.
  10. ^Lei, Jing; Rinaldo, Alessandro (February 2015). "Consistency of spectral clustering in stochastic block models".The Annals of Statistics.43 (1):215–237.arXiv:1312.2050.doi:10.1214/14-AOS1274.ISSN 0090-5364.S2CID 88519551.
  11. ^Mossel, Elchanan; Neeman, Joe; Sly, Allan (September 2013). "Belief Propagation, Robust Reconstruction, and Optimal Recovery of Block Models".The Annals of Applied Probability.26 (4):2211–2256.arXiv:1309.1380.Bibcode:2013arXiv1309.1380M.doi:10.1214/15-AAP1145.S2CID 184446.
  12. ^Fathi, Reza (April 2019). "Efficient Distributed Community Detection in the Stochastic Block Model".arXiv:1904.07494 [cs.DC].
  13. ^Karrer, Brian; Newman, Mark E J (2011)."Stochastic blockmodels and community structure in networks".Physical Review E.83 (1) 016107.arXiv:1008.3926.Bibcode:2011PhRvE..83a6107K.doi:10.1103/PhysRevE.83.016107.PMID 21405744.S2CID 9068097.Archived from the original on 2023-02-04. Retrieved2021-06-16.
  14. ^Peixoto, Tiago (2014)."Hierarchical block structures and high-resolution model selection in large networks".Physical Review X.4 (1) 011047.arXiv:1310.4377.Bibcode:2014PhRvX...4a1047P.doi:10.1103/PhysRevX.4.011047.S2CID 5841379.Archived from the original on 2021-06-24. Retrieved2021-06-16.
  15. ^Galhotra, Sainyam; Mazumdar, Arya; Pal, Soumyabrata;Saha, Barna (February 2018). "The Geometric Block Model".AAAI.32.arXiv:1709.05510.doi:10.1609/aaai.v32i1.11905.S2CID 19152144.
  16. ^Airoldi, Edoardo; Blei, David; Feinberg, Stephen; Xing, Eric (May 2007)."Mixed membership stochastic blockmodels".Journal of Machine Learning Research.9:1981–2014.arXiv:0705.4485.Bibcode:2007arXiv0705.4485A.PMC 3119541.PMID 21701698.
  17. ^Martin Gerlach; Tiago Peixoto; Eduardo Altmann (2018)."A network approach to topic models".Science Advances.4 (7) eaaq1360.arXiv:1708.01677.Bibcode:2018SciA....4.1360G.doi:10.1126/sciadv.aaq1360.PMC 6051742.PMID 30035215.
  18. ^Alyson Fox; Geoffrey Sanders; Andrew Knyazev (2018). "Investigation of Spectral Clustering for Signed Graph Matrix Representations".2018 IEEE High Performance extreme Computing Conference (HPEC). pp. 1–7.doi:10.1109/HPEC.2018.8547575.ISBN 978-1-5386-5989-2.OSTI 1476177.S2CID 54443034.
  19. ^[1]Archived 2023-02-04 at theWayback Machine DARPA/MIT/AWS Graph Challenge
  20. ^[2]Archived 2023-02-04 at theWayback Machine DARPA/MIT/AWS Graph Challenge Champions
  21. ^A. J. Uppal; J. Choi; T. B. Rolinger; H. Howie Huang (2021). "Faster Stochastic Block Partition Using Aggressive Initial Merging, Compressed Representation, and Parallelism Control".2021 IEEE High Performance Extreme Computing Conference (HPEC). pp. 1–7.doi:10.1109/HPEC49654.2021.9622836.ISBN 978-1-6654-2369-4.S2CID 244780210.
  22. ^David Zhuzhunashvili; Andrew Knyazev (2017). "Preconditioned spectral clustering for stochastic block partition streaming graph challenge (Preliminary version at arXiv.)".2017 IEEE High Performance Extreme Computing Conference (HPEC). pp. 1–6.arXiv:1708.07481.doi:10.1109/HPEC.2017.8091045.ISBN 978-1-5386-3472-1.S2CID 19781504.
  23. ^Lisa Durbeck; Peter Athanas (2020). "Incremental Streaming Graph Partitioning".2020 IEEE High Performance Extreme Computing Conference (HPEC). pp. 1–8.doi:10.1109/HPEC43674.2020.9286181.ISBN 978-1-7281-9219-2.S2CID 229376193.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Stochastic_block_model&oldid=1333382776"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2026 Movatter.jp