NYC Data Science Academy| Blog
Bootcamps
Lifetime Job Support Available Financing Available
Bootcamps
Data Science with Machine Learning Flagship 🏆 Data Analytics Bootcamp Artificial Intelligence Bootcamp New Release 🎉
Free Lesson
Intro to Data Science New Release 🎉
Find Inspiration
Find Alumni with Similar Background
Job Outlook
Occupational Outlook Graduate Outcomes Must See 🔥
Alumni
Success Stories Testimonials Alumni Directory Alumni Exclusive Study Program
Courses
View Bundled Courses
Financing Available
Bootcamp Prep Popular 🔥 Data Science Mastery Data Science Launchpad with Python View AI Courses Generative AI for Everyone New 🎉 Generative AI for Finance New 🎉 Generative AI for Marketing New 🎉
Bundle Up
Learn More and Save More
Combination of data science courses.
View Data Science Courses
Beginner
Introductory Python
Intermediate
Data Science Python: Data Analysis and Visualization Popular 🔥 Data Science R: Data Analysis and Visualization
Advanced
Data Science Python: Machine Learning Popular 🔥 Data Science R: Machine Learning Designing and Implementing Production MLOps New 🎉 Natural Language Processing for Production (NLP) New 🎉
Find Inspiration
Get Course Recommendation Must Try 💎 An Ultimate Guide to Become a Data Scientist
For Companies
For Companies
Corporate Offerings Hiring Partners Candidate Portfolio Hire Our Graduates
Students Work
Students Work
All Posts Capstone Data Visualization Machine Learning Python Projects R Projects
Tutorials
About
About
About Us Accreditation Contact Us Join Us FAQ Webinars Subscription An Ultimate Guide to
Become a Data Scientist
    Login
NYC Data Science Acedemy
Bootcamps
Courses
Students Work
About
Bootcamps
Bootcamps
Data Science with Machine Learning Flagship
Data Analytics Bootcamp
Artificial Intelligence Bootcamp New Release 🎉
Free Lessons
Intro to Data Science New Release 🎉
Find Inspiration
Find Alumni with Similar Background
Job Outlook
Occupational Outlook
Graduate Outcomes Must See 🔥
Alumni
Success Stories
Testimonials
Alumni Directory
Alumni Exclusive Study Program
Courses
Bundles
financing available
View All Bundles
Bootcamp Prep
Data Science Mastery
Data Science Launchpad with Python NEW!
View AI Courses
Generative AI for Everyone
Generative AI for Finance
Generative AI for Marketing
View Data Science Courses
View All Professional Development Courses
Beginner
Introductory Python
Intermediate
Python: Data Analysis and Visualization
R: Data Analysis and Visualization
Advanced
Python: Machine Learning
R: Machine Learning
Designing and Implementing Production MLOps
Natural Language Processing for Production (NLP)
For Companies
Corporate Offerings
Hiring Partners
Candidate Portfolio
Hire Our Graduates
Students Work
All Posts
Capstone
Data Visualization
Machine Learning
Python Projects
R Projects
About
Accreditation
About Us
Contact Us
Join Us
FAQ
Webinars
Subscription
An Ultimate Guide to Become a Data Scientist
Tutorials
Data Analytics
  • Learn Pandas
  • Learn NumPy
  • Learn SciPy
  • Learn Matplotlib
Machine Learning
  • Boosting
  • Random Forest
  • Linear Regression
  • Decision Tree
  • PCA
Interview by Companies
  • JPMC
  • Google
  • Facebook
Artificial Intelligence
  • Learn Generative AI
  • Learn ChatGPT-3.5
  • Learn ChatGPT-4
  • Learn Google Bard
Coding
  • Learn Python
  • Learn SQL
  • Learn MySQL
  • Learn NoSQL
  • Learn PySpark
  • Learn PyTorch
Interview Questions
  • Python Hard
  • R Easy
  • R Hard
  • SQL Easy
  • SQL Hard
  • Python Easy
Data Science Blog > Student Works > An analysis of motorcyclist habits through online sales

An analysis of motorcyclist habits through online sales

Thomas Boulenger
Posted on Feb 25, 2016

Contributed by Thomas Boulenger. He is currently in the NYC Data Science Academy 12 week full time Data Science Bootcamp program taking place between January 11th to April 1st, 2016. This post is based on his third class project – Web scraping (due on the 6th week of the program).

sportster-customized

Overview and assumptions

