Confidence distribution

In

In recent years, there has been a surge of renewed interest in confidence distributions.

), but it uses a sample-dependent distribution function on the parameter space (instead of a point or an interval) to estimate the parameter of interest.

A simple example of a confidence distribution, that has been broadly used in statistical practice, is a

Just as a Bayesian posterior distribution contains a wealth of information for any type of

History

Neyman (1937)^[8] introduced the idea of "confidence" in his seminal paper on confidence intervals which clarified the frequentist repetition property. According to Fraser,^[9] the seed (idea) of confidence distribution can even be traced back to Bayes (1763)^[10] and Fisher (1930).^[1] Although the phrase seems to first be used in Cox (1958).^[11] Some researchers view the confidence distribution as "the Neymanian interpretation of Fisher's fiducial distributions",^[12] which was "furiously disputed by Fisher".^[13] It is also believed that these "unproductive disputes" and Fisher's "stubborn insistence"^[13] might be the reason that the concept of confidence distribution has been long misconstrued as a fiducial concept and not been fully developed under the frequentist framework.^[6][14] Indeed, the confidence distribution is a purely frequentist concept with a purely frequentist interpretation, although it also has ties to Bayesian and fiducial inference concepts.

Definition

Classical definition

Classically, a confidence distribution is defined by inverting the upper limits of a series of lower-sided confidence intervals.^{[15][16][page needed]} In particular,

For every α in (0, 1), let (−∞, ξ_n(α)] be a 100α% lower-side confidence interval for θ, where ξ_n(α) = ξ_n(X_n,α) is continuous and increasing in α for each sample X_n. Then, H_n(•) = ξ_n⁻¹(•) is a confidence distribution for θ.

Efron stated that this distribution "assigns probability 0.05 to θ lying between the upper endpoints of the 0.90 and 0.95 confidence interval, etc." and "it has powerful intuitive appeal".^[16] In the classical literature,^[3] the confidence distribution function is interpreted as a distribution function of the parameter θ, which is impossible unless fiducial reasoning is involved since, in a frequentist setting, the parameters are fixed and nonrandom.

To interpret the CD function entirely from a frequentist viewpoint and not interpret it as a distribution function of a (fixed/nonrandom) parameter is one of the major departures of recent development relative to the classical approach. The nice thing about treating confidence distributions as a purely frequentist concept (similar to a point estimator) is that it is now free from those restrictive, if not controversial, constraints set forth by Fisher on fiducial distributions.^[6][14]

The modern definition

The following definition applies;^[12][17][18] Θ is the parameter space of the unknown parameter of interest θ, and χ is the sample space corresponding to data X_n={X₁, ..., X_n}:

A function H_n(•) = H_n(X_n, •) on χ × Θ → [0, 1] is called a confidence distribution (CD) for a parameter θ, if it follows two requirements:

• (R1) For each given X_n ∈ χ, H_n(•) = H_n(X_n, •) is a continuous cumulative distribution function on Θ;
• (R2) At the true parameter value θ = θ₀, H_n(θ₀) ≡ H_n(X_n, θ₀), as a function of the sample X_n, follows the uniform distribution U[0, 1].

Also, the function H is an asymptotic CD (aCD), if the U[0, 1] requirement is true only asymptotically and the continuity requirement on H_n(•) is dropped.

In nontechnical terms, a confidence distribution is a function of both the parameter and the random sample, with two requirements. The first requirement (R1) simply requires that a CD should be a distribution on the parameter space. The second requirement (R2) sets a restriction on the function so that inferences (point estimators, confidence intervals and hypothesis testing, etc.) based on the confidence distribution have desired frequentist properties. This is similar to the restrictions in point estimation to ensure certain desired properties, such as unbiasedness, consistency, efficiency, etc.^[6][19]

A confidence distribution derived by inverting the upper limits of confidence intervals (classical definition) also satisfies the requirements in the above definition and this version of the definition is consistent with the classical definition.^[18]

Unlike the classical fiducial inference, more than one confidence distributions may be available to estimate a parameter under any specific setting. Also, unlike the classical fiducial inference, optimality is not a part of requirement. Depending on the setting and the criterion used, sometimes there is a unique "best" (in terms of optimality) confidence distribution. But sometimes there is no optimal confidence distribution available or, in some extreme cases, we may not even be able to find a meaningful confidence distribution. This is not different from the practice of point estimation.

A definition with measurable spaces

