Asymptotic expansions for the distribution of quadratic forms in normal variables

Higher order asymptotic expansions for the distribution of quadratic forms in normal variables are obtained. The Cornish-Fisher inverse expansions for the percentiles of the distribution are also given. The resulting formula for a definite quadratic form guarantees accuracy almost up to fourth decimal place if the distribution is not very skew. The normalizing transformation investigated by Jensen and Solomon (1972, J. Amer. Statist. Assoc., 67, 898–902) is reconsidered based on the rate of convergence to the normal distribution.Faculty of Science, Kyushu University, 6-10-1 Hakozaki, Higashi-ku, Fukuoka 812 Japan     

Weighted bootstrap for U-statistics

In this paper we investigate the weighted bootstrap for U-statistics and its properties. Under very general choices of random weights and certain regularity conditions, we show that the weighted bootstrap method with U-statistics provides second-order accurate approximations to the distribution of U-statistics. We shall prove this via one-term Edgeworth expansions of weighted U-statistics.     

Edgeworth expansion in censored linear regression model

For the censored simple linear regression model, we establish a oneterm Edgeworth expansion for the Koul, Susarla and Van Ryzin type estimator of the regression coefficient. Our approach is to represent the estimator of the regression coefficient as an asymptoticU-statistic plus some ignorable terms and hence apply the known results on the Edgeworth expansions for asymptoticU-statistic. The counting process and martingale techniques are used to provide the proof of the main results.     

Some distribution theory relating to confidence regions in multivariate calibration

Summary In the problem of multivariate calibration, Williams (1959,Regression Analysis, Wiley) and Wood (1982, to appear inProc. 11th Internat. Bio. Conf.) have proposed a decomposition of the usual Hotelling'sT ² statistic into the sum of two statistics for use in constructing confidence regions. This paper presents general results for the moment terms basic to Fujikoshi and Nishii's (1984,Hiroshima J. Math.,14, 215–225) approach to the distributions of these statistics, and presents simple alternative approximations to their percentiles.     

Error bounds for asymptotic expansion of the scale mixtures of the normal distribution

Summary LetX be a standard normal random variable and let σ be a positive random variable independent ofX. The distribution of η=σX is expanded around that ofN(0, 1) and its error bounds are obtained. Bounds are given in terms of E(σ ²V−σ ²−1) ^k whereσ ²Vσ ⁻² denotes the maximum of the two quantitiesσ ² andσ ⁻², andk is a positive integer, and of E(σ ²−1) ^k , ifk is even. The Institute of Statistical Mathematics     

Edgeworth Expansions for Studentized Versions of Symmetric Finite Population Statistics

We show the validity of the one-term Edgeworth expansion for Studentized asymptotically linear statistics based on samples drawn without replacement from finite populations. Replacing the moments defining the expansion by their estimators, we obtain an empirical Edgeworth expansion. We show the validity of the empirical Edgeworth expansion in probability.     

An Asymptotic Expansion for the Distribution of Hotelling'sT-Statistic under Nonnormality

In this paper we obtain an asymptotic expansion for the distribution of Hotelling'sT²-statisticT²under nonnormality when the sample size is large. In the derivation we find an explicit Edgeworth expansion of the multivariatet-statistic. Our method is to use the Edgeworth expansion and to expand the characteristic function ofT².     

Closer asymptotic approximations for the distributions of the power divergence goodness-of-fit statistics

Summary The members of the power divergence family of statistics all have an asymptotically equivalent χ² distribution (Cressie and Read [1]). An asymptotic expansion for the distribution function is derived which shows that the speed of convergence to this asymptotic limit is dependent on λ. Known results for Pearson'sX ² statistic and the log-likelihood ratio statistic then appear as special cases in a continuum rather than as separate (unrelated) expansions.     

Higher order asymptotics in estimation for two-sided Weibull type distributions