The motorcycle industry is one that deals with passion, its History although parallel to that of the automotive industry, has long been restricted to a small circle of purists. Motorcycle buyers quite often remain loyal to a brand they like and they identify to. This implies that one often expects a given brand owner to show some behaviour typically associated to the brand's image. On average one would assume that BMW owners care more about security than others, drive their motorcycles more often and for longer distances than others, including holiday road-trips. On the other hand, Ducati owners are expected to use their motorcycle less often, but to care more about the look, the sound and general appeal of their vehicle.

Within a given brand, the same sort of differences in habits are expected from one particular model to another. Although maybe not as a strong distinction as the one by brand, we still expect to see something here.

An online sales listing website: Cycle-Crunch.com

Screen Shot 2016-02-24 at 09.50.21

Unlike other Western motorcycle markets, in the US the sales have long been dominated by one Historic emblematic brand: Harley-Davidson. However, in recent years the company have grown worried as their historical pool of customers, Caucasian 35+ men, ages and gets thinner over time. On the company's website, one may find information about the Demographics of their customers and shows that not only is the company aware of the matter, but also willing to address it.

Screen Shot 2016-02-24 at 11.39.38

Aim of the analysis

Here we wish to get a glimpse at the US market sales for motorcycles by scraping the website we have selected to obtain a dataset that will summarise the listed motorcycles' information (price, year of production, brand, model, city, color, listed date). The hope is then that the data we obtain allows us to say something about some brands characteristics, the owners habits. If a listed motorcycle has had a high average mileage over its lifetime, how likely is it to be an Harley-Davidson? Or a Yamaha? Can we find better suited indicators to decide wether a given listed motorcycle belongs to a given brand, or to be a given model?

Scraping with BeautifulSoup

To ease the scraping, we use the BeautifulSoup package in Python. Running a loop, we sift through all the results front pages by updating our current page url after we are done retrieving all the information we needed. More specifically, for all front page of the search results we find every box's url link to go look at the page dedicated to each listed model and get all the available details (such as mileage and descriptions which are for instance missing from the front page). See the Appendix A below for the detailed code.

The dataset

After a few hours we were able to get the information contained in 1,282 front pages, giving us the details of 12,812 listed motorcycles. We thus obtain the following dataset

CycCrDat

Some of the variables are categorical (color, condition, listed, make i.e. brand, model and state), so in order to compare them with the numerical variables, we first add some new variables representing an index for any of these categorical variables. It turns out that most of the above variables are quite poorly correlated to one another.

Screen Shot 2016-02-24 at 16.49.05

In order to get a better sense of what the data we have obtained looks like, we can plot the distribution of the motorcycle listed per brand

Screen Shot 2016-02-24 at 16.48.19

This seems a clear sign of the long time domination of Harley-davidson over the US market, but the magnitude by which Harley dominates over all the other brands is still a little surprising. When restricting at heavyweight motorcycles (those with engines of 601 cubic centimeters or greater), Harley Davidson represents roughly 50% of the US market share. However, in recent years signs of a decline have started to show, with their market share for the first semester of 2015 falling at 47.5%, below the psychological mark of the 50%. Yet here we see no sign of such a decrease in popularity of Harley Davidson models. If anything, the trend even seems to go in the opposite direction, with an even greater proportion of Harley Davidson among models to be sold as new than among models to be sold as used.

Screen Shot 2016-02-24 at 17.14.30

This kind of distribution seems a very particular feature of the US market where motorcycles owners often buy in the "chopper" category, that is low, heavy,  very stable and relatively slow motorbikes but with a lot of torque. Among the heavyweight motorcycle category, the best selling model in the US would thus looks like the following Street Glide Special

hd-street-glide-special

Characteristics by brand

Let us now look a little more closely at our data and see wether we can find patterns to characterise brands or models.

- Evolution of prices

Though few data are available for years of production < 1980 (only 35 motorcycles listed), we see a funny bouncing showing how valuable old-style models have become, and that the rebound in devaluation seems to happen some time in the 1980s.

Screen Shot 2016-02-24 at 17.25.29

Based on this data it thus seems one should wait 35 years on average before a model becomes a classic, and that provided it is in good condition, its listed price starts increasing again as it ages. Zooming in recent years we see the slopes for Harley-Davidson and Indians appear steeper than the others, probably indicating a worse depreciation rate than for the most common competitor Japanese brands plotted here (Yamaha, Honda, Suzuki, Kawasaki).

