A saddlepoint approximation to the distribution of the sum of independent non‐identically beta random variables

The exact distribution of the sum of more than two independent beta random variables has not been known. Even in terms of approximations, only the normal approximation is known for the sum. Motivated by Murakami [Statistica Neerlandica, 2014, doi:10.1111/stan.12032], we derive here a saddlepoint approximation for the distribution of sum. An extensive simulation study shows that it always performs better than the normal approximation.     

A linear approximation to the power function of a test

In this paper we obtain a linear approximation to the power function of a test that is very accurate for small sample sizes. This is especially useful for robust tests where not many power functions are available. The approximation is based on the von Mises expansion of the tail probability functional and on the Tail Area Influence Function (TAIF). The goals of the paper are, first to extend the definition of the TAIF to the case of non identically distributed random variables, defining the Partial Tail Area Influence Functions and the Vectorial Tail Area Influence Function; second, to obtain exact expressions for computing these new influence functions; and, finally, to find accurate approximations to the power function, that can be used in the case of non identically distributed random variables. We include some examples of the application of this linear approximation to tests that involve the Huber statistic and also saddlepoint tests, so proving that the approximations apply not only to simple problems but also to complex ones.     

Testing for the bivariate Poisson distribution

This paper studies goodness-of-fit tests for the bivariate Poisson distribution. Specifically, we propose and study several Cramér–von Mises type tests based on the empirical probability generating function. They are consistent against fixed alternatives for adequate choices of the weight function involved in their definition. They are also able to detect local alternatives converging to the null at a certain rate. The bootstrap can be used to consistently estimate the null distribution of the test statistics. A simulation study investigates the goodness of the bootstrap approximation and compares their powers for finite sample sizes. Extensions for testing goodness-of-fit for the multivariate Poisson distribution are also discussed.     

Saddlepoint approximations to the mean and variance of the extended hypergeometric distribution

Conditional inference on 2 x 2 tables with fixed margins and unequal probabilities is based on the extended hypergeometric distribution. If the support of the distribution is large, exact calculation of the conditional mean and variance of the table entry may be computationally demanding. This paper proposes a single‐saddlepoint approximation to the mean and variance. While the approximation achieves acceptable accuracy for ordinary practical purposes, an alternative saddlepoint approximation is provided that gives much closer to exact results. It improves the accuracy of current approximations to up to more than four powers of ten.     

Saddlepoint approximations for the sum of independent non‐identically distributed binomial random variables

We discuss saddlepoint approximations to the distribution of the sum of independent non‐identically distributed binomial random variables. We examine the accuracy of the saddlepoint methods for a sum of 10 binomials with different sets of parameter values. The numerical results indicate that the saddlepoint approximations provide very accurate estimates for the probability mass function and the right‐tail probabilities for the cumulative distribution function of the sum.     

Another look at the Tail Area Influence Function

A new way of computing the Tail Area Influence Function (TAIF) exactly is proposed and a new finite sample robustness measure, based on the TAIF, is introduced. The main properties of this robustness measure are also studied, for both finite and asymptotic sample sizes. Next, a very accurate approximation to the finite sample power function of a test is obtained; this is based on the TAIF plus an iterative procedure. The results are valid when there are no nuisance parameters.     

Spline Regression in the Presence of Categorical Predictors

We consider the problem of estimating a relationship nonparametrically using regression splines when there exist both continuous and categorical predictors. We combine the global properties of regression splines with the local properties of categorical kernel functions to handle the presence of categorical predictors rather than resorting to sample splitting as is typically done to accommodate their presence. The resulting estimator possesses substantially better finite‐sample performance than either its frequency‐based peer or cross‐validated local linear kernel regression or even additive regression splines (when additivity does not hold). Theoretical underpinnings are provided and Monte Carlo simulations are undertaken to assess finite‐sample behavior; and two illustrative applications are provided. An implementation in R is available; see the R package ‘crs’ for details. Copyright © 2014 John Wiley & Sons, Ltd.     

Nonlinear regression and curved exponential families. Improvement of the approximation to the asymptotic distribution

