Laplace record data

Basic properties of upper record values _XT(1),_XT(2),…,_XT(n)$X_{T(1)},X_{T(2)},\dots ,X_{T(n)}$ from a symmetric two-parameter Laplace distribution are established. In particular, unimodality of the density function and the exact expression of the mode are derived. Moreover, we obtain approximations of the first and second moment and the variance of _XT(k)$X_{T(k)}$ which provide close approximations even for moderate k. Additionally, limit laws and simulation of Laplace records are considered. Finally, we discuss maximum likelihood estimation in a location-scale family of Laplace distributions. We obtain nice representations of the estimators provided that the location parameter is unknown and present interesting properties of the established estimators. Some illustrative examples complete the presentation.     

Representing a function of two variables as a superposition over angles of ridge functions with a given shape via linear programming

A ridge function with shape function g   in the horizontal direction is a function of the form g(x)h(y,0)$g(x)h(y,0)$. Along each horizontal line it has the shape g(x)$g(x)$, multiplied by a function h(y,0)$h(y,0)$ which depends on the y-value of the horizontal line. Similarly a ridge function with shape function g   in the vertical direction has the form g(y)h(x,π/2)$g(y)h(x,\pi /2)$. For a given shape function g it may or may not be possible to represent an arbitrary   function f(x,y)$f(x,y)$ as a superposition over all angles of a ridge function with shape g   in each direction, where h=_hf=_hf,g$h=h_f=h_{f,g}$ depends on the functions f and g   and also on the direction, θ:h=_hf,g(·,θ)$\theta :h=h_{f,g}(·,\theta )$. We show that if g   is Gaussian centered at zero then this is always possible and we give the function _hf,g$h_{f,g}$ for a given f(x,y)$f(x,y)$. For highpass or for odd shapes g  , we show it is impossible to represent an arbitrary f(x,y)$f(x,y)$, i.e. in general there is no _hf,g$h_{f,g}$. Note that our problem is similar to tomography, where the problem is to invert the Radon transform, except that the use of the word inversion is here somewhat “inverted”: in tomography f(x,y)$f(x,y)$ is unknown and we find it by inverting the projections of f  ; here, f(x,y)$f(x,y)$ is known, g(z)$g(z)$ is known, and _hf(·,θ)=_hf,g(·,θ)$h_f(·,\theta )=h_{f,g}(·,\theta )$ is the unknown.     

The q-cluster distribution

The paper studies the three-parameter generalization of the logarithmic distribution that is obtained as the cluster distribution for the generalized Euler distribution. The diagnostic statistic, R(x)=_xpx/[(x-1)_px-1]$R(x)={\mathit{xp}}_x/[(x-1)p_{x-1}]$, is constant for the logarithmic distribution. For the new distribution it can decrease, stay constant, or increase as x increases. The relative values of the extra parameters determine the flatness/hollowness of the distribution and its tail behaviour. Kemp's q-logarithmic distribution and the Euler cluster distribution are special cases. Fitted data sets illustrate the properties of the distribution and its limiting forms.     

Truncated linear estimation of a bounded multivariate normal mean

We consider the problem of estimating the mean θ$\theta$ of an _Np(θ,_Ip)${\mathrm{N}}_{\mathrm{p}}(\theta ,I_p)$ distribution with squared error loss ∥δ−θ∥²$\parallel \delta -\theta {\parallel }^2$ and under the constraint ∥θ∥≤m$\parallel \theta \parallel \le m$, for some constant m>0$m>0$. Using Stein's identity to obtain unbiased estimates of risk, Karlin's sign change arguments, and conditional risk analysis, we compare the risk performance of truncated linear estimators with that of the maximum likelihood estimator _δmle${\delta }_{\mathrm{mle}}$. We obtain for fixed (m,p)$(m,p)$ sufficient conditions for dominance. An asymptotic framework is developed, where we demonstrate that the truncated linear minimax estimator dominates _δmle${\delta }_{\mathrm{mle}}$, and where we obtain simple and accurate measures of relative improvement in risk. Numerical evaluations illustrate the effectiveness of the asymptotic framework for approximating the risks for moderate or large values of p.     

Prediction in moving average processes

