How to read this notebook

This notebook reproduces, in a self-contained and runnable form, the experiment behind the thesis "Factors Influencing the Evolution of Culture in Population Encounters"

The structure mirrors the scientific argument:

  1. Experiment design - what is simulated, the metrics, and the experiment plan.
  2. Hypotheses - the four predictions tested.
  3. Analysis - one section per hypothesis: the result, the figure, and a short interpretation.

Experiment design

Evaluates the influence of population size, network connectivity, and training set overlap on assimilation or hybridization of behavior and knowledge

Date: 2026-03-20 (Edoardo Cecchini (UGA))

Hypothesis: Larger population assimilates smaller ones; lower connectivity leads to hybridization; larger training set overlap increazes hybridization

10 runs; 150000 games

Experimental setting: Two populations stabilize independently on the same environment and then interact together

The question

This experiment asks: what happens when two populations, each having independently developed its own culture, meet and have to live together? Two outcomes are conceivable:

  • Assimilation - one population's knowledge is adopted by the other.
  • Hybridization - the two combine into something new that belongs to neither original culture.

The goal is to find out which of these occurs, and what influence it.

The model in brief

Agents represent the world with a formal ontology and use it to make a decision about each object they encounter. Through repeated interaction they adapt: whenever two agents disagree, the less successful one adjusts its ontology toward the more successful one (failure-driven adaptation). A population that has reached internal agreement on every object is said to have a stable culture.

The experiment extends this baseline with a two-phase protocol:

  • Phase 1 — Stabilization. Two populations each stabilize in isolation, developing distinct cultures. Their decision functions ($h_1$, $h_2$) and ontologies sets ($\mathcal{O}_1$, $\mathcal{O}_2$) are recorded as baselines.
  • Phase 2 — Encounter. The two populations merge into one and interact for a fixed budget of iterations. The merged decision function ($h'$) and ontologies ($\mathcal{O}'$) are recorded at the end.

Variables

Three variables are manipulated to see what shapes the outcome:

  • Population distribution ($\alpha$) - the relative size of the two populations (population 1 will be the minority in unbalanced configurations).
  • Network connectivity ($\eta$) - how densely the two groups are allowed to interact.
  • Environment overlap ($\omega$) - how much of the world the two populations shared during training (their prior "common ground").

The four hypotheses

  • H1 - Stabilization. The encounter drives the merged population to behavioral consensus (final behavioral consensus → 1).
  • H2 - Population asymmetry. The larger population assimilates the smaller one.
  • H3 - Network connectivity. Lower connectivity:
    • H3.1: Reduces majority assimilation / increases hybridization.
    • H3.2: Yields a more semantically diverse final culture.
  • H4 - Environment overlap. Higher overlap:
    • H4.1: Shifts the outcome toward hybridization.
    • H4.2: Lowers the final semantic diversity.

Experiment plan

Fixed parameters:

  • 2 populations totalizing
  • 30 agents,
  • 10 runs per configuration,
  • 150000 encounter iterations.

Variable grid:

  • population distribution $\alpha \in$ 0.12,0.25,0.5,
  • network connectivity $\eta \in$ 0.2,0.5,1.0,
  • environment overlap $\omega \in$ 0.0,0.33,0.66,1.0

Explored as a full factorial design.

Measuring cultural change: two complementary distances

Because an agent's knowledge lives on two levels, the outcome is measured with two distances:

  • Decision distance $d_d(h_x, h_y)$ — the behavioral level: the fraction of objects on which two cultures decide differently. It answers "do they act alike?"
  • Ontology distance $d_o(\mathcal{O}_x, \mathcal{O}_y)$ — the semantic level: the average structural dissimilarity between the agents' ontologies (average-linkage over all agent pairs). It answers "do they think alike?"

Turning distances into outcomes

Geometric classification scheme. Let $A$ and $B$ be the two pre-encounter cultures and $C$ the merged one. Place a sphere of radius $\rho = \tfrac{d(A,B)}{2}$ on each of $A$ and $B$. Then:

  • $C$ inside $A$'s sphere only → assimilation by Population 1
  • $C$ inside $B$'s sphere only → assimilation by Population 2
  • $C$ outside both → hybridization

