- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
UPDATE2:一个更好的标题(现在我明白了这个问题)应该是:在 scipy optimize.fmin 中输入的正确语法是什么?
更新:已请求可运行代码,因此函数定义已替换为可运行代码。示例输入数据已被硬编码为 numpy 数组“数据”。
我正在尝试用 scipy 优化一个函数,但我真的卡住了,必须寻求帮助。一个零长度数组被传递给优化器中的一个方法,我不明白为什么,也不知道如何解决这个问题。
这段代码试图做什么的简要概述:
运行此代码时出现的错误是: ValueError:大小为零的数组到没有身份的最大缩减操作
这是一段可运行的代码片段,它在我的机器上产生了错误。
#!/usr/bin/env python2.7
import numpy as np
from scipy import optimize
def p_of_r(m, r): ## this calculates p(r|m) for each datum r
r_range = np.arange(0, r+1, 1, dtype='int')
p_r = []
p_r = np.array([0.0 for a in r_range])
for x in r_range:
if x == 0:
p_r[x] = np.exp(-1 * m)
else:
total = 0.0
for y in np.arange(0, x, 1, dtype='int'):
current = ( p_r[y] ) / (x - y + 1)
total = current + total
p_r[x] = ( m / x ) * total
return p_r
def likelihood_function(m, *data): # calculates P(m|data) using entire data set
p_r = p_of_r(m, np.ma.max(data))
p_r_m = np.array([p_r[y] for y in data])
bigP = np.prod(p_r_m)
return bigP
def main():
data = np.array( [10, 10, 7, 19, 9, 23, 26, 7, 164, 16 ] )
median_r = np.median(data)
def Drake(m):
return median_r / m - np.log(m)
m_initial = optimize.broyden1(Drake, 1)
def helper(x, *args):
helper_value = -1 * likelihood_function(x, *args)
return helper_value
# here is the actual optimize.fmin
fmin_result = optimize.fmin(helper, x0=[m_initial], args=data)
print fmin_result
# for i in np.arange(0.0, 25.0, 0.1):
# print i, helper(i, data)
if __name__ == "__main__" : main()
错误本身:ValueError:大小为零的数组到没有身份的最大缩减操作
下面提供了回溯。
ValueError Traceback (most recent call last)
/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/IPython/utils/py3compat.pyc in execfile(fname, *where)
176 else:
177 filename = fname
--> 178 __builtin__.execfile(filename, *where)
/Users/deyler/bin/MSS-likelihood-minimal.py in <module>()
43 print fmin_result
44
---> 45 if __name__ == "__main__" : main()
/Users/deyler/bin/MSS-likelihood-minimal.py in main()
40
41
---> 42 fmin_result = optimize.fmin(helper, x0=[m_initial], args=data)
43 print fmin_result
44
/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/optimize/optimize.pyc in fmin(func, x0, args, xtol, ftol, maxiter, maxfun, full_output, disp, retall, callback)
371 'return_all': retall}
372
--> 373 res = _minimize_neldermead(func, x0, args, callback=callback, **opts)
374 if full_output:
375 retlist = res['x'], res['fun'], res['nit'], res['nfev'], res['status']
/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/optimize/optimize.pyc in _minimize_neldermead(func, x0, args, callback, xtol, ftol, maxiter, maxfev, disp, return_all, **unknown_options)
436 if retall:
437 allvecs = [sim[0]]
--> 438 fsim[0] = func(x0)
439 nonzdelt = 0.05
440 zdelt = 0.00025
/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/scipy/optimize/optimize.pyc in function_wrapper(*wrapper_args)
279 def function_wrapper(*wrapper_args):
280 ncalls[0] += 1
--> 281 return function(*(wrapper_args + args))
282
283 return ncalls, function_wrapper
/Users/deyler/bin/MSS-likelihood-minimal.py in helper(x, *args)
33 m_initial = optimize.broyden1(Drake, 1)
34 def helper(x, *args):
---> 35 helper_value = -1 * likelihood_function(x, *args)
36 return helper_value
37
/Users/deyler/bin/MSS-likelihood-minimal.py in likelihood_function(m, *data)
21
22 def likelihood_function(m, *data):
---> 23 p_r = p_of_r(m, np.ma.max(data))
24 p_r_m = np.array([p_r[y] for y in data])
25 bigP = np.prod(p_r_m)
/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/ma/core.pyc in max(obj, axis, out, fill_value)
5899 # If obj doesn't have a max method,
5900 # ...or if the method doesn't accept a fill_value argument
-> 5901 return asanyarray(obj).max(axis=axis, fill_value=fill_value, out=out)
5902 max.__doc__ = MaskedArray.max.__doc__
5903
/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/ma/core.pyc in max(self, axis, out, fill_value)
5159 # No explicit output
5160 if out is None:
-> 5161 result = self.filled(fill_value).max(axis=axis, out=out).view(type(self))
5162 if result.ndim:
5163 # Set the mask
/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/numpy/core/_methods.pyc in _amax(a, axis, out, keepdims)
8 def _amax(a, axis=None, out=None, keepdims=False):
9 return um.maximum.reduce(a, axis=axis,
---> 10 out=out, keepdims=keepdims)
11
12 def _amin(a, axis=None, out=None, keepdims=False):
ValueError: zero-size array to reduction operation maximum which has no identity
最佳答案
正确 fmin syntax是:
args: tuple, optional
Extra arguments passed to func, i.e. f(x,*args).
fmin_result = optimize.fmin(helper, x0=[m_initial], args=(data,))
下一个结果是预期的吗?
Optimization terminated successfully.
Current function value: -0.000000
Iterations: 16
Function evaluations: 32
[ 5.53610656]
关于python - SciPy optimize.fmin ValueError : zero-size array to reduction operation maximum which has no identity,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/20494719/
考虑一下,我们在具有整数弧容量的有向网络中有一个非整数最大流。 是否有算法可以将这个流量转化为整数最大流量? 它的运行时间是多少? 这不是作业问题。 最佳答案 如果您正在寻找具有积分弧容量的最大 s-
我有一个可以接受 byte[] 的 WCF 服务。我正在使用 HttpClient 创建客户端并收到以下错误。我在网上看到你必须在服务器和客户端上设置 readerQuotas,但是如何在 HttpC
我有一个如下所示的数据框: df = pd.DataFrame({'A':[100,300,500,600], 'B':[100,200,300,400],
我开始使用 cocoapi评估使用 Object Detection API 训练的模型。在阅读了解释平均精度 (mAP) 和召回率的各种来源后,我对 cocoapi 中使用的“最大检测”参数感到困惑
这个问题已经有答案了: 已关闭10 年前。 Possible Duplicate: Size-limited queue that holds last N elements in Java Java
您好,我需要帮助制作以下算法: 假设一个二维空间域,具有 xmax、xmin、ymin、ymax,空间中有“n~10,000”个点。 浏览点位置列表。 当框中的点数达到最大数量(假设为 10 个)时,
我在 Android 应用程序中有一个类,它包含一个字节数组作为纹理的数据源。我使用帧缓冲区将一些东西渲染到该纹理上,然后在屏幕上渲染该纹理。这非常有效。 但是,我只能使用 151 个纹理来执行此操作
系统调用中可以传递多少个参数?我检查了内核文件 /asm/unistd.h,没有看到包含超过 4 个参数的系统调用。 最佳答案 这取决于您使用的架构。对于 i386,系统调用号旁边最多有 6 个参数。
题目地址:https://leetcode-cn.com/problems/maximum-average-subtree/ 题目描述 Given the root of a binary tre
题目地址:https://leetcode.com/problems/sliding-window-maximum/ 题目描述 Given an array nums, there is a sl
题目地址:https://leetcode-cn.com/problems/maximum-distance-in-arrays/ 题目描述 Given m arrays, and each ar
题目地址:https://leetcode.com/problems/maximum-average-subarray-i/description/open in new window 题目描述
题目地址:https://leetcode.com/problems/third-maximum-number/description/ 题目描述 Given a non-empty array
我有一个很难重现的错误。另外,有人告诉我写日志文件是一种安全责任。所以我想在异常中尽可能多地捕捉。 我找不到任何说明 C# 异常的最大长度是多少的地方。 我想粘贴一条 XML 消息(1 或 2K),也
是否可以限制 QML 应用的最大 FPS? 我在低端 iten atom 硬件中获得 60FPS 和 30% 的 CPU 使用率使用 win32 Angle 驱动程序(openGL 软件无法使用),我
我已启用 SeriLog(最新版本)自记录功能,并且看到数百条消息说 Maximum destructuring depth reached 不知道这意味着什么以及这是否是我需要担心的问题。 有谁知道
我有一个具有两个属性的模型(我只显示两个,因为只需要这两列) 我的模型 place_id ---------------------- user_id 1 ----------------------
考虑以下情况:要设置一个群集,其中每台计算机都具有32G GB的RAM。和16个CPU核心 如何根据信息(32G GB RAM和16 CPU CORE)确定以下参数 ) yarn.scheduler.
完成如下函数定义: -- maxi x y returns the maximum of x and y 我应该不使用 Haskell 中的“最大”函数来完成这个任务。 maxi :: Int
Haskell 标准库中是否有与 maximum 等效的安全值? *Main Control.Monad.State Data.List> maximum [] *** Exception: Prel
我是一名优秀的程序员,十分优秀!