From Machine Learning - Multiple Linear Regression - Using Python 3, how can I create a scatter plot given 2 independent variables (year, size)? And what would the linear regression formula be?

Hi Casey,You know that if you have a single predictor, then you'd have a 2D plot, right? There's 1 independent X and the a dependent Y (say 'size' and 'price' in that case). To plot 2 independent variables, you will need yet another dimension. In that case you will have a 3D plot. To achieve that you can use the mplot3d package (built on top of matplotlib). Here's the documentation: https://matplotlib.org/mpl_toolkits/mplot3d/So, assuming that your 2 independent variables are 'size' and 'year', while the dependent is 'price', for the specific exercise you are referring to, you can use code in the lines of: from mpl_toolkits.mplot3d import axes3d, Axes3Dfig = plt.figure()ax = plt.axes(projection="3d")ax.scatter3D(data['size'], data['year'], data['price'], cmap='hsv')plt.show() This will result into a 3D scatter plot:Now what if you want to plot a regression line? First, it is important to note that it is a regression line, in 2D, however, in 3D that would be a plane.To find that plane with statsmodels (and within the given example), we can write: results.params[0] + results.params[1]*data['size'] + results.params[2]*data['year'] Note that here, results.params contains the coefficients of the regression. There are 3 of them: constant, coef. for size and coef. for year. So the regression line becomes the abovementioned equation. Finally, we can incorporate that in the graph, using a 3D line (plane). from mpl_toolkits.mplot3d import axes3d, Axes3Dfig = plt.figure()ax = plt.axes(projection="3d") x_line = data['size']y_line = data['year']z_line = results.params[0] + results.params[1]*data['size'] + results.params[2]*data['year']ax.plot3D(x_line, y_line, z_line, 'orange')ax.scatter3D(data['size'], data['year'], data['price'], cmap='hsv')plt.show() And what you'll see is:Best,The 365 Team

Ok, I am almost there, thank you, but the following code displayed the same looking plane with no dots. What's a scatter plot with no dots?! Any ideas on why the dots are missing? fig = plt.figureax = plt.axes(projection="3d") x_line = data['size']y_line = data['year']z_line = results.params[0] + results.params[1]*data['size'] + results.params[2]*data['year']ax.plot3D(x_line, y_line, z_line, 'orange')plt.show()

Hi Casey,The scatter itself is created with the code: ax.scatter3D(data['size'], data['year'], data['price'], cmap='hsv') I believe you are missing this bit. This scatter represents the points. The rest of the code (the one you've implemented) draws the regression line. To get only the scatter, please refer to the first code cell in the original answer: fig = plt.figure()ax = plt.axes(projection="3d")ax.scatter3D(data['size'], data['year'], data['price'], cmap='hsv')plt.show() Hope this helps,The 365 Team

I had to add this from mpl_toolkits.mplot3d import axes3d, Axes3D

So I tried this because the code following the # leads to an IndexError: Index is out of bounds. Can you explain why that is? Also, the code below yields a different plane than yours, though I don't know how to post it fig = plt.figure() ax = plt.axes(projection="3d") x_line = data['size'] y_line = data['year'] z_line = data['price'] #z_line = results.params[0] + results.params[1]*data['size'] + results.params[2]*data['year'] ax.plot3D(x_line, y_line, z_line, 'orange') ax.scatter3D(data['size'], data['year'], data['price'], cmap='hsv') plt.show()