This is applied both behaviorally (using $d_d$ with $h_1,h_2,h'$) and semantically (using $d_o$ with the ontology sets). The scheme is intuitive but discrete; where it proves too coarse, the analysis falls back on a continuous distance-difference metric introduced where it is first needed.

Experiment

Performed 2026-05-04 by Edoardo Cecchini (UGA).

Lazy lavender hash: 2bf30aa31d3ce78a4ff1f39b665588093fec056c

Parameter file: params.sh

Executed command (script.sh):

#!/bin/bash

. params.sh

EOUTPUT=${OUTPUT}${LABEL}
mkdir -p ${EOUTPUT}

SVDIR="${EOUTPUT}/end-metrics"
git --git-dir=lazylav/.git rev-parse --verify HEAD > ${EOUTPUT}/lazylav.txt

java -Dlog.level=INFO -cp ${JPATH}:logback.xml \
    fr.inria.exmo.lazylavender.engine.ExperimentalPlan \
    ${OPT} \
    -DsaveDir=${SVDIR} \
#    -Dexperiment=fr.inria.exmo.lazylavender.expe.PopulationEncounter \
    -DresultDir=${EOUTPUT} \
#    -DrevisionModality=addjoin

Raw results

Full results are available on Zenodo: DOI:10.5281/zenodo.20735019

Analysis

The remaining sections load the run's outputs and then work through the four hypotheses, one at a time. Each section states the prediction, shows the result, and interprets it.

First, the end-of-run metrics and the post-encounter interaction logs are loaded, and each run is labelled by its $(\alpha, \eta, \omega)$ configuration.

Each run is then classified independently under the geometric scheme — once on the behavioral distance and once on the ontology distance. Classifying per run (rather than averaging distances first) matters: the scheme's radius depends on the pre-encounter distance of that run, so averaging across runs would produce a threshold never actually observed.

Hypothesis 1 - The encounter leads to stabilization

Prediction. As long as the two groups can interact, repairing disagreements should push the merged population to a single behavioral consensus, the cumulative success rate should climb toward 1.

The success rate is cumulative ($\tfrac{\text{successful games}}{\text{games played}}$), so it trends upward as the population converges.

No description has been provided for this image

Because the success rate carries the weight of every early failed game, it approaches but never quite reaches 1 within the fixed budget. A cleaner check is the final consensus: the average share of agents agreeing on the majority decision for each object at the end of the run.

Average Final Consensus across all runs: 0.9996296296296296
Found 9 runs where final consensus is not 1.0:
 network_connectivity  population_distribution  environment_overlap  run  finalConsensus
                  0.5                     0.12                 0.66    4        0.989583
                  0.5                     0.25                 0.66    3        0.975000
                  0.5                     0.25                 0.66   10        0.983333
                  1.0                     0.50                 0.33    6        0.981250
                  0.2                     0.50                 0.00   10        0.997917
                  0.2                     0.25                 0.66    1        0.970833
                  1.0                     0.50                 0.00    8        0.989583
                  1.0                     0.50                 0.00   10        0.989583
                  0.5                     0.50                 0.66   10        0.989583

Finding: H1 is supported. Almost every run reaches total final consensus (mean ≈ 0.9996); the few runs that stop just short reflect the fixed iteration budget meeting the simulator's randomness, not a failure of convergence. The success-rate curve is a good picture of the trajectory; the final-consensus value is the precise verdict.

Hypothesis 2 - The larger population assimilates the smaller

Prediction. With skewed sizes, minority agents are drawn into cross-cultural interactions more often and have fewer peers to reinforce their knowledge, so the majority's behavior should prevail.

Below, runs are grouped by population distribution $\alpha$ and counted by outcome (behavioral, then semantic), with a $\chi^2$ test and Cramér's V effect size. As stated before, by convention, for $\alpha < 0.5$ Population 1 is the minority.

Decision Outcome (%)
decision_outcome assimilation_pop1 assimilation_pop2 hybridization
population_distribution
0.12 0.0 99.2 0.8
0.25 0.8 89.2 10.0
0.50 28.3 40.0 31.7
Ontology Outcome (%)
ontology_outcome assimilation_pop2 hybridization
population_distribution
0.12 12.5 87.5
0.25 1.7 98.3
0.50 0.8 99.2
Chi-square (decision outcome vs population_distribution):
chi2=138.2698, p=6.6228e-29, dof=4, Cramér's V=0.4382

Chi-square (ontology outcome vs population_distribution):
chi2=21.4035, p=2.2505e-05, dof=2, Cramér's V=0.2438
No description has been provided for this image

Reading the tables. Behaviorally, assimilation by the majority (Population 2) rises from a near-even split at $\alpha=0.50$ to near-total dominance ($\approx 99\%$) at $\alpha=0.12$, a highly significant shift with a medium-high effect size. The semantic tables are dominated by "hybridization" almost everywhere: this is an artifact, not a finding. Ontologies in this model only ever grow (classes split, never merge), so post-encounter ontologies sit far from both baselines and land outside both spheres by construction. This is exactly the limitation that motivates the continuous metric below.

A continuous view: the distance difference

To get past the three-bin scheme, define the normalized distance difference

$$\Delta d = \frac{d(A,C) - d(B,C)}{d(A,B)}$$

so that:

  • $\Delta d > 0$ means the final culture leans toward Population 2.
  • $\Delta d < 0$ means the final culture leans toward Population 1.
  • $\Delta d \to 0$ means a balanced (hybrid) outcome.

Values near $\pm 1$ mean near-complete dominance.

Normalising by $d(A,B)$ makes conditions comparable.

Kruskal-Wallis on dd_diff_strength across population distributions:
H=138.1268, p=1.0142e-30
popDist=0.12: mean dd_diff_strength = 0.8588
popDist=0.25: mean dd_diff_strength = 0.6308
popDist=0.5: mean dd_diff_strength = 0.0451
Eta-squared: 0.3923
Kruskal-Wallis on do_diff_strength across population distributions:
H=203.9578, p=5.1420e-45
popDist=0.12: mean do_diff_strength = 0.3464
popDist=0.25: mean do_diff_strength = 0.1956
popDist=0.5: mean do_diff_strength = 0.0207
Eta-squared: 0.5657
No description has been provided for this image
No description has been provided for this image

Finding: H2 is supported. Assimilation strength toward the majority decreases monotonically as the minority grows, on both the behavioral ($\Delta d_d$) and semantic ($\Delta d_o$) dimensions, with large effect sizes. The continuous metric adds what the bins could not: the change from $\alpha=0.12$ to $0.25$ is not just fewer assimilation runs but weaker assimilation in each. At $\alpha=0.50$ the distribution is centred on zero, confirming the slight "Population 2 advantage" seen in the table was stochastic drift, not a real tendency. Notably the semantic effect is the stronger of the two: because ontologies only grow, every agent keeps its origin's structural fingerprint, anchoring the merged ontology toward whichever population is more numerous.

Hypothesis 3 - Network connectivity and information bottlenecks

Two claims are tested: lower connectivity should (H3.1) reduce majority assimilation / increase hybridization, and (H3.2) leave a more semantically diverse final culture.

H3.1 - Connectivity and assimilation

Since H2 showed population size dominates, connectivity is first isolated in the balanced case ($\alpha = 0.5$), where the size effect is neutralised.

Decision Outcome (%)
decision_outcome assimilation_pop1 assimilation_pop2 hybridization
network_connectivity
0.2 27.5 42.5 30.0
0.5 32.5 37.5 30.0
1.0 25.0 40.0 35.0
Ontology Outcome (%)
ontology_outcome assimilation_pop2 hybridization
network_connectivity
0.2 0.0 100.0
0.5 0.0 100.0
1.0 2.5 97.5
Chi-square (decision outcome vs network_connectivity):
chi2=0.7473, p=9.4537e-01, dof=4, Cramér's V=0.0558

Chi-square (ontology outcome vs network_connectivity):
chi2=2.4561, p=2.9286e-01, dof=2, Cramér's V=0.0826
Kruskal-Wallis on dd_diff_strength across network connectivity (symmetric cases):
H=0.5066, p=7.7625e-01
netCon=0.2: mean dd_diff_strength = 0.1043
netCon=0.5: mean dd_diff_strength = 0.0073
netCon=1.0: mean dd_diff_strength = 0.0260
Eta-squared: -0.0138

Kruskal-Wallis on do_diff_strength across network connectivity (symmetric cases):
H=0.4311, p=8.0611e-01
netCon=0.2: mean do_diff_strength = 0.0257
netCon=0.5: mean do_diff_strength = 0.0277
netCon=1.0: mean do_diff_strength = 0.0085
Eta-squared: -0.0134

In the balanced case there is no effect of connectivity. The hypothesis expects a stronger effect in the unbalanced configurations, examined next.

Decision Outcome (%)
decision_outcome assimilation_pop1 assimilation_pop2 hybridization
network_connectivity
0.2 1.2 90.0 8.8
0.5 0.0 96.2 3.8
1.0 0.0 96.2 3.8
Ontology Outcome (%)
ontology_outcome assimilation_pop2 hybridization
network_connectivity
0.2 5.0 95.0
0.5 11.2 88.8
1.0 5.0 95.0
Chi-square (decision outcome vs network_connectivity):
chi2=4.6828, p=3.2142e-01, dof=4, Cramér's V=0.0988

Chi-square (ontology outcome vs network_connectivity):
chi2=2.4561, p=2.9286e-01, dof=2, Cramér's V=0.0826
Kruskal-Wallis on dd_diff_strength across network connectivity (asymmetric cases):
H=0.5760, p=7.4978e-01
netCon=0.2: mean dd_diff_strength = 0.7093
netCon=0.5: mean dd_diff_strength = 0.7686
netCon=1.0: mean dd_diff_strength = 0.7582
Eta-squared: -0.0060

Kruskal-Wallis on do_diff_strength across network connectivity (asymmetric cases):
H=1.5657, p=4.5709e-01
netCon=0.2: mean do_diff_strength = 0.2569
netCon=0.5: mean do_diff_strength = 0.2881
netCon=1.0: mean do_diff_strength = 0.2679
Eta-squared: -0.0018

The heatmaps below give the clearest view: mean distance difference across every combination of $\alpha$ (rows) and $\eta$ (columns), behaviorally and then semantically.

No description has been provided for this image

The interesting metric here is the ontology distance. Behavioral convergence is already guaranteed by H1 (everyone ends up deciding like the majority), so a behavioral measure can only restate that. The semantic measure can in principle reveal whether restricted connectivity let the minority's structure survive even after its decisions aligned.

No description has been provided for this image
Kruskal-Wallis H-statistic: 1.0121
P-value: 6.0287e-01

Finding: H3.1 is not supported. The heatmaps show clean horizontal bands: colour barely changes along a row (varying $\eta$) but shifts sharply between rows (varying $\alpha$). Population size drives the outcome; connectivity does not, in either dimension, at any size. A connectivity signal hinted at in an earlier 5-run pilot did not survive the larger 10-run design. The reason is temporal: a sparse network may slow convergence, but given a long-enough budget the pressure to agree (H1) overrides any buffering, and because ontologies only grow the minority ends up accumulating the same adaptations regardless of topology.

H3.2 - Connectivity and final semantic diversity

Here "diversity" is the internal heterogeneity of the merged population: the average pairwise ontology distance among all its agents at the end of the run.

Kruskal-Wallis H-statistic: 3.8200
P-value: 1.4808e-01
netCon=0.2: mean final distance = 0.4894
netCon=0.5: mean final distance = 0.4760
netCon=1.0: mean final distance = 0.4800

Finding: H3.2 is not supported. Mean final diversity is essentially flat across connectivity levels (≈ 0.476–0.489) and the difference is not significant. This null result runs counter to the single-population literature, where network fragmentation preserves diversity; here the hard constraint to reach agreement appears strong enough to wash that effect out. It may also signal that the decision task is simply not complex enough for connectivity to matter — a candidate for future work.

Hypothesis 4 - Environment overlap and cultural dynamics

Two claims: higher overlap should (H4.1) shift outcomes toward hybridization, and (H4.2) lower the final semantic diversity. Both are analysed in the balanced case ($\alpha = 0.5$).

H4.1 - Overlap and assimilation direction

Decision Outcome (%)
decision_outcome assimilation_pop1 assimilation_pop2 hybridization
environment_overlap
0.00 30.0 43.3 26.7
0.33 26.7 46.7 26.7
0.66 26.7 50.0 23.3
1.00 30.0 20.0 50.0
Ontology Outcome (%)
ontology_outcome assimilation_pop2 hybridization
environment_overlap
0.00 3.3 96.7
0.33 0.0 100.0
0.66 0.0 100.0
1.00 0.0 100.0
Chi-square (decision outcome vs environment_overlap):
chi2=8.6001, p=1.9735e-01, dof=6, Cramér's V=0.1893

Chi-square (ontology outcome vs environment_overlap):
chi2=3.0252, p=3.8775e-01, dof=3, Cramér's V=0.1588

There is a visible jump in hybridization at full overlap ($\omega = 1$); the extreme conditions (0.0 vs 1.0) are compared directly to test it.

chi2=4.7094, p=9.4923e-02, dof=2, Cramér's V=0.2802
No description has been provided for this image

The continuous metric clarifies what is happening at $\omega = 1$.

Kruskal-Wallis on dd_diff_strength across overlap values (popDist=0.5):
H=3.7845, p=2.8569e-01
overlap=0.0: mean dd_diff_strength = 0.0285
overlap=0.33: mean dd_diff_strength = 0.0485
overlap=0.66: mean dd_diff_strength = 0.1956
overlap=1.0: mean dd_diff_strength = -0.1508
Kruskal-Wallis on do_diff_strength across overlap values (popDist=0.5):
H=4.0623, p=2.5482e-01
overlap=0.0: mean do_diff_strength = 0.0548
overlap=0.33: mean do_diff_strength = 0.0094
overlap=0.66: mean do_diff_strength = 0.0170
overlap=1.0: mean do_diff_strength = 0.0014
No description has been provided for this image

The behavioral band explodes at high overlap while the semantic band stays narrow.

This is a denominator artifact: with full overlap the two populations are already almost identical before they meet, so the pre-encounter distance $d_d(h_1,h_2)$ shrinks toward zero and the normalised metric amplifies any tiny drift. The raw (non-normalised) difference confirms this, its spread shrinks as overlap grows, the exact opposite of the normalised version.

Kruskal-Wallis (popDist=0.5): H=79.8080, p=3.3746e-17
Eta-squared: 0.6621
overlap=0.0: dd_h1_h2 mean = 0.4813
overlap=0.33: dd_h1_h2 mean = 0.3229
overlap=0.66: dd_h1_h2 mean = 0.1625
overlap=1.0: dd_h1_h2 mean = 0.1042
Kruskal-Wallis (Raw Difference, popDist=0.5): H=4.3347, p=2.2752e-01
Eta-squared: 0.0115

--- Raw Behavioral Distance Difference (dd_diff_nn) ---
Overlap = 0.0: Mean = 0.0063 | Std Dev = 0.1693
Overlap = 0.33: Mean = 0.0187 | Std Dev = 0.1127
Overlap = 0.66: Mean = 0.0333 | Std Dev = 0.0995
Overlap = 1.0: Mean = -0.0250 | Std Dev = 0.0847
No description has been provided for this image

Finding: H4.1 is not supported. Neither the geometric scheme nor the raw distance difference shows any genuine, directional effect of overlap on the outcome at balanced sizes. The apparent spike at $\omega = 1$ is an artifact of the near-zero baseline distance — and conceptually, when two populations are already identical, "hybridization" loses its meaning: the encounter only settles a few leftover disagreements.

H4.2 - Overlap and final semantic diversity

Kruskal-Wallis (popDist=0.5): H=8.0023, p=4.5964e-02
Eta-squared: 0.0431
overlap=0.0: mean = 0.4562
overlap=0.33: mean = 0.4901
overlap=0.66: mean = 0.4917
overlap=1.0: mean = 0.4879

Finding: H4.2 is not supported. Higher overlap did not lower final diversity.

Instead, it was observed that the completely disjoint condition ($\omega = 0$) produced the least diverse final culture (mean ≈ 0.456 vs ≈ 0.49 for all overlapping conditions), a threshold effect separating only the disjoint case from the rest, significant with a small effect size.

The hypothesis assumed: more overlap → more similar before contact → fewer adaptations → less diverse result. The first link does hold — let's verify it before interpreting.

Kruskal-Wallis (popDist=0.5): H=64.3425, p=6.9340e-14
Eta-squared: 0.5288
overlap=0.0: do_O1_O2 mean = 0.7404
overlap=0.33: do_O1_O2 mean = 0.6851
overlap=0.66: do_O1_O2 mean = 0.5572
overlap=1.0: do_O1_O2 mean = 0.5171
No description has been provided for this image

Interpretation. The pre-encounter distance does fall monotonically with overlap (dashed line): more shared training really does make populations more alike before contact. But final diversity (solid line) is lowest exactly where pre-encounter distance was highest, the paradox.

The resolution is in what kind of adaptation occurs, not how much. With disjoint training ($\omega = 0$), the two populations decide differently on a large share of objects, so the encounter forces a deep, system-wide renegotiation: agents rebuild compatible knowledge almost from scratch and converge onto a uniform shared structure, uniform, hence less internally diverse.

With any overlap, they already agree on much, so the encounter only needs targeted touch-ups, and each population largely keeps the structural diversity it built during stabilization.

However, this finding is only marginally significant and should be treated as an observation requiring further investigation to be conclusively verified.

Summary

Hypothesis Prediction Outcome
H1 Encounter → behavioral consensus Supported
H2 Larger population assimilates smaller Supported
H3.1 Lower connectivity → less assimilation Not supported
H3.2 Lower connectivity → more diversity Not supported
H4.1 Higher overlap → more hybridization Not supported
H4.2 Higher overlap → less diversity Not supported

Bottom line. The outcome of a cultural encounter in this model is governed almost entirely by relative population size: the larger group assimilates the smaller, more strongly the more lopsided the sizes. Neither the structure of contact nor the amount of shared prior knowledge had a comparable effect, contrary to what single-population studies would suggest.

Appendix — Version check

A reproducibility safeguard: the simulator version (Lazy Lavender hash) recorded with the run is compared against the version requested in the parameter file.

Experiment performed with requested version: 2bf30aa31d3ce78a4ff1f39b665588093fec056c