Machine Learning Project -House Prices: Advanced Regression Techniques

Posted on Mar 25, 2018

In 2016, Kaggle released a competition called House Prices: Advanced Regression Techniques. The goal of the competition was to predict the final sale price of homes in Ames, Iowa. The dataset itself came with 79 explanatory variables describing just about every aspect of a residential home.


To tackle the competition, I used various techniques I learned from the Bootcamp, as well as any insight gained from investigating Kaggler’s kernels.  Using the workflow depicted below, I was able to make predictions of my own.

Most of the work done dealt with transforming, cleaning, imputing, and aggregating data. Once complete, models were made using the following machine learning algorithms.

  • Lasso Regularization
  • Ridge Regularization
  • ElasticNet
  • Gradient Boosting Regression

Data Structure and Data Preprocessing


A histogram plot shows the distribution of the target variable ‘SalePrice’ as being was right-skewed. Before moving any further, I decided to obtain a normal distribution by way of log-transformation.

To get a better understanding of missingness with the given data sets (train and test), I determined the percentage of missing values and built bar plots.

In the train dataset, 19 out of 81 features had missing values. As for the test dataset, 33 out of 80 features had missing values.

Some of the features had outliers.  Looking at the graph below we can see that there are two homes that are very spacious yet are extremely low in price. When making my model I must ensure that they are robust to outliers.

To check for multicollinearity, I created a correlation matrix with respect to the target variable ‘SalePrice’.

Some variables were highly correlated such as 'GarageArea' and 'GarageCars'.  This makes sense since the size of a garage determines('GarageArea') the number of cars that fit in it ('GarageCars'). Other variables that were highly correlated show the same type of dependency.


As you can see the categorical features ‘Alley, ‘PoolQC’, and ‘MscFeature’ have over 90% of its column missing. However, looking at the data description, the NA’s in categorical variables actually means “not present”.  For example, NA in the ‘Alley’ variable meant that the home has no alley.

To correct the issue, I first concatenated the train and test set. I then replaced the Na’s with “None”. Researching further, I realized that two of the numerical features (‘MSSubClass’, ‘Mosold’) were actually categorical.  ‘MoSold’ used numerical values that represented months. For example, “1” meant January, and “2” meant February. In ‘MSSubClass’, numerical values identified the type of dwelling involved in the sale. Both features were converted into strings. For the remaining numerical features, NA’s were replaced by the value zero.

Feature Engineering

Since the data has many categorical variables, I converted them into dummy/indicator variables. As for numerical features, quite a few had a highly skewed distribution which can lead to poor models. To combat this, these features were log-transformed.

Ridge Regression

Regularization reduces the magnitude of the features’ coefficients to improve the accuracy of a model. By shrinking parameters, we prevent multicollinearity. With ridge regression, an L2 penalty is added to the loss functions to implement the reduction. Its’ tuning parameter lambda controls the relative impact of the coefficient estimates. Using a grid of values, I determined the best choice for lambda (alpha in the ridge function) was 0.0005.

Lasso Regression

Lasso has an L1 penalty which allows coefficients to shrink to zero one by one to reduce their variance. Hence, lasso can be used to select the most important features. In this case,101 features were chosen and the other 119 features were eliminated. “GrLivArea” which represents the above-ground living area seems to be the central feature in determining the sale price of a home.


ElasticNet combines L1 with L2 regularization to create a linear regression model to correct overfitting. Here the alpha tuning parameter sets the ratio between them. It is also known to be robust to outliers. This performed slightly better than lasso and ridge when submitting predictions to lasso

Gradient Boosting Regression

Gradient boosting uses weak prediction models and combines them to make a strong “learner”. The loss was set to “huber” for it to be robust to outliers. My model received an R squared core of .98 which means that 98% of the variance in the target can be explained by the explanatory variables.


The advantage of ridge and lasso is that regularization helps to balance between overfitting and underfitting. It's easy to assume that elasticnet would the top choice for modeling.  However, after many submissions and fine tuning, xgboost had the best predictions.

As the final prediction, I decided to average the predictions of lasso, elasticnet, and gradientboostingregressor. Model ensembling is known to increase accuracy on a variety of machine learning tasks. As a result, I was able to land in the top 16% of 4074 Kaggle teams. Meticulously prepping the data through imputation and feature engineering prior to modeling ensured good results. Furthermore, fine-tuning parameters of the model help prevent overfitting.

If you would like more information, visit my GitHub for more details.

About Author

Related Articles

Leave a Comment

No comments found.

View Posts by Categories

Our Recent Popular Posts

View Posts by Tags

#python #trainwithnycdsa 2019 airbnb Alex Baransky alumni Alumni Interview Alumni Reviews Alumni Spotlight alumni story Alumnus API Application artist aws 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 Bundles California Cancer Research capstone Career Career Day citibike clustering Coding Course Demo Course Report D3.js data Data Analyst data science Data Science Academy Data Science Bootcamp Data science jobs Data Science Reviews Data Scientist Data Scientist Jobs data visualization Deep Learning Demo Day Discount dplyr employer networking feature engineering Finance Financial Data Science Flask gbm Get Hired ggplot2 googleVis Hadoop higgs boson Hiring hiring partner events Hiring Partners Industry Experts Instructor Blog Instructor Interview Job Job Placement Jobs Jon Krohn JP Morgan Chase Kaggle Kickstarter lasso regression Lead Data Scienctist Lead Data Scientist leaflet linear regression Logistic Regression machine learning Maps matplotlib Medical Research Meet the team meetup Networking neural network Neural networks New Courses nlp NYC NYC Data Science nyc data science academy NYC Open Data NYCDSA NYCDSA Alumni Online Online Bootcamp Online Training Open Data painter pandas Part-time Portfolio Development prediction Prework Programming PwC python python machine learning python scrapy python web scraping python webscraping Python Workshop R 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 Selenium sentiment analysis Shiny Shiny Dashboard Spark Special Special Summer Sports statistics streaming Student Interview Student Showcase SVM Switchup Tableau team TensorFlow Testimonial tf-idf Top Data Science Bootcamp twitter visualization web scraping Weekend Course What to expect word cloud word2vec XGBoost yelp