# Simulation der sprint Daten #
# =========================== #

set.seed(123)

n1 = 10
n2 = n1
n = n1 + n2
M = 3
b = rep(rnorm(n1 + n2, 0, 0.5), each = M)
sex = c(rep(0, n1), rep(1, n2))
beta = c(rep(10.6, n1 * M), rep(11.9, n2 * M))
eps = rnorm(n, 0, 0.2)
times = beta + b + eps
athlete = rep(1:n, each = M)

sprint <- data.frame(athlete, sex = rep(sex, each = M), times)