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