A SADDLEPOINT APPROXIMATION TO THE DISTRIBUTION FUNCTION OF THE ANDERSON-DARLING TEST STATISTIC

The Anderson-Darling goodness-of-fit test has a highly skewed and non-standard limit distribution. Various attempts have been made to tabulate the associated critical points, using both theoretical approximations and simulation methods. We show that a standard saddlepoint approximation performs well in both tails of the distribution. It is markedly superior to other theoretical approximations in the lower tail of the distribution.     

Saddlepoint Approximations for a Product and Its Application

A technique of the saddlepoint approximation with double exponential base, SPA_D is developed to evaluate the probability of a product of two random variables, which may be independent or dependent, normal or contaminated normal random variables. The SPA_D shows a slightly better approximation as compared to the saddlepoint approximation with Lagannani–Rice formula. However, both methods get remarkable results when applied to evaluate the tail probabilities of the Reynolds stress for soil erosion prediction.     

Saddlepoint Approximations for the Difference of Order Statistics and Studentized Sample Quantiles

For a sample from a given distribution the difference of two order statistics and the Studentized quantile are statistics whose distribution is needed to obtain tests and confidence intervals for quantiles and quantile differences. This paper gives saddlepoint approximations for densities and saddlepoint approximations of the Lugannani–Rice form for tail probabilities of these statistics. The relative errors of the approximations are n ⁻¹ uniformly in a neighbourhood of the parameters and this uniformity is global if the densities are log-concave.     

Saddlepoint Approximations for Stopped-Sum Distributions

One of the common used classes of distributions is the stopped-sum class. This class includes Hermite distribution, Polya–Aeppli distribution, Poisson-Gamma distribution, and Neyman type A. This article introduces the saddlepoint approximations to the stopped-sum class in continuous and discrete settings. We discuss approximations for mass/density and cumulative distribution functions of stopped-sum distributions. Examples of continuous and discrete distributions from the Poisson stopped-sum class are presented. Comparisons between saddlepoint approximations and the exact calculations show the great accuracy of the saddlepoint methods.     

Improved Measures of the Spread of Data for Some Unknown Complex Distributions Using Saddlepoint Approximations

Measures of the spread of data for random sums arise frequently in many problems and have a wide range of applications in real life, such as in the insurance field (e.g., the total claim size in a portfolio). The exact distribution of random sums is extremely difficult to determine, and normal approximation usually performs very badly for this complex distributions. A better method of approximating a random-sum distribution involves the use of saddlepoint approximations.

Saddlepoint approximations are powerful tools for providing accurate expressions for distribution functions that are not known in closed form. This method not only yields an accurate approximation near the center of the distribution but also controls the relative error in the far tail of the distribution.

In this article, we discuss approximations to the unknown complex random-sum Poisson–Erlang random variable, which has a continuous distribution, and the random-sum Poisson-negative binomial random variable, which has a discrete distribution. We show that the saddlepoint approximation method is not only quick, dependable, stable, and accurate enough for general statistical inference but is also applicable without deep knowledge of probability theory. Numerical examples of application of the saddlepoint approximation method to continuous and discrete random-sum Poisson distributions are presented.     

Saddlepoint p-values for group of linear rank tests

Many nonparametric tests in one sample problem, matched pairs, and competingrisks under censoring have the same underlying permutation distribution. This article proposes a saddlepoint approximation to the exact p-values of these tests instead of the asymptotic approximations. The performance of the saddlepoint approximation is assessed by using simulation studies that show the superiority of the saddlepoint methods over the asymptotic approximations in several settings. The use of the saddlepoint to approximate the p-values of class of two sample tests under complete randomized design is also discussed.     

Saddlepoint Approximations to the Density and the Distribution Functions of Linear Combinations of Ratios of Partial Sums

A class of ratios of partial sums, including Normal, Weibull, Gamma, and Exponential distributions, is considered. The distribution of a linear combination of ratios of partial sums from this class is characterized by the distribution of a linear combination of Dirichlet components. This article presents two saddlepoint approaches to calculate the density and the distribution function for such a class of linear combinations. A simulation study is conducted to assess the performance of the saddlepoint methods and shows the great accuracy of the approximations over the usual asymptotic approximation. Applications of the presented approximations in statistical inferences are discussed.     

