Data Study on Factors affecting marathon results

Posted on Aug 19, 2016
The skills the author demoed here can be learned through taking Data Science with Machine Learning bootcamp with NYC Data Science Academy.

We run faster than you: an data analysis of marathon results from Switzerland

The Jungfrau is one of the main summits in the Bernese Alps in Switzerland. Each year, data shows thousands of people choose to run a marathon through the valley...and up the mountain. In 2010, I was one of them.

On one of the steeper sections, an elderly gentleman flew past me. What had just happened?! He was not the only person with grey or white hair running in the marathon. In addition, Swiss citizens were distinguished by a different color number start number. I noticed the Swiss cheered a bit louder for their own. So where do people come from who compete, and do the Swiss have a home-field advantage? This blog article will answer these questions.

The data was scraped using Python from the Jungfrau Marathon Website. Data for the years 2015, 2014, 2013, 2012, and 2010 were collected. Other years were not collected for this analysis because they used a completely different website format. The data analysis was done in R. Click here to see the Python and R code.

Data on External and Personal Factors

All factors which influence your marathon time can be easily split into external factors and personal (internal) factors. In this blog article we will investigate one external factor, the number of participants in the race, and three personal factors: gender, age, and nationality. There are myriad additional factors which could be considered but were not available in the scraped data, so we will limit our analysis to these four.

Distribution of Race Times

Data Study on Factors affecting marathon results

Data Analysis

There were 15% fewer race participants in 2012 as compared to the average: around 3200 compared to an average from the other years of nearly 3800. Interestingly, 2012 had a higher median race time than the other four years. The drop in the marathon times in 2012 could have been due from many faster runners deciding to not participate or harsher conditions.

The mean of the race times of one of the years is significantly different from the others. This difference might appear significant to anyone looking at the boxplot, but in order for it to be statistically significant, a statistical test called a one-way ANOVA must be used to determine if the mean of race times of one year significantly deviate from the mean of other years.

The test calculates the probability that the values from one year will be observed had they been drawn from the distribution of race times represented by the other years. In this case, the probability is less than the commonly-used cut-off value of 5% (p-value 0.0003). This means there is not just a practical difference seen by plain observation of the boxplot but also a statistical difference.

Mean Results

The mean race times of 2012 differ significantly from the mean race times of 2015, the year with the next closest median value. In this case, the practical observation of the boxplot was used to inform which years to compare. In order to determine this, two steps must be taken, each with its own statistical tool:

  • 1. calculate the probability of observing the distribution of race times from 2012 under the assumption that they came from the distribution represented by the race times from 2015 using a two-way T-test, a special version of the one-way ANOVA mentioned above
  • 2. perform a test to see if the shape of the distributions from the two years is similar, also called having similar variances, using an F-test.

The results: 1. there almost no chance of seeing the race times from 2012 if they came from the distribution of times represented by the times from 2015 (p-value 2.2e-16), and 2. the variances of the the two distributions are not significantly different (p-value 0.427).

Data on Gender

Data Study on Factors affecting marathon results

Your gender appears to affect your race time, but statistically the two cannot be compared. Why not? The assumption that the two distributions are similar, meaning they have the same variance, is not true. To determine if the means of two distributions are significantly different, a statistical test called a two-way T-test is performed.

The result shows that indeed the mean man's time is significantly less than the mean woman's time (p-value 2.2e-16). The test which checks if the variances of the two distributions are similar, the F-test, shows that the variances are indeed significantly different (p-value 2.2e-8). The assumption needed to claim the result of the T-test is not true, invalidating the T-test and not allowing a direct statistical comparison of the two genders.

Data on Age

03_Your_age_impacts

The subjective experience of being passed by an elderly person was shown to not extend to the whole race population because the age group bins cannot be directly statistically compared. The F-test, being in this case a ratio of variances for the age groups [17,35) and [57,77], indicated that the variances are significantly different (p-value 2.2e-16).  The median time visually increases with age when the ages were binned into seven equal quantile bins, but the statistics are not able to confirm this observation.

Nationality

05_Proximity_to_Switzerland

The race roster showed that Switzerland had the most participants and four of its neighboring countries were in the top ten: Germany, Italy, France, and Austria. Non-bordering countries rounding out the top ten were Great Britain, Netherlands, Belgium, Luxembourg, and USA. The US was the only non-European country in the top ten.

Data Findings

