Empirical characteristic function approach to goodness-of-fit tests for the Cauchy distribution with parameters estimated by MLE or EISE

We consider goodness-of-fit tests of the Cauchy distribution based on weighted integrals of the squared distance between the empirical characteristic function of the standardized data and the characteristic function of the standard Cauchy distribution. For standardization of data Gürtler and Henze (2000,Annals of the Institute of Statistical Mathematics,52, 267–286) used the median and the interquartile range. In this paper we use the maximum likelihood estimator (MLE) and an equivariant integrated squared error estimator (EISE), which minimizes the weighted integral. We derive an explicit form of the asymptotic covariance function of the characteristic function process with parameters estimated by the MLE or the EISE. The eigenvalues of the covariance function are numerically evaluated and the asymptotic distributions of the test statistics are obtained by the residue theorem. A simulation study shows that the proposed tests compare well to tests proposed by Gürtler and Henze and more traditional tests based on the empirical distribution function.     

Tail indices for $$\mathbf{A}\mathbf{X}+\mathbf{B}$$ A X + B Recursion with Triangular Matrices

Journal of Theoretical Probability - We study multivariate stochastic recurrence equations (SREs) with triangular matrices. If coefficient matrices of SREs have strictly positive entries, the...     

Prediction of Components in Random Sums

We consider predictions of the random number and the magnitude of each iid component in a random sum based on its distributional structure, where only a total value of the sum is available and where iid random components are non-negative. The problem is motivated by prediction problems in a Poisson shot noise process. In the context, although conditional moments are best possible predictors under the mean square error, only a few special cases have been investigated because of numerical difficulties. We replace the prediction problem of the process with that of a random sum, which is more general, and establish effective numerical procedures. The methods are based on conditional technique together with the Panjer recursion and the Fourier transform. In view of numerical experiments, procedures work reasonably. An application in the compound mixed Poisson process is also suggested.     

Integral representations of one-dimensional projections for multivariate stable densities

We consider the numerical evaluation of one-dimensional projections of general multivariate stable densities introduced by Abdul-Hamid and Nolan [H. Abdul-Hamid, J.P. Nolan, Multivariate stable densities as functions of one dimensional projections, J. Multivariate Anal. 67 (1998) 80-89]. In their approach higher order derivatives of one-dimensional densities are used, which seems to be cumbersome in practice. Furthermore there are some difficulties for even dimensions. In order to overcome these difficulties we obtain the explicit finite-interval integral representation of one-dimensional projections for all dimensions. For this purpose we utilize the imaginary part of complex integration, whose real part corresponds to the derivative of the one-dimensional inversion formula. We also give summaries on relations between various parametrizations of stable multivariate density and its one-dimensional projection.     

Estimation of the tail index for lattice-valued sequences

If one applies the Hill, Pickands or Dekkers–Einmahl–de Haan estimators of the tail index of a distribution to data which are rounded off one often observes that these estimators oscillate strongly as a function of the number k of order statistics involved. We study this phenomenon in the case of a Pareto distribution. We provide formulas for the expected value and variance of the Hill estimator and give bounds on k when the central limit theorem is still applicable. We illustrate the theory by using simulated and real-life data.