# Constructed data: 16 WT and 16 KO mice, 4 cages per genotype, 4 mice per cage. The knockout raises taxon 1 four-fold;
# no other taxon is changed in absolute terms. Python 3.13, NumPy, pandas. Seed 1, the first one tried.
import numpy as np, pandas as pd
from model import study, NTAXA
rng = np.random.default_rng(1)
rows = study(rng)
d = pd.DataFrame([(m, g, c, *x) for m, g, c, x in rows], columns=["mouse", "genotype", "cage"] + [f"taxon_{i + 1:02d}" for i in range(NTAXA)])
d.to_csv("counts.csv", index=False)