Inference in nonlinear models is usually based on the asymptotic normal distribution, based on linearizing the model. The accuracy of this approximation can in many cases be improved by a reparametrization. Systematic methods for doing this will be described. Sometimes a saddlepoint approximation can be used, and this offers several advantages compared to the asymptotic distribution and the Edgeworth expansion. The improved methods are unfortunately not commonly used. It will be discussed why this is so. The methods will be illustrated by a series of examples.     

Estimation of nonlinear models with mismeasured regressors using marginal information

We consider the estimation of nonlinear models with mismeasured explanatory variables, when information on the marginal distribution of the true values of these variables is available. We derive a semi‐parametric MLE that is shown to be $\sqrt{n}$ consistent and asymptotically normally distributed. In a simulation experiment we find that the finite sample distribution of the estimator is close to the asymptotic approximation. The semi‐parametric MLE is applied to a duration model for AFDC welfare spells with misreported welfare benefits. The marginal distribution of the correctly measured welfare benefits is obtained from an administrative source. Copyright © 2010 John Wiley & Sons, Ltd.     

Forecast encompassing tests and probability forecasts

We consider tests of forecast encompassing for probability forecasts, for both quadratic and logarithmic scoring rules. We propose test statistics for the null of forecast encompassing, present the limiting distributions of the test statistics, and investigate the impact of estimating the forecasting models' parameters on these distributions. The small‐sample performance is investigated, in terms of small numbers of forecasts and model estimation sample sizes. We show the usefulness of the tests for the evaluation of recession probability forecasts from logit models with different leading indicators as explanatory variables, and for evaluating survey‐based probability forecasts. Copyright © 2009 John Wiley & Sons, Ltd.     

The Contribution of Structural Break Models to Forecasting Macroeconomic Series

This paper compares the forecasting performance of models that have been proposed for forecasting in the presence of structural breaks. They differ in their treatment of the break process, the model applied in each regime and the out‐of‐sample probability of a break. In an extensive empirical evaluation, we demonstrate the presence of breaks and their importance for forecasting. We find no single model that consistently works best in the presence of breaks. In many cases, the formal modeling of the break process is important in achieving a good forecast performance. However, there are also many cases where rolling window forecasts perform well. Copyright © 2014 John Wiley & Sons, Ltd.     

Approximating the distribution of Greenwood's statistic

Statistics based on uniform spacings are often used in goodness-of-fit problems. In this paper special attention is paid to the distribution of Greenwood's statistic. Although its asymptotic distribution is normal, the normal approximation is extremely bad, even for large sample sizes. It is shown that the Edge-worth expansion yields a considerably better approximation for the distribution of this statistic. Furthermore, an overview is given of the higher order asympto-tics for the sum of functions of uniform spacings, of which Greenwood's statistic is a special case.     

Replacing Sample Trimming with Boundary Correction in Nonparametric Estimation of First‐Price Auctions

Two‐step nonparametric estimators have become standard in empirical auctions. A drawback concerns boundary effects which cause inconsistencies near the endpoints of the support and bias in finite samples. To cope, sample trimming is typically used, which leads to non‐random data loss. Monte Carlo experiments show this leads to poor performance near the support boundaries and on the interior due to bandwidth selection issues. We propose a modification that employs boundary correction techniques, and we demonstrate substantial improvement in finite‐sample performance. We implement the new estimator using oil lease auctions data and find that trimming masks a substantial degree of bidder asymmetry and inefficiency in allocations. Copyright © 2014 John Wiley & Sons, Ltd.     

Exploring the relation between the r* approximation and the Edgeworth expansion

In this paper we study the relation between the r* saddlepoint approximation and the Edgeworth expansion when quite general assumptions for the statistic under consideration are fulfilled. We will show that the two term Edgeworth expansion approximates the r* formula up to an O(n ^?3/2) remainder; this provides a new way of looking at the order of the error of the r* approximation. This finding will be used to inspect the close connection between the r* formula and the Edgeworth B adjustment introduced in Phillips (Biometrika 65:91–98, 1978). We will show that, whenever an Edgeworth expansion exists, this adjustment approximates both the distribution function of the statistic and the r* formula to the same order degree as the Edgeworth expansion. Some numerical examples for the sample mean and U-statistics are given in order to shed light on the theoretical discussion.     

