Random spectral measure for non Gaussian moving averages

We study the distribution of phases and amplitudes for the spectral representation of weighted moving averages of a general noise measure. The simple independent structure, known for the Gaussian case, and involving Rayleigh amplitude and uniform phase distributions, is lost for the non Gaussian noise case. We show that the amplitude/phase distributions exhibit a rich and more complex structure depending not just on the covariance of the process but specifically on the form of the kernel and the noise distribution. We present a theoretical tool for studying these distributions that follows from a proof of the spectral theorem that yields an explicit expression for the spectral measure. The main interest is in noise measures based on second-order Lévy motions since such measures are easily available through independent sampling. We approximate the spectral stochastic measure by independent noise increments which allows us to obtain amplitude/phase distributions that is of fundamental interest for analyzing processes in the frequency domain. For the purpose of approximating the moving average process through sums of trigonometric functions, we assess the mean square error of discretization of the spectral representation. For a specified accuracy, the approximation is explicitly given. We illustrate the method for the moving averages driven by the Laplace motion.     

Bayesian analysis of semiparametric Bernstein polynomial regression models for data with sample selection

In regression analysis, a sample selection scheme often applies to the response variable, which results in missing not at random observations on the variable. In this case, a regression analysis using only the selected cases would lead to biased results. This paper proposes a Bayesian methodology to correct this bias based on a semiparametric Bernstein polynomial regression model that incorporates the sample selection scheme into a stochastic monotone trend constraint, variable selection, and robustness against departures from the normality assumption. We present the basic theoretical properties of the proposed model that include its stochastic representation, sample selection bias quantification, and hierarchical model specification to deal with the stochastic monotone trend constraint in the nonparametric component, simple bias corrected estimation, and variable selection for the linear components. We then develop computationally feasible Markov chain Monte Carlo methods for semiparametric Bernstein polynomial functions with stochastically constrained parameter estimation and variable selection procedures. We demonstrate the finite-sample performance of the proposed model compared to existing methods using simulation studies and illustrate its use based on two real data applications.     

Modeling of semi-competing risks by means of first passage times of a stochastic process

In semi-competing risks one considers a terminal event, such as death of a person, and a non-terminal event, such as disease recurrence. We present a model where the time to the terminal event is the first passage time to a fixed level c in a stochastic process, while the time to the non-terminal event is represented by the first passage time of the same process to a stochastic threshold S, assumed to be independent of the stochastic process. In order to be explicit, we let the stochastic process be a gamma process, but other processes with independent increments may alternatively be used. For semi-competing risks this appears to be a new modeling approach, being an alternative to traditional approaches based on illness-death models and copula models. In this paper we consider a fully parametric approach. The likelihood function is derived and statistical inference in the model is illustrated on both simulated and real data.     

Small Sample Tests for Shape Parameters of Gamma Distributions

The introduction of shape parameters into statistical distributions provided flexible models that produced better fit to experimental data. The Weibull and gamma families are prime examples wherein shape parameters produce more reliable statistical models than standard exponential models in lifetime studies. In the presence of many independent gamma populations, one may test equality (or homogeneity) of shape parameters. In this article, we develop two tests for testing shape parameters of gamma distributions using chi-square distributions, stochastic majorization, and Schur convexity. The first one tests hypotheses on the shape parameter of a single gamma distribution. We numerically examine the performance of this test and find that it controls Type I error rate for small samples. To compare shape parameters of a set of independent gamma populations, we develop a test that is unbiased in the sense of Schur convexity. These tests are motivated by the need to have simple, easy to use tests and accurate procedures in case of small samples. We illustrate the new tests using three real datasets taken from engineering and environmental science. In addition, we investigate the Bayes’ factor in this context and conclude that for small samples, the frequentist approach performs better than the Bayesian approach.     

Filtration of ASTA: a weak convergence approach

Consider a stochastic process (X,A), where X represents the evolution of a system over time, and A is an associated point process that has stationary independent increments. Suppose we are interested in estimating the time average frequency of the process X being in a set of states. Often it is more convenient to have a sampling procedure for estimating the time average based on averaging the observed values of X(T_n) (T_n being a point of A) over a long period of time: the event average of the process. In this paper we examine the situation when the two procedures—event averaging and time averaging—produce the same estimate (the ASTA property: Arrivals See Time Averages). We prove a result stronger than ASTA. Under a lack-of-anticipation assumption we prove that the point process, A, restricted to any set of states, has the same probabilistic structure as the original point process. In particular, if the original point process is Poisson the new point process is still Poisson with the same parameter as the original point process. We develop our results in the more general setting of a stochastic process (X,A), that is, a process with an imbedded cumulative process, A={A(t),t0}, which is assumed to be a Levy process with non-decreasing sample paths. This framework allows for modeling fluid processes, as well as compound Poisson processes with non-integer increments. First, we state the result in discrete time; the discrete-time result is then extended to the continuous-time case using limiting arguments and weak-convergence theory. As a corollary we give a proof of ASTA under weak conditions and a simple, intuitive proof of (Poisson Arrivals See Time Averages) under the standard conditions. The results are useful in queueing and statistical sampling theory.     

