Я пытаюсь вычислить бета-pert с помощью scipy.stats.beta.cdf, но всякий раз, когда запускается приведенный ниже код, он возвращает [нан нан нан нан нан нан нан нан нан]. Мне нужна помощь в расчете бета-pert.
import scipy
from scipy import stats
from scipy.stats import beta
import numpy as np
import random
random_values = np.random.randint(low=5000, high=625000, size=10)
print ('impart random' + ' ' + str(random_values) )
minimum = 5000
mode = 215000
maximum = 625000
d = (minimum + 4*mode + maximum)/6
print ("intermediate value" + ' ' + str(d) )
alpha = 6*((d - minimum)/(maximum - minimum))
beta = 6*((maximum - d)/(maximum - minimum))
location = minimum
scale = maximum - minimum
pert = scipy.stats.beta.cdf(random_values, alpha, beta, location, scale)
print ('pert' + ' ' + str(pert))
Спасибо за помощь.