How to make user defined functions for binned_statistic(如何为 binned_statistic 制作用户定义的函数)
问题描述
我正在使用 scipy stats 包沿轴获取统计信息,但我无法使用 binned_statistic 获取百分位数统计信息.我已经概括了下面的代码,我试图在一系列 x 箱中使用 x、y 值的数据集的第 10 个百分位数,但它失败了.
I am using scipy stats package to take statistics along the an axis, but I am having trouble taking the percentile statistic using binned_statistic
. I have generalized the code below, where I am trying taking the 10th percentile of a dataset with x, y values within a series of x bins, and it fails.
我当然可以使用 np.std
进行函数选项,例如中值,甚至是 numpy 标准差.但是,我无法弄清楚如何使用 np.percentile
因为它需要 2 个参数(例如 np.percentile(y, 10)
),但是它给了我一个 <代码>ValueError: 统计不理解 错误.
I can of course do function options, like median, and even the numpy standard deviation using np.std
. However, I cannot figure out how to use np.percentile
because it requires 2 arguments (e.g. np.percentile(y, 10)
), but then it gives me a ValueError: statistic not understood
error.
import numpy as np
import scipy.stats as scist
y_median = scist.binned_statistic(x,y,statistic='median',bins=20,range=[(0,5)])[0]
y_std = scist.binned_statistic(x,y,statistic=np.std,bins=20,range=[(0,5)])[0]
y_10 = scist.binned_statistic(x,y,statistic=np.percentile(10),bins=20,range=[(0,5)])[0]
print y_median
print y_std
print y_10
我不知所措,甚至尝试过像这样的用户定义函数,但没有运气:
I am at a loss and have even played around with user defined functions like this, but with no luck:
def percentile10():
return(np.percentile(y,10))
任何帮助,不胜感激.
谢谢.
推荐答案
你定义的函数的问题是它根本没有参数!它需要一个 y
与您的样本相对应的参数,如下所示:
The problem with the function you defined is that it takes no arguments at all! It needs to take a y
argument that corresponds to your sample, like this:
def percentile10(y):
return(np.percentile(y,10))
为了简洁起见,您也可以使用 lambda
函数:
You could also use a lambda
function for brevity:
scist.binned_statistic(x, y, statistic=lambda y: np.percentile(y, 10), bins=20,
range=[(0, 5)])[0]
这篇关于如何为 binned_statistic 制作用户定义的函数的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持编程学习网!
本文标题为:如何为 binned_statistic 制作用户定义的函数


- CTR 中的 AES 如何用于 Python 和 PyCrypto? 2022-01-01
- 计算测试数量的Python单元测试 2022-01-01
- 我如何透明地重定向一个Python导入? 2022-01-01
- YouTube API v3 返回截断的观看记录 2022-01-01
- 使用公司代理使Python3.x Slack(松弛客户端) 2022-01-01
- 检查具有纬度和经度的地理点是否在 shapefile 中 2022-01-01
- ";find_element_by_name(';name';)";和&QOOT;FIND_ELEMENT(BY NAME,';NAME';)";之间有什么区别? 2022-01-01
- 如何使用PYSPARK从Spark获得批次行 2022-01-01
- 我如何卸载 PyTorch? 2022-01-01
- 使用 Cython 将 Python 链接到共享库 2022-01-01