Nicht-asymptotische Anpassungsgütetests und Modellauswahl in bewerteten stochastischen Blockmodellen

Grundlegende Informationen

• Paper-ID: 2510.13636
• Titel: Non-asymptotic goodness-of-fit tests and model selection in valued stochastic blockmodels
• Autoren: Félix Almendra-Hernández, Miles Bakenhus, Vishesh Karwa, Mitsunori Ogawa, Sonja Petrović
• Klassifizierung: stat.ME cs.SI math.ST stat.TH
• Veröffentlichungsdatum: 16. Oktober 2025
• Paper-Link: https://arxiv.org/abs/2510.13636

Zusammenfassung

Diese Arbeit untersucht Anpassungsgütetests für bewertete stochastische Blockmodelle (valued stochastic blockmodel, SBM). Bewertete SBMs sind eine universelle Modellierungsmethode für Netzwerkdaten, die Knoten in Blöcke gruppiert und Verbindungen zwischen Knoten durch Zählungen oder Etiketten misst. Diese Modellfamilie ermöglicht verschiedene Dyaden-Stichprobenschemata, einschließlich klassischer SBMs, Poisson-SBMs und gekennzeichneter SBMs sowie bestimmter zensierter Kantenbeobachtungen. Der Artikel konzentriert sich auf endliche Stichprobentests für nicht-Bernoulli-SBMs, leitet explizite Markov-Basis-Bewegungen ab, die zur Erzeugung von Referenzverteilungsstichproben erforderlich sind, und definiert Anpassungsgütestatistiken zur Bestimmung der Modellanpassung. Für gekennzeichnete SBMs (einschließlich zensierter Kantenmodelle) wird das asymptotische Verhalten dieser Statistiken untersucht.

Forschungshintergrund und Motivation

Problemdefinition

1. Kernproblem: Wie führt man Anpassungsgütetests für nicht-Bernoulli-bewertete stochastische Blockmodelle durch, besonders in endlichen Stichprobensituationen?
2. Bedeutung:
   • Bei der Netzwerkdatenanalyse ist die Bestimmung, ob die Blockzugehörigkeit von Knoten die Netzwerkbildung beeinflusst, eine Schlüsselfrage
   • Bestehende Methoden konzentrieren sich hauptsächlich auf einfache Graphen (Bernoulli-Dyaden), während reale Netzwerkdaten häufig Zählungen oder mehrere Arten von Verbindungen enthalten
   • Endliche Stichprobentests haben praktischen Wert bei kleinen Datensätzen

Einschränkungen bestehender Methoden

1. Klassische SBM-Einschränkungen: Die meisten bestehenden Rahmenwerke verwenden einfache Graphen und modellieren Dyaden als Bernoulli-Zufallsvariablen
2. Probleme asymptotischer Methoden: Großstichprobenkriterien wie BIC versagen bei Netzwerkmodellen, wenn die Parameterdimension wächst
3. Mangel an theoretischen Garantien: Bestehende Methoden bieten keine theoretischen Garantien für die Nullhypothesenverteilung und asymptotische Macht

Forschungsmotivation

• Erweiterung von Markov-Basis-Methoden aus der algebraischen Statistik auf nicht-Bernoulli-Netzwerkmodelle
• Konstruktion eines teilweise Bayesschen Testrahmenwerks, der Parameterunsicherheit berücksichtigt
• Bereitstellung theoretischer Anleitung für Modellauswahl, besonders für die Wahl der Blockanzahl

Kernbeiträge

1. Theoretische Beiträge:
   • Ableitung expliziter Markov-Basen für Poisson-SBMs und gekennzeichnete SBMs
   • Beweis der Konsistenz von Interpolationsschätzern
   • Etablierung asymptotischer Theorie für Anpassungsgütestatistiken