A confidence distribution^[20] ${\displaystyle C}$ for a parameter ${\displaystyle \gamma }$ in a measurable space is a distribution estimator with ${\displaystyle C(A_{p})=p}$ for a family of confidence regions ${\displaystyle A_{p}}$ for ${\displaystyle \gamma }$ with level ${\displaystyle p}$ for all levels ${\displaystyle 0<p<1}$. The family of confidence regions is not unique.^[21] If ${\displaystyle A_{p}}$ only exists for ${\displaystyle p\in I\subset (0,1)}$, then ${\displaystyle C}$ is a confidence distribution with level set ${\displaystyle I}$. Both ${\displaystyle C}$ and all ${\displaystyle A_{p}}$ are measurable functions of the data. This implies that ${\displaystyle C}$ is a random measure and ${\displaystyle A_{p}}$ is a random set. If the defining requirement ${\displaystyle P(\gamma \in A_{p})\geq p}$ holds with equality, then the confidence distribution is by definition exact. If, additionally, ${\displaystyle \gamma }$ is a real parameter, then the measure theoretic definition coincides with the above classical definition.

Examples

Example 1: Normal mean and variance

Suppose a normal sample X_i ~ N(μ, σ²), i = 1, 2, ..., n is given.

(1) Variance σ² is known

Let, Φ be the cumulative distribution function of the standard normal distribution, and ${\displaystyle F_{t_{n-1}}}$ the cumulative distribution function of the Student ${\displaystyle t_{n-1}}$ distribution. Both the functions ${\displaystyle H_{\mathit {\Phi }}(\mu )}$ and ${\displaystyle H_{t}(\mu )}$ given by

${\displaystyle H_{\Phi }(\mu )={\mathit {\Phi }}\left({\frac {{\sqrt {n}}(\mu -{\bar {X}})}{\sigma }}\right),\quad {\text{and}}\quad H_{t}(\mu )=F_{t_{n-1}}\left({\frac {{\sqrt {n}}(\mu -{\bar {X}})}{s}}\right),}$

satisfy the two requirements in the CD definition, and they are confidence distribution functions for μ.^[3] Furthermore,

${\displaystyle H_{A}(\mu )={\mathit {\Phi }}\left({\frac {{\sqrt {n}}(\mu -{\bar {X}})}{s}}\right)}$

satisfies the definition of an asymptotic confidence distribution when n→∞, and it is an asymptotic confidence distribution for μ. The uses of ${\displaystyle H_{t}(\mu )}$ and ${\displaystyle H_{A}(\mu )}$ are equivalent to state that we use ${\displaystyle N({\bar {X}},\sigma ^{2})}$ and ${\displaystyle N({\bar {X}},s^{2})}$ to estimate ${\displaystyle \mu }$, respectively.

(2) Variance σ² is unknown

For the parameter μ, since ${\displaystyle H_{\mathit {\Phi }}(\mu )}$ involves the unknown parameter σ and it violates the two requirements in the CD definition, it is no longer a "distribution estimator" or a confidence distribution for μ.^[3] However, ${\displaystyle H_{t}(\mu )}$ is still a CD for μ and ${\displaystyle H_{A}(\mu )}$ is an aCD for μ.

For the parameter σ², the sample-dependent cumulative distribution function

${\displaystyle H_{\chi ^{2}}(\theta )=1-F_{\chi _{n-1}^{2}}((n-1)s^{2}/\theta )}$

is a confidence distribution function for σ².^[6] Here, ${\displaystyle F_{\chi _{n-1}^{2}}}$ is the cumulative distribution function of the ${\displaystyle \chi _{n-1}^{2}}$ distribution .

In the case when the variance σ² is known, ${\displaystyle H_{\mathit {\Phi }}(\mu )={\mathit {\Phi }}\left({\frac {{\sqrt {n}}(\mu -{\bar {X}})}{\sigma }}\right)}$ is optimal in terms of producing the shortest confidence intervals at any given level. In the case when the variance σ² is unknown, ${\displaystyle H_{t}(\mu )=F_{t_{n-1}}\left({\frac {{\sqrt {n}}(\mu -{\bar {X}})}{s}}\right)}$ is an optimal confidence distribution for μ.

Example 2: Bivariate normal correlation

Let ρ denotes the

:

${\displaystyle z={1 \over 2}\ln {1+r \over 1-r}}$

