A new multivariate model involving geometric sums and maxima of exponentials

We study the joint distribution of (X, Y, N), where N has a geometric distribution while X and Y are, respectively, the sum and the maximum of N i.i.d. exponential random variables. We present fundamental properties of this class of mixed trivariate distributions, and discuss their applications. Our results include explicit formulas for the marginal and conditional distributions, joint integral transforms, moments and related parameters, stability properties, and stochastic representations. We also derive maximum likelihood estimators for the parameters of this distribution, along with their asymptotic properties, and briefly discuss certain generalizations of this model. An example from finance, where N represents the duration of the growth period of the daily log-returns of currency exchange rates, illustrates the modeling potential of this model.     

A mixed bivariate distribution with exponential and geometric marginals

We study the joint distribution of X and N, where N has a geometric distribution and X is the sum of N i.i.d. exponential variables, independent of N. We present basic properties of this class of mixed bivariate distributions, and discuss their possible applications. Our results include marginal and conditional distributions, joint integral transforms, infinite divisibility, and stability with respect to geometric summation. We also discuss maximum likelihood estimation connected with this distribution. An example from finance, where N represents the number of consecutive positive daily log-returns of currency exchange rates, illustrates the modeling potential of these laws.     

On Censored Bivariate Random Variables: Copula,Characterization, and Estimation

Let (X, Y) be a bivariate random vector with joint distribution function F_{X, Y}(x, y) = C(F(x), G(y)), where C is a copula and F and G are marginal distributions of X and Y, respectively. Suppose that (X_i, Y_i), i = 1, 2, …, n is a random sample from (X, Y) but we are able to observe only the data consisting of those pairs (X_i, Y_i) for which X_i ? Y_i. We denote such pairs as (X*_i, Y_i*), i = 1, 2, …, ν, where ν is a random variable. The main problem of interest is to express the distribution function F_{X, Y}(x, y) and marginal distributions F and G with the distribution function of observed random variables X* and Y*. It is shown that if X and Y are exchangeable with marginal distribution function F, then F can be uniquely determined by the distributions of X* and Y*. It is also shown that if X and Y are independent and absolutely continuous, then F and G can be expressed through the distribution functions of X* and Y* and the stress–strength reliability P{X ? Y}. This allows also to estimate P{X ? Y} with the truncated observations (X*_i, Y_i*). The copula of bivariate random vector (X*, Y*) is also derived.     

On the Distribution of Matrix Quadratic Forms

A characterization of the distribution of the multivariate quadratic form given by X A X′, where X is a p × n normally distributed matrix and A is an n × n symmetric real matrix, is presented. We show that the distribution of the quadratic form is the same as the distribution of a weighted sum of non central Wishart distributed matrices. This is applied to derive the distribution of the sample covariance between the rows of X when the expectation is the same for every column and is estimated with the regular mean.     

Two-Sample Tests Based on the Integrated Empirical Process

Abstract

Let X ₁, …, X _m and Y ₁, …, Y _n be independent random variables, where X ₁, …, X _m are i.i.d. with continuous distribution function (df) F, and Y ₁, …, Y _n are i.i.d. with continuous df G. For testing the hypothesis H ₀: F = G, we introduce and study analogues of the celebrated Kolmogorov–Smirnov and one- and two-sided Cramér-von Mises statistics that are functionals of a suitably integrated two-sample empirical process. Furthermore, we characterize those distributions for which the new tests are locally Bahadur optimal within the setting of shift alternatives.     

CHARACTERIZATIONS OF SHIFT-INVARIANT DISTRIBUTIONS BASED ON SUMMATION MODULO ONE

Let X₁Y_1,…, Y_n be independent random variables. We characterize the distributions of X and Y_j satisfying the equation {X+Y₁++Y_n}=_dX, where {Z} denotes the fractional part of a random variable Z. In the case of full generality, either X is uniformly distributed on [0,1), or Y_j has.a shifted lattice distribution and X is shift-invariant. We also give a characterization of shift-invariant distributions. Finally, we consider some special cases of this equation.     

Exact distribution of the convolution of negative binomial random variables