2. Methodische Beiträge:
   • Vorschlag bedingter Anpassungsgütetests für feste und unbekannte Blockzuweisungen
   • Konstruktion eines teilweise Bayesschen p-Wert-Berechnungsrahmens
   • Entwicklung eines MCMC-basierten Fasersampling-Algorithmus
3. Anwendungsbeiträge:
   • Anwendung der Methode auf die Analyse von Wirt-Parasit-Wechselwirkungen in Ökologischen Netzwerken
   • Validierung der Methoden-Macht und Fehlerrate erster Art auf simulierten Daten
   • Bereitstellung praktischer Richtlinien für Modellauswahl

Methodische Details

Aufgabendefinition

Gegeben ein bewertetes Netzwerk $G = (G_{uv})_{1≤u<v≤n}$, wobei $G_{uv}$ die Verbindungsstärke (Zählung oder Etikett) zwischen dem Knotenpaar $\{u,v\}$ darstellt, ist das Ziel:

1. Zu testen, ob das Netzwerk einem gegebenen bewerteten SBM entspricht
2. Durchführung von Modellanpassungstests, wenn Blockzuweisungen unbekannt sind
3. Auswahl einer geeigneten Blockanzahl $k$

Modellarchitektur

1. Definition bewerteter SBMs

Für $n$ Knoten und $k$ Blöcke nimmt ein bewertetes SBM an:

• Bedingte Unabhängigkeit: $G_{uv} \perp\!\!\!\perp G_{u'v'} | Z$ für beliebige zwei Dyaden
• Exponentialfamilienform: $P_θ(G = g | Z = z) = \prod_{1≤u<v≤n} \frac{h(g_{uv})\exp⟨T_z(g_{uv}), θ_{z_uz_v}⟩}{ψ(θ_{z_uz_v})}$

wobei $h$ das Basismass ist, $T_z$ die suffiziente Statistik, und $θ$ der Parametervektor.

2. Spezialfälle

• Klassisches SBM: $G_{uv} | Z = z \sim \text{Bernoulli}(θ_{z_uz_v})$
• Poisson-SBM: $G_{uv} | Z = z \sim \text{Poisson}(θ_{z_uz_v})$
• Gekennzeichnetes SBM: $G_{uv} | Z = z \sim \text{Multinomial}(N, (θ^{(ℓ)}_{z_uz_v})^L_{ℓ=1})$

3. Markov-Basis-Konstruktion

Markov-Basis für Poisson-SBM: $B = \{ε_{uv} - ε_{u'v'} : z_u = z_{u'}, z_v = z_{v'}\}$

