alphaFeedback

WorksMathematics

An Asymptotic Framework for Closed-Form Estimators from Surrogate Estimating Equations

AI contributions
  • ChatGPT
View details
Ancillary data
1 code link
View

Abstract

Many statistical estimators, including maximum likelihood and method of moments estimators, are defined as roots of estimating equations. In complex models, these equations often cannot be solved explicitly; although numerical methods can provide approximate solutions, closed-form approximations can reduce computation and facilitate theoretical analysis by making an estimator's dependence on the data more transparent. In this paper, we develop a general approach to constructing such approximations and a framework for quantifying their accuracy. Our theory accommodates sequences of parameters that approach the boundary of the parameter space or diverge to infinity, and provides conditions under which a closed-form approximate estimator has the same first-order asymptotic distribution as its exact counterpart. We focus on approximations to maximum likelihood estimators, although the theory applies to general estimating equations. We apply the framework to iid models from the gamma, zeta, von Mises, and generalized Pareto families, as well as to gamma regression. Together, these examples cover scalar and vector parameters, iid and non-iid settings, and normal and nonnormal asymptotic regimes. Each example is accompanied by a simulation study, and the gamma regression example also includes a real-world data application.

Introduction

Many statistical estimators are defined implicitly as roots of data-dependent estimating equations. Given a parameter space Θ⊆Rd\Theta\subseteq\mathbb{R}^{d} and an estimating equation qn:Θ→Rd\boldsymbol{q}_{n}: \Theta\to\mathbb{R}^{d} of the form

qn(θ∣Xn)=0,(1)\boldsymbol{q}_{n}(\theta\mid\boldsymbol{X}_{n}) = 0, \tag*{(1)}

