Question

我正在尝试使用statsmodels模块执行多个样本比较和Tukey HSD，但我不断收到此错误消息，“当p> = .9时，ValueError：v必须大于1。”我曾尝试在互联网上查找可能的解决方案，但无济于事。任何熟悉此模块的人都有可能帮助我破译我做错了什么以提示此错误。我使用Python 2.7x和spyder。以下是我的数据和打印语句的示例。谢谢！

import numpy as np
from statsmodels.stats.multicomp import (pairwise_tukeyhsd,MultiComparison)

###--- Here are the data I am using:
data1 = np.array([ 1,     1,     1,     1,   976,    24,     1,     1,    15, 15780])
data2 = np.array(['lau15', 'gr17', 'fri26', 'bays29', 'dantzig4', 'KAT38','HARV50', 'HARV10', 'HARV20', 'HARV41'], dtype='|S8')

####--- Here's my print statement code: 
print pairwise_tukeyhsd(data1, data2, alpha=0.05)

Answer 1

为了使测试有效，您似乎需要提供比每组单个观察更多的数据。

最小示例：

from statsmodels.stats.multicomp import pairwise_tukeyhsd,MultiComparison

data=[1,2,3]
groups=['a','b','c']

print("1st try:")
try:
        print(pairwise_tukeyhsd(data,groups, alpha=0.05))
except ValueError as ve:
        print("whoops!", ve)


data.append(2)
groups.append('a')
print("2nd try:")
try:
        print( pairwise_tukeyhsd(data, groups, alpha=0.05))
except ValueError as ve:
        print("whoops!", ve)

输出：

1st try:
/home/user/.local/lib/python3.7/site-packages/numpy/core/fromnumeric.py:3367: RuntimeWarning: Degrees of freedom <= 0 for slice
  **kwargs)
/home/user/.local/lib/python3.7/site-packages/numpy/core/_methods.py:132: RuntimeWarning: invalid value encountered in double_scalars
  ret = ret.dtype.type(ret / rcount)
whoops! v must be > 1 when p >= .9
2nd try:
Multiple Comparison of Means - Tukey HSD, FWER=0.05 
====================================================
group1 group2 meandiff p-adj  lower    upper  reject
----------------------------------------------------
     a      b      0.5   0.1  -16.045  17.045  False
     a      c      1.5   0.1  -15.045  18.045  False
     b      c      1.0   0.1 -18.1046 20.1046  False
----------------------------------------------------

statsmodels错误消息：“ ValueError：当p> = .9时v必须大于1”

1 个答案: