# Constructed data: 12 heterozygous x heterozygous litters genotyped at weaning; knockout pups survive to weaning at 60%
# of the rate of their littermates. 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, 12, viability=0.6), columns=["litter", "pup", "sex", "genotype"]).to_csv("pups.csv", index=False)