We may next attempt to roughly average the amount by which motorcycles depreciate every year among every brand. To do so, for every brand we subset the used motorcycles and keep only the models for which at least one new such model is also listed on sale. For every model we then average the price over the ones listed as new and use it as our "Model Price as New". Within the dataset for used motorcycles, we may then calculate the yearly depreciation ratio. We then obtain the following graph

Screen Shot 2016-02-24 at 17.25.47

We see what is probably a big difference in Harley-Davidson owners habits: the difference between the price of used models and their price as new tend to be negative, as most riders love to customize their motorcycle, often spending half the price as new in accessories/painting jobs/exhaust pipes etc.. Of course this means that overall, they lose more money than others every year, but there is no simple way to find this out for sure. A possibility would be to use the pictures available in the description, compare them with original models and quantify how much the two differ.

Although we don't have enough observations for the other brands to make a definitive claim, it seems the depreciation per year is pretty low and similar for Honda, Suzuki and Yamaha. Indian on the other hand seems to be showing a negative feature, similar to Harley-Davidson. While one could assume owners of Honda, Suzuki or Yamaha tend to use their vehicle for more practical reason such as commuting, Indians owners may sow a love for customisation similar to Harley-Davidson owners.

-Yearly Mileage

How often do owners use their motorcycle is another pretty obvious characterising feature. The only way we have to measure this here is through the yearly mileage that we estimate by the ratio of the kilometres the motorcycle yet totalizes with the time since it has been produced

Screen Shot 2016-02-25 at 11.21.54

Note that given there are so many more Harley-Davidson listed than any other brand, there is an overfit to the distribution that doesn't allow us to initially compare it to the others. In order to produce the graph above we have reduced the resolution on the "ydist" (yearly mileage) variable in the following fashion

Surprisingly, the only obvious peak of yearly mileage appears for Honda, when the others seem scattered all across the spectrum. We don't even see an obvious peak for Harley-Davidson, although the number of motorbikes listed tends to decrease as the yearly mileage increases, it appears almost evenly distributed. The lowest peaks appear for Honda then Yamaha, when the highest peaks appear for Suzuki. Once again, Indian seems to closely follow Harley-Davidson's distribution. We start to clearly identify, at least within our dataset, two different kind of observations. One is of the like of Harley-Davidson and Indian, with a large distribution of different yearly mileage, probably corresponding to different sort of usage for different kind of customers. The other most obvious one appears for Suzuki and show a relatively large peak around 3.500 km a year on average and is likely to correspond to commuters. The other Japanese brands show much lower peaks, suggesting most owners do not drive many miles a year.

Weather information

One thing that most riders dread is a rainy day and cold temperatures, as one might assume such things make the ride unpleasant to most riders. Owning a motorcycle in a region where cold temperatures and rainy conditions are common probably says a lot about an owner's behaviour. A regular commuter will have to drive almost every day and often face driving in the rain, while more occasional users will only get on a ride when the weather is nice.

In order to say something about the correlation between the local weather and an owner's habits, we must first obtain the weather's data. To do so we use the cities listed in our Dataframe from Cycle-Crunch.com and use it to sift through the online database provided by weatherbase.com:

Screen Shot 2016-03-31 at 21.33.42

 

Once again, we use the BeautifulSoup package in python to scrape the website. The procedure turns out a little more tricky than previously, as the website expects a response when sending cookies. One has to use the Mechanize package in Python to simulate a browser so that weatherbase.com does not load an error page. Furthermore there are some subtleties occurring when the initial query leads to multiple results or if the given city does not match any recorded location of the database. In that case the query's response is often a new search page with multiple results in a small radius (often 5-20 miles) of the initial query. Using the first closest location, we may thus obtain a pretty accurate approximation of the location's weather conditions. For the sake of simplicity we retain only annual figures for the following features: Average Precipitation, Average Number of Days With Precipitation, Average Number of Days Below 32F/0C, Average Number of Days Above 90F/32C, Average Number of Rainy Days.

Despite using the closest location when there is no result for the listed city, we still obtain some rows in the initial motorcycles DataFrame left with only partial weather information (one of the selected feature may simply be missing for a given location). Anyhow, we end up with additional valuable information to work with. See Appendix B for more details about the code for scraping the weather data.

The idea now is to try to looks at how the different brands and the new weather information behave with respect to one another. We for instance try to answer questions such as:

  • What is the distribution of every brand with respect to the weather features? In other words, what is the conditional probability P(feature split|brand)?
  • Given that the average number of days with precipitation is between x and y, what is our sample probability to obtain such and such brands when a listed sale is picked at random?

