BOUNTY GOAL
To get bounty I am asking to build a function that serves as an analog of FindDistribution
. ItYou can also simply re-implement FindDistribution
if you know how it works. Your solution can work in same or different ways, have same or different syntax and settingsoptions, but it needs to give a similar output: best fitted distributions to a data, desirably ranked by fitness. Minimal working function demoed on ay data example is OKacceptable, even in a quite raw state.
ORIGINAL POST
Friends, especially those with math and stats background, perhaps you could enlighten. We face the following general problem quite often. Perhaps as a community we should find solution working especially well inside Wolfram ecosystem. Moreover - I think some of you might have such strategies - would be nice to share with community.
FindDistribution
is an excellent automation for search of analytic statistical models fitting data. But it is sophisticated large machine learning basedimplementation and we cannot control the algorithm there unless we explicitly rewrite or reinvent it. I wonder if we can build a simple but exhaustive algorithm taking in account all available in Wolfram Language statistical distributions (see related) that acts similar to FindDistribution
.
For the data consider any that is listed in the APPLICATIONS section of docs on FindDistribution. Or any data you like.
I suggest the following very simple idea for a start:
Take all set of analytical distributions in Wolfram Language (suitable to your data - e.g. all continuous or all discrete) -- also related
Use
EstimatedDistribution
orFindDistributionParameters
to fit all of themRun some tests - like
PearsonChiSquareTest
orDistributionFitTest
or anything similar and find top models according to some criteriaIf even top models are still far from a good fit use things like
MixtureDistribution
andTransformedDistribution
etc. to derive better models. See Derived Statistical DistributionsWe could also think of aadditional criteria that chooses models of smaller complexity - in terms of fitted parameters number e e.g. number of fitted parameters.
This is very blueprint-ish :-) I lack the deeper vision that takes in combination both: stats knowledge and Wolfram ecosystem structure. I hope some of you got the insight.
Feel free to demonstrate any strategy on any simple data I described above. Thsi might give some scope:
http://reference.wolfram.com/language/guide/ProbabilityAndStatistics.html