Markov-Basis für gekennzeichnetes SBM: $B = \{ε^{(ℓ)}_{uv} + ε^{(ℓ')}_{u'v'} - ε^{(ℓ')}_{uv} - ε^{(ℓ)}_{u'v'} : ℓ,ℓ' ∈ [L], z_u = z_{u'}, z_v = z_{v'}\}$

Technische Innovationen

1. Bedingter Testrahmen

• Faserdefinition: $F_{z,t} := \{g ∈ G : T_z(g) = t\}$
• Bedingte Verteilung: $P(G = g | T_z(g) = t) = \frac{h(g)}{\sum_{g' ∈ F_{z,t}} h(g')}$
• Exakter p-Wert: $p(z,g) = P(\text{GoF}_z(G) ≥ \text{GoF}_z(g) | T_z(g))$

2. Teilweise Bayessche Methode

Für unbekannte Blockzuweisungen wird der teilweise Bayessche p-Wert definiert als: $p_b(g) = \sum_{z ∈ Z_{n,k}} p(z,g)P(Z = z | g)$

Diese Methode behandelt die Unsicherheit der Blockzuweisung durch Mittelung über die Posteriori-Verteilung.

3. Anpassungsgütestatistiken

Poisson-SBM: $\text{GoF}_z(g) = \sum_{u=1}^n \sum_{i=1}^k \frac{(m_{ui} - n_iθ̂_{z_ui})^2}{n_iθ̂_{z_ui}}$

Gekennzeichnetes SBM: $\text{GoF}_z(g) = χ^2_{BC}(g,z) = \sum_{u=1}^n \sum_{i=1}^k \sum_{ℓ=1}^{L-1} \frac{(m^{(ℓ)}_{ui} - n_iθ̂^{(ℓ)}_{z_ui})^2}{n_iθ̂^{(ℓ)}_{z_ui}}$

Experimentelle Einrichtung

Datensätze

1. Simulierte Daten:
   • Knotenzahl: $n = 50, 100$
   • 4 verschiedene Verbindungsmatrizen $θ^{(1)}, θ^{(2)}, θ^{(3)}, θ^{(4)}$
   • 100 Graphen pro Einstellung generiert
2. Reale Daten:
   • Parasitärer Pilzartenetzwerk (154 Knoten)
   • Baumartennetzwerk (51 Knoten)
   • Kantengewichte stellen die Anzahl gemeinsamer Wirt-/Parasitenarten dar

Bewertungsmetriken

1. Fehlerrate erster Art: Ablehnungsrate der Nullhypothese auf Signifikanzniveau 0,05
2. Testmacht: Modellablehnungsrate bei verschiedenen Blockanzahlen
3. Modellauswahlgenauigkeit: Vergleich mit ICL-Kriterium

Vergleichsmethoden

• ICL-Kriterium (Integrated Complete Likelihood)
• Variationelles EM-Algorithmus für Blockzuweisungsschätzung
• sbm R-Paket-Implementierung

Implementierungsdetails

• MCMC-Kettenlänge: durch numGraphs-Parameter gesteuert
• Blockzuweisungsschätzung: Verwendung des variationellen EM-Algorithmus
• Posteriori-Wahrscheinlichkeitsschwelle: $ε = 1/m$, wobei $m$ die Trägergröße ist

Experimentelle Ergebnisse

Hauptergebnisse

1. Validierung von Macht und Fehlerrate erster Art

Ergebnisse bei $n=50$:

Ergebnisse bei $n=100$:

2. Anwendung auf reale Daten

Baumartennetzwerk:

• Blöcke 3-7: p-Wert = 0
• Blöcke 8-9: p-Wert = 0,01
• Block 10: p-Wert = 0,19
• Blöcke 11-15: p-Wert ≥ 0,68

Pilzartennetzwerk:

• Blöcke 3-17: p-Wert = 0
• Blöcke 18-21: p-Wert = 0,01
• Block 22: p-Wert = 0,07

Experimentelle Erkenntnisse

1. Methodeneffektivität: Ablehnungsraten nahe 1 bei Verwendung von 2 oder 3 Blöcken, nahe 0 bei 4 oder 5 Blöcken, wie erwartet
2. Stichprobengrösseneffekt: Grössere Stichprobenumfänge ($n=100$) bieten stärkere statistische Macht
3. Unterschiede zu bestehenden Methoden:
   • Diese Methode wählt 10 Blöcke für Baumartennetzwerk, 22 für Pilznetzwerk
   • ICL-Kriterium wählt 7 Blöcke für Baumartennetzwerk, 9 für Pilznetzwerk
   • Unterschiede können auf die Konservativität der Methode und strenge Anforderungen an Modellanpassung zurückgeführt werden

Verwandte Arbeiten

Netzwerk-Anpassungsgütetests

1. Spektrale Methoden: Lei (2016) spektrale Anpassungsgütetests
2. ERGM-Methoden: Hunter et al. (2008) Referenzverteilungsvergleichsmethoden
3. Verbesserte Tests: Hu et al. (2021) direkte Behandlung von Rechenkosten und theoretischen Garantien

Stochastische Blockmodelle

1. Klassische SBMs: Holland et al. (1983) latente Blockerweiterung
2. Bewertete Netzwerke: Krivitsky (2012) ERGM-Erweiterung auf Zählnetzwerke
3. Modellauswahl: Wang und Bickel (2017) Likelihood-Modellauswahl

Algebraische Statistik

1. Markov-Basen: Diaconis und Sturmfels (1998) Grundlagentheorie
2. Netzwerkanwendungen: Karwa et al. (2023) endliche Stichprobentests für Bernoulli-SBMs
3. Dynamische Konstruktion: Gross et al. (2014) dynamische Markov-Basis-Methoden

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Theoretischer Beitrag: Erfolgreiche Erweiterung algebraischer Statistikmethoden auf nicht-Bernoulli-Netzwerkmodelle mit strengem theoretischem Fundament
2. Methodeneffektivität: Vorgeschlagene Testmethode zeigt gute statistische Eigenschaften auf simulierten und realen Daten
3. Praktischer Wert: Bietet praktikable Lösung für Modellauswahl in bewerteten Netzwerken

Einschränkungen

1. Rechenkomplexität: MCMC-Methoden können bei grossen Netzwerken auf Konvergenzprobleme stossen
2. Konservativität: Methode kann zu konservativ sein, was zur Auswahl mehr Blöcke führt
3. Blockzuweisungsabhängigkeit: Methode hängt von der Qualität der Blockzuweisungsschätzung ab

Zukünftige Richtungen

1. Zusammengesetzte Markov-Ketten: Entwicklung von Methoden zum Sampling mehrerer Fasern
2. Rechnerische Optimierung: Verbesserung der MCMC-Konvergenz und Recheneffizienz
3. Erweiterte Anwendungen: Integration mit dynamischen und mehrschichtigen Netzwerken

Tiefgreifende Bewertung

Stärken

1. Theoretische Strenge: Vollständiger theoretischer Rahmen mit Konsistenzbeweisen und asymptotischer Analyse
2. Methodische Neuheit: Erste Anwendung von Markov-Basis-Methoden auf nicht-Bernoulli-SBMs
3. Umfassende Experimente: Einschliesslich Machtanalyse, Fehlerrate-erster-Art-Validierung und reale Datenanwendungen
4. Klare Darstellung: Gut strukturiertes Papier mit präzisen technischen Details

Mängel

1. Rechnerische Herausforderungen: Rechenkomplexität kann bei grossen Netzwerken zum Engpass werden
2. Parametersensitivität: Methode ist relativ empfindlich gegenüber der Qualität der Blockzuweisungsschätzung
3. Begrenzte Vergleiche: Vergleiche mit anderen nicht-asymptotischen Methoden sind nicht ausreichend

Auswirkungen

1. Akademischer Wert: Eröffnet neue Richtungen für Schnittstellen zwischen Netzwerkstatistik und algebraischer Statistik
2. Praktischer Wert: Bietet Werkzeuge für Netzwerkanalyse in Ökologie, Sozialwissenschaften und anderen Bereichen
3. Reproduzierbarkeit: Detaillierte Algorithmusbeschreibungen ermöglichen einfache Implementierung und Reproduktion

Anwendungsszenarien

1. Kleine bis mittlere Netzwerke: Methode funktioniert gut bei Netzwerken mit bis zu einigen hundert Knoten
2. Bewertete Netzwerkdaten: Besonders geeignet für Zähl- oder Multi-Label-Netzwerkdaten
3. Strenge statistische Inferenz: Anwendungsszenarien, die präzise statistische Tests erfordern

Literaturverzeichnis

Wichtige Referenzen umfassen:

• Diaconis, P. and Sturmfels, B. (1998). Algebraic algorithms for sampling from conditional distributions
• Holland, P. W., Laskey, K. B., and Leinhardt, S. (1983). Stochastic blockmodels: First steps
• Karwa, V. et al. (2023). Monte Carlo goodness-of-fit tests for degree corrected and related stochastic blockmodels
• Wang, Y. X. R. and Bickel, P. J. (2017). Likelihood-based model selection for stochastic block models