Consider for instance the weather feature "average number of days with Temperatures below 0°C". We split the feature in two, the number of such days is either less or more than 140 days a year. In the code below, note we can either split automatically in an even number of splits or provide a manual one. After trying a few cut-offs, here we decided to select 140 days for our cut-off value as it revealed an interesting contrast. We obtain the following set of proportions for the brands with at least 10 listed models:

Screen Shot 2016-04-02 at 18.03.57

 

We observe that 84.7% of the 59 yamaha listed where in "cold" zones with at least 140 days with temperatures below 0°C. This seems a common phenomenon among Japanese brands, as Honda, Suzuki and Kawasaki respectively have 68.7%, 92% and 72.2% of their observations listed in the "cold" zone (with 32, 25 and 18 total observations listed). This is in sheer contrast to Harley-Davidson, Buell, Big Dog or Victory for which the majority of observations are listed in the "hot" zones with fewer than 140 days with temperatures below 0°C, or respectively 76.6%, 72.7%, 88% and 87.5% of their total observations.

This suggests that this particular weather feature might be a good classifier. Considering now each of those splits, we may look at the conditional probability of belonging to a brand knowing that the observation lies within one of the previous two splits.

Screen Shot 2016-04-02 at 18.05.28

For both splits, the vast majority of observations are Harley-Davidson. However they vary from 97.5% in the "cold" region to 90.6% in the "hot" region, quite a significant drop given the skewness of the motorcycles data. We may also note the following interesting fact: given that an observation is NOT an Harley-Davidson, in the "cold" region, we have a 31.25% chance that it is a Yamaha, whereas in the "hot" it is only a 5.8% chance. For the brand Yamaha, there is a factor 6 between cold and hot regions, when a priori assuming we do not have an Harley-Davidson. Although Harley-Davidson massively predominates, we can thus recover some ways to hone in on Japanese brands by splitting the data between these regions with different weather conditions.

 

To obtain the previous tables and more given any of the weather feature we selected, we used the following generic code:

 

Clustering

Let us finally make a final attempt to differentiate the different brands as best as we can by clustering the observations with both the yearly mileage and the weather  features for the average number of days with precipitation and below 0°C/above 32°C. We use the K-means algorithm implemented by the SciKit library in Python. Note we first normalise our features to make them more comparable to one another and prevent the apparition of artificial clusters. Given the large number of observations and how they are vastly clouded by the majority class of Harley-Davidson, we may try to use different numbers of clusters and see if a large number of them could pinpoint a given brand other than Harley-Davidson. We try this with Yamaha and its 59 observations for which we have a complete record of average number of days below 0°C. To determine what number of clusters we should select, we test a range of possibilities from 2 clusters up to 250. In order to decide which number of cluster is best we record the the probability of an observation being a Yamaha among every cluster. Surprisingly, the Yamaha's family of observation remains stable and mostly contained within a single cluster of 48 observations so long as we select less than 101 clusters.

The "purest" cluster is achieved the ridicule number of 246 clusters, in which case one of them contains 15 Yamaha, representing 35.7% of its total 42 observations. Another run led to a different 30.7% (in K-means, the centroids being initialised randomly, the final result may vary from one run to another). In the 25th cluster, there were then 20 Yamahas while only 22 Harley-Davidsons, 12 Suzuki and 11 Honda. In other word that cluster contained 66% of Japanese brands.

Screen Shot 2016-04-02 at 20.40.39

Now a more stable clustering using only 92 clusters achieved 15.5% of Yamahas and 45% of Japanese brands and 55% of Harley-Davidson. That cluster contained 81% of the Yamahas observations and was much less likely to be an Harley-Davidson than in the rest of the clusters. Of course, repeating the process a certain number of times and recording the best results we may refine these results a little. But the hurdle of clearing the clusters of any Harley-Davidson to differentiate it from other brands is still difficult.

 

Appendix A: code for scraping Cycle-Crunch.com

 

 

Appendix B: scraping weatherbase.com

 

 

About Author

Thomas Boulenger

Thomas holds a PhD in pure Mathematics from Paris, and started to pick up an interest for Machine Learning while still a Postdoc at Basel University. After spending a few years in academia, he decided to join NYC...
View all posts by Thomas Boulenger >

Leave a Comment

Cancel reply

You must be logged in to post a comment.

