时间戳

首页 > 解决方案 > 使用 matplotlib 将数据绘制为 3D 线性模型

python - 使用 matplotlib 将数据绘制为 3D 线性模型

问题描述

我想在下面的链接中绘制一些数据。

当我在 OLS 公式中有多个“电视和广播”并且我只想用“销售”来绘制这两个时，我该怎么办？因为如果我喜欢下面的代码（链接），它会向我显示其他未定义的错误（电视和广播除外）..

感谢帮助！

https://stackoverflow.com/a/26434204/14208684

这是链接的代码：

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
import statsmodels.formula.api as sm
from matplotlib import cm

csv = pd.read_csv('http://www-bcf.usc.edu/~gareth/ISL/Advertising.csv', index_col=0)
model = sm.ols(formula='Sales ~ TV + Radio', data = csv)
fit = model.fit()

fit.summary()

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

x_surf = np.arange(0, 350, 20)                # generate a mesh
y_surf = np.arange(0, 60, 4)
x_surf, y_surf = np.meshgrid(x_surf, y_surf)

exog = pd.core.frame.DataFrame({'TV': x_surf.ravel(), 'Radio': y_surf.ravel()})
out = fit.predict(exog = exog)
ax.plot_surface(x_surf, y_surf,
                out.reshape(x_surf.shape),
                rstride=1,
                cstride=1,
                color='None',
                alpha = 0.4)

ax.scatter(csv['TV'], csv['Radio'], csv['Sales'],
           c='blue',
           marker='o',
           alpha=1)

ax.set_xlabel('TV')
ax.set_ylabel('Radio')
ax.set_zlabel('Sales')

plt.show()

标签： pythonpandasmatplotliblinear-regression

解决方案

推荐阅读