Your situation is called an “unbalanced panel”. In PyMC, you can handle this using indicator variables. Here’s a similar situation. In your case, you will have two indicators, one for time and one for runner. You’ll make a single runner_mean
with shape 100, and a single year_mean
with shape 4, then you’ll use the indices to broadcast them together.