For the stationary invertible moving average process of order one with unknown innovation distribution F, we construct root-n   consistent plug-in estimators of conditional expectations E(h(_Xn+1)|_X1,…,_Xn)$E(h(X_{n+1})|X_1,\dots ,X_n)$. More specifically, we give weak conditions under which such estimators admit Bahadur-type representations, assuming some smoothness of h or of F. For fixed h it suffices that h   is locally of bounded variation and locally Lipschitz in _L2(F)$L_2(F)$, and that the convolution of h and F   is continuously differentiable. A uniform representation for the plug-in estimator of the conditional distribution function P(_Xn+1?·|_X1,…,_Xn)$P(X_{n+1}?·|X_1,\dots ,X_n)$ holds if F has a uniformly continuous density. For a smoothed version of our estimator, the Bahadur representation holds uniformly over each class of functions h that have an appropriate envelope and whose shifts are F-Donsker, assuming some smoothness of F. The proofs use empirical process arguments.     

On asymptotic properties of the parameter estimator for a type of SPDE

In this paper, we study a random field _U?(t,x)$U_{?}(t,x)$ governed by some type of stochastic partial differential equations with an unknown parameter θ$\theta$ and a small noise ?$?$. We construct an estimator of θ$\theta$ based on the continuous observation of N   Fourier coefficients of _U?(t,x)$U_{?}(t,x)$, and prove the strong convergence and asymptotic normality of the estimator when the noise ?$?$ tends to zero.     

On the maximum of periodic integer-valued sequences with exponential type tails via max-semistable laws

In this work we study the limiting distribution of the maximum term of periodic integer-valued sequences with marginal distribution belonging to a particular class where the tail decays exponentially. This class does not belong to the domain of attraction of any max-stable distribution. Nevertheless, we prove that the limiting distribution is max-semistable when we consider the maximum of the first k_n observations, for a suitable sequence {_kn}$\{k_n\}$ increasing to infinity. We obtain an expression for calculating the extremal index of sequences satisfying certain local conditions similar to conditions D^(m)(_un)$D^{(m)}(u_n)$, m∈N$m\in \mathbb{N}$, defined by Chernick et al. (1991). We apply the results to a class of max-autoregressive sequences and a class of moving average models. The results generalize the ones obtained for the stationary case.     

Distributions associated with (k1,k2) events on semi-Markov binary trials

We derive neat expressions for the probability generating functions of relevant waiting times associated with (_k1,_k2)$(k_1,k_2)$ events on semi-Markov binary trials. These lead to evaluation of relevant probabilities associated with numbers of occurrence of such events on a string of a fixed length. Our methodology is general enough and provides a template for treating more general events than those of type (_k1,_k2)$(k_1,k_2)$. Also, the same template is extendable to semi-Markov trials with more than two outcomes.     

On the extremal behavior of a nonstationary normal random field

Let X={_Xn}n?1$\mathbf{X}=\{X_{\mathbf{n}}{\}}_{\mathbf{n}?\mathbf{1}}$ be a nonstationary random field satisfying a long range weak dependence for each coordinate at a time and a local dependence condition that avoids clustering of exceedances of high values. For these random fields, the probability of no exceedances of high values can be approximated by exp(−τ)$\exp (-\tau )$, where τ$\tau$ is the limiting mean number of exceedances. We present a class of nonstationary normal random fields for which this result can be applied.     

Empirical Bayes regression analysis with many regressors but fewer observations

In this paper, we consider the prediction problem in multiple linear regression model in which the number of predictor variables, p, is extremely large compared to the number of available observations, n  . The least-squares predictor based on a generalized inverse is not efficient. We propose six empirical Bayes estimators of the regression parameters. Three of them are shown to have uniformly lower prediction error than the least-squares predictors when the vector of regressor variables are assumed to be random with mean vector zero and the covariance matrix (1/n)X^tX$(1/n){\mathit{X}}^{\mathrm{t}}\mathit{X}$ where X^t=(_x1,…,_xn)${\mathit{X}}^{\mathrm{t}}=({\mathit{x}}_1,\dots ,{\mathit{x}}_n)$ is the p×n$p\times n$ matrix of observations on the regressor vector centered from their sample means. For other estimators, we use simulation to show its superiority over the least-squares predictor.     

Local asymptotic minimax theory for block-decreasing densities

In this paper, we study Lebesgue densities on (0,∞)^d$(0,\infty {)}^d$ that are non-increasing in each coordinate, while keeping all other coordinates fixed, from the perspective of local asymptotic minimax lower bound theory. In particular, we establish a local optimal rate of convergence of the order n^−1/(d+2)$n^{-1/(d+2)}$.