Likelihood-based inference for a class of multivariate diffusions with unobserved paths

This paper presents a Markov chain Monte Carlo algorithm for a class of multivariate diffusion models with unobserved paths. This class is of high practical interest as it includes most diffusion driven stochastic volatility models. The algorithm is based on a data augmentation scheme where the paths are treated as missing data. However, unless these paths are transformed so that the dominating measure is independent of any parameters, the algorithm becomes reducible. The methodology developed in Roberts and Stramer [2001a. On inference for partial observed nonlinear diffusion models using the metropolis-hastings algorithm. Biometrika 88(3); 603–621] circumvents the problem for scalar diffusions. We extend this framework to the class of models of this paper by introducing an appropriate reparametrisation of the likelihood that can be used to construct an irreducible data augmentation scheme. Practical implementation issues are considered and the methodology is applied to simulated data from the Heston model.     

Some observations on the KMT dyadic scheme

We point out some useful properties of the KMT dyadic scheme, which was used to construct a Brownian bridge approximation to the uniform empirical process.     

Minimax optimality of the Shiryayev–Roberts change-point detection rule

The result of Pollak [1985. Optimal detection of a change in distribution. Ann. Statist. 13, 206–227] proving the asymptotic optimality in sequential change-point detection of a suitable Shirayayev–Roberts stopping rule up to terms that vanish in the limit is generalized from the case of two completely specified distributions to that of a composite alternative hypothesis in a multidimensional exponential family. An explicit asymptotic lower bound on the expected Kullback–Leibler information required to detect a change-point is derived and is shown to be attained by a Shirayayev–Roberts stopping rule.     

Multivariate spatial U-quantiles: A Bahadur–Kiefer representation,a Theil–Sen estimator for multiple regression,and a robust dispersion estimator

A leading multivariate extension of the univariate quantiles is the so-called “spatial” or “geometric” notion, for which sample versions are highly robust and conveniently satisfy a Bahadur–Kiefer representation. Another extension of univariate quantiles has been to univariate U-quantiles, on the basis of which, for example, the well-known Hodges–Lehmann location estimator has a natural formulation. Generalizing both extensions, we introduce multivariate spatial U-quantiles and develop a corresponding Bahadur–Kiefer representation. New statistics based on spatial U-quantiles are presented for nonparametric estimation of multiple regression coefficients, extending the classical Theil–Sen nonparametric simple linear regression slope estimator, and for robust estimation of multivariate dispersion. Some other applications are mentioned as well.     

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.     

Union–intersection principle and constrained statistical inference

Most statistical models arising in real life applications as well as in interdisciplinary research are complex in their designs, sampling plans, and associated probability laws, which in turn are often constrained by inequality, order, functional, shape or other restraints. Optimality of conventional likelihood ratio based statistical inference may not be tenable here, although the use of restricted or quasi-likelihood has spurred in such environments. S.N. Roy's ingenious union–intersection principle provides an alternative avenue, often having some computational advantages, increased scope of adaptability, and flexibility beyond conventional likelihood paradigms. This scenario is appraised here with some illustrative examples, and with some interesting problems of inference on stochastic ordering (dominance) in parametric as well as beyond parametric setups.     

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.     

Truncated sequential change-point detection based on renewal counting processes II

The standard approach in change-point theory is to base the statistical analysis on a sample of fixed size. Alternatively, one observes some random phenomenon sequentially and takes action as soon as one observes some statistically significant deviation from the “normal” behaviour. The present paper is a continuation of Gut and Steinebach [2002. Truncated sequential change-point detection based on renewal counting processes. Scand. J. Statist. 29, 693–719] the main point being that here we look in more detail into the behaviour of the relevant stopping times, in particular the time it takes from the actual change-point until the change is detected, more precisely, we prove asymptotics for stopping times under alternatives.     

A probablistic analysis for an optimal screening problem

Summary We consider a lotL formed byN apparently similar unitsW ₁,…,W _N, where each of theW _i may come from one of two different populationsP ₁ andP ₂;T ₁,…,T _N denote the corresponding lifetimes. The units fromP _i undergo a failure of kindi and their survival function isS _i (t). We assume that the failure rate function are known and that the units fromP ₁ are ?substandard?: λ ₁ (t)≥λ ₂ (t), ∀t≥0. We want to putW ₁,…,W _N under a pre-operational test (burn-in test) in order to eliminate at least a great part of the substandard units and we face the problem of obtaining a rule for stopping the test under the assumption that, with the failure of a unit, it is possible to recognize the population from which the unit comes. Such a problem will be formalized as an optimal stopping problem for a suitably defined Markov process. Our study shall evidentiate some fundamental aspects of the problem and the role of the prior distribution of the (random) numberM ₀ of those units inL coming fromP ₁ (substandard). The latter distribution has a great influence on the form of the solution. This research was supported by the C.N.R. Project ?Statistica Bayesiana e Simulazione in Affidalità e Modellistica Biologica?.     

