Moments of Special Normally Distributed Matrices

In solving systems of simultaneous random linear algebraic equations some approximating methods lead to the problem of determinating moments of special random matrices and vectors. In this article corresponding formulas are provided for moments of some normally distributed matrices. The deduced relations can be considered as a generalization of the known formulas for the central moments of normally distributed random variables.     

The tube method for the moment index in projection pursuit

The projection pursuit index defined by a sum of squares of the third and the fourth sample cumulants is known as the moment index proposed by Jones and Sibson [1987. What is projection pursuit? J. Roy. Statist. Soc. Ser. A 150, 1–36]. The limiting distribution of the maximum of the moment index under the null hypothesis that the population is multivariate normal is shown to be the maximum of a Gaussian random field with a finite Karhunen–Loève expansion. An approximate formula for tail probability of the maximum, which corresponds to the p-value, is given by virtue of the tube method through determining Weyl's invariants of all degrees and the critical radius of the index manifold of the Gaussian random field.     

Occupancy numbers for dynamic heterogeneous populations: Estimate of particles lifetimes

A dynamic model of a heterogeneous population is studied. Particles belonging to a population are divided, at every time t, into a finite number of classes according to their types and the partition changes over time. The role of the occupancy numbers, namely the cardinality of each class, is highlighted. The relationship between the stochastic process of occupancy numbers and the process of particle types is analyzed. The main goal of this paper is the estimation of the lifetime of each particle at a given time t, when the observed data are the history of the process of the number of dead particles up to t. Furthermore, a discrete time approximation of the filter is given.     

A generalization of the multivariate slash distribution

In this paper, we propose a generalization of the multivariate slash distribution and investigate some of its properties. We show that the new distribution belongs to the elliptically contoured distributions family, and can have heavier tails than the multivariate slash distribution. Therefore, this generalization of the multivariate slash distribution can be considered as an alternative heavy-tailed distribution for modeling data sets in a variety of settings. We apply the generalized multivariate slash distribution to two real data sets to provide some illustrative examples.     

Finite mixtures of multivariate Poisson distributions with application

In the present paper we examine finite mixtures of multivariate Poisson distributions as an alternative class of models for multivariate count data. The proposed models allow for both overdispersion in the marginal distributions and negative correlation, while they are computationally tractable using standard ideas from finite mixture modelling. An EM type algorithm for maximum likelihood (ML) estimation of the parameters is developed. The identifiability of this class of mixtures is proved. Properties of ML estimators are derived. A real data application concerning model based clustering for multivariate count data related to different types of crime is presented to illustrate the practical potential of the proposed class of models.     

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.     

Stochastic orderings related to ranking in an nth order record process

Let X₁, X₂, …, X_m be successive observations on m objects, numbered 1,2, …, m. If X₁ belongs to the n largest observations among X₁,X₂,…,X_i, object i is called a record, i = 1,2,…,m; n?1.In investigating the influence of the ranking of the objects on the expected number of records, a hierarchy of stochastic order relations between random variables arises.It is these order relations and their relationship with known stochastic orderings that are studied in this paper.     

Remarks on an integral functional driven by sub-fractional Brownian motion

This paper studies the functionals where is a one-dimension sub-fractional Brownian motion with index H∈(0,1). It shows that there exists a constant p_H∈(1,2) such that p-variation of the process (j=1,2) is equal to 0 if p>p_H, where ?_j, j=1,2, are the local time and weighted local time of S^H, respectively. This extends the classical results for Brownian motion.     

Multivariate Matsumoto–Yor property is rather restrictive

