Three and four parameter distributions of events and intervals based on the Kummer and Tricomi functions are shown to be statistically tractable. The Kummer distributions cover events less as well as more clustered than the Poisson while the Tricomi distributions have two and three parameter negative binomials as special cases. This flexibility and the moment relations between the parameters make them attractive in cases where the Poisson distribution proves inadequate. Fitzgerald (2000) has shown that statistically interesting results involving Tricomi functions arise when the rate parameter of a Poisson process is randomised by an inverted beta or a Tricomi exponential variate. These results are now developed into a flexible method for treating series of events more clustered than the Poisson and the intervals between them. The method offered there for computing Tricomi functions by statistical means is improved and a method of generating Tricomi exponential variates is also given. Then, using truncated beta distributions as randomisers, waiting-time distributions based on Kummer functions are derived; a statistical method of computing the Kummer function emerges and leads to another accurate method of computing the Tricomi function.
IntroductionMoment relations between the parameters of the distributions are readily obtained as are the probabilities of a given number of occurrences and the statistical properties of the inter-occurrence times.A hyper-Poisson distribution based on Kummer functions and having four parameters is also treated as it covers series of events whether more or less clustered than the Poisson. In a very similar vein, but having a negative binomial as a special case, a four-parameter process based on the Tricomi function can be found. Mathematically and statistically both are generally as tractable as the three-parameter randomised versions but moment solutions in the four-parameter case present problems.With all these processes it is shown that, when the distribution for the exceedances possesses an inverse, the cumulative distribution functions of the return periods have the same analytic form as that of the occurrences and the values can be estimated by a method such as bisection.