python - Getting a pdf from scipy.stats in a generic way -
i running goodness of fit tests using scipy.stats in python 2.7.10.
for distrname in distrnamelist:     distr = getattr(distributions, distrname)     param = distr.fit(sample)     pdf   = distr.pdf(???)   what pass distr.pdf() values of best-fit pdf on list of sample points of interest, called abscissas?
to evaluate pdf @ abscissas, pass abcissas first argument pdf. specify parameters, use * operator unpack param tuple , pass values distr.pdf:
pdf = distr.pdf(abscissas, *param)   for example,
import numpy np import scipy.stats stats  distrnamelist = ['beta', 'expon', 'gamma'] sample = stats.norm(0, 1).rvs(1000) abscissas = np.linspace(0,1, 10) distrname in distrnamelist:     distr = getattr(stats.distributions, distrname)     param = distr.fit(sample)     pdf = distr.pdf(abscissas, *param)     print(pdf)      
Comments
Post a Comment