胖尾分布的Levene检验

时间:2019-09-25 18:01:27

标签: python-3.x scipy statistics

我正在尝试对表示经验数据频率分布的两个值数组进行Levene检验。我使用scipy.stats levene package,并想添加'trimmed'参数,因为我的数据是胖尾的,并且提到这对于此类数据是最好的。为了记录,我也将数据标准化。分布看起来像this,您可以从here下载数据

从我尝试使用'trimmed'参数的文档中:

levene(list(group),list(group2),'trimmed')[1]

当不添加'trimmed'时可以使用,但是我对结果不满意:

levene(list(group),list(group2))[1]

我收到以下错误消息

> Traceback (most recent call last):   File "C:\Program
> Files\JetBrains\PyCharm Community Edition
> 2019.2.1\helpers\pydev\pydevd.py", line 2060, in <module>
>     main()   File "C:\Program Files\JetBrains\PyCharm Community Edition 2019.2.1\helpers\pydev\pydevd.py", line 2054, in main
>     globals = debugger.run(setup['file'], None, None, is_module)   File "C:\Program Files\JetBrains\PyCharm Community Edition
> 2019.2.1\helpers\pydev\pydevd.py", line 1405, in run
>     return self._exec(is_module, entry_point_fn, module_name, file, globals, locals)   File "C:\Program Files\JetBrains\PyCharm Community
> Edition 2019.2.1\helpers\pydev\pydevd.py", line 1412, in _exec
>     pydev_imports.execfile(file, globals, locals)  # execute the script   File "C:\Program Files\JetBrains\PyCharm Community Edition
> 2019.2.1\helpers\pydev\_pydev_imps\_pydev_execfile.py", line 18, in execfile
>     exec(compile(contents+"\n", file, 'exec'), glob, loc)   File "D:/1. Study/DELFT/Thesis/Thesis!/07 -
> Workspace/Kick_off/Pycharm/Efficiency T-Tests.py", line 108, in
> <module>
>     row2.append(levene(list(group),list(group2),'trimmed')[1])   File "C:\Users\Acer\Anaconda3\envs\MSc_Thesis\lib\site-packages\scipy\stats\morestats.py",
> line 2320, in levene
>     Yci[j] = func(args[j])   File "C:\Users\Acer\Anaconda3\envs\MSc_Thesis\lib\site-packages\scipy\stats\morestats.py",
> line 2310, in <lambda>
>     func = lambda x: np.median(x, axis=0)   File "C:\Users\Acer\Anaconda3\envs\MSc_Thesis\lib\site-packages\numpy\lib\function_base.py",
> line 3497, in median
>     overwrite_input=overwrite_input)   File "C:\Users\Acer\Anaconda3\envs\MSc_Thesis\lib\site-packages\numpy\lib\function_base.py",
> line 3385, in _ureduce
>     axis = _nx.normalize_axis_tuple(axis, nd)   File "C:\Users\Acer\Anaconda3\envs\MSc_Thesis\lib\site-packages\numpy\core\numeric.py",
> line 1607, in normalize_axis_tuple
>     axis = tuple([normalize_axis_index(ax, ndim, argname) for ax in axis])   File
> "C:\Users\Acer\Anaconda3\envs\MSc_Thesis\lib\site-packages\numpy\core\numeric.py",
> line 1607, in <listcomp>
>     axis = tuple([normalize_axis_index(ax, ndim, argname) for ax in axis]) numpy.AxisError: axis 0 is out of bounds for array of dimension
> 0

似乎由于某种原因,阵列被展平了。谁能提出不同的建议?谢谢。

0 个答案:

没有答案