Modified Normal-based Approximation to the Percentiles of Linear Combination of Independent Random Variables with Applications

A modified normal-based approximation for calculating the percentiles of a linear combination of independent random variables is proposed. This approximation is applicable in situations where expectations and percentiles of the individual random variables can be readily obtained. The merits of the approximation are evaluated for the chi-square and beta distributions using Monte Carlo simulation. An approximation to the percentiles of the ratio of two independent random variables is also given. Solutions based on the approximations are given for some classical problems such as interval estimation of the normal coefficient of variation, survival probability, the difference between or the ratio of two binomial proportions, and for some other problems. Furthermore, approximation to the percentiles of a doubly noncentral F distribution is also given. For all the problems considered, the approximation provides simple satisfactory solutions. Two examples are given to show applications of the approximation.     

An F Approximation to the Distribution of a Linear Combination of Chi-squared Variables.

A three-parameter F approximation to the distribution of a positive linear combination of central chi-squared variables is described. It is about as easy to implement as the Satterthwaite-Welsh and Hall-Buckley-Eagleson approximations. Some reassuring properties of the F approximation are derived, and numerical results are presented. The numerical results indicate that the new approximation is superior to the Satterthwaite approximation and, for some purposes, better than the Hall-Buckley-Eagleson approximation. It is not quite as good as the Gamma-Weibull approximation due to Solomon and Stephens, but is easier to implement because iterative methods are not required.     

A new approximation to the distribution function of the studentized maximum root

The Studentized maximum root (SMR) distribution is useful for constructing simultaneous confidence intervals around product interaction contrasts in replicated two-way ANOVA. A three-moment approximation to the SMR distribution is proposed. The approximation requires the first three moments of the maximum root of a central Wishart matrix. These values are obtained by means of numerical integration. The accuracy of the approximation is compared to the accuracy of a two-moment approximation for selected two-way table sizes. Both approximations are reasonably accurate. The three-moment approximation is generally superior.     

A comparison of tests of homogeneity for sparse contingency tables

Five tests of homogeneity for a 2x(k+l) contingency table are compared using Monte Carlo techniques. For these studiesit is assumed that k becomes large in such a way that thecontingency table is sparse for 2xk of the cells, but the sample size in two of the cells remains large. The test statistics studied are: the chi-square approximation to the Pearson test statistic, the chi-square approximation to the likelihood ratio statistic, the normal approximation to Zelterman's (1984)the normal approximation to Pearson's chi-square, and the normal approximation to the likelihood ratio statistic. For the range of parameters studied the chi-square approximation to Pearson's statistic performs consistently well with regard to its size and power.     

Saddlepoint approximations to sensitivities of tail probabilities of random sums and comparisons with Monte Carlo estimators

This article proposes computing sensitivities of upper tail probabilities of random sums by the saddlepoint approximation. The considered sensitivity is the derivative of the upper tail probability with respect to the parameter of the summation index distribution. Random sums with Poisson or Geometric distributed summation indices and Gamma or Weibull distributed summands are considered. The score method with importance sampling is considered as an alternative approximation. Numerical studies show that the saddlepoint approximation and the method of score with importance sampling are very accurate. But the saddlepoint approximation is substantially faster than the score method with importance sampling. Thus, the suggested saddlepoint approximation can be conveniently used in various scientific problems.     

The null distribution of multivariate kurtosis

The paper introduces a x²-approximation to multivariate kurtosis b_2,punder normality. It requires calculating the third moment of b_2,pwhich is obtained. We compare the approximation with simulated percentage points and the normal approximation, and find it to be adequate for p=l and 2. For p=3, the simple average of this estimate and the normal approximation is found to be generally superior to either approximation on its own. For p=4, the normal approximation is best for non-extreme values of ∝     

Saddlepoint approach to inference for response probabilities under the logistic response model

This paper studies lower confidence limits of response probabilities based on sensitivity testing data set. The saddlepoint approximation to a conditional distribution is developed. Based on it we give a modified algorithm to find approximate confidence limits for the parameter of interest. A simulation study shows that the saddlepoint approximation with proper corrections gives better coverage probability than the direct saddlepoint approximation and the asymptotic normality approximation. Finally, we apply the proposed approximation to a real data set.     

A new approximation to the incomplete beta function

A new approximation to the incomplete beta function is proposed. This approximation compares favourably well with the Mudholkar-Chaubey (1976) approximation and it is infact superior for equal degrees of freedom. Moreover’ it is simpler to evaluate.     

Bivariate symmetry tests for complete and competing risks data: a saddlepoint approach

A class of bivariate symmetry tests for complete data and competing risks data is considered. Saddlepoint approximation for the exact p-values of the underlying permutation distribution of these tests is derived. Several simulation studies are conducted to evaluate the performance of the saddlepoint approximation and the asymptotic approximation. The saddlepoint approximation was found to be highly accurate and superior to the asymptotic approximations in replicating the exact permutation significance.     

Sequential determination of the number of bootstrap samples

The primary purpose of this paper is that of developing a sequential Monte Carlo approximation to an ideal bootstrap estimate of the parameter of interest. Using the concept of fixed-precision approximation, we construct a sequential stopping rule for determining the number of bootstrap samples to be taken in order to achieve a specified precision of the Monte Carlo approximation. It is shown that the sequential Monte Carlo approximation is asymptotically efficient in the problems of estimation of the bias and standard error of a given statistic. Efficient bootstrap resampling is discussed and a numerical study is carried out for illustrating the obtained theoretical results.     