Second‐order analysis of anisotropic spatiotemporal point process data

Second‐order orientation methods provide a natural tool for the analysis of spatial point process data. In this paper, we extend to the spatiotemporal setting the spatial point pair orientation distribution function. The new space–time orientation distribution function is used to detect space–time anisotropic configurations. An edge‐corrected estimator is defined and illustrated through a simulation study. We apply the resulting estimator to data on the spatiotemporal distribution of fire ignition events caused by humans in a square area of 30 × 30 km² for 4 years. Our results confirm that our approach is able to detect directional components at distinct spatiotemporal scales. © 2014 The Authors. Statistica Neerlandica © 2014 VVS.     

De omvang van de steekproef bij een enkelvoudig steekproefsysteem

Summary  (Sample size for a single sampling scheme).
The operating characteristic of a sampling scheme may be specified by the producers 1 in 20 risk point ( p ₁), at which the probability of rejecting a batch is 0.05, and the consumers 1 in 20 risk point ( p ₂) at which the probability of accepting a batch of that quality is also 0.05.
A nomogram is given (fig. 2) to determine for single sampling schemes and for given values of p₁ and p ₂ the necessary sample size ( n ) and the allowable number of defectives in the sample ( c ).
The nomogram may reversedly be used to determine the producers and consumers 1 in 20 risk points for a given single sampling scheme.
The curves in this nomogram were computed from a table of percentage points of the χ² distribution. For v > 30 Wilson and Hilferty's approximation to the χ² distribution was used.     

A weighted Kuiper statistic for goodness of fit

A variance-weighted Kuiper statistic for goodness of fit is studied. The exact finite sample distribution can be obtained through modification of Noé's (1972) algorithm. Asymptotic distribution theory for the statistic is available from Jaeschke (1979) and Eicker (1979), but this theory does not lead to useful approximations with finite sample sizes less than 100. Monte Carlo power studies demonstrate that the weighted Kuiper statistic is especially sensitive to alternatives that are not stochastically ordered relative to the postulated null distribution.     

Bayesian semiparametric estimation of discrete duration models: an application of the dirichlet process prior

This paper proposes a Bayesian estimator for a discrete time duration model which incorporates a non‐parametric specification of the unobserved heterogeneity distribution, through the use of a Dirichlet process prior. This estimator offers distinct advantages over the Nonparametric Maximum Likelihood estimator of this model. First, it allows for exact finite sample inference. Second, it is easily estimated and mixed with flexible specifications of the baseline hazard. An application of the model to employment duration data from the Canadian province of New Brunswick is provided. Copyright © 2001 John Wiley & Sons, Ltd.     

Soft Questions,Hard Answers: Jacob Bernoulli's Probability in Historical Context

On the 300th anniversary of the publication of Jacob Bernoulli's Ars Conjectandi, the book is revisited in historical context and seen in a new light. Bernoulli's announced goal of extending the application of probability from games to civil life was frustrated by the a priori difficulty of specifying equally likely cases. The great approximation theorem Bernoulli proved and suggested could lead to empirical applications was impractical for all but unrealistically large sample sizes, but it was immensely influential in the development of the mathematics of probability.     

Estimating the returns to schooling: a likelihood approach based on normal mixtures

In this paper we develop likelihood‐based methods for statistical inference in a joint system of equations for the choice of length of schooling and earnings. The model for schooling choice is assumed to be an ordered probit model, whereas the earnings equation contains variables that are flexible transformations of schooling and experience, with corresponding coefficients that are allowed to be heterogeneous across individuals. Under the assumption that the distribution of the random terms of the model can be expressed as a finite mixture of multinormal distributions, we show that the joint probability distribution for schooling and earnings can be expressed on closed form. In an application of our method on Norwegian data, we find that the mixed Gaussian model offers a substantial improvement in fit to the (heavy‐tailed) empirical distribution of log‐earnings compared to a multinormal benchmark model. Copyright © 2010 John Wiley & Sons, Ltd.