A Saddlepoint Approximation to the Distribution of the Half-Life Estimator in a Stationary Autoregressive Model

We derive saddlepoint approximations for the distribution and density functions of the half-life estimated by OLS from autoregressive time-series models. Our results are used to prove that none of the integer-order moments of these half-life estimators exist. This provides an explanation for the very large estimates of persistency, and the extremely wide confidence intervals, that have been reported by various authors, i.e., in the empirical economics literature relating to purchasing power parity.     

Saddlepoint Approximation for Sample Quantiles with Some Applications

In this article, we use the integral form of the binomial distribution to derive saddlepoint approximations for sample quantiles. As an application, we present the calculation of the tail probability of the empirical log-likelihood ratio statistic for quantiles. Simulation results are also given to show that our approximations are extremely accurate.     

Partial Saddlepoint Approximations for Transformed Means

The full saddlepoint approximation for real valued smooth functions of means requires the existence of the joint cumulant generating function for the entire vector of random variables which are being transformed. We propose a mixed saddlepoint-Edgeworth approximation requiring the existence of a cumulant generating function for only part of the random vector considered, while retaining partially the relative nature of the errors. Tail probability approximations are obtained under conditions which enable the approximations to be used in resampling situations and hence to obtain a result on the relative error of coverage in the case of the bootstrap approximation to the confidence interval for the Studentized mean.     

Kolmogorov-Smirnov Test Statistic and Critical Values for the Erlang-3 and Erlang-4 Distributions

Following a procedure applied to the Erlang-2 distribution in a recent paper, an adjusted Kolmogorov-Smirnov statistic and critical values are developed for the Erlang-3 and -4 cases using data from Monte Carlo simulations. The test statistic produced features of compactness and ease of implementation. It is quite accurate for sample sizes as low as ten.     

Tail Exactness of Multivariate Saddlepoint Approximations

We consider a log-concave density f in R ^m satisfying certain weak conditions, particularly on the Hessian matrix of log f . For such a density, we prove tail exactness of the multivariate saddlepoint approximation. The proof is based on a local limit theorem for the exponential family generated by f . However, the result refers not to asymptotic behaviour under repeated sampling, but to a limiting property at the boundary of the domain of f . Our approach does not apply any complex analysis, but relies totally on convex analysis and exponential models arguments.     

Testing the Difference of Two Binomial Proportions: Comparison of Continuity Corrections for Saddlepoint Approximation

We carried out a simulation study based on the methodology of Newcombe (1998 Newcombe , R. G. ( 1998 ). Interval estimation for the difference between independent proportions: comparison of eleven methods . Statist. Med. 17 : 873 – 890 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) to compare tests for the difference of two binomial proportions by applying different continuity corrections on saddlepoint approximation to tail probabilities. In this article, we proposed a new continuity correction based on the least common multiple of two sample sizes. We evaluated that the best test should have the actual Type I error rates that are, on the whole, closest to α, but not exceeding α, where α is nominal level of significance.     

Saddlepoint approximations to tail expectations under non-Gaussian base distributions: option pricing applications

The saddlepoint approximation formulas provide versatile tools for analytic approximation of the tail expectation of a random variable by approximating the complex Laplace integral of the tail expectation expressed in terms of the cumulant generating function of the random variable. We generalize the saddlepoint approximation formulas for calculating tail expectations from the usual Gaussian base distribution to an arbitrary base distribution. Specific discussion is presented on the criteria of choosing the base distribution that fits better the underlying distribution. Numerical performance and comparison of accuracy are made among different saddlepoint approximation formulas. Improved accuracy of the saddlepoint approximations to tail expectations is revealed when proper base distributions are chosen. We also demonstrate enhanced accuracy of the generalized saddlepoint approximation formulas under non-Gaussian base distributions in pricing European options on continuous integrated variance under the Heston stochastic volatility model.     

