# Constructed data: 8 mice per group; control growth 0.15 per day on the log scale (doubling about 4.6 days),
# drug 0.10 per day (doubling about 6.9 days). The seed (1) was the first one run; it was not selected.
import numpy as np, pandas as pd
from model import simulate, DAYS
rng = np.random.default_rng(1)
rows = []
for grp, rate in [("vehicle", 0.15), ("drug", 0.10)]:
    v = simulate(rng, 8, rate)
    for i in range(8):
        for j, day in enumerate(DAYS):
            if not np.isnan(v[i, j]):
                rows.append((f"{grp[0]}{i+1}", grp, int(day), round(v[i, j], 1)))
pd.DataFrame(rows, columns=["mouse", "group", "day", "volume_mm3"]).to_csv("tumours.csv", index=False)