We consider the estimation problem of a location parameter on a sample of size n from a two-sided Weibull type density % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% GaamOzaiaacIcacaWG4bGaeyOeI0IaeqiUdeNaaiykaiabg2da9iaa% doeacaGGOaGaeqySdeMaaiykaiGacwgacaGG4bGaaiiCaiaacIcacq% GHsislcaGG8bGaamiEaiabgkHiTiabeI7aXjaacYhadaahaaWcbeqa% aiabeg7aHbaakiaacMcaaaa!52AD!\[f(x - \theta ) = C(\alpha )\exp ( - |x - \theta |^\alpha )\] for –<x<, –<< and 1<a<3/2, where % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gaam4qaiaacIcacqaHXoqycaGGPaGaeyypa0JaeqySdeMaai4laiaa% cUhacaaIYaGaeu4KdCKaaiikaiaaigdacaGGVaGaeqySdeMaaiykai% aac2haaaa!4B0E!\[C(\alpha ) = \alpha /\{ 2\Gamma (1/\alpha )\} \]. Then the bound for the distribution of asymptotically median unbiased estimators is obtained up to the 2a-th order, i.e., the order % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% GaamOBamaaCaaaleqabaGaeyOeI0IaaiikaiaaikdacqaHXoqycqGH% sislcaaIXaGaaiykaiaac+cacaaIYaaaaaaa!4444!\[n^{ - (2\alpha - 1)/2} \]. The asymptotic distribution of a maximum likelihood estimator (MLE) is also calculated up to the 2a-th order. It is shown that the MLE is not 2a-th order asymptotically efficient. The amount of the loss of asymptotic information of the MLE is given.     

Existence and nonexistence results for reaction-diffusion equations in product of cones

Problems of existence and nonexistence of global nontrivial solutions to quasilinear evolution differential inequalities in a product of cones are investigated. The proofs of the nonexistence results are based on the test-function method developed, for the case of the whole space, by Mitidieri, Pohozaev, Tesei and Véron. The existence result is established using the method of supersolutions.     

Normalizing transformations of some statistics of Gaussian ARMA processes

In this paper, the authors investigate Edgeworth type expansions of certain transformations of some statistics of Gaussian ARMA processes. They also investigate transformations which will make the second order part of the Edgeworth expansions vanish. Some numerical studies are made and they show that the above transformations give better approximations than the usual approximation.This work is supported by Contract N00014-K-0292 of the Office of Naval Research and Contract F49620-85-C-0008 of the Air Force Office of Scientific Research. The United States Government is authorized to reproduce and distribute reprints for governmental purposes notwithstanding any copyright notation hereon.The work of this author was done at the Center for Multivariate Analysis. His permanent address is Department of Mathematics, Hiroshima University, Naka-ku, Hiroshima 730, Japan.     

Probabilistic Choice in Games: Properties of Rosenthal’s t-Solutions

The t-solutions introduced in R. W. Rosenthal (1989, Int J Game Theory 18:273–292) are quantal response equilibria based on the linear probability model. Choice probabilities in t-solutions are related to the determination of leveling taxes in taxation problems. The set of t-solutions coincides with the set of Nash equilibria of a game with quadratic control costs. Evaluating the set of t-solutions for increasing values of t yields that players become increasingly capable of iteratively eliminating never-best replies and eventually only play rationalizable actions with positive probability. These features are not shared by logit quantal response equilibria. Moreover, there exists a path of t-solutions linking uniform randomization to Nash equilibrium     

A local breakdown property of robust tests in linear regression

The breakdown slope, as a useful summary measure of local stability for estimators and test statistics, has been studied recently by He, Simpson, and Protnoy (1990, J. Amer. Statist. Assoc., 85). It is shown here that all regression estimates based on residuals alone in linear models have zero breakdown slopes in contamination neighborhoods, even though they can have breakdown points as high as one-half. The breakdown functions of tests based on the S-estimation are investigated. It is also shown that the Generalized M-estimators can have better local breakdown robustness. One way to obtain regression estimators with desirable local and global breakdown properties is discussed.     

Certain Families of Elliptic Functions Defined by <Emphasis Type="Italic">q</Emphasis>-Series

We introduce certain families of elliptic functions involving the Weierstrass ℘-function via q-series systematically and investigate their analytic and algebraic properties. Especially, we give examples of non-linear algebraic ordinary differential equations satisfied by some such elliptic functions. 2000 Mathematics Subject Classification: Primary–11M36, 33E05     

On the behavio r of complete g raphs in $${\mathbb{R}\sp{n+1}}$$ with Lp-finite r-cu rvatu re