has the limiting distribution ${\displaystyle N({1 \over 2}\ln {{1+\rho } \over {1-\rho }},{1 \over n-3})}$ with a fast rate of convergence, where r is the sample correlation and n is the sample size.

The function

${\displaystyle H_{n}(\rho )=1-{\mathit {\Phi }}\left({\sqrt {n-3}}\left({1 \over 2}\ln {1+r \over 1-r}-{1 \over 2}\ln {{1+\rho } \over {1-\rho }}\right)\right)}$

is an asymptotic confidence distribution for ρ.^[22]

An exact confidence density for ρ is^[23][24]

${\displaystyle \pi (\rho |r)={\frac {\nu (\nu -1)\Gamma (\nu -1)}{{\sqrt {2\pi }}\Gamma (\nu +{\frac {1}{2}})}}(1-r^{2})^{\frac {\nu -1}{2}}\cdot (1-\rho ^{2})^{\frac {\nu -2}{2}}\cdot (1-r\rho )^{\frac {1-2\nu }{2}}F({\frac {3}{2}},-{\frac {1}{2}};\nu +{\frac {1}{2}};{\frac {1+r\rho }{2}})}$

where ${\displaystyle F}$ is the Gaussian hypergeometric function and ${\displaystyle \nu =n-1>1}$ . This is also the posterior density of a Bayes matching prior for the five parameters in the binormal distribution.^[25]

The very last formula in the classical book by Fisher gives

${\displaystyle \pi (\rho |r)={\frac {(1-r^{2})^{\frac {\nu -1}{2}}\cdot (1-\rho ^{2})^{\frac {\nu -2}{2}}}{\pi (\nu -2)!}}\partial _{\rho r}^{\nu -2}\left\{{\frac {\theta -{\frac {1}{2}}\sin 2\theta }{\sin ^{3}\theta }}\right\}}$

where ${\displaystyle \cos \theta =-\rho r}$ and ${\displaystyle 0<\theta <\pi }$. This formula was derived by C. R. Rao.^[26]

Example 3: Binormal mean

Let data be generated by ${\displaystyle Y=\gamma +U}$ where ${\displaystyle \gamma }$ is an unknown vector in the

and ${\displaystyle U}$ has a binormal and known distribution in the plane. The distribution of ${\displaystyle \Gamma ^{y}=y-U}$ defines a confidence distribution for ${\displaystyle \gamma }$. The confidence regions ${\displaystyle A_{p}}$ can be chosen as the interior of ellipses centered at ${\displaystyle \gamma }$ and axes given by the eigenvectors of the covariance matrix of ${\displaystyle \Gamma ^{y}}$. The confidence distribution is in this case binormal with mean ${\displaystyle \gamma }$, and the confidence regions can be chosen in many other ways.^[21] The confidence distribution coincides in this case with the Bayesian posterior using the right Haar prior.^[27] The argument generalizes to the case of an unknown mean ${\displaystyle \gamma }$ in an infinite-dimensional Hilbert space, but in this case the confidence distribution is not a Bayesian posterior.^[28]

Using confidence distributions for inference

Confidence interval

From the CD definition, it is evident that the interval ${\displaystyle (-\infty ,H_{n}^{-1}(1-\alpha )],[H_{n}^{-1}(\alpha ),\infty )}$ and ${\displaystyle [H_{n}^{-1}(\alpha /2),H_{n}^{-1}(1-\alpha /2)]}$ provide 100(1 − α)%-level confidence intervals of different kinds, for θ, for any α ∈ (0, 1). Also ${\displaystyle [H_{n}^{-1}(\alpha _{1}),H_{n}^{-1}(1-\alpha _{2})]}$ is a level 100(1 − α₁ − α₂)% confidence interval for the parameter θ for any α₁ > 0, α₂ > 0 and α₁ + α₂ < 1. Here, ${\displaystyle H_{n}^{-1}(\beta )}$ is the 100β% quantile of ${\displaystyle H_{n}(\theta )}$ or it solves for θ in equation ${\displaystyle H_{n}(\theta )=\beta }$. The same holds for a CD, where the confidence level is achieved in limit. Some authors have proposed using them for graphically viewing what parameter values are consistent with the data, instead of coverage or performance purposes.^[29][30]

Point estimation

