Рост хоккеистов: анализируем данные всех чемпионатов мира в текущем веке
На днях завершился очередной чемпионат мира по хоккею.
За просмотром матчей родилась идея. Когда в перерывах телевизионная камера показывает уходящих в раздевалку игроков, трудно не заметить, насколько они огромные. На фоне тренеров, функционеров команд, сотрудников ледовой арены, журналистов или просто фанатов они, как правило, выглядят очень внушительно.
И я задался вопросами. Действительно ли хоккеисты выше обычных людей? Как изменяется рост хоккеистов со временем в сравнении с обычными людьми? Есть ли устойчивые межстрановые различия?
ДанныеIIHF, организация, проводящая чемпионаты мира по хоккею, каждый год публикует составы участвующих команд с информацией о росте и весе каждого игрока. Архив этих данных тут.
Я собрал вместе данные всех чемпионатов мира с 2001 по 2016 годы. От года к году формат предоставления данных слегка меняется, что требует некоторых усилий по их очистке. Не представляя, как грамотно автоматизировать процесс, все данные копировал вручную, что заняло чуть больше 3 часов. Объединенный датасет выложил в открытый доступ.
Растут ли хоккеисты? Грубое (периодное) сравнениеДля начала сравним средний рост игроков на всех 16 чемпионатах мира.
Положительный тренд очевиден. За полтора десятилетия средний рост хоккеиста на чемпионате мира увеличился почти на 2 сантиметра (левая панель). Как будто бы незначительный прирост на фоне довольно большой вариации (правая панель). Много это или мало? Чтобы ответить на вопрос, надо корректно сравнить с населением (но об этом ближе к концу статьи).
Когортный анализБолее корректный способ изучения изменения в росте подразумевает сравнение по когортам рождения. Тут мы сталкиваемся с любопытным нюансом — некоторые хоккеисты участвовали не в одном чемпионате мира. Вопрос: вычищать ли повторные записи для одних и тех же людей? Если нам интересен средний рост хоккеиста на чемпионате (как на картинке выше), пожалуй, не имеет смысла зачищать. Но если мы хотим проследить изменение роста хоккеистов как таковое, на мой взгляд, было бы неправильно присваивать больший вес тем игрокам, которые регулярнее попадали на чемпионаты мира. Поэтому для дальнейшего анализа я очистил данные от повторных записей одних и тех же игроков.
Общее количество наблюдений сократилось с 6292 до 3333. Если хоккеист участвовал более чем в одном чемпионате мира, данные о росте и весе я усреднял, поскольку рост и (в особенности) вес отдельно взятого хоккеиста мог меняться со временем. Сколько же раз хоккеисты удостаиваются чести сыграть за национальные сборные на чемпионатах мира? В среднем чуть менее 2 раз.
Но есть и уникумы. Посмотрим, кто из игроков принял участие как минимум в 10 чемпионатах мира. Таких игроков оказалось 14.
name country position birth cohort av.height times_participated av.weight bmi 1 ovechkin alexander RUS F 1985-09-17 1985 188.45 11 98.36 27.70 2 nielsen daniel DEN D 1980-10-31 1980 182.27 11 79.73 24.00 3 staal kim DEN F 1978-03-10 1978 182.00 10 87.80 26.51 4 green morten DEN F 1981-03-19 1981 183.00 12 85.83 25.63 5 masalskis edgars LAT G 1980-03-31 1980 176.00 12 79.17 25.56 6 ambuhl andres SUI F 1983-09-14 1983 176.80 10 83.70 26.78 7 granak dominik SVK D 1983-06-11 1983 182.00 10 79.50 24.00 8 madsen morten DEN F 1987-01-16 1987 189.82 11 86.00 23.87 9 redlihs mikelis LAT F 1984-07-01 1984 180.00 10 80.40 24.81 10 cipulis martins LAT F 1980-11-29 1980 180.70 10 82.10 25.14 11 holos jonas NOR D 1987-08-27 1987 180.18 11 91.36 28.14 12 bastiansen anders NOR F 1980-10-31 1980 190.00 11 93.64 25.94 13 ask morten NOR F 1980-05-14 1980 185.00 10 88.30 25.80 14 forsberg kristian NOR F 1986-05-05 1986 184.50 10 87.50 25.70Александр Овечкин, 11 раз! Но тут надо отметить, что не для всех хоккеистов в принципе возможно было поучаствовать во всех 16 чемпионатах: зависит когорты рождения (насколько игровая карьера пересеклась именно с этим периодом наблюдения), от того, участвовала ли сборная игрока во всех чемпионатах мира (см. рисунок 3) и попадал ли игрок стабильно в сборную; наконец есть еще НХЛ, стабильно отвлекающий лучших из лучших от участия в чемпионатах мира.
Растут ли хоккеисты? Регрессионный анализРегрессионный анализ позволяет более корректно ответить на вопрос об изменении роста игроков. В данном случаем с помощью мультиноминальной линейной регрессии предсказыватся рост хоккеиста в зависимость от когорты рождения. Включая в спецификацию регресиионной модели различные дополнительные (контрольные) переменные, мы получаем значение наиболее интересующего нас коэффициента "при прочих равных". Например, добавляя к объясняющим переменным помимо когорты рождения позицию игрока на поле, мы получаем взаимосвязь роста и когорты, очищенную от эффекта различий в зависимости от позиции; добавляя в контрольны переменные страны, получаем результат, очищенный от межстрановых различий. Разумеется, если контрольные переменные сами оказываются значимыми, на это тоже стоит обратить внимание. Регрессионные модели (особенно линейные регрессии) очень чувствительны к выбросам (см., например, эту статью). Не вдаваясь глубоко в эту обширную тему, я лишь убрал из анализа когорты, для которых мы имеем слишком небольшое количество представителей.
Не желая резать данные сильно, я убрал только когорты 1963, 1997 и 1998 годов рождения, для которых у нас есть менее 10 игроков.
Итак, результаты рагрессионного анализа. В каждой следующей модели я добавляю одну переменную. Зависимая переменная: рост хоккеиста. Объясняющие перемеенные: 1) когорта рождения; 2) + позиция на поле (сравнение с защитниками); 3) + страна (сравнение с Россией).
Statistical models Model 1 Model 2 Model 3 (Intercept) -10.17 (27.67) -18.64 (27.01) 32.59 (27.00) cohort 0.10 (0.01) *** 0.10 (0.01) *** 0.08 (0.01) *** positionF -2.59 (0.20) *** -2.59 (0.20) *** positionG -1.96 (0.31) *** -1.93 (0.30) *** countryAUT -0.94 (0.55) countryBLR -0.95 (0.53) countryCAN 1.13 (0.46) * countryCZE 0.56 (0.49) countryDEN -0.10 (0.56) countryFIN 0.20 (0.50) countryFRA -2.19 (0.69) ** countryGER -0.61 (0.51) countryHUN -0.61 (0.86) countryITA -3.58 (0.61) *** countryJPN -5.24 (0.71) *** countryKAZ -1.16 (0.57) * countryLAT -1.38 (0.55) * countryNOR -1.61 (0.62) ** countryPOL 0.06 (1.12) countrySLO -1.55 (0.58) ** countrySUI -1.80 (0.53) *** countrySVK 1.44 (0.50) ** countrySWE 1.18 (0.48) * countryUKR -1.82 (0.59) ** countryUSA 0.54 (0.45) R 2 0.01 0.06 0.13 Adj. R 2 0.01 0.06 0.12 Num. obs. 3319 3319 3319 RMSE 5.40 5.27 5.10 *** p < 0.001, ** p < 0.01, * p < 0.05 Интерпретация моделейМодель 1. Увеличение когорты на один год соответсвует увеличению роста хоккеистов на 0.1 см. Коэффициент статистически значим, но при этом модель объясняет лишь 1% вариации зависимой переменной. В принципе это не проблема, поскольку моделирование носит объясняющий характер, задача предсказания не ставится. Тем не менее, низкий коэффициент детерминации показывает, что должны быть другие переменные, гораздо лучше объясняющие различия между хоккеистами в росте.
Модель 2. Защитники — самые высокие игроки в хоккее. Вратари ниже на 2 см, нападающие — на 2.6 см. Все коэффициенты статистически значимы. Объясненная вариация зависимой переменной возрастает до 6%. При этом коэффициент при переменной когорта рождения не изменяется.
Модель 3. Добавление контрольных переменных для стран любопытно по двум причинам. Во-первых, некоторые различия статистически значимы и интересны сами по себе. Так, например, шведы, словаки и канадцы статистически значимо выше наших игроков. Большинство же наций значительно ниже нас, японцы аж на 5.2 см, итальянцы — на 3.6 см, французы — на 2.2 см (см. также рисунок 4). Во-вторых, введение контрольных переменных для стран значительно уменьшает коэффициент при переменной когорта рождения — до 0.08. Это значит, что межстрановые различия объясняют часть различий по когортам рождения. Коэффициент детерминации модели возрастает до 13%.
Наиболее полная модель показывает, что увеличение роста хоккеистов происходит со скоростью 0.08 см в год. Это означает прирост 0.8 см за десятилетие или на 2.56 см за 32 года с 1964 по 1996. Обратите внимание, что при учете контрольных переменных скорость увеличения роста хоккеистов оказывается примерно в полтора раза ниже, чем при более грубом анализе средних значений (рисунок 1): 0.8 см за десятилетие против примерно 1.2 см.
Прежде чем мы, наконец, постараемся понять, насколько значительным оказывается увеличение роста, хочу обратить внимание еще на один любопытный момент. Введение контрольных переменных подразумевает фиксацию различий между категориями при едином наклоне регрессионной линии (единый коэффициент при главной объясняющей переменной). Это не всегда хорошо и может замаскировать значительные различия в тесноте связи между исследуемыми переменными в подвыборках. Так, например, раздельное моделирование зависимости роста игроков от амплуа (рисунок 5) показывает, что взаимосвязь наиболее ярко выражена для вратарей и наименее заметна для защитников.
Statistical models Model 3 D Model 3 F Model 3 G (Intercept) 108.45 (46.46) * 49.32 (36.73) -295.76 (74.61) *** cohort 0.04 (0.02) 0.07 (0.02) *** 0.24 (0.04) *** countryAUT 0.14 (0.96) -2.01 (0.75) ** 0.47 (1.47) countryBLR 0.30 (0.87) -1.53 (0.73) * -2.73 (1.55) countryCAN 1.55 (0.78) * 0.39 (0.62) 3.45 (1.26) ** countryCZE 0.87 (0.84) 0.30 (0.67) 0.63 (1.36) countryDEN -0.60 (0.95) 0.10 (0.75) -0.19 (1.62) countryFIN -0.55 (0.89) -0.04 (0.67) 2.40 (1.32) countryFRA -3.34 (1.15) ** -2.06 (0.93) * 1.39 (2.07) countryGER 0.48 (0.85) -1.40 (0.72) -0.65 (1.33) countryHUN -1.32 (1.47) -0.70 (1.16) 0.65 (2.39) countryITA -2.08 (1.08) -4.78 (0.82) *** -2.02 (1.62) countryJPN -4.13 (1.26) ** -6.52 (0.94) *** -2.27 (1.98) countryKAZ -1.23 (0.95) -1.82 (0.79) * 1.79 (1.58) countryLAT -0.73 (0.95) -1.39 (0.75) -3.42 (1.49) * countryNOR -3.25 (1.07) ** -1.06 (0.85) -0.10 (1.66) countryPOL 0.82 (1.89) -0.58 (1.55) 0.37 (2.97) countrySLO -1.57 (0.99) -1.54 (0.79) -2.25 (1.66) countrySUI -1.98 (0.91) * -2.36 (0.71) *** 1.12 (1.47) countrySVK 2.94 (0.87) *** 0.81 (0.67) -0.70 (1.50) countrySWE 0.75 (0.81) 1.24 (0.65) 1.37 (1.33) countryUKR -1.37 (1.01) -1.77 (0.80) * -3.71 (1.66) * countryUSA 0.76 (0.78) -0.08 (0.62) 2.58 (1.26) * R 2 0.09 0.10 0.24 Adj. R 2 0.07 0.09 0.20 Num. obs. 1094 1824 401 RMSE 5.08 5.08 4.87 *** p < 0.001, ** p < 0.01, * p < 0.05Раздельное моделирование показывает, что в когортах 1964-1996 годов рождения, средний рост хоккеистов, участвовавших в чемпионатах мира в 2001-2016 годах, увеличивался со скоростью 0.4 см за десятиление для защитников, 0.7 см — для нападающих и (!) 2.4 см — для вратарей. За три десятиления средний рост вратарей увеличился на 7 см!
Пришло время сравнить эти изменения со средними значениями для населения.
Сравнение с населениемРезультаты регрессионного анализа фиксируют значительные межстрановые различия. Поэтому сравнивать имеет смысл по странам: хоккеистов определенной страны с мужским населением этой же страны.
Для сравнения роста хоккеистов со средними показателями мужского населения я использовал данные из релевантной научной статьи (PDF). Данные я скопировал из статьи (использовав замечательную программку tabula) и тоже разместил в открытом доступе.
К сожалению, данные о динамике роста населения пересекаются лишь с 8 странами из моего хоккейного датасета: Австрия, Дания, Финляндия, Франция, Германия, Италия, Норвегия, Швеция.
Во всех проанализировнных странах хоккеисты выше стеднестатистических мужчин на 2-5 см. Но это не удивительно — в спорте значительная селекция. Примечательно другое. В развитых странах мира особенно бурное увеличение роста мужского населения происходило в первой середине 20 века. В когортах примерно 1960-х годов рождения рост мужчин приблизился к плато и пеерстал бурно увеличиваться. Тренд среднего роста хоккеистов во всех странах (кроме почему-то Дании) как будто бы продолжил приостановившийся многолетний тренд всего мужского населения. Для когорт европейцев, родившихся в первой половине 20 века, темпы увеличения среднего роста варьировались от 1.18 до 1.74 см за десятилетие в зависимости от страны (рисунок 7). Начиная с 1960-х годов этот показатель опустился до уровня 0.15-0.80 за 10 лет.
На фоне стагнирующего тренда в населении увеличение роста хоккеистов выглядит весьма внушительным. А акселерация среди вратарей вообще беспрецедентна. Не стоит забывать и про селекцию. Расхождение трендов в населении и среди хоккеистов, вероятно, свидетельствует об усиливающейся селекции — хоккей требует все большего роста для успешной карьеры.
Селекция в спортеПроглядывая научную литературу по теме я наткнулся на примечательный результат. Оказывается, в профессиональном спорте преобладают люди, рожденные в первой половине года. Объясняется это тем, что спортивные секции, как правило, формируют детские команды по когортам рождения. Таким образом, рожденные в начале года, всегда имеют чуть больше прожитого времени за плечами, что зачастую прямо выражается в физическом превосходстве над сверстниками, рожденными под конец года. Нетрудно проверить этот результат на нашем датасете.
Действительно, респределение довольно сильно смещено в сторону ранних месяцев. Если разбить данные по декадам рождения, то невооруженным глазом видно, что эффект усиливается со временем (рисунок 9). Косвенно это свидетельствует о том, что селекция в хоккее становится жестче.
На будущееЛюбопытно будет посмотреть, влияют ли физические данные на игровую статистику хоккеистов. Наткнулся на занимательную статью, опубликованную в очень приличном научном журнале, в которой авторы нашли корреляцию между соотношением пропорций лица хоккеиста и средним количеством штрафных минут за игру.
ReproducibilityПолный R скрипт, воспроизводящий результаты моей статьи, тут. Использована версия R-3.2.4 Все пакеты по состоянию на 2016-03-14. В случае пакетных несовместимостей, данный код будет гарантированно воспроизведен при использовании пакета checkpoint с указанием соответствующей даты.