We study the behavior of complete graphs in with L ^p -finite r-curvature, that is, whose length of the r-th Newton transformation |P _r | is in L ^p , for some p ≥ 1. Moreover, we use a monotonicity formulae to establish an L ^p -lower bound for |P _r | in balls. As application, we prove some new Bernstein-type results.      

The generalized de Rham-Hodge theory aspects of Delsarte-Darboux type transformations in multidimension

The differential-geometric and topological structure of Delsarte transmutation operators and their associated Gelfand-Levitan-Marchenko type eqautions are studied along with classical Dirac type operator and its multidimensional affine extension, related with selfdual Yang-Mills eqautions. The construction of soliton-like solutions to the related set of nonlinear dynamical system is discussed.     

Optimal designs for estimating the coefficients of the lower frequencies in trigonometric regression models

In the common Fourier regression model we determine the optimal designs for estimating the coefficients corresponding to the lower frequencies. An analytical solution is provided which is found by an alternative characterization of c-optimal designs. Several examples are provided and the performance of the D-optimal design with respect to the estimation of the lower order coefficients is investigated. The results give a complete answer to an open question which was recently raised in the literature.     

On the Wl q-Regularity of Solutions to Systems of Differential Equations in the Case When the Equations Are Constructed from Discontinuous Functions

Some solution, final in a sense from the standpoint of the theory of Sobolev spaces, is obtained to the problem of regularity of solutions to a system of (generally) nonlinear partial differential equations in the case when the system is locally close to elliptic systems of linear equations with constant coefficients. The main consequences of this result are Theorems 5 and 8. According to the first of them, the higher derivatives of an elliptic C ^l -smooth solution to a system of lth-order nonlinear partial differential equations constructed from C ^l -smooth functions meet the local Hoelder condition with every exponent , 0<<1. Theorem 8 claims that if a system of linear partial differential equations of order l with measurable coefficients and right-hand sides is uniformly elliptic then, under the hypothesis of a (sufficiently) slow variation of its leading coefficients, the degree of local integrability of lth-order partial derivatives of every W ^l _q,loc-solution, q>1, to the system coincides with the degree of local integrability of lower coefficients and right-hand sides.     

Fuzzy Relation Equations for Compression/Decompression Processes of Colour Images in the RGB and YUV Colour Spaces

We use particular fuzzy relation equations for compression/decompression of colour images in the RGB and YUV spaces, by comparing the results of the reconstructed images obtained in both cases. Our tests are made over well known images of 256×256 pixels (8 bits per pixel in each band) extracted from Corel Gallery. After the decomposition of each image in the three bands of the RGB and YUV colour spaces, the compression is performed using fuzzy relation equations of “min - →_t” type, where “t” is the Lukasiewicz t-norm and “→_t” is its residuum. Any image is subdivided in blocks and each block is compressed by optimizing a parameter inserted in the Gaussian membership functions of the fuzzy sets, used as coders in the fuzzy equations. The decompression process is realized via a fuzzy relation equation of max-t type. In both RGB and YUV spaces we evaluate and compare the root means square error (RMSE) and the consequentpeak signal to noise ratio (PSNR) on the decompressed images with respect to the original image under several compression rates.     

On the area swept under the occupation process of an M/M/1 queue in a busy period

We compute in this paper the distribution of the area swept under the occupation process of an M/M/1 queue during a busy period. For this purpose, we use the expression of the Laplace transform of the random variable established in earlier studies as a fraction of Bessel functions. To get information on the poles and the residues of , we take benefit of the fact that this function can be represented by a continued fraction. We then show that this continued fraction is the even part of an S fraction and we identify its successive denominators by means of Lommel polynomials. This allows us to numerically evaluate the poles and the residues. Numerical evidence shows that the poles are very close to the numbers as . This motivated us to formulate some conjectures, which lead to the derivation of the asymptotic behaviour of the poles and the residues. This is finally used to derive the asymptotic behaviour of the probability survivor function . The outstanding property of the random variable is that the poles accumulate at 0 and its tail does not exhibit a nice exponential decay but a decay of the form for some positive constants c and , which indicates that the random variable has a Weibull-like tail. This revised version was published online in June 2006 with corrections to the Cover Date.