Я пытаюсь выполнить тест Левена на двух массивах значений, представляющих частотное распределение эмпирических данных.Я использую scipy.stats levene package и хотел бы добавить параметр 'trimmed', так как мои данные толстые, и упоминается, что это было бы лучше для таких данных.Для записи я также нормализовал данные.Распределения выглядят как this , и вы можете загрузить данные с здесь
Из документов, которые я пытался использовать с помощью параметра trimmed:
levene(list(group),list(group2),'trimmed')[1]
Когда не добавляется «обрезано», это работает, но я не уверен в результатах:
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
Похоже, что по какой-то причине массив уплощен.Кто-нибудь может подсказать, что делать по-другому?Спасибо.