Python ValueError: alpha <= 0 Dirichlet 分布

Question

我有一个选举的投票结果。我想为每次投票从 Dirichlet 分布中模拟 50000 次抽签。但是当我运行这个单元格时。我收到此错误，但无法弄清楚我错在哪里。任何帮助，将不胜感激。

我的代码如下：

number= 50000

for i in range(len(dAll)):
    vector =list(dAll[i][3:7])

    random=np.random.dirichlet(vector,number)

    postn=[]
    for i in range(number): postn.append(random[i].argmax(axis=0))
    postn=np.array(postn)
    
    PNPD=100*len(np.where(postn ==0)[0])/number
    ProbPNPD.append(PNPD)
    
    PL=100*len(np.where(postn ==1)[0])/number
    ProbPL.append(PL)
    
    AD=100*len(np.where(postn ==2)[0])/number
    ProbAD.append(AD)
    
    Others=100*len(np.where(postn ==3)[0])/number
    ProbOthers.append(Others)

#Probability of victory for each party
    setprint=False
    if setprint == True:
        print('PNPD',PNPD,'%')
        print('PL',PL,'%')
        print('AD',AD,'%')
        print('Others',Others,'%')

错误代码：

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-7-16ba25b5e62c> in <module>
      4     vector =list(dAll[i][3:7])
      5 
----> 6     random=np.random.dirichlet(vector,number)
      7 
      8     postn=[]

mtrand.pyx in numpy.random.mtrand.RandomState.dirichlet()

ValueError: alpha <= 0

我的数据集中没有任何负数。 check in here

Answer 1

好吧，当我改变这部分时：

for i in range(len(dAll)):
    vector =list(dAll[i][3:7]) --> vector =list(dAll[i][3:**6**])

我解决这个问题。但是我仍然必须添加此列进行分析。