I've been doing some data analysis with Scipy. So far I accomplished this with continuous distributions:
- I can fit a probability distribution to a set of data points using a maximum likelihood fit. For example using
stats.chi2.fit(data_points)
. - I can test if my samples comes from a specific distribution. For example using
stats.kstest(data_points, 'norm')
I tried to do the same analysis with discrete distributions (like Poisson, Binomial and Geometric), but it seems like Scipy doesn't offer similar functionality for discrete distributions. What's the standard approach for fitting and testing for discrete distributions?