1

我在 statsmodels 中使用 Tukeys 成对比较来可视化组之间的差异。当我像这样运行代码时，代码运行良好：

from statsmodels.stats.multicomp import pairwise_tukeyhsd
from matplotlib import pyplot as plt 
import numpy as np 
import matplotlib

%matplotlib notebook
fig,ax = plt.subplots()
T = np.asanyarray([1,2,0,1,2,3,4,5,10])
N = np.asanyarray(['a','a','a','b','b','b','c','c','c'])

tk = pairwise_tukeyhsd(T,N)
tk.plot_simultaneous(ax=ax)

它输出 在此处输入图像描述

但是，当我尝试使用 compare_name ='c' 时：

tk.plot_simultaneous(comparison_name ='c',ax=ax)

它失败了，我得到以下回溯：

TypeError                                 Traceback (most recent call last)
<ipython-input-10-9d84068180f9> in <module>()
     19 
     20 tk = pairwise_tukeyhsd(T,N)
---> 21 tk.plot_simultaneous(comparison_name ='c',ax=ax)

C:\Anaconda3\lib\site-packages\statsmodels\sandbox\stats\multicomp.py in plot_simultaneous(self, comparison_name, ax, figsize, xlabel, ylabel)
    733                 if self.groupsunique[i] == comparison_name:
    734                     continue
--> 735                 if (min(maxrange[i], maxrange[midx]) -
    736                                          max(minrange[i], minrange[midx]) < 0):
    737                     sigidx.append(i)

TypeError: only integer scalar arrays can be converted to a scalar index

我怎样才能解决这个问题？

4

0 回答 0