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 > Predicting Higgs Boson Signal

Predicting Higgs Boson Signal

Danli, Diego De Lazzari and Linlin Cheng
Posted on Sep 2, 2016

The Higgs boson was the last hold-out particle remaining hidden during the quest to check the accuracy of the Standard Model of Physics. After decades of search, on March 14, 2013 CERN officially announced the confirmation of the Higgs boson. This Kaggle competition is about building a machine learning algorithm to filter out signal from background noise using selected & simulated ATLAS experiment data.

The code can be found here.


The Missingness Pattern

The first thing caught our eyes while reading the data was the enormous amount of -999s in both the training and the test datasets. So we did a quick check using the aggr function inside the VIM package.

The graph below shows the missingness pattern of the training dataset. To the left is the ratio of missingness for each column, from .15 to more than .70 (The x ticks might be a little bit hard to identify but don't worry we'll talk about them later).  And to the right, each row represents a possible combination of missingness in terms of variables.  From this graph alone we can tell that there are several variables acting weird. There's a group of 7 variables that always be missing together, and sometimes a group of other 3 variables follows. The only variable that is potentially MCAR ( Missing completely at random) is DER_mass_MMC (The narrow red bar at the very left of the combination graph).

Screen Shot 2016-09-04 at 12.52.03 AM

Since the ticks generated in previous graph was intentionally minimized to fit the space thus hard to read, we have double checked the missingness and discovered an interesting fact, The missingness pattern has a strong correlation with the value of variable PRI_jet_num! When it's equal to 0 or 1, there are bunch of variables be missing completely, and for the variable DER_mass_MMC, which we suspected to be MCAR, it has small amount of missingness in all of the cases.

Picture1


Data Distribution

Next we decided to take a look at some of the EDAs, namely data distribution and correlation, to further explore the structure of our data. What been shown below is the box plots for all the variables, we found that the means of some them are substantially different for background and signal, also many of the features have outliers, however, we finally decided only to take out the extreme ones (those ones which are in the red circles).

Screen Shot 2016-09-04 at 7.04.14 AM

As of correlation, we kicked out all the NAs (i.e. -999s) before drawing the graph, otherwise many of the features would appear highly correlated. And the final graph shows some multicollinearity, which we think is acceptable since a lot of the features are actually derived from others (those ones with names starting with a "DER"), but we kept this graph just in case we would use any linear models in the future.

Screen Shot 2016-09-04 at 7.16.27 AM


Feature Selection

We also did a PCA (Principal Component Analysis) to help us bringing down some dimensions, or aid our feature selection.

Since PCA does not take NAs as input, we chose to do a random imputation before the analysis, simply because the box plots told us many of the variables are highly skewed and mean imputation might not be a good choice.

The scree plot suggests 11 PCs, but looking at the result to the right we found it's impossible to use eigenvectors to predict since the 11 of them can only explain 70% of the variance, and in terms of feature selection, all of the variables are contributing at least something to at least one of the PCs, so we can not dump any of them frivolously. Thus we moved forward quickly to other solutions.

Screen Shot 2016-09-04 at 7.28.09 AM Screen Shot 2016-09-04 at 7.28.25 AM

There's another method to examine feature importance which is quite popular, the tree algorithms. We fitted a random forest for complete cases and plotted out the importance graph(the little yellow thunderbolt indicates that variable has missing values). We can see the DER_mass_MMC really popped out in terms of importance, however, what we consider even more valuable is that, many of the missing-in-bulk variables are not trivial at all. So our data cleaning strategy finally came down to:

  • Kick out the extreme outliers
  • For small amount of missingness, we do random imputation
  • Keep as many variables as we can, unless it really doesn't make sense to use them

Screen Shot 2016-09-04 at 7.39.16 AM


Modeling - First & Second Try 

We subsetted the training data into 3 pieces by the value of PRI_jet_num. When it's 0 and 1, deleted those variables that are completely missing, when it's 2 and 3, imputed the small portion of NAs, and for each subsets of data, we tuned 3 models, namely random forest, XGBoost and AdaBoost (we used only tree based models considering the multicollinearity mentioned above and lack of computing time for complex algorithms like SVM and neural network), ensemble them using caretStack.

When we test, we also need to subset the data and plug into 3 different models.

Screen Shot 2016-09-04 at 7.52.44 AM

We actually did a simple ensemble with only AdaBoost and XGBoost before stacking everything together.

  • ADABOOST (iter=100, nu =.03, maxdepth = 10. Accuracy after threshold = 80%, Best threshold at p = 0.75 ) on the full raw dataset => AMS 3.45
  • XGBOOST (usual parameters, accuracy after threshold 80%, ) => AMS 3.55
  • ENSEMBLE: (ADABOOST + XGBOOST)/2 (average on predictions (unknown data) , averaged and then rounded. Threshold at p = 0.8) => AMS 3.57

By simply taking an average of the 2 models' predictions, we achieved an AMS (approximate median significance, which is the official evaluation metric of this competition) score of 3.57. However, after we stacked the models using a meta-model GBM (Gradient Boosting Machine), the score bounced back to 3.26.


Modeling - Third Try

We noticed that something might be wrong with our stacking method, thus we changed into another one. Manually stacking meta features.

Screen Shot 2016-09-04 at 8.51.12 AM

About Authors

Danli

Danli Zeng is a young professional with 5 years' experience in MARKETING and MEDIA. She was specialized in integrated media planning and ROMI analysis for FMCG industry. Having worked on all sides of agency, media and client, she...
View all posts by Danli >

Diego De Lazzari

Researcher, developer and data scientist. Diego De Lazzari is an applied physicist with a rather diverse background. He spent 8 years in applied research, developing computational models in the field of Plasma Physics (Nuclear Fusion) and Geophysics. As...
View all posts by Diego De Lazzari >

Linlin Cheng

View all posts by Linlin Cheng >

Leave a Comment

Cancel reply

You must be logged in to post a comment.

No comments found.

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