Applied Mathematics
Vol.5 No.6(2014), Article ID:44594,6 pages DOI:10.4236/am.2014.56091

Statistically Dual Distributions and Estimation

Sergey Bityukov1,2, Nikolai Krasnikov2,3, Saralees Nadarajah4, Vera Smirnova1

1Institute for High Energy Physics, Protvino, Russia

2Institute for Nuclear Research RAS, Moscow, Russia

3Joint İnstitute for Nuclear Research, Dubna, Russia

4School of Mathematics, University of Manchester, Manchester, UK

Email: Vera.Smirnova @ihep.ru, Serguei.Bitioukov@cern.ch

Copyright © 2014 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 24 January 2014; revised 24 February 2014; accepted 4 March 2014

ABSTRACT

The reconstruction of a parameter by the measurement of a random variable depending on the parameter is one of the main tasks in statistics. In statistical inference, the concept of a confidence distribution and, correspondingly, confidence density has often been loosely referred to as a distribution function on the parameter space that can represent confidence intervals of all levels for a parameter of interest. In this short note, the notion of statistically dual distributions is discussed. Based on properties of statistically dual distributions, a method for reconstructing the confidence density of a parameter is proposed.

Keywords:Distribution Theory, Confidence Distribution, Measurement, Error Theory, Data Analysis: Algorithms and Implementation, Data Management

1. Introduction

Let denote the observed number of events in a simple Poisson process. Its distribution can be described by a gamma random variable, with the probability density function (pdf) that looks like a Poisson probability:

where is a variable and is a parameter (in the case of the Poisson distribution, is a parameter and is a variable). It means, as shown below, that we can estimate the value and error of Poisson distribution parameter by the mean of a Poisson random variable and by using the corresponding gamma distribution. This approach is also correct in the case of the normal distribution, Cauchy distribution, Laplace distribution, the inverse gamma distribution.

Let us name such distributions, which allow one to exchange the parameter and the variable, conserving the same formula for the distribution of probabilities, as statistically dual distributions [1] . In many cases statistical duality of such type can be used for construction of confidence intervals for parameters.

In the next section, we show that Poisson and gamma distributions are statistically dual distributions and that the normal and Cauchy distributions are statistically self-dual distributions. An application of statistical duality for estimation of parameters is discussed in Section 3.

2. Statistically Dual Distributions

Definition: If a function can be expressed as a family of pdfs for variable given parameter, and a family of pdfs for variable given parameter, so that , then and are said to be statistically dual.

This definition is a purely probabilistic (and, in this sense, a frequentist) definition. Nevertheless, statistically dual distributions considered also belong to conjugate families defined in the Bayesian framework (see, for example, [2] ).

The statistical duality of Poisson and gamma distributions follows by a simple example. Let us consider the gamma distribution with pdf

Replacing, and by a, and, respectively, we obtain the following formula for the pdf

where a is a scale parameter and n + 1> 0 is a shape parameter. If a = 1 then the pdf of is

(1)

The Poisson distribution is a popular model for counts. For instance, if there are n events of a certain kind then it is reasonable to say

(2)

One can see that the parameter and the variable in Equations (1) and (2) are exchanged. In other aspects the formulae are identical. As a result these distributions (gamma and Poisson) are statistically dual distributions. These distributions are connected by the identity [3] (this identity arises in other forms in [4] -[6] )

(3)

for any and integer.

Another example of statistically dual distributions is the normal distribution with mean a and variance:

where x is a real variable, and are parameters. Here, we can exchange the parameter and the variable x without changing the formula for the pdf. It allows one to estimate the parameter by the mean value of x. In this case, the new pdf with variable and parameters and is

.

Hence, the normal distribution can be named as a statistically self-dual distribution. The identity analogous to (3) is

or, simply

(4)

for any real b, c and d.

The Cauchy distribution also has statistical self-duality like the normal distribution. The pdf of the Cauchy distribution is

where is a real variable, and is a real parameter. Here, we can exchange the parameter and the variable without altering the pdf. In this case,

Hence, the Cauchy distribution is also a statistically self-dual distribution. The identity analogous to (4) is

for any real, c and d.

The same property applies to several other distributions, for example, the Laplace distribution.

3. Estimation of Parameters

The identity (3) can be written in form [7]

(5)

that is,

for any, and and non negative integer.

The definition of the confidence interval for a Poisson parameter, , is [7]

(6)

where

This definition is consistent with the identity (5). It contrasts with other frequentist definitions of confidence intervals. The right hand side of (6) represents the frequentist definition.

Let us suppose that is the pdf of the Poisson parameter1 if number of observed events is equal to. It is a conditional pdf. It follows from the formulae (1), (5) that is a gamma pdf.

On the other hand: if is not equal to this pdf and the pdf of the Poisson parameter is some other function then

(7)

This identity is correct for any, and. The sum on the left hand side determines the boundary conditions on the confidence interval.

If we subtract Equation (7) from Equation (5) then we have

We can choose and arbitrarily. Let us choose them so that is not equal in the interval. For example, and. In this case, we have

and hence a contradiction (i.e.) everywhere except for a finite set of points). As a result we can mix Bayesian and frequentist probabilities without logical inconsistencies. The identity (5) leaves no place for any prior except for the uniform prior const2. Actually, it shows that the pdf is the “confidence density” [9] 3.

Statistically dual distributions allow one to exchange the parameter and the random variable. It means that one can construct the confidence interval for the parameter of the gamma distribution (Equation (1)) because

For the normal distribution the identity (4) can be written as

(8)

for any and.