In this note, we derive the exact distribution of S by using the method of generating function and BELL polynomials, where S = X₁ + X₂ + ??? + X_n, and each X_i follows the negative binomial distribution with arbitrary parameters. As a particular case, we also obtain the exact distribution of the convolution of geometric random variables.     

A new class of bivariate lifetime distributions

This paper introduces a new class of bivariate lifetime distributions. Let {X_i}_{i ? 1} and {Y_i}_{i ? 1} be two independent sequences of independent and identically distributed positive valued random variables. Define T₁ = min?(X₁, …, X_M) and T₂ = min?(Y₁, …, Y_N), where (M, N) has a discrete bivariate phase-type distribution, independent of {X_i}_{i ? 1} and {Y_i}_{i ? 1}. The joint survival function of (T₁, T₂) is studied.     

Global convergence of the log-concave MLE when the true distribution is geometric

Let X₁, …, X_n be i.i.d. from a discrete probability mass function (pmf) p. In Balabdaoui et al. [(2013), ‘Asymptotic Distribution of the Discrete Log-Concave mle and Some Applications’, JRSS-B, in press], the pointwise limit distribution of the log-concave maximum-likelihood estimator (MLE) was derived in both the well- and misspecified settings. In the well-specified setting, the geometric distribution was excluded, classified as being degenerate. In this article, we establish the global asymptotic theory of the log-concave MLE of a geometric pmf in all ?_q distances for q∈{1, 2, …}∪{∞}. We also show how these asymptotic results could be used in testing whether a pmf is geometric.     

Limiting distributions of Kolmogorov-Lévy-type statistics under the alternative

Let X_i, 1 ≤ i ≤ n, be independent identically distributed random variables with a common distribution function F, and let G be a smooth distribution function. We derive the limit distribution of □{ρ_α(F_n, G) - _α(F, G)}, where F_n is the empirical distribution function based on X₁,…,X_n and _α is a Kolmogorov-Lévy-type metric between distribution functions. For α ≤ 0 and two distribution functions F and G the metric p_α is given by p_α(F, G) = inf {? ≤ 0: G(x - α?) - ? F(x) ≤ G(x + α?) + ? for all x ?}.     

Asymptotic tail behavior of a random sum with conditionally dependent subexponential summands

In this paper, we obtain some results for the asymptotic behavior of the tail probability of a random sum S_τ = ∑^τ_{k = 1}X_k, where the summands X_k, k = 1, 2, …, are conditionally dependent random variables with a common subexponential distribution F, and the random number τ is a non negative integer-valued random variable, independent of {X_k: k ? 1}.     

Samples of iid Lognormals: Approximations for Characteristics of Minima

We concentrate on characteristics of minima X_N from samples of iid lognormals X_1i ～ 2pLND. We demonstrate that the distribution of X_N for 2 ≤ N ≤ 1,000 may be fitted more accurately by 2pLND than by the limiting Gumbel distribution. An extended power model is established to represent the quotients C_N/C₁, where C_N is the mean, standard deviation, or p-quantile of X_N and C₁ is the corresponding characteristic of X_1i. Our empirical comparisons show that this model provides not only more accurate estimates than alternating approximations but it is also much simpler than its competitors.     

Commuting Matrices in the Queue Length and Sojourn Time Analysis of MAP/MAP/1 Queues

Queues with Markovian arrival and service processes, i.e., MAP/MAP/1 queues, have been useful in the analysis of computer and communication systems and different representations for their stationary sojourn time and queue length distribution have been derived. More specifically, the class of MAP/MAP/1 queues lies at the intersection of the class of QBD queues and the class of semi-Markovian queues. While QBD queues have a matrix exponential representation for their queue length and sojourn time distribution of order N and N², respectively, where N is the size of the background continuous time Markov chain, the reverse is true for a semi-Markovian queue. As the class of MAP/MAP/1 queues lies at the intersection, both the queue length and sojourn time distribution of a MAP/MAP/1 queue has an order N matrix exponential representation. The aim of this article is to understand why the order N² distributions of the sojourn time of a QBD queue and the queue length of a semi-Markovian queue can be reduced to an order N distribution in the specific case of a MAP/MAP/1 queue. We show that the key observation exists in establishing the commutativity of some fundamental matrices involved in the analysis of the MAP/MAP/1 queue.     