General Saddlepoint Approximation for Testing Separate Location-Scale Families of Hypotheses

For testing separate families of hypotheses, the likelihood ratio test does not have the usual asymptotic properties. This paper considers the asymptotic distribution of the ratio of maximized likelihoods (RML) statistic in the special case of testing separate scale or location-scale families of distributions. We derive saddlepoint approximations to the density and tail probabilities of the log of the RML statistic. These approximations are based on the expansion of the log of the RML statistic up to the second order, which is shown not to depend on the location and scale parameters. The resulting approximations are applied in several cases, including normal versus Laplace, normal versus Cauchy, and Weibull versus log-normal. Our results show that the saddlepoint approximations are satisfactory, even for fairly small sample sizes, and are more accurate than normal approximations and Edgeworth approximations, especially for tail probabilities that are the values of main interest in hypothesis testing problems.     

Saddlepoint Approximation for the Distribution of the Modified Signed Root of Likelihood Ratio Statistics Near the Mean

The Tracy-Singh product, which can be viewed as generalized Kronecher product, is studied and used widely in matrix theory and statistics. Therefore, some well-known Kantorovich inequalities are generalized to Tracy-Singh product in this article. Furthermore, some applications to statistics are presented.     

A Proof of the Asymptotic Equivalence of Two-Tail Probability Approximations

A modified efficient jump algorithm is proposed for the Markov Chain Monte Carlo draws of the exponential power distribution. Bayesian inference based on the exponential power error term and that on the normal error term are compared. Unbiasedness of the LAD estimator is proven.     

Monte Carlo Approximations of the Quantiles of a Sample Statistic

We consider in this article the problem of numerically approximating the quantiles of a sample statistic for a given population, a problem of interest in many applications, such as bootstrap confidence intervals. The proposed Monte Carlo method can be routinely applied to handle complex problems that lack analytical results. Furthermore, the method yields estimates of the quantiles of a sample statistic of any sample size though Monte Carlo simulations for only two optimally selected sample sizes are needed. An analysis of the Monte Carlo design is performed to obtain the optimal choices of these two sample sizes and the number of simulated samples required for each sample size. Theoretical results are presented for the bias and variance of the numerical method proposed. The results developed are illustrated via simulation studies for the classical problem of estimating a bivariate linear structural relationship. It is seen that the size of the simulated samples used in the Monte Carlo method does not have to be very large and the method provides a better approximation to quantiles than those based on an asymptotic normal theory for skewed sampling distributions.     

Alternative Approximations to Value-At-Risk: A Comparison

This article compares three value-at-risk (VaR) approximation methods suggested in the literature: Cornish and Fisher (1937 Cornish, E.A., Fisher, R.A. (1937). Moments and cumulants in the specification of distributions. Revue de l’Institut International de Statistique 5:307–320.[Crossref] , [Google Scholar]), Sillitto (1969 Sillitto, G.P. (1969). Derivation of approximants to the inverse distribution function of a continuous univariate population from the order statistics of a sample. Biometrika 56:641–650.[Crossref], [Web of Science ®] , [Google Scholar]), and Liu (2010 Liu, W.-H. (2010). Estimation and testing of portfolio value-at-risk based on L-comoment matrices. Journal of Futures Markets 30:897–908.[Crossref], [Web of Science ®] , [Google Scholar]). Simulation results are obtained for three families of distributions: student-t, skewed-normal, and skewed-t. We recommend the Sillitto approximation as the best method to evaluate the VaR when the financial return has an unknown, skewed, and heavy-tailed distribution.     

Saddlepoint density and distribution functions for the ratio of two linear functions and the product of generalized gamma variates

The generalized gamma distribution is a flexible and attractive distribution because it incorporates several well-known distributions, i.e., gamma, Weibull, Rayleigh, and Maxwell. This article derives saddlepoint density and distribution functions for the ratio of two linear functions of generalized gamma variables and the product of n independent generalized gamma variables. Simulation studies are used to evaluate the accuracy of the saddlepoint approximations. The saddlepoint approximations are fast, easy, and very accurate.