Joint modeling of degradation and failure time data

This paper surveys some approaches to model the relationship between failure time data and covariate data like internal degradation and external environmental processes. These models which reflect the dependency between system state and system reliability include threshold models and hazard-based models. In particular, we consider the class of degradation–threshold–shock models (DTS models) in which failure is due to the competing causes of degradation and trauma. For this class of reliability models we express the failure time in terms of degradation and covariates. We compute the survival function of the resulting failure time and derive the likelihood function for the joint observation of failure times and degradation data at discrete times. We consider a special class of DTS models where degradation is modeled by a process with stationary independent increments and related to external covariates through a random time scale and extend this model class to repairable items by a marked point process approach. The proposed model class provides a rich conceptual framework for the study of degradation–failure issues.     

Characterizations of exponentiality within the HNBUE class and related tests

A harmonic new better than used in expectation (HNBUE) variable is a random variable which is dominated by an exponential distribution in the convex stochastic order. We use a recently obtained condition on stochastic equality under convex domination to derive characterizations of the exponential distribution and bounds for HNBUE variables based on the mean values of the order statistics of the variable. We apply the results to generate discrepancy measures to test if a random variable is exponential against the alternative that is HNBUE, but not exponential.     

Parameter estimation of state space models for univariate observations

This paper contributes to the problem of estimation of state space model parameters by proposing estimators for the mean, the autoregressive parameters and the noise variances which, contrarily to maximum likelihood, may be calculated without assuming any specific distribution for the errors. The estimators suggested widen the scope of the application of the generalized method of moments to some heteroscedastic models, as in the case of state-space models with varying coefficients, and give sufficient conditions for their consistency. The paper includes a simulation study comparing the proposed estimators with maximum likelihood estimators. Finally, these methods are applied to the calibration of the meteorological radar and estimation of area rainfall.     

A consistent system of conditional probabilities which is not compatible with any random field

A consistent, Markovian family of conditional densities is constructed which is not compatible with any random field.     

Asymptotic distributions of statistics and parameter estimates for mixed Poisson processes

Mixed Poisson processes have been used as natural models for events occurring in continuous or discrete time. Our main result is the derivation of the joint asymptotic distributions of statistics, including parameter estimators, computed in different time intervals from data generated by mixed Poisson processes. These distributions can be used, for example, to test the hypothesis about the adequacy of the mixed Poisson process against data. We provide some simulation results and test the model on actual market research data.     

On quadratic prediction problems in finite populations

For a finite population and its linear model, Liu and Rong proposed a notion of optimal invariant quadratic unbiased prediction (OIQUP) and offered two methods for studying this notion, in which the first is incomplete. In this note, we mainly aim at fulfilling the first approach used by Liu and Rong by considering a transformed matrix equation set through permutation matrix techniques. Solvability of the matrix equation set, optimality of the resulting predictor, and equivalence of the representations of OIQUP, derived in this note and by Liu and Rong, are investigated in detail. In addition, an application to predicting population variance is conducted based on a simulated population.     

Stopping times of one-sample rank order sequential probability ratio tests

Weed, Bradley and Grovindarajulu (1974) propose one-sample probability ratio tests based on Lehmann alternatives. They also study the finite sure termination of the stopping times. Motivated by Stein's proof of (1946) of the termination of a sequential probability ratio test (SPRT) in the case of independent and identically distributed (i.i.d.) random variables and the work of Sethuraman (1970) for the two- sample rank order SPRT, we obtain a very mild condition (namely, that a certain random variable U(Z) is not identically zero) for the finite sure termination of the existence of the moment generating function (m.g.f.) for the stopping time of one- sample rank order SPRT's.     

Some aspects of decision making under uncertainty

Two discrete-time insurance models are studied in the framework of cost approach. The models being non-deterministic one deals with decision making under uncertainty. Three different situations are investigated: (1) underlying processes are stochastic however their probability distributions are given; (2) information concerning the distribution laws is incomplete; (3) nothing is known about the processes under consideration. Mathematical methods useful for establishing the (asymptotically) optimal control are demonstrated in each case. Algorithms for calculation of critical levels are proposed. Numerical results are presented as well.     

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.     

Convex rearrangements,generalized Lorenz curves,and correlated Gaussian data

We propose a statistical index for measuring the fluctuations of a stochastic process ξ$\xi$. This index is based on the generalized Lorenz curves and (modified) Gini indices of econometric theory.