This identity (8) also shows that the conditional distribution (if observed value is) of the true parameter obeys the normal distribution with mean and variance (here, in contrast to the previous example, is an unbiased estimator of the parameter). As a result we can construct the distribution of error and confidence intervals for 4, taking into account systematics and statistical uncertainties in accordance with standard analysis of errors [14] .

In case of the Cauchy distribution

for any and. This implies that the “confidence density” of the Cauchy distribution parameters is itself a Cauchy pdf.

4. Conclusions

We have discussed the notion of statistically dual distributions. The relation between the measurement of a casual variable and estimation of the given distribution parameter is discussed for three pairs of statistically dual distributions.

The proposed approach allows one to construct the distribution of the estimator of the distribution parameter by using statistically dual distributions. For example, the confidence density of the Poisson distribution parameter can be built by Monte Carlo by using properties of statistically dual distributions [15] . This notion is used for the construction of unified approach to measurement error and missing data [16] . The numerical example of the use of statistical duality for the confidence intervals construction is considered in the paper [17] .

In summary, statistical duality gives a clear frequentist interpretation of the “confidence density” of a parameter. It allows one to construct confidence intervals easily.

Acknowledgements

The authors are grateful to V. A. Kachanov, Louis Lyons, V. A. Matveev and V. F. Obraztsov for useful comments, R. D. Cousins, Yu. P. Gouz, G. Kahrimanis, V. Taperechkina and C. Wulz for fruitful discussions. This work has been particularly supported by the grant RFBR 13-02-00363.

References

  1. Bityukov, S.I., Krasnikov, N.V., Smirnova, V.V. and Taperechkina, V.A. (2006) Statistically Dual Distributions in Statisticalinference. In: Lyons, L. and Unel, M.K., Eds., Proceedings of Statistical Problems in Particle Physics, Astrophysics and Cosmology, Imperial College Press, Oxford, England, 102. http://dx.doi.org/10.1142/9781860948985_0023
  2. Bernardo, J.M. and Smith, A.F.M. (1994) Bayesian Theory. John Wiley and Sons, Chichester.
  3. Bityukov, S.I. (2002) On the Signal Significance in the Presence of Systematic and Statistical Uncertainties. Journal of High Energy Physics, 09, 060. http://dx.doi.org/10.1088/1126-6708/2002/09/060
  4. Jaynes, E.T. (1983) Confidenceintervals vs Bayesian Intervals. In: Rosenkrantz, R.D., Ed., E. T. Jaynes: Papers on Probability, Statistics and Statistical Physics, D. Reidel Publishing Company, Dordrecht, 165.
  5. Frodesen, A.G., Skjeggestad, O. and Toft, H. (1979) Probability and Statistics in Particle Physics. Universitetsforlaget, Bergen-Oslo-Tromso, 97.
  6. Cousins, R.D (1995) Why Isn’t Every Physicist a Bayesian? American Journal of Physics, 63, 398.
  7. Bityukov, S.I., Krasnikov, N.V. and Taperechkina, V.A. (2000) Confidence Intervals for Poisson Distribution Parameter. Preprint IFVE 2000-61, Protvino.
  8. Casella, G. and Berger, R.L. (2001) Statistical Inference. 2nd Edition, Duxbury Press.
  9. Efron, B. (1998) R. A. Fisher in the 21st Century. Statistical Science, 13, 95.
  10. Xie, M. and Singh, K. (2013) Confidence Distribution, the Frequentist Distribution Estimator of a Parameter: A Review. International Statistical Review, 81, 3-39.
  11. Schweder, T. and Hjort, N.L. (2003) Frequentist Analogies of Priors and Posteriors, Econometrics and the Philosophy of Economics. Princeton University Press, Princeton, 285.
  12. Bityukov, S., Krasnikov, N., Nadarajah, S. and Smirnova, V. (2010) Confidence Distributions in Statistical Inference. AIP Conference Proceedings, 130, 346. http://dx.doi.org/10.1063/1.3573637
  13. Fisher, R.A. (1930) Inverse Probability. Proceedings of the Cambridge Philosophical Society, 26, 528.
  14. Eadie, W.T., Drijard, D., James, F.E., Roos, M. and Sadoulet, B. (1971) Statistical Methods in Experimental Physics. North Holland, Amsterdam.
  15. Bityukov, S.I., Medvedev, V.A., Smirnova, V.V. and Zernii, Yu.V. (2004) Experimental Test of the Probability Densityfunction of True Value of Poisson Distribution Parameter by Single Observation of Number of Events. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 534, 228-231. http://dx.doi.org/10.1016/j.nima.2004.07.092
  16. Blackwell, M., Honaker, J. and King, G. (2012) Multiple Overimputation: A Unified Approach to Measurement Error and Missing Data. Political Methodology, Committee on Concepts and Methods, Working Paper Series, Paper 36.

NOTES

1The definition of conjugate families as stated in is: given a family F of pdf's indexed by a parameter, a family, , of prior distributions is said to be conjugate for the family F if the posterior distribution of is in the family for all, all priors and all possible data x.

2Bayesian methods suppose that, where is the prior pdf for.

3The confidence density is a natural notion in the frequentist concept of confidence distributions -. Confidence distributions are a way to represent all possible confidence intervals. The area under the confidence distribution between any two points gives the confidence that the parameter value will lie between those points. In Bayesian notation, the confidence density can be considered as a posteriori pdf with the assumption that we have a uniform prior.

4In this case, the definition of confidence intervals as

where is the observed value of 𝑥, and, coincides with the definition of fiducial intervals introduced by R. A. Fisher .