【问题标题】:How to plot the predicted value against all features of a dataset如何针对数据集的所有特征绘制预测值
【发布时间】:2021-08-05 11:28:46
【问题描述】:

我想可视化波士顿数据集的数据点,并绘制线性回归平面。但是,我收到一个值错误。我正在使用 colab。以下是我运行的代码。

import pandas as pd
import matplotlib.pyplot as plt
from sklearn.datasets import load_boston
from sklearn.model_selection import train_test_split
from sklearn.linear_model import LinearRegression

data1 = load_boston()
df1 = pd.DataFrame(data1.data)
df1.columns = data1.feature_names
df1['price']=data1.target

X = df1.iloc[:,0:13]
Y = df1.iloc[:,13]

xtrain,xtest,ytrain,ytest = train_test_split(X,Y,test_size = 0.2,random_state =0)

lin_reg = LinearRegression()
lin_reg.fit(xtrain,ytrain)
ypredict = lin_reg.predict(xtest)

plt.scatter(xtrain, ytrain, color = 'red')
plt.plot(xtrain, lin_reg.predict(xtrain), color = 'blue')

最后两行出现错误。

ValueError                                Traceback (most recent call last)

<ipython-input-54-4f89e0554222> in <module>()
----> 1 plt.scatter(xtrain, ytrain, color = 'red')
      2 plt.plot(xtrain, lin_reg.predict(xtrain), color = 'blue')

3 frames

/usr/local/lib/python3.7/dist-packages/matplotlib/axes/_axes.py in scatter(self, x, y, s, c, marker, cmap, norm, vmin, vmax, alpha, linewidths, verts, edgecolors, plotnonfinite, **kwargs)
   4389         y = np.ma.ravel(y)
   4390         if x.size != y.size:
-> 4391             raise ValueError("x and y must be the same size")
   4392 
   4393         if s is None:

ValueError: x and y must be the same size

我知道 X 有 13 列,而 Y 有 1 列。这就是错误显示的原因。但我不知道如何纠正它。

有人可以帮忙吗?

【问题讨论】:

    标签: python pandas matplotlib scikit-learn linear-regression


    【解决方案1】:
    • 必须单独绘制每个要素。
    • 请记住,'price' 是目标,因变量,lin_reg.predict(xtrain) 是训练数据的预测价格。
    # predicted price from xtrain
    ypred_train = lin_reg.predict(xtrain)
    
    # create the figure
    fig, axes = plt.subplots(ncols=4, nrows=4, figsize=(20, 20))
    
    # flatten the axes to make it easier to index
    axes = axes.flatten()
    
    # iterate through the column values, and use i to index the axes
    for i, v in enumerate(xtrain.columns):
        
        # seclect the column to be plotted
        data = xtrain[v]
        
        # plot the actual price against the features
        axes[i].scatter(x=data, y=ytrain, s=35, ec='white', label='actual')
        
        # plot predicted prices against the features
        axes[i].scatter(x=data, y=ypred_train, c='pink', s=20, ec='white', alpha=0.5, label='predicted')
    
        # set the title and ylabel
        axes[i].set(title=f'Feature: {v}', ylabel='price')
    
    # set a single legend
    axes[12].legend(title='Price', bbox_to_anchor=(1.05, 1), loc='upper left')
    
    # delete the last 3 unused axes
    for v in range(13, 16):
        fig.delaxes(axes[v])
    

    • 如果您要将所有内容都绘制成一个图,那么它会过于拥挤且毫无用处

    • 您还可以通过将df1 从宽格式融合到长格式来绘制带有seaborn.relplot 的所有数据。
      • 但是,在图形级图的顶部添加预测值更加困难。
    import seaborn as sns
    
    dfm = df1.melt(id_vars='price', value_vars=df1.columns[:-1], var_name='Feature')
    
    p = sns.relplot(kind='scatter', data=dfm, x='value', y='price', height=3,
                    col='Feature', col_wrap=4, facet_kws={'sharex': False})
    

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 1970-01-01
      • 2021-12-30
      • 1970-01-01
      • 2018-08-19
      • 2014-03-25
      • 2019-02-03
      • 1970-01-01
      • 1970-01-01
      相关资源
      最近更新 更多