Point estimators can also be constructed given a confidence distribution estimator for the parameter of interest. For example, given H_n(θ) the CD for a parameter θ, natural choices of point estimators include the median M_n = H_n⁻¹(1/2), the mean ${\displaystyle {\bar {\theta }}_{n}=\int _{-\infty }^{\infty }t\,\mathrm {d} H_{n}(t)}$, and the maximum point of the CD density

${\displaystyle {\widehat {\theta }}_{n}=\arg \max _{\theta }h_{n}(\theta ),h_{n}(\theta )=H'_{n}(\theta ).}$

Under some modest conditions, among other properties, one can prove that these point estimators are all consistent.^[6][22] Certain confidence distributions can give optimal frequentist estimators.^[28]

Hypothesis testing

One can derive a p-value for a test, either one-sided or two-sided, concerning the parameter θ, from its confidence distribution H_n(θ).^[6][22] Denote by the probability mass of a set C under the confidence distribution function ${\displaystyle p_{s}(C)=H_{n}(C)=\int _{C}\mathrm {d} H(\theta ).}$ This p_s(C) is called "support" in the CD inference and also known as "belief" in the fiducial literature.^[31] We have

(1) For the one-sided test K₀: θ ∈ C vs. K₁: θ ∈ C^c, where C is of the type of (−∞, b] or [b, ∞), one can show from the CD definition that sup_θ ∈ CP_θ(p_s(C) ≤ α) = α. Thus, p_s(C) = H_n(C) is the corresponding p-value of the test.

(2) For the singleton test K₀: θ = b vs. K₁: θ ≠ b, P_{{K0: θ = b}}(2 min{p_s(C_lo), one can show from the CD definition that p_s(C_up)} ≤ α) = α. Thus, 2 min{p_s(C_lo), p_s(C_up)} = 2 min{H_n(b), 1 − H_n(b)} is the corresponding p-value of the test. Here, C_lo = (−∞, b] and C_up = [b, ∞).

See Figure 1 from Xie and Singh (2011)^[6] for a graphical illustration of the CD inference.

Implementations

A few statistical programs have implemented the ability to construct and graph confidence distributions.

R, via the concurve,^[32][33] pvaluefunctions,^[34] and episheet^[35] packages

Excel, via episheet^[36]

Stata, via concurve^[32]

See also

• Coverage probability

References

Bibliography

• Xie, M. and Singh, K. (2013). [1] "Confidence Distribution, the Frequentist Distribution Estimator of a Parameter: A Review". International Statistical Review, 81, 3–39.
• Schweder, T and Hjort, N L (2016). [2]Confidence, Likelihood, Probability: Statistical Inference with Confidence Distributions. London: Cambridge University Press.
  ISBN 9781139046671
• Fisher, R A (1956). Statistical Methods and Scientific Inference. New York: Hafner.
  ISBN 0-02-844740-9
  .
• Fisher, R. A. (1955). "Statistical methods and scientific induction" J. Roy. Statist. Soc. Ser. B. 17, 69—78. (criticism of statistical theories of Jerzy Neyman and Abraham Wald from a fiducial perspective)
• Hannig, J. (2009). "On generalized fiducial inference". Statistica Sinica, 19, 491–544.
• Lawless, F. and Fredette, M. (2005). "Frequentist prediction intervals and predictive distributions." Biometrika. 92(3) 529–542.
• Lehmann, E.L. (1993). "The Fisher, Neyman–Pearson theories of testing hypotheses: one theory or two?" Journal of the American Statistical Association 88 1242–1249.
• Neyman, Jerzy (1956). "Note on an Article by Sir Ronald Fisher". Journal of the Royal Statistical Society. Series B (Methodological) 18 (2): 288–294.
  JSTOR 2983716
  . (reply to Fisher 1955, which diagnoses a fallacy of "fiducial inference")
• Schweder T., Sadykova D., Rugh D. and Koski W. (2010) "Population Estimates From Aerial Photographic Surveys of Naturally and Variably Marked Bowhead Whales" Journal of Agricultural Biological and Environmental Statistics 2010 15: 1–19
• Bityukov S., Krasnikov N., Nadarajah S. and Smirnova V. (2010) "Confidence distributions in statistical inference". AIP Conference Proceedings, 1305, 346-353.
• Singh, K. and Xie, M. (2012). "CD-posterior --- combining prior and data through confidence distributions." Contemporary Developments in Bayesian Analysis and Statistical Decision Theory: A Festschrift for William E. Strawderman. (D. Fourdrinier, et al., Eds.). IMS Collection, Volume 8, 200 -214.