First of all thanks a lot for the detailed answer. Although I think I am either missing something or I haven’t explained it clearly.
The following line is fine. One mean of racetimes for each year. 7 years, 7 mean race times for each year.
year_mean = pm.Normal(‘year_mean’, dims=[‘year’])
I have a question concerning the following line
runner_mean = pm.Normal(‘runner_mean’, dims=[‘runner’])
This means one mean race time for each runner. But i don’t want this. If a runner1 has run 3 races in year1 , 2 races in year2, 5 races in year3 i want at the end 3 mean race times for each one of the years. If runner2 has run in two years i need two means for him, one for each year.
So each runner should have one mean for every of the years he has run but runners in general have different number of means as they have raced in different number of years.