Linear Increments with Non‐monotone Missing Data and Measurement Error

Linear increments (LI) are used to analyse repeated outcome data with missing values. Previously, two LI methods have been proposed, one allowing non‐monotone missingness but not independent measurement error and one allowing independent measurement error but only monotone missingness. In both, it was suggested that the expected increment could depend on current outcome. We show that LI can allow non‐monotone missingness and either independent measurement error of unknown variance or dependence of expected increment on current outcome but not both. A popular alternative to LI is a multivariate normal model ignoring the missingness pattern. This gives consistent estimation when data are normally distributed and missing at random (MAR). We clarify the relation between MAR and the assumptions of LI and show that for continuous outcomes multivariate normal estimators are also consistent under (non‐MAR and non‐normal) assumptions not much stronger than those of LI. Moreover, when missingness is non‐monotone, they are typically more efficient.     

The Ornstein-Uhlenbeck Dirichlet process and other time-varying processes for Bayesian nonparametric inference

This paper introduces a new class of time-varying, measure-valued stochastic processes for Bayesian nonparametric inference. The class of priors is constructed by normalising a stochastic process derived from non-Gaussian Ornstein-Uhlenbeck processes and generalises the class of normalised random measures with independent increments from static problems. Some properties of the normalised measure are investigated. A particle filter and MCMC schemes are described for inference. The methods are applied to an example in the modelling of financial data.     

Threshold Estimation for Stochastic Processes with Small Noise

Consider a process satisfying a stochastic differential equation with unknown drift parameter, and suppose that discrete observations are given. It is known that a simple least squares estimator (LSE) can be consistent but numerically unstable in the sense of large standard deviations under finite samples when the noise process has jumps. We propose a filter to cut large shocks from data and construct the same LSE from data selected by the filter. The proposed estimator can be asymptotically equivalent to the usual LSE, whose asymptotic distribution strongly depends on the noise process. However, in numerical study, it looked asymptotically normal in an example where filter was chosen suitably, and the noise was a Lévy process. We will try to justify this phenomenon mathematically, under certain restricted assumptions.     

APPROXIMATION OF EXCEEDANCE PROCESSES IN LARGE POPULATIONS

We consider a population of n individuals. Each of these individuals generates a discrete time branching stochastic process. We study the number of ancestors S(n,t) whose offspring at time t exceeds level θ(t), where θ(t) is some positive valued function. It is proved that S(n,t) may be approximated as t → ∞ and n → ∞ by some stochastic processes with independent increments.

     

Uniform Convergence of the Non-Weighted Poverty Measures

We consider in this note the weak convergence, in the frame of the empirical processes theory, of the nonweighted poverty measures viewed as stochastic processes defined on some space of bounded functions and indexed by real numbers or monotone functions. The results include the asymptotic behavior of the Foster–Greer–Thorbecke process of poverty indices. We use them to follow up the poverty evolution in poor countries between two periods with appropriate curves.     

Computation of the Two-Sample Smirnov Statistics

A simple method is given for evaluating the one- or two-sided Smirnov statistic for comparing two independent samples from continuous populations.     

Property of bivariate poisson distribution and its application to stochastic processes

Assuming that the random vectors X ₁ and X ₂ have independent bivariate Poisson distributions, the conditional distribution of X ₁ given X ₁?+?X ₂?=?n is obtained. The conditional distribution turns out to be a finite mixture of distributions involving univariate binomial distributions and the mixing proportions are based on a bivariate Poisson (BVP) distribution. The result is used to establish two properties of a bivariate Poisson stochastic process which are the bivariate extensions of the properties for a Poisson process given by Karlin, S. and Taylor, H. M. (1975). A First Course in Stochastic Processes, Academic Press, New York.     

Bathtub distributions: a review