Approximations for the doubly noncentral-F distribution

The approximate normality of the cube root of the noncentral chi-square observed by Aty (1954) and an Edgeworth-series expansion are used to derive an approximation for the doubly noncentral-F distribution. Another approximation in terms of a noncentral-F distribution is also proposed. Both these approximations are seen to compare favorably with some earlier approximations due to Das Gupta (1968) and Tiku (1972). The problem of approximating the cumulants of the doubly noncentral-F variable, which is pivotal in Tiku’s approximation, is examined and use of a noncentral-F distribution is seen to provide a good solution for it. A FORTRAN routine for the Edgeworth-series approximation is given.     

A simple approximation for the doubly noncentral t-distribution

A simple approximation for the doubly noncentral t-distribution, based upon the Fieller-Geary Theorem (1930) and approximate normality of the square root of the noncentral chi-square variable observed by Patnaik (1949), is developed, This approximation and an Edgeworth series expansion associated with it are evaluated. The simple approximation is seen to be reasonably accurate for most practical purposes.     

FITTING SPATIAL CORRELATION MODELS: APPROXIMATING A LIKELIHOOD APPROXIMATION

The paper considers a class of spatial correlation models (stationary Gaussian processes) which includes (spatial) conditional autoregressive, simultaneous autoregressive, moving average and direct covariance models. Given observations on a finite rectangular lattice, a likelihood approximation for estimating the parameters in the spectral density of the model is discussed. The approximation consists of applying the trapezoidal rule, with a her grid of frequencies than the usual Fourier frequencies, to compute the integral in an appraximation due to Whittle (1954) and later modified by Guyon (1984). With this approximation, a Fisher scoring type algorithm has a simple form and in some casea reduces to iteratively reweighted least squares. Methods for computing the unbiased two-dimensional periodogram required by the method are presented and the accuracy of the approximation is discussed. The asymptotic distribution of the parameter estimates computed from the likelihood approximation is also given.     

Computationally tractable approximate and smoothed Polya trees

A discrete approximation to the Polya tree prior suitable for latent data is proposed that enjoys surprisingly simple and efficient conjugate updating. This approximation is illustrated in two applied contexts: the implementation of a nonparametric meta-analysis involving studies on the relationship between alcohol consumption and breast cancer, and random intercept Poisson regression for Ache armadillo hunting treks. The discrete approximation is then smoothed with Gaussian kernels to provide a smooth density for use with continuous data; the smoothed approximation is illustrated on a classic dataset on galaxy velocities and on recent data involving breast cancer survival in Louisiana.     

An approximation to the exact distribution of the wilcoxon-mann-whitney rank sum test statistic

An approximation to the exact distribution of the Wilcoxon rank sum test (Mann-Whitney U-test) and the Siegel-Tukey test based on a linear combination of the two-sample t-test applied to ranks and the normal approximation is compared with the usual normal approximation. The normal approximation results in a conservative test in the tails while the linear combination of the test statistics provides a test that has a very high percentage of agreement with tables of the exact distribution. Sample sizes 3≤m, n≤50 were considered.     

SMALL SAMPLE INFERENCE IN LOCATION SCALE MODELS WITH STUDENT(λ) ERRORS

ABSTRACT

A third order accurate approximation to the p value in testing either the location or scale parameter in a location scale model with Student(λ) errors is introduced. The third order approximation is developed via an asymptotic method, based on exponential models and the saddlepoint approximation. Techniques are presented for the numerical computation of all quantities required for the third order approximation. To compare the accuracy of various asymptotic methods a numerical example and simulation study are included. The numerical example and simulation study illustrate that the third order method presented leads to a more accurate p value approximation compared to first order methods in Student(λ) models with small samples.     

Moment closure based parameter inference of stochastic kinetic models

Parameter inference for stochastic kinetic models is a topic that spans many disciplines. Although it is possible to carry out exact inference using partial observations of a stochastic process, it is often computationally impractical. In this paper we use the moment closure approximation of the underlying stochastic process as a fast approximation of the likelihood. We show that this approximation is fast and accurate, even when the population numbers are small.     

负二项分布的两种近似分布及其比较

负二项分布是一个重要的离散型随机变量的分布,可以用泊松分布和正态分布作为其近似分布,通过对两种近似分布进行比较分析,结果表明：在参数q很小时,泊松近似的精度好于正态近似,而且在参数q很小时,即便r不是很大,用泊松分布也能获得负二项分布较好的近似;当参数q较大时,泊松近似效果不好,相比之下,正态近似的结果不错。     

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.     

An error-bounded polynomial approximation for bivariate normal probabilities

Although the bivariate normal distribution is frequently employed in the development of screening models, the formulae for computing bivariate normal probabilities are quite complicated. A simple and accurate error-bounded, noniterative approximation for bivariate normal probabilities based on a simple univariate normal quadratic or cubic approximation is developed for use in screening applications. The approximation, which is most accurate for large absolute correlation coefficients, is especially suitable for screening applications (e.g., in quality control), where large absolute correlations between performance and screening variables are desired. A special approximation for conditional bivariate normal probabilities is also provided which in quality control screening applications improves the accuracy of estimating the average outgoing product quality. Some anomalies in computing conditional bivariate normal probabilities using BNRDF and NORDF in IMSL are also discussed.