Minimum-Distance Estimator for Stable Exponent

Assume that X ₁, X ₂,…, X _n is a sequence of i.i.d. random variables with α-stable distribution (α ∈ (0,2], the stable exponent, is the unknown parameter). We construct minimum distance estimators for α by minimizing the Kolmogorov distance or the Cramér–von-Mises distance between the empirical distribution function G _n , and a class of distributions defined based on the sum-preserving property of stable random variables. The minimum distance estimators can also be obtained by minimizing a U-statistic estimate of an empirical distribution function involving the stable exponent. They share the same invariance property with the maximum likelihood estimates. In this article, we prove the strong consistency of the minimum distance estimators. We prove the asymptotic normality of our estimators. Simulation study shows that the new estimators are competitive to the existing ones and perform very closely even to the maximum likelihood estimator.     

BAYESIAN SMOOTHING WITH OCCASIONAL JUMPS

Consider a sequence of independent random variables X ₁, X ₂,…,X _n observed at n equally spaced time points where X _i has a probability distribution which is known apart from the values of a parameter θ _i ∈ R which may change from observation to observation. We consider the problem of estimating θ = (θ₁, θ₂,…,θ _n ) given the observed values of X ₁, X ₂,…,X _n . The paper proposes a prior distribution for the parameters θ for which sets of parameter values exhibiting no change, or no change apart from a few sudden large changes, or lots of small changes, all have positive prior probability. Markov chain sampling may be used to calculate Bayes estimates of the parameters. We report the results of a Monte Carlo study based on Poisson distributed data which compares the Bayes estimator with estimators obtained using cubic splines and with estimators derived from the Schwarz criterion. We conclude that the Bayes method is preferable in a minimax sense since it never produces the disastrously large errors of the other methods and pays only a modest price for this degree of safety. All three methods are used to smooth mortality rates for oesophageal cancer in Irish males aged 65–69 over the period 1955 through 1994.     

A bivariate distribution with Lomax and geometric margins

We develop a stochastic model describing the joint distribution of $(X,N)$, where $N$ has a geometric distribution while $X$ is the sum of $N$ dependent, heavy-tail Pareto components. Models of this form arise in many applications, ranging from hydro-climatology to finance and insurance. We present fundamental properties of this vector, including marginal and conditional distributions, moments, representations, and parameter estimation. We also include an example from finance, illustrating modeling potential of this new bivariate distribution.     

Prediction based on linear combinations of order statistics and bivariate concomitants in the case of multivariate elliptical distributions

In this paper, by considering a (3n+1) -dimensional random vector (X₀, X^T, Y^T, Z^T)^T having a multivariate elliptical distribution, we derive the exact joint distribution of (X₀, a^TX_(n), b^TY_[n], c^TZ_[n])^T, where a, b, c∈?ⁿ, X_(n)=(X₍₁₎, …, X_(n))^T, X₍₁₎<···<X_(n), is the vector of order statistics arising from X, and Y_[n]=(Y_[1], …, Y_[n])^T and Z_[n]=(Z_[1], …, Z_[n])^T denote the vectors of concomitants corresponding to X_(n) ((Y_[r], Z_[r])^T, for r=1, …, n, is the vector of bivariate concomitants corresponding to X_(r)). We then present an alternate approach for the derivation of the exact joint distribution of (X₀, X_(r), Y_[r], Z_[r])^T, for r=1, …, n. We show that these joint distributions can be expressed as mixtures of four-variate unified skew-elliptical distributions and these mixture forms facilitate the prediction of X_(r), say, based on the concomitants Y_[r] and Z_[r]. Finally, we illustrate the usefulness of our results by a real data.     

Upper tolerance limit for the largest observation from censored samples

This paper deals with obtaining an upper tolerance limit for a largest observation X_(n) in an ordered sample of size n from a continuous distribution where the first m observations X₍₁₎ < X₍₂₎ < … < X_(m), l ≤ m < n, have been observed. A criterion of “goodness” of tolerance limit is developed, and a method is given to obtain the best tolerance limit. This method is applied to exponential and Pareto distributions.