Matsumoto and Yor [2001. An analogue of Pitman's 2M-X$2M-X$ theorem for exponential Wiener functionals. Part II: the role of the GIG laws. Nagoya Math. J. 162, 65–86] discovered an interesting invariance property of a product of the generalized inverse Gaussian (GIG) and the gamma distributions. For univariate random variables or symmetric positive definite random matrices it is a characteristic property for this pair of distributions. It appears that for random vectors the Matsumoto–Yor property characterizes only very special families of multivariate GIG and gamma distributions: components of the respective random vectors are grouped into independent subvectors, each subvector having linearly dependent components. This complements the version of the multivariate Matsumoto–Yor property on trees and related characterization obtained in Massam and Weso?owski [2004. The Matsumoto–Yor property on trees. Bernoulli 10, 685–700].     

Geostochastic calculus

A stochastic calculus for a family of continuous measure-valued Markov processes is developed. Such processes arise naturally in the construction of stochastic models of spatially distributed populations. The stochastic calculus is a tool whereby a class of density-dependent models can be studied in terms of the multiplicative measure diffusion process. In this paper the stochastic integral is introduced in the space-time setting and a Cameron-Martin-Girsanov theorem is established.     

Bayesian variable window scan statistics

In this article two-stage hierarchical Bayesian models are used for the observed occurrences of events in a rectangular region. Two Bayesian variable window scan statistics are introduced to test the null hypothesis that the observed events follow a specified two-stage hierarchical model vs an alternative that indicates a local increase in the average number of observed events in a subregion (clustering). Both procedures are based on a sequence of Bayes factors and their p$p$-values that have been generated via simulation of posterior samples of the parameters, under the null and alternative hypotheses. The posterior samples of the parameters have been generated by employing Gibbs sampling via introduction of auxiliary variables. Numerical results are presented to evaluate the performance of these variable window scan statistics.     

Birth-multiple catastrophe processes

Birth-multiple catastrophe processes are analyzed where the birth transition rates are assumed to be constant while catastrophes are distinguished by having possibly different destinations and possibly different transition rates. The transient probability functions of such birth-multiple catastrophe systems are determined. The solution method uses dual processes, randomization, and sample path counting. Solutions are explicit in terms of being a finite linear combination of products of exponential functions of time, t, and nonnegative integer powers of t. The coefficients within this expansion follow a pattern of rational functions of the transition rates.     

Limiting behavior of the maximum of the partial sum for asymptotically negatively associated random variables under residual Cesàro alpha-integrability assumption

Asymptotically negative association is a special dependence structure. By relating such dependence condition to residual Cesàro alpha-integrability and to strongly residual Cesàro alpha-integrability, some L_p-convergence and complete convergence results of the maximum of the partial sum are derived, respectively. In addition, some of these conclusions are based on a new Rosenthal type inequality concerning asymptotically negatively associated random variables, which is of independent interest.     

A note on replacement policy comparisons from NBUC lifetime of the unit

Age and block replacement policies are commonly used in order to reduce the number of in-service failures when the systems are functioning indefinitely. In reliability theory, the lifetime of a system can be modeled by means of the NBUC aging class that is characterized throughout comparisons of the residual lives in the sense of the icx order. The purpose of this paper is to establish stochastic comparisons between the age (block) replacement policy and a renewal process with no planned replacements when the lifetime of the unit is NBUC. Supported by Ministerio de Ciencia y Tecnología under grant BFM2000-0362     

Extremal limit theorems for observations separated by random power law waiting times

This paper develops extreme value theory for random observations separated by random waiting times whose exceedence probability falls off like a power law. In the case where the waiting times between observations have an infinite mean, a limit theorem is established, where the limit is comprised of an extremal process whose time index is randomized according to the non-Markovian hitting time process for a stable subordinator. The resulting limit distributions are shown to be solutions of fractional differential equations, where the order of the fractional time derivative coincides with the power law index of the waiting time. The probability that the limit process remains below a threshold is also computed. For waiting times with finite mean but infinite variance, a two-scale argument yields a fundamentally different limit process. The resulting limit is an extremal process whose time index is randomized according to the first passage time of a positively skewed stable Lévy motion with positive drift. This two-scale limit provides a second-order correction to the usual limit behavior.     

On Vervaat and Vervaat-error-type processes for partial sums and renewals

We study the asymptotic behaviour of stochastic processes that are generated by sums of partial sums of i.i.d. random variables and their renewals. We conclude that these processes cannot converge weakly to any nondegenerate random element of the space D[0,1]$D[0,1]$. On the other hand, we show that their properly normalized integrals as Vervaat-type stochastic processes converge weakly to a squared Wiener process. Moreover, we also deal with the asymptotic behaviour of the deviations of these processes, the so-called Vervaat-error-type processes.     

Construction of non-Gaussian random fields with any given correlation structure

A positive definite function can be thought of as the covariance function of a Gaussian random field, according to the celebrated Kolmogorov existence theorem. A question of great theoretical and practical interest is: how could one construct a non-Gaussian random field with the given positive definite function as its covariance function? In this paper we demonstrate a novel and simple method for constructing many such non-Gaussian random fields, with the corresponding finite-dimensional distributions identified. Also, we show how to construct a non-Gaussian random field with a given negative definite function as its variogram.