diff --git a/examples/regression.py b/examples/regression.py
index ee06ade..50ed0dc 100644
--- a/examples/regression.py
+++ b/examples/regression.py
@@ -1,4 +1,3 @@
-from sklearn.datasets import load_boston
 from sklearn.metrics import mean_squared_error
 from sklearn.model_selection import train_test_split
 
@@ -7,7 +6,12 @@
 
 if __name__ == "__main__":
 
-    X, Y = load_boston(return_X_y=True)
+    #Load Boston housing dataset
+    data_url = "http://lib.stat.cmu.edu/datasets/boston"
+    raw_df = pd.read_csv(data_url, sep="\s+", skiprows=22, header=None)
+    X = np.hstack([raw_df.values[::2, :], raw_df.values[1::2, :2]])
+    Y = raw_df.values[1::2, 2]
+    
     X_train, X_test, Y_train, Y_test = train_test_split(X, Y, test_size=0.2)
 
     ngb = NGBRegressor(Dist=Normal).fit(X_train, Y_train)