a root — when it exists — is known as a Z-estimator. When qn\boldsymbol{q}_{n} is a score function for some statistical model, for example, any maximum likelihood estimator (MLE) that is an interior point of Θ\Theta is a root of qn\boldsymbol{q}_{n}. While score equations are perhaps the most prominent sources of Z-estimators, method of moments and generalized method of moments estimators under suitable regularity conditions [e.g. 17] also fit into this framework. The study of such estimators dates back to Huber [ [23]; modern treatments are given by van der Vaart [ [38] and van der Vaart and Wellner [ [39] who develop asymptotic theory of these estimators.

In the simplest models, (1) can be solved explicitly, yielding a closed-form1 function of the data Xn\boldsymbol{X}_{n}. More commonly, however, the estimating equation contains a nonlinear — often transcendental — function with no explicit inverse, couples several parameters together, or has a nearly singular Jacobian near the boundary of Θ\Theta. Classical examples include the beta and Dirichlet distributions [6, 29], the Student’s tt distribution [25, 19], the von Mises distribution [34, 21], and families used in extreme value analysis [32, 16]. Although modern root-finding methods and software often make solving such equations routine, an explicit approximation avoids the need for an initial value, a stopping criterion, and repeated computations; it can also be useful in regimes where numerical root-finding is unstable, and can facilitate theoretical analysis by making the estimator’s dependence on the data more transparent. Work on specific models has therefore pursued closed-form or rapidly computable alternatives when exact likelihood estimation requires iteration; examples include Ye and Chen [2017] for the gamma distribution, Clauset et al. [ [9] for the zeta distribution, and Ho Chang et al. [ [20] for the Dirichlet distribution.

In this paper, we present an approach to construct closed-form approximations to roots of intractable estimating equations of the form (1) by replacing the original estimating equation with a tractable surrogate that admits a closed-form solution, and we develop a general framework for determining when such approximations are accurate and for quantifying the resulting error. The framework considers sequences of data-generating parameters that approach a limiting value; importantly, this limiting value may lie on the boundary of the parameter space or be a formal value such as infinity.

The gamma model provides a simple illustration. Suppose that X1,X2,…,XnX_{1}, X_{2}, \ldots, X_{n} constitute a random sample from the Gamma(α,β)(\alpha, \beta) distribution in the shape-rate parameterization, where the rate parameter β\beta is known. The MLE α^n\hat{\alpha}_n maximizes the log-likelihood function, and differentiating the log-likelihood shows that it is the unique solution of

ψ(α^n)=log⁡(βGn),(2)\psi(\hat{\alpha}_n) = \log(\beta G_n), \tag*{(2)}

where ψ:=(log⁡(Γ))′\psi:= (\log(\Gamma))' is the digamma function and Gn:=(∏i=1nXi)1/nG_n := \left(\prod_{i=1}^{n} X_i\right)^{1/n} is the geometric mean of the sample. Because ψ\psi has no standard inverse in closed form, α^n\hat{\alpha}_n is usually obtained numerically, using Newton’s method, for example; see Greenwood and Durand [ [15]; Choi and Wette [ [8]; Minka [ [30] for computational treatments. However, (2) is equivalent to

exp⁡(−ψ(α^n))=1βGn,(3)\exp(-\psi(\hat{\alpha}_n)) = \frac{1}{\beta G_n}, \tag*{(3)}

and, as α→∞\alpha\to\infty, one can show that

exp⁡(−ψ(α))=1α−1/2+O(1α3).\exp(-\psi(\alpha)) = \frac{1}{\alpha- 1/2} + O\left(\frac{1}{\alpha^3}\right).

Replacing the left side of (3) by its leading approximation and solving yields the explicit estimator

α~n:=12+β(∏i=1nXi)1/n.(4)\tilde{\alpha}_n := \frac{1}{2} + \beta\left(\prod_{i=1}^{n} X_i\right)^{1/n}. \tag*{(4)}

We show in § 3 that, for data-generating sequences {αn}n≥1\{\alpha_n\}_{n\geq1} satisfying αn→∞\alpha_n \to\infty,

∣α~n−α^n∣=oPαn(1),\left|\tilde{\alpha}_n-\hat{\alpha}_n\right|=o_{P_{\alpha_n}}(1),

so the two estimators are equivalent on the original parameter scale. If, in addition, αn3/n→∞\alpha_n^3/n \to\infty, then

nαn∣α~n−α^n∣=oPαn(1),\sqrt{\frac{n}{\alpha_n}}\left|\tilde{\alpha}_n-\hat{\alpha}_n\right|=o_{P_{\alpha_n}}(1),

and under the same condition, we further show that

nαn(α~n−αn)→dN(0,1)andnαn(α^n−αn)→dN(0,1).\sqrt{\frac{n}{\alpha_n}}(\tilde{\alpha}_n-\alpha_n)\xrightarrow{d}\mathcal{N}(0,1) \quad\text{and}\quad \sqrt{\frac{n}{\alpha_n}}(\hat{\alpha}_n-\alpha_n)\xrightarrow{d}\mathcal{N}(0,1).

Thus, the approximate estimator and the MLE have the same first-order asymptotic distribution.

The preceding calculation suggests a general strategy. First, rewrite the exact estimating equation qn(θ∣Xn)=0\boldsymbol{q}_n(\theta\mid X_n)=\boldsymbol{0} in a form that exposes the behavior of its intractable component near the limiting parameter value. Next, replace that component by a finite truncation of its asymptotic expansion, retaining sufficiently many expansion terms to control the approximation error while keeping the resulting surrogate estimating equation q~n(θ∣Xn)=0\widetilde{\boldsymbol{q}}_n(\theta\mid X_n)=\boldsymbol{0} explicitly solvable. The root θ~n\tilde{\theta}_n of the surrogate equation is then a candidate approximation to the selected exact root θ^n\hat{\theta}_n. Whether it inherits the first-order asymptotic distribution of θ^n\hat{\theta}_n depends not only on convergence of the parameter to its limiting value, but on its rate of approach relative to the sample size as well.

Although deriving θ~n\tilde{\theta}_n can be straightforward, as in the gamma example above, establishing that it is close to θ^n\hat{\theta}_n requires a careful analysis. The essential problems are that both the equations and their selected roots are random, the sample size and the data-generating parameter move jointly, and the accuracy of the approximation depends on the scale at which the estimators will be compared. Accordingly, even if ∥qn(θ∣Xn)−q~n(θ∣Xn)∥→p0\|\boldsymbol{q}_n(\theta\mid X_n)-\widetilde{\boldsymbol{q}}_n(\theta\mid X_n)\|\xrightarrow{p}0 for each fixed θ∈Θ\theta\in\Theta, it does not follow that ∥θ^n−θ~n∥→p0\|\hat{\theta}_n-\tilde{\theta}_n\|\xrightarrow{p}0. Pointwise control need not extend uniformly to a neighborhood containing both θ^n\hat{\theta}_n and θ~n\tilde{\theta}_n, and even a small perturbation can produce a large change, especially when the Jacobian of the estimating equation is nearly singular as in the general framework of Rice [ [33]; a deterministic antecedent appears in de Bruijn [ [12].

To address these difficulties, we introduce one-to-one transformations of the parameter, which we call vanishing charts, under which the transformed data-generating parameter converges to zero along the regime under consideration. Such regimes can differ substantially from standard fixed-parameter asymptotics, exhibiting altered convergence rates, non-Gaussian limits, or even failure of MLE existence [37, 36, 14]. After localizing the exact and surrogate roots on a common deterministic set in vanishing chart coordinates, our analysis separates three factors governing their discrepancy: the error incurred when qnq_n is replaced by q~n\tilde{q}_n, the sensitivity of the surrogate root to perturbations, and the magnification that occurs when errors in chart coordinates are translated back to the original parameter scale. Our main result combines uniform control of the equation error with a lower bound on the smallest singular value of the surrogate Jacobian to bound the distance between the roots in chart coordinates. Further results accommodate vector parameters with different componentwise rates, transfer these bounds to the original parameter scale, and identify when θ~n\tilde{\theta}_n inherits the first-order asymptotic distribution of θ^n\hat{\theta}_n.

We study the gamma family in detail: the known-rate case is scalar, joint estimation of the shape and rate introduces different componentwise rates, and gamma regression shows that the observations need not be identically distributed (a real-world data example compares the exact and approximate gamma regression fits for measurements of black cherry trees). Our remaining examples highlight other features of the construction, with mathematical details left to the appendix. The zeta distribution is a discrete example where the MLE can fail to be finite. For the von Mises model, we preserve the explicit location estimator while approximating a ratio of Bessel functions for the concentration parameter. The generalized Pareto model has nonregular, parameter-dependent support, requiring the likelihood to be profiled prior to truncation. For all of these models, we use R simulations to assess the finite-sample accuracy of the approximate estimators and to illustrate the distinct features of each model.

The remainder of the paper develops the general theory in § 2, presents the worked examples in § 3 and § 4, and concludes with a discussion in § 5; technical details are deferred to the appendices.

Theory and Practical Guidance

Setup

For each n≥1n \ge1, let Xn:=(X1,…,Xn)X_n := (X_1,\ldots,X_n) be the random vector of data and consider the sequence of parametric models {Pn,θ:θ∈Θ⊆Rd}n≥1\{\mathbb{P}_{n,\theta}:\theta\in\Theta\subseteq\mathbb{R}^d\}_{n\ge1}, where Pn,θ\mathbb{P}_{n,\theta} denotes the joint law of XnX_n. In the iid setting, Xn=(X1,…,Xn)X_n = (X_1,\ldots,X_n) and Pn,θ=Pθ⊗n\mathbb{P}_{n,\theta} = \mathbb{P}_{\theta}^{\otimes n}. When a joint density or mass function pn,θp_{n,\theta} exists, an interior MLE at which log⁡(pn,θ(Xn))\log(p_{n,\theta}(X_n)) is differentiable is a root of the joint score function ∇log⁡(pn,θ(Xn))\nabla\log(p_{n,\theta}(X_n)). For a differentiable map F:Rd→RmF:\mathbb{R}^d \to\mathbb{R}^m, let JF(x)JF(x) denote its Jacobian at xx, and for a square matrix AA, let σmin⁡(A)\sigma_{\min}(A) denote its smallest singular value. ∥⋅∥\|\cdot\| denotes the Euclidean norm on vectors and its induced operator norm on matrices. For z∈Rmz \in\mathbb{R}^m and r>0r > 0, we write Br(z):={x∈Rm:∥x−z∥≤r}B_r(z) := \{x \in\mathbb{R}^m:\|x-z\| \le r\} for the closed ball of radius rr centered at zz. When it is useful to display the governing law explicitly, we write

Zn→Pn,θnzZ_n \xrightarrow{\mathbb{P}_{n,\theta_n}} z

References

  1. [1]Abeyasekera, S. and Collett, D. (1982). On the estimation of the parameters of the von mises distribution. Communications in Statistics - Theory and Methods, 11(18):2083–2090.DOI
  2. [2]Abramowitz, M. and Stegun, I. A. (1965). Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables. Dover Publications, New York, NY.
  3. [3]Apostol, T. M. (1976). Introduction to Analytic Number Theory. Undergraduate Texts in Mathematics. Springer, New York, NY.
  4. [4]Banerjee, A., Dhillon, I. S., Ghosh, J., and Sra, S. (2005). Clustering on the Unit Hypersphere using von Mises-Fisher Distributions. Journal of Machine Learning Research, 6(46):1345–1382.DOI
  5. [5]Barndorff-Nielsen, O. and Cox, D. R. (1979). Edgeworth and Saddle-Point Approximations with Statistical Applications. Journal of the Royal Statistical Society. Series B (Methodological), 41(3):279–312.
  6. [6]Beckman, R. J. and Tietjen, G. L. (1978). Maximum likelihood estimation for the beta distribution. J. Stat. Comput. Simul.; (United Kingdom), 7.DOI
  7. [7]Bertsekas, D. (2009). Convex Optimization Theory. Athena Scientific.
  8. [8]Choi, S. C. and Wette, R. (1969). Maximum Likelihood Estimation of the Parameters of the Gamma Distribution and Their Bias. Technometrics, 11(4):683–690.
  9. [9]Clauset, A., Shalizi, C. R., and Newman, M. E. J. (2009). Power-Law Distributions in Empirical Data. SIAM Review, 51(4):661–703.
  10. [10]Comtet, L. (1974). Advanced Combinatorics: The Art of Finite and Infinite Expansions. Springer Netherlands, Dordrecht.DOI
  11. [11]Daniels, H. E. (1983). Saddlepoint approximations for estimating equations. Biometrika, 70(1):89–96.DOI
  12. [12]de Bruijn, N. G. (1981). Asymptotic Methods in Analysis. Dover Publications, New York, NY.
  13. [13]Field, C. (1982). Small Sample Asymptotic Expansions for Multivariate M-Estimates. The Annals of Statistics, 10(3):672–689.DOI
  14. [14]Geyer, C. J. (1994). On the Asymptotics of Constrained M-Estimation. The Annals of Statistics, 22(4):1993–2010.DOI
  15. [15]Greenwood, J. A. and Durand, D. (1960). Aids for Fitting the Gamma Distribution by Maximum Likelihood. Technometrics, 2(1):55–65.DOI
  16. [16]Grimshaw, S. D. (1993). Computing Maximum Likelihood Estimates for the Generalized Pareto Distribution. Technometrics, 35(2):185–191.DOI
  17. [17]Hansen, L. P. (1982). Large Sample Properties of Generalized Method of Moments Estimators. Econometrica, 50(4):1029–1054.DOI
  18. [18]Hardy, G. H. and Wright, E. . M. (2008). An Introduction To The Theory Of Numbers. Oxford University Press.
  19. [19]Hasannasab, M., Hertrich, J., Laus, F., and Steidl, G. (2021). Alternatives to the EM algorithm for ML estimation of location, scatter matrix, and degree of freedom of the Student t distribution. Numerical Algorithms, 87(1):77–118.
  20. [20]Ho Chang, J., Kyu Lee, S., and Kim, H.-M. (2023). An asymptotically efficient closed-form estimator for the Dirichlet distribution. Stat, 12(1):e640.DOI
  21. [21]Hornik, K. and Grün, B. (2014). On maximum likelihood estimation of the concentration parameter of von Mises–Fisher distributions. Computational Statistics, 29(5):945–957.
  22. [22]Hosking, J. R. M. and Wallis, J. R. (1987). Parameter and Quantile Estimation for the Generalized Pareto Distribution. Technometrics, 29(3):339–349.
  23. [23]Huber, P. J. (1967). The behavior of maximum likelihood estimates under nonstandard conditions. In Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability, Volume 1: Statistics, volume 5.1, pages 221–234. University of California Press.
  24. [24]Kotz, S., Balakrishnan, N., and Johnson, N. L. (2000). Continuous Multivariate Distributions, Volume 1: Models and Applications. Wiley-Interscience, New York.
  25. [25]Liu, C. and Rubin, D. B. (1995). ML Estimation of the t Distribution Using EM and Its Extensions, ECM and ECME. Statistica Sinica, 5(1):19–39.
  26. [26]Mardia, K. V. and Jupp, P. E. (1999). Directional Statistics. Wiley Series in Probability and Statistics. John Wiley & Sons, Ltd.DOI
  27. [27]McCullagh, P. and Nelder, J. A. (1989). Generalized Linear Models. Routledge, New York, 2 edition.
  28. [28]Meyer, H. A. (1953). Forest Mensuration. Penns Valley Publishers, State College, PA.
  29. [29]Minka, T. (2000). Estimating a Dirichlet Distribution.
  30. [30]Minka, T. (2002). Estimating a Gamma Distribution.
  31. [31]Ortega, J. M. and Rheinboldt, W. C. (1970). Iterative Solution of Nonlinear Equations in Several Variables. SIAM.
  32. [32]Prescott, P. and Walden, A. T. (1980). Maximum likelihood estimation of the parameters of the generalized extreme-value distribution. Biometrika, 67(3):723–724.DOI
  33. [33]Rice, J. R. (1966). A Theory of Condition. SIAM Journal on Numerical Analysis, 3(2):287–310.DOI
  34. [34]Schou, G. (1978). Estimation of the Concentration Parameter in von Mises-Fisher Distributions. Biometrika, 65(2):369–377.DOI
  35. [35]Seber, G. A. F. and Lee, A. J. (2012). Linear Regression Analysis. John Wiley & Sons.
  36. [36]Self, S. G. and Liang, K.-Y. (1987). Asymptotic Properties of Maximum Likelihood Estimators and Likelihood Ratio Tests under Nonstandard Conditions. Journal of the American Statistical Association, 82(398):605–610.DOI
  37. [37]Smith, R. (1985). Maximum likelihood estimation in a class of nonregular cases. Biometrika, 72(1):67–90.DOI
  38. [38]van der Vaart, A. W. (1998). Asymptotic Statistics. Cambridge Series in Statistical and Probabilistic Mathematics. Cambridge University Press, Cambridge.
  39. [39]van der Vaart, A. W. and Wellner, J. A. (2023). Weak Convergence and Empirical Processes: With Applications to Statistics. Springer Series in Statistics. Springer International Publishing, Cham.
  40. [40]Willink, R. (2003). Relationships Between Central Moments and Cumulants, with Formulae for the Central Moments of Gamma Distributions. Communications in Statistics: Theory & Methods, 32(4):701.DOI
  41. [41]Wood, S. N. (2017). Generalized Additive Models: An Introduction with R, Second Edition. Chapman and Hall/CRC, New York, 2 edition.

Paper details

vlnc:W_z7wb7trv
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

ChatGPT · View statement →

Looking for related papers…
Contents