Bathtub distributions are characterized by bathtub failure rate functions . These are possibly more realisitic models than the monotone failure rate models . A systematic account of such distributions is not available and this review aims to give such an account . We give some easily verifiable conditions to check the bathtub property of a distribution along with methods to construct such distributions . We also discuss some stochastic and reliablity mechanisms which lead to bathtub distributions. These include mixtures ( stochastic failure rate models ) , series system , stochastic differential equation models and so on. We also review inference on bathtub distributions. The paper concludes with a rather exhaustive list of bathtub distributions.     

Conjugacy as a Distinctive Feature of the Dirichlet Process

Abstract.  Recently the class of normalized random measures with independent increments, which contains the Dirichlet process as a particular case, has been introduced. Here a new technique for deriving moments of these random probability measures is proposed. It is shown that, a priori , most of the appealing properties featured by the Dirichlet process are preserved. When passing to posterior computations, we obtain a characterization of the Dirichlet process as the only conjugate member of the whole class of normalized random measures with independent increments.     

The finite-time ruin probability of a discrete-time risk model with GARCH discounted factors and dependent risks

The finite-time ruin probability of a discrete-time risk model with dependent stochastic discount factors and dependent insurance and financial risks is investigated in this paper. Assume that the stochastic discount factors follow a GARCH process and the one-period insurance and financial risks form a sequence of independent and identically distributed random pairs, which are the copies of a random pair with a bivariate Sarmanov dependent distribution. When the common distribution of claim-sizes is heavy-tailed, we establish an asymptotic estimate for the finite-time ruin probability. Applying the result to a special case, we also get conservative asymptotic bounds. A numerical simulation is given at the end of the paper.     

Robust Model Selection for Stochastic Processes

We address the problem of robust model selection for finite memory stochastic processes. Consider m independent samples, with most of them being realizations of the same stochastic process with law Q, which is the one we want to retrieve. We define the asymptotic breakdown point γ for a model selection procedure and also we devise a model selection procedure. We compute the value of γ which is 0.5, when all the processes are Markovian. This result is valid for any family of finite order Markov models but for simplicity we will focus on the family of variable length Markov chains.     

A Robust Test for Weak Instruments

This article deals with the estimation of continuous-time stochastic volatility models of option pricing. We argue that option prices are much more informative about the parameters than are asset prices. This is confirmed in a Monte Carlo experiment that compares two very simple strategies based on the different information sets. Both approaches are based on indirect inference and avoid any discretization bias by simulating the continuous-time model. We assume an Ornstein-Uhlenbeck process for the log of the volatility, a zero-volatility risk premium, and no leverage effect. We do not pursue asymptotic efficiency or specification issues; rather, we stick to a framework with no overidentifying restrictions and show that, given our option-pricing model, estimation based on option prices is much more precise in samples of typical size, without increasing the computational burden.     

On the Bootstrap and Monotone Likelihood in the Cox Proportional Hazards Regression Model

Recent literature has provided encouragement for using the bootstrap for inference on regression parameters in the Cox proportional hazards (PH) model. However, generating and performing the necessary partial likelihood computations on multitudinous bootstrap samples greatly increases the chances of incurring problems with monotone likelihood at some point in the analysis. The only symptom of monotone likelihood may be a failure to converge in the numerical maximization procedure, and so the problem might naively be dismissed by deleting the offending data set and replacing it with a new one. This strategy is shown to lead to potentially high selection biases in the subsequent summary statistics. This note discusses the importance of keeping track of these monotone likelihood cases and provides recommendations for their use in interpreting bootstrap findings, and for avoiding unwanted biases that may result from high rates of occurrence. In many cases, high monotone likelihood rates indicate that a more highly-specified model may be preferred. Special consideration is given to the problem of high monotone likelihood incidence in Monte Carlo studies of the bootstrap.     

Goodness-of-fit test for exponentiality based on spacings for general progressive Type-II censored data

There have been numerous tests proposed to determine whether or not the exponential model is suitable for a given data set. In this article, we propose a new test statistic based on spacings to test whether the general progressive Type-II censored samples are from exponential distribution. The null distribution of the test statistic is discussed and it could be approximated by the standard normal distribution. Meanwhile, we propose an approximate method for calculating the expectation and variance of samples under null hypothesis and corresponding power function is also given. Then, a simulation study is conducted. We calculate the approximation of the power based on normality and compare the results with those obtained by Monte Carlo simulation under different alternatives with distinct types of hazard function. Results of simulation study disclose that the power properties of this statistic by using Monte Carlo simulation are better for the alternatives with monotone increasing hazard function, and otherwise, normal approximation simulation results are relatively better. Finally, two illustrative examples are presented.