Google July 28, 2020
Google Although web sites we backlink to beneath are considerably not related to ours, we feel they may be truly really worth a go by means of, so possess a look.
www.Youtube.com November 17, 2016
Some truly nice and utilitarian information on this web site, also I believe the layout contains wonderful features.
wood rug November 12, 2016
This is the right blog for anybody who really wants to understand this topic. You know so much its almost hard to argue with you (not that I actually will need to...HaHa). You definitely put a brand new spin on a topic which has been written about for decades. Great stuff, just wonderful!
Elisa October 8, 2016
Great info it is definitely. We've been looking for this information.
Deon October 4, 2016
I wanted to thank you once more for this amazing web site you have developed here. It's full of useful tips for those who are definitely interested in this specific subject, primarily this very post. Your all so sweet as well as thoughtful of others plus reading your website posts is an excellent delight in my experience. And such a generous treat! Jeff and I usually have enjoyment making use of your tips in what we need to do next week. Our listing is a kilometer long and tips is going to be put to fine use.
www.youtube.com September 28, 2016
Great ? I should definitely pronounce, impressed with your site. I had no trouble navigating through all the tabs as well as related info ended up being truly easy to do to access. I recently found what I hoped for before you know it in the least. Quite unusual. Is likely to appreciate it for those who add forums or something, site theme . a tones way for your client to communicate. Nice task.
Steven O'Toole May 11, 2016
Awesome job with this BeautifulSoup and analysis. Had a great time looking this over!

View Posts by Categories

All Posts 2399 posts
AI 7 posts
AI Agent 2 posts
AI-based hotel recommendation 1 posts
AIForGood 1 posts
Alumni 60 posts
Animated Maps 1 posts
APIs 41 posts
Artificial Intelligence 2 posts
Artificial Intelligence 2 posts
AWS 13 posts
Banking 1 posts
Big Data 50 posts
Branch Analysis 1 posts
Capstone 206 posts
Career Education 7 posts
CLIP 1 posts
Community 72 posts
Congestion Zone 1 posts
Content Recommendation 1 posts
Cosine SImilarity 1 posts
Data Analysis 5 posts
Data Engineering 1 posts
Data Engineering 3 posts
Data Science 7 posts
Data Science News and Sharing 73 posts
Data Visualization 324 posts
Events 5 posts
Featured 37 posts
Function calling 1 posts
FutureTech 1 posts
Generative AI 5 posts
Hadoop 13 posts
Image Classification 1 posts
Innovation 2 posts
Kmeans Cluster 1 posts
LLM 6 posts
Machine Learning 364 posts
Marketing 1 posts
Meetup 144 posts
MLOPs 1 posts
Model Deployment 1 posts
Nagamas69 1 posts
NLP 1 posts
OpenAI 5 posts
OpenNYC Data 1 posts
pySpark 1 posts
Python 16 posts
Python 458 posts
Python data analysis 4 posts
Python Shiny 2 posts
R 404 posts
R Data Analysis 1 posts
R Shiny 560 posts
R Visualization 445 posts
RAG 1 posts
RoBERTa 1 posts
semantic rearch 2 posts
Spark 17 posts
SQL 1 posts
Streamlit 2 posts
Student Works 1687 posts
Tableau 12 posts
TensorFlow 3 posts
Traffic 1 posts
User Preference Modeling 1 posts
Vector database 2 posts
Web Scraping 483 posts
wukong138 1 posts

Our Recent Popular Posts

AI 4 AI: ChatGPT Unifies My Blog Posts
by Vinod Chugani
Dec 18, 2022
Meet Your Machine Learning Mentors: Kyle Gallatin
by Vivian Zhang
Nov 4, 2020
NICU Admissions and CCHD: Predicting Based on Data Analysis
by Paul Lee, Aron Berke, Bee Kim, Bettina Meier and Ira Villar
Jan 7, 2020

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 ChatGPT 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 football 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 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

NYC Data Science Academy

NYC Data Science Academy teaches data science, trains companies and their employees to better profit from data, excels at big data project consulting, and connects trained Data Scientists to our industry.

NYC Data Science Academy is licensed by New York State Education Department.

Get detailed curriculum information about our
amazing bootcamp!

Please enter a valid email address
Sign up completed. Thank you!

Offerings

  • HOME
  • DATA SCIENCE BOOTCAMP
  • ONLINE DATA SCIENCE BOOTCAMP
  • Professional Development Courses
  • CORPORATE OFFERINGS
  • HIRING PARTNERS
  • About

  • About Us
  • Alumni
  • Blog
  • FAQ
  • Contact Us
  • Refund Policy
  • Join Us
  • SOCIAL MEDIA

    © 2025 NYC Data Science Academy
    All rights reserved. | Site Map
    Privacy Policy | Terms of Service
    Bootcamp Application