# Constructed data: 5 WT and 5 KO mice, 40 neurons each, 20 stimulus trials per neuron. By design WT and KO have the same
# average response (ΔF/F 0.05). 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)
rows = [(m, g, n, t + 1, round(float(v), 4)) for m, g, n, r in study(rng) for t, v in enumerate(r)]
pd.DataFrame(rows, columns=["mouse", "genotype", "neuron", "trial", "dff"]).to_csv("responses.csv", index=False)