怎样在python中做多元回归_如何在python中绘制多元回归三维图
我不是科学家,所以请假设我不知道有经验的程序员的行话,也不知道复杂的科学绘图技术。Python是我唯一知道的语言(初学者+或者中级)。
任务:将多元回归(z=f(x,y))的结果绘制为三维图形上的二维平面(例如,我可以使用OSX的绘图实用程序,也可以使用R在这里实现)。
经过一周的搜索和阅读matplotlib、seaborn和mayavi的各种文档,我终于找到了Simplest way to plot 3d surface given 3d points,这听起来很有前途。这是我的数据和代码:
首先尝试matplotlib:shape: (80, 3)
type:
zmul:
[[ 0.00000000e+00 0.00000000e+00 5.52720000e+00]
[ 5.00000000e+02 5.00000000e-01 5.59220000e+00]
[ 1.00000000e+03 1.00000000e+00 5.65720000e+00]
[ 1.50000000e+03 1.50000000e+00 5.72220000e+00]
[ 2.00000000e+03 2.00000000e+00 5.78720000e+00]
[ 2.50000000e+03 2.50000000e+00 5.85220000e+00]
……]
import matplotlib
from matplotlib.ticker import MaxNLocator
from matplotlib import cm
from numpy.random import randn
from scipy import array, newaxis
Xs = zmul[:,0]
Ys = zmul[:,1]
Zs = zmul[:,2]
surf = ax.plot_trisurf(Xs, Ys, Zs, cmap=cm.jet, linewidth=0)
fig.colorbar(surf)
ax.xaxis.set_major_locator(MaxNLocator(5))
ax.yaxis.set_major_locator(MaxNLocator(6))
ax.zaxis.set_major_locator(MaxNLocator(5))
fig.tight_layout()
plt.show()
我只得到一个空的三维坐标框架,其中包含以下错误消息:
运行时错误:qhull Delaunay三角剖分计算中的错误:奇异输入数据(exitcode=2);使用python verbose选项(-v)查看原始qhull错误。
使用mayavi进行第二次尝试:
相同的数据,分成3个numpy数组:type:
X: [ 0 500 1000 1500 2000 2500 3000 ….]
type:
Y: [ 0. 0.5 1. 1.5 2. 2.5 3. ….]
type:
Z: [ 5.5272 5.5922 5.6572 5.7222 5.7872 5.8522 5.9172 ….]
代码:from mayavi import mlab
def multiple3_triple(tpl_lst):
X = xs
Y = ys
Z = zs
# Define the points in 3D space
# including color code based on Z coordinate.
pts = mlab.points3d(X, Y, Z, Z)
# Triangulate based on X, Y with Delaunay 2D algorithm.
# Save resulting triangulation.
mesh = mlab.pipeline.delaunay2d(pts)
# Remove the point representation from the plot
pts.remove()
# Draw a surface based on the triangulation
surf = mlab.pipeline.surface(mesh)
# Simple plot.
mlab.xlabel("x")
mlab.ylabel("y")
mlab.zlabel("z")
mlab.show()
我只得到这个:
如果这很重要,我将在OSX 10.9.3上使用64位版本的Enthough's Canopy
对于我做错了什么会很感激。
编辑:张贴最终的代码,以防有人帮助。'''After the usual imports'''
def multiple3(tpl_lst):
mul = []
for tpl in tpl_lst:
calc = (.0001*tpl[0]) + (.017*tpl[1])+ 6.166
mul.append(calc)
return mul
fig = plt.figure()
ax = fig.gca(projection='3d')
'''some skipped code for the scatterplot'''
X = np.arange(0, 40000, 500)
Y = np.arange(0, 40, .5)
X, Y = np.meshgrid(X, Y)
Z = multiple3(zip(X,Y))
surf = ax.plot_surface(X, Y, Z, rstride=1, cstride=1,cmap=cm.autumn,
linewidth=0, antialiased=False, alpha =.1)
ax.set_zlim(1.01, 11.01)
ax.set_xlabel(' x = IPP')
ax.set_ylabel('y = UNRP20')
ax.set_zlabel('z = DI')
ax.zaxis.set_major_locator(LinearLocator(10))
ax.zaxis.set_major_formatter(FormatStrFormatter('%.02f'))
fig.colorbar(surf, shrink=0.5, aspect=5)
plt.show()
怎样在python中做多元回归_如何在python中绘制多元回归三维图相关推荐
- python如何做动画_如何在Python中动画动态系统的运动?
假设你有一个4个微分方程,它代表了一个动态系统(推车上的钟摆)的运动,你用来求解这些方程scipy.integrate.odeint持续10秒,间隔0.01秒.在 最后得到尺寸为(1000,4)的解矩 ...
- python怎么做回归分析_如何在Python中进行二维回归分析?
这是使用scipy.optimize.curve_fit拟合表面的示例Python代码,它使原始数据生成3D散点图,对错误进行3D散点图绘制,绘制表面图和轮廓图.更改它以使用您自己的数据和功能,您应该 ...
- python3提取字符串中的数字_如何在Python中从字符串中提取数字?
14 回复 | 直到 1 年前 1 430 3 年前 如果只想提取正整数,请尝试以下操作: >>> str = "h3110 23 cat 444.4 rabbit 11 ...
- python 指定证书验证_如何在python中验证SSL证书?
我需要验证我的自定义CA签署了证书.使用OpenSSL命令行实用程序很容易做到: # Custom CA file: ca-cert.pem # Cert signed by above CA: bo ...
- python定义一个列表_如何在Python中创建用户定义的列表?
你需要将怪物的数量发送给使用它的函数:def welcome(): number = monsters() print('Alright, ' + str(number) + ' monsters.' ...
- python隐藏启动台_如何在Python中启动后台进程?
如何在Python中启动后台进程? 我正在尝试将shell脚本移植到更易读的python版本. 原始shell脚本在后台使用"&"启动多个进程(实用程序,监视器等). 如何 ...
- python 运行r语言_如何在R中运行Python
python 运行r语言 尽管我很喜欢R,但很显然Python还是一种很棒的语言-既适用于数据科学又适用于通用计算. R用户想要在Python中做一些事情可能有充分的理由. 也许这是一个很棒的库,还没 ...
- python列表去空值_如何在Python列表中的列表中删除nan / null值? - python
因此,我有一个具有NaN值的数据框,并将该数据框中的所有行转换为一个列表,然后将其添加到另一个列表中. Index 1 2 3 4 5 6 7 8 9 10 ... 71 72 73 74 75 76 ...
- python的loc函数_如何在pandas中使用loc、iloc函数进行数据索引(入门篇)
在数据分析过程中,很多时候我们需要从数据表中提取出我们需要的部分,而这么做的前提是我们需要先索引出这一部分数据.今天我们就来探索一下,如何在pandas中使用loc函数和iloc函数索引数据. 今天我 ...
最新文章
- 理解C#中的string类型
- LiveVideoStack线上分享第五季(七):开源流媒体服务器:为何一定得再撸个新的...
- c语言 pow优化_c程序代码优化的一些方法
- 国内首家,腾讯云云开发“全家桶”来了
- MongoDB入门简介
- 电子邮件乱码的解决方法
- [Linux]gocron定时任务平台的部署
- 单片机ADC采样算法----限幅平均滤波法
- Python Cheat Sheet 中文版
- js 实现文件上传 php,JS+php后台实现文件上传功能详解
- 计算机显示器刷新率怎么调,显示器刷新率怎么超频?电脑显示器提高屏幕刷新率超频教程...
- ASP入门教程 1小时ASP入门,非常简单
- 帕斯卡计算机的控制原理,液体静力学基本方程式和帕斯卡原理
- 记一次很坑很坑的报错java.lang.Exception: The class is not public.
- php 美团配送 生成签名
- Python+OpenCV实用案例应用教程:建立自定义物体检测器
- 文件系统测试工具整理
- python3.0如何画表格_怎么用python画表格?
- Linux查看pcie槽位
- 方方格子access_Quick Access-Quick Access(项目路径管理AE脚本)下载 v1.21官方版--pc6下载站...
热门文章
- kswapd0进程对于CPU占有率高的情况下排查到黑客植入脚本,与黑客斗智斗勇的三个回合
- linux 命令行 过滤,Linux常用文本过滤命令(find grep wc awk sed sort uniq split)
- 智慧农业项目建设体系之数字化育种体系建设
- 双模sa_为什么NSA/SA双模的说法是错误的?为什么5G出来后你觉得4G慢了?
- 2011MacbookPro win7安装
- iframe嵌套页面
- FAT文件系统原理的详细分析
- 解决aspera 下载ENA 数据报错Session Stop (Error: Client unable to connect to server (check UDP port and fire
- 智慧城市安全体系实例——智能门锁安全吗?
- 反爬虫兵法演绎04 _ 爬虫的首轮攻势:如何低调地拿到自己想要的数据?