在熊猫中将`scipy.stats`函数与`apply()`一起使用时出错

时间:2019-10-15 21:32:25

标签: python pandas dataframe scipy pandas-apply

我正在尝试计算与数据框中各个条目相关的百分位数(使用列中的值分布)。我确定我缺少了一些基本的内容,但是无法弄清为什么在运行以下代码时出现错误,

from scipy.stats import percentileofscore as pctl
import pandas as pd
import numpy as np

data = np.arange(100).reshape(20,5)
df = pd.DataFrame(data)

def f(series):
    r= series.index
    return pctl(series.values, series.iloc[r])

df.apply(f)

这是我得到的错误,

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-7-4d3ad4c6f441> in <module>
----> 1 df.apply(f)

C:\Python\Miniconda\envs\leiap\lib\site-packages\pandas\core\frame.py in apply(self, func, axis, broadcast, raw, reduce, result_type, args, **kwds)
   6012                          args=args,
   6013                          kwds=kwds)
-> 6014         return op.get_result()
   6015 
   6016     def applymap(self, func):

C:\Python\Miniconda\envs\leiap\lib\site-packages\pandas\core\apply.py in get_result(self)
    316                                       *self.args, **self.kwds)
    317 
--> 318         return super(FrameRowApply, self).get_result()
    319 
    320     def apply_broadcast(self):

C:\Python\Miniconda\envs\leiap\lib\site-packages\pandas\core\apply.py in get_result(self)
    140             return self.apply_raw()
    141 
--> 142         return self.apply_standard()
    143 
    144     def apply_empty_result(self):

C:\Python\Miniconda\envs\leiap\lib\site-packages\pandas\core\apply.py in apply_standard(self)
    246 
    247         # compute the result using the series generator
--> 248         self.apply_series_generator()
    249 
    250         # wrap results

C:\Python\Miniconda\envs\leiap\lib\site-packages\pandas\core\apply.py in apply_series_generator(self)
    275             try:
    276                 for i, v in enumerate(series_gen):
--> 277                     results[i] = self.f(v)
    278                     keys.append(v.name)
    279             except Exception as e:

<ipython-input-6-347aa35ccd44> in f(series)
      1 def f(series):
      2     r= series.index
----> 3     return pctl(series.values, series.iloc[r])

C:\Python\Miniconda\envs\leiap\lib\site-packages\scipy\stats\stats.py in percentileofscore(a, score, kind)
   1785 
   1786     """
-> 1787     if np.isnan(score):
   1788         return np.nan
   1789     a = np.asarray(a)

C:\Python\Miniconda\envs\leiap\lib\site-packages\pandas\core\generic.py in __nonzero__(self)
   1574         raise ValueError("The truth value of a {0} is ambiguous. "
   1575                          "Use a.empty, a.bool(), a.item(), a.any() or a.all()."
-> 1576                          .format(self.__class__.__name__))
   1577 
   1578     __bool__ = __nonzero__

ValueError: ('The truth value of a Series is ambiguous. Use a.empty, a.bool(), a.item(), a.any() or a.all().', 'occurred at index 0')

2 个答案:

答案 0 :(得分:3)

问题是scipy.stats.percentileofscore(a, score, kind='rank')的值不是score参数的数组。

(正在发布NumPy解决方案,但这是不对的)该怎么做:

def f(series):

    return series.apply(lambda x:pctl(series,x) )

df.apply(f)

答案 1 :(得分:1)

根据文档,score参数应为int或float。参见https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.percentileofscore.html

series.iloc [r]是一个系列。您需要为分数使用标量值。