# Constructed data: 8 control and 8 treated mice; the treatment lowers the true burden by 1.5 log10 (3.0 to 1.5).
# Python 3.13, NumPy, pandas. Seed 1, the first one tried.
import numpy as np, pandas as pd
from model import study
rng = np.random.default_rng(1)
pd.DataFrame(study(rng, mean_trt=1.5), columns=["group", "mouse", "dilution", "colonies"]).to_csv("cfu.csv", index=False)