06_The_Italians_are_fastest

The Italians proved to be the fastest out of the top five participating countries, edging out host country Switzerland in median race time.

The one-way ANOVA indicates one of the group mean race times is different than the others. The F-test shows, however, that, like age and gender analyses, the mean race time of Italians also cannot be compared to that of Switzerland due to differing variances (p-value 1.42e-06).

Interestingly enough, the race times of the Germans and the Dutch are significantly different from one another (t-test p-value: 0.0007, F-test p-value: 0.092).

Investigating the boxplot closely shows that the range between the bottom and top of the box (interquartile range) of the Italians appears larger than the other nationalities. This invalidates statistical comparisons of their distributions with the two-way T-test. That of the Germans and Dutch are more similar.

Gender does not change time distribution by nationality

07_Gender_does_not

So if the Germans are faster than the Dutch, does this apply to both genders? Surprisingly, the nationalities were the same when considering median race time across genders. Perhaps partners train together and therefore have a similar race profile? That, however, is an analysis for another time requiring additional data.

Takeaway

- The Italians have the lowest median time
- Gender does not change nation rankings
- Age is running still matters: younger is better

About Author

David Richard Steinmetz

David became a data scientist for two reasons: it was the part of his previous jobs he loved, and he saw the need companies have at interpreting their data to meet business goals. With a PhD in Materials...
View all posts by David Richard Steinmetz >

Leave a Comment

No comments found.

View Posts by Categories


Our Recent Popular Posts


View Posts by Tags

#python #trainwithnycdsa 2019 2020 Revenue 3-points agriculture air quality airbnb airline alcohol Alex Baransky algorithm alumni Alumni Interview Alumni Reviews Alumni Spotlight alumni story Alumnus ames dataset ames housing dataset apartment rent API Application artist aws bank loans beautiful soup Best Bootcamp Best Data Science 2019 Best Data Science Bootcamp Best Data Science Bootcamp 2020 Best Ranked Big Data Book Launch Book-Signing bootcamp Bootcamp Alumni Bootcamp Prep boston safety Bundles cake recipe California Cancer Research capstone car price Career Career Day citibike classic cars classpass clustering Coding Course Demo Course Report covid 19 credit credit card crime frequency crops D3.js data data analysis Data Analyst data analytics data for tripadvisor reviews data science Data Science Academy Data Science Bootcamp Data science jobs Data Science Reviews Data Scientist Data Scientist Jobs data visualization database Deep Learning Demo Day Discount disney dplyr drug data e-commerce economy employee employee burnout employer networking environment feature engineering Finance Financial Data Science fitness studio Flask flight delay gbm Get Hired ggplot2 googleVis H20 Hadoop hallmark holiday movie happiness healthcare frauds higgs boson Hiring hiring partner events Hiring Partners hotels housing housing data housing predictions housing price hy-vee Income Industry Experts Injuries Instructor Blog Instructor Interview insurance italki Job Job Placement Jobs Jon Krohn JP Morgan Chase Kaggle Kickstarter las vegas airport lasso regression Lead Data Scienctist Lead Data Scientist leaflet league linear regression Logistic Regression machine learning Maps market matplotlib Medical Research Meet the team meetup methal health miami beach movie music Napoli NBA netflix Networking neural network Neural networks New Courses NHL nlp NYC NYC Data Science nyc data science academy NYC Open Data nyc property NYCDSA NYCDSA Alumni Online Online Bootcamp Online Training Open Data painter pandas Part-time performance phoenix pollutants Portfolio Development precision measurement prediction Prework Programming public safety PwC python Python Data Analysis python machine learning python scrapy python web scraping python webscraping Python Workshop R R Data Analysis R language R Programming R Shiny r studio R Visualization R Workshop R-bloggers random forest Ranking recommendation recommendation system regression Remote remote data science bootcamp Scrapy scrapy visualization seaborn seafood type Selenium sentiment analysis sentiment classification Shiny Shiny Dashboard Spark Special Special Summer Sports statistics streaming Student Interview Student Showcase SVM Switchup Tableau teachers team team performance TensorFlow Testimonial tf-idf Top Data Science Bootcamp Top manufacturing companies Transfers tweets twitter videos visualization wallstreet wallstreetbets web scraping Weekend Course What to expect whiskey whiskeyadvocate wildfire word cloud word2